Dear Shri, Thank you. For each bus the jacobian is stored in the way you mentioned? But how would I extract it, for example in the 9 bus system if I want the 4th element of each jacobian matrix for all 9 busses which elements will this be in the full jacobian matrix? Hiranya Sent from my iPhone
> On 3 Jun 2015, at 02:26, "Abhyankar, Shrirang G." <[email protected]> wrote: > > The Jacobian is stored in a sparse matrix format. The Jacobian values not > displayed are all zeros. > > MATPOWER orders the Jacobian as follows: > > J = [dP_dVa dP_dVm > dQ_dVa dQ_dVm] > > where P and Q are the real and reactive power mismatch equations and Va, > Vm are the voltage angles and magnitudes of the buses. The size of the > Jacobian is 2*nbus X 2*nbus where nbus is the number of buses. > > If you only want the sensitivities w.r.t. voltage magnitudes then extract > the dQ_dVm sub-matrix from the Jacobian. > > > Shri > > -----Original Message----- > From: Hiranya Seneviratne <[email protected]> > Reply-To: MATPOWER discussion forum <[email protected]> > Date: Tuesday, June 2, 2015 at 4:54 PM > To: MATPOWER discussion forum <[email protected]> > Subject: Re: Jacobian matrix values > >> 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.htm >>>>> l >>>>> >>>>> >>>>> -----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 > > >
