Ok, to make it a bit easier, since it seems the Jacobian is something people
want access to, I extracted the relevant code into a function makeJac that you
can now call directly using a MATPOWER case struct with consecutive bus numbers
(internal ordering). So you should be able to do something like ...
>> r = runpf('case4gs');
>> J = makeJac(r);
>> full(J)
ans =
41.1962 -26.1025 -15.0937 -4.1183 -2.1655
-25.6508 44.3270 0 7.0968 0
-14.7397 0 39.6093 0 5.8753
6.3048 -10.3721 0 42.9754 0
3.8683 0 -9.6932 0 38.3182
I'll include this function in MATPOWER 4.
--
Ray Zimmerman
Senior Research Associate
211 Warren Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645
makeJac.m
Description: makeJac.m
On Jun 24, 2010, at 7:40 AM, RAVIKUMAR V wrote: Hai , This is the answer for displaying the jacobian matrix from the load flow, as given by Dr Ray Zimmerman, from the previous post. You can try this. ---------------------------------------------------------
The easiest way is probably to set a breakpoint in the Matlab debugger
at line 126 in newtonpf.m, as both Ybus and J (the Jacobian)
are sitting right there and can be viewed or saved or whatever. For
example,
>> dbstop 126 newtonpf
>> runpf('case4gs')
126 converged = 1;
K>> full(Ybus)
ans =
8.9852 -44.8360i -3.8156 +19.0781i -5.1696 +25.8478i 0
-3.8156 +19.0781i 8.9852 -44.8360i 0 -5.1696
+25.8478i
-5.1696 +25.8478i 0 8.1933 -40.8638i -3.0237
+15.1185i
0 -5.1696 +25.8478i -3.0237 +15.1185i 8.1933
-40.8638i
K>> full(J)
ans =
41.1963 -26.1026 -15.0938 -4.1183 -2.1655
-25.6508 44.3270 0 7.0969 0
-14.7398 0 39.6096 0 5.8755
6.3048 -10.3721 0 42.9755 0
3.8683 0 -9.6932 0 38.3186
K>> save mymatrices Ybus J
K>> dbcont
Newton's method power flow converged in 3 iterations.-------------------------------------------------- On Thu, Jun 24, 2010 at 3:05 AM, indra g.c. <[email protected]> wrote:
Dear all,
I am new in the list and also new user of MATPOWER .
I require to print the base case Jacobian Matrix in MATPOWER
Can anybody help me how can I do ?
Best regards,
Indra Bahadur GC.
-- ------------------------------------------------------------- WITH REGARDS, V.RAVIKUMAR PANDI, -------------------------------------------------------------
HAVE A NICE DAY -------------------------------------------------------------
|
