Dear Shri,

Thanks a lot. I have tried finding the jacobian matrix when i ran the power 
flow on case 9 in the example cases given in matlab.
I assigned fullJac to be = 1.
However, I obtained the following values for J:
J =

   (1,1)      17.1204
   (4,1)     -17.1204
  (10,1)       0.7195
  (13,1)       0.7195
   (2,2)      15.8554
   (8,2)     -15.8554
  (11,2)       1.6300
  (17,2)       1.6300
   (3,3)      17.1013
   (6,3)     -17.1013
  (12,3)       0.8500
  (15,3)       0.8500
   (1,4)     -17.1204
   (4,4)      38.2940
   (5,4)     -10.0631
   (9,4)     -10.9179
  (10,4)      -0.7195
  (13,4)      -3.2220
  (14,4)       2.1536
  (18,4)       1.6615
   (4,5)     -10.1682
   (5,5)      15.3734
   (6,5)      -5.5701
  (13,5)       1.5848
  (14,5)      -3.9680
  (15,5)       0.6817
   (3,6)     -17.1013
   (5,6)      -5.3103
   (6,6)      32.3713
   (7,6)      -9.6479
  (12,6)      -0.8500
  (14,6)       1.8143
  (15,6)      -2.4536
  (16,6)       1.3623
   (6,7)      -9.6999
   (7,7)      22.9890
   (8,7)     -13.5172
  (15,7)       0.9218
  (16,7)      -3.6932
  (17,7)       0.8399
   (2,8)     -15.8554
   (7,8)     -13.3411
   (8,8)      35.1757
   (9,8)      -5.4819
  (11,8)      -1.6300
  (16,8)       2.3309
  (17,8)      -2.7834
  (18,8)       1.9295
   (4,9)     -11.0054
   (8,9)      -5.8031
   (9,9)      16.3998
  (13,9)       0.9177
  (17,9)       0.3135
  (18,9)      -3.5910
   (1,10)      0.7195
   (4,10)     -0.7195
  (10,10)     17.6018
  (13,10)    -17.1204
   (2,11)      1.6300
   (8,11)     -1.6300
  (11,11)     16.1446
  (17,11)    -15.8554
   (3,12)      0.8500
   (6,12)     -0.8500
  (12,12)     17.0284
  (15,12)    -17.1013
   (1,13)      0.7290
   (4,13)      3.2644
   (5,13)     -2.1820
   (9,13)     -1.6834
  (10,13)    -17.3458
  (13,13)     38.7981
  (14,13)    -10.1956
  (18,13)    -11.0616
   (4,14)     -1.6246
   (5,14)      2.2225
   (6,14)     -0.6989
  (13,14)    -10.4239
  (14,14)     15.1448
  (15,14)     -5.7102
   (3,15)      0.8471
   (5,15)     -1.8082
   (6,15)      2.4453
   (7,15)     -1.3577
  (12,15)    -17.0438
  (14,15)     -5.2924
  (15,15)     32.2624
  (16,15)     -9.6154
   (6,16)     -0.9353
   (7,16)      1.7179
   (8,16)     -0.8521
  (15,16)     -9.8411
  (16,16)     22.6136
  (17,16)    -13.7140
   (2,17)      1.6362
   (7,17)     -2.3398
   (8,17)      2.7940
   (9,17)     -1.9369
  (11,17)    -15.9161
  (16,17)    -13.3922
  (17,17)     35.3104
  (18,17)     -5.5029
   (4,18)     -0.9583
   (8,18)     -0.3274
   (9,18)      1.1393
  (13,18)    -11.4924
  (17,18)     -6.0599
  (18,18)     16.0813



I was wondering why some of the values are missing and why it is an 18*18 
matrix.

What I want is J4 for each bus, in order to find the voltage sensitivity.
I will be extremely grateful if you could please help me out as my project 
deadline is next week.

Thank you
Hiranya
On 2 Jun 2015, at 17:08, Abhyankar, Shrirang G. <[email protected]> wrote:

> fullJac is a flag that you can set to either access the full Jacobian
> matrix or a partial Jacobian. With fullJac = 1, the routine returns the
> full Jacobian that has partial derivatives of the power flow equations
> w.r.t to the voltage magnitudes and angles at all buses. If fullJac = 0 or
> not given, the Jacobian as used for the power flow calculation, i.e.,
> Jacobian ignoring reactive power equations and voltage magntidues at PV
> buses, real, reactive power equations, voltage magnitudes and angles at
> swing buses. 
> 
> Shri
> 
> -----Original Message-----
> From: Hiranya Seneviratne <[email protected]>
> Reply-To: MATPOWER discussion forum <[email protected]>
> Date: Tuesday, June 2, 2015 at 10:44 AM
> To: MATPOWER discussion forum <[email protected]>
> Subject: Re: Jacobian matrix values
> 
>> Thank you. 
>> However what is fulljac? As I am unable to find the value for that.
>> Thank you
>> 
>> Sent from my iPhone
>> 
>>> On 2 Jun 2015, at 16:29, "Abhyankar, Shrirang G." <[email protected]>
>>> wrote:
>>> 
>>> http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.1/makeJac.html
>>> 
>>> 
>>> -----Original Message-----
>>> From: Hiranya Seneviratne <[email protected]>
>>> Reply-To: MATPOWER discussion forum <[email protected]>
>>> Date: Tuesday, June 2, 2015 at 4:04 AM
>>> To: MATPOWER discussion forum <[email protected]>
>>> Subject: Jacobian matrix values
>>> 
>>>> Dear Ray,
>>>> Is it possible to find out the values of the jacobian matrix once the
>>>> power flow is calculated in Matpower?
>>>> Thank you
>>>> Hiranya 
>>>> 
>>>> Sent from my iPhone
>>> 
>>> 
>>> 
>> 
> 
> 
> 



Reply via email to