Scribd
Upload
45 views15 pages

Formation of Jacobian Matrix

The document discusses the formation of the Jacobian matrix for the Newton-Raphson power flow method. It involves computing the diagonal and off-diagonal elements of the J1-J4 submatrices that make up the full Jacobian matrix. The diagonal elements are computed by iterating through each bus and summing the voltage and admittance products. The off-diagonal elements are similarly computed by iterating through bus pairs, with the exception of the swing bus which is held fixed. The submatrices are then assembled into the full Jacobian matrix for the power flow problem.

Uploaded by

hananwar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PPT, PDF, TXT or read online on Scribd
45 views15 pages

Formation of Jacobian Matrix

The document discusses the formation of the Jacobian matrix for the Newton-Raphson power flow method. It involves computing the diagonal and off-diagonal elements of the J1-J4 submatrices that make up the full Jacobian matrix. The diagonal elements are computed by iterating through each bus and summing the voltage and admittance products. The off-diagonal elements are similarly computed by iterating through bus pairs, with the exception of the swing bus which is held fixed. The submatrices are then assembled into the full Jacobian matrix for the power flow problem.

Uploaded by

hananwar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PPT, PDF, TXT or read online on Scribd
You are on page 1/ 15

Formation of Jacobian matrix

Newton- Raphson power flow solution

For large system ,more efficient and practical .


No iteration is independent of system size
For typical bus of the power system current
entering bus i is given by

clc
nbus=3; sb=1;
v=[1.0 1.0 1.0];delta=[0 0 0]
yb=[65-j*45 -40+j*20 -25+j*25
-40+j*20 60-j*60 -20+j*40
-25+j*25 -20+j*40 45-j*65]

clear J
% computation of diagonal elements of J1 to J4
for i=1:nbus
J1(i ,i)=0.0;J3 (i, i)=0.0;J22=0;J44=0;
for j=1:nbus
if j~=i
J1(i, i)=J1(i ,i)+v(i)*v(j)*ym( i ,j)*sin(angle( i,j)-delta(i)+delta(j));
J3(i ,i)=J3(i ,i)+v(i)*v(j)*ym( i,j)*cos(angle( i ,j)-delta(i)+delta(j));
J22=J22+v(j)*ym(i,j)*cos(angle(i,j)-delta(i)+delta(j));
J44=J44+v(j)*ym(i,j)*sin(angle(i,j)-delta(i)+delta(j));
end
end
if i~=sb
J1(i, i)=J1(i ,i);
J3(i, i)=J3(i ,i);
J2(i,i)=2*v(i)*ym(i,i)*cos(angle(i,i))+J22;
J4(i,i)=-2*v(i)*ym(i,i)*sin(angle(i,i))-J44;
end
end

% computation of OFF diagonal elements of J1 to J4


for i=1:nbus

for j=1:nbus

if (j~=i) & (i~=sb) & (j~=sb)

J1(i, j)=-v(i)*v(j)*ym(i, j)*sin(angle(i, j)-delta(i)+delta(j));

J2(i ,j)=v(i)*ym(i ,j )*cos(angle(i, j )-delta(i)+delta(j));

J3(i , j)=-v(i)*v(j)*ym(i, j)*cos(angle(i, j)-delta(i)+delta(j));

J4(i,j)=-v(i)*ym(i,j)*sin(angle(i,j)-delta(i)+delta(j));

end

end
end

%arrangement of jacobian matrix


for i=1:nbus
for j=1:nbus
if i~=sb & j~=sb
J(i-1,j-1)=J1(i, j);
J(i-1,j+nbus-2)=J2(i, j);
J(i+nbus-2,j-1)=J3(i , j);
J(i+nbus-2,j+nbus-2)=J4(i ,j);
end
end
end
J

function [ym,angle]=matrixrec2polar(a)
a=[65-j*45 -40+j*20 -25+j*25
-40+j*20 60-j*60 -20+j*40
-25+j*25 -20+j*40 45-j*65]
ym=sqrt(real(a).^2+imag(a).^2)
angle=atan2(imag(a),real(a))

yb =

65.0000 -45.0000i -40.0000 +20.0000i -25.0000 +25.0000i


-40.0000 +20.0000i 60.0000 -60.0000i -20.0000 +40.0000i
-25.0000 +25.0000i -20.0000 +40.0000i 45.0000 -65.0000i

ym =

79.0569 44.7214 35.3553


44.7214 84.8528 44.7214
35.3553 44.7214 79.0569

angle =

-0.6055 2.6779 2.3562


2.6779 -0.7854 2.0344
2.3562 2.0344 -0.9653

J1 =

J2 =

45.0000
0
0
0
60.0000 -40.0000
0
-40.0000 65.0000

0
0
0

0
0
60.0000 -20.0000
-20.0000 45.0000

J4 =

J3 =

0
0
0

-65 0 0
0 -60 20
0 20 -45

J=
60.0000
-40.0000
-60.0000
20.0000

-40.0000
65.0000
20.0000
-45.0000

60.0000
-20.0000
60.0000
-40.0000

-20.0000
45.0000
-40.0000
65.0000

0
0
60.0000 -40.0000
-40.0000 65.0000

You might also like