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
-------------------------------------------------------------
