《電力系統(tǒng)分析》（附錄Matlab程序）公開課

附錄Matlab程序

第四章牛頓?拉夫遜法潮流計(jì)算：

%牛頓-拉夫遜法潮流計(jì)算

yO=O.25j;%接地導(dǎo)納都為此值

z24=O.OI5j;

zl2=0.04+0.25j;

z23=0.()8+0.3j;

zl3=O.l+O.35j;

z35=0.03j;

k=1.05;

%求導(dǎo)納矩陣

Yll=yO+l/zl3+l/zl2;

Y12=-l/zl2;

Y21=Y12;

Y13=-l/zl3;

Y31=Y13;

Y14=0;

Y41=0;

Y15=0;

Y51=0;

Y22=y()+I/ZI2+1/z23+y()+1/z24/kA2;

Y23=-l/z23;

Y32=Y23;

Y24=-l/z24/k;

Y42=Y24;

Y25=0；

Y52=0;

Y33=yO+I/ZI3+1/z23+1Zz35/kA2;

Y35=-l/z35/k;

Y53=Y35;

Y34=0;

Y43=0;

Y44=1/Z24;

Y45=0;

Y54=0;

Y55=l/z35;

YY=[Y!1,Y12,Y13,YI4,Y15;Y21,Y22,Y23,Y24,Y25;Y31,Y32,Y33,Y34,Y35;Y41,Y42,Y43,Y4

4.Y45;Y51,Y52,Y53,Y54,Y55]

%形成Bl矩陣(導(dǎo)納矩隴去掉平衡節(jié)點(diǎn)形成)

%B1=[YI1,Y12,YI3,Y14;Y21,Y22,Y23,Y24;Y31,Y32,Y33,Y34;Y41,Y42,Y43,Y44];

%B1l=imag(Bl)

%形成B2矩陣(導(dǎo)納矩除去掉平衡節(jié)點(diǎn)和PV節(jié)點(diǎn)形成)

%B2=[Y1LY12,Y13;Y21,Y22,Y23;Y31,Y32,Y33];

%B22=imag(B2)

G=real(YY);

B=imag(YY);

%賦電壓初值

U=[1.0,1.0,1.Oj.05,1.05]';

Dlt=[0,0,0,0,0],;

P=[-1.6,-2,-375]';

Q=[-0.8,-L-1.3]*；

Ep=le-4;%迭代精度為0.0001

fork=1:10

forii=l:4

aa=0;

forjj=l:5

aa=aa+U(ii)*U(jj)*(G(ii,jj)*cos(Dlt(ii)-Dlt(jj))+B(ii,jj)*sin(Dlt(ii)-Dlt(jj)));

end

dP(ii)=P(ii)-aa;

end

forii=l:3

bb=0;

forjj=l:5

bb=bb+U(ii)*U(jj)*(G(ii,jj)*sin(Dlt(ii)-Dlt(jj))-B(ii,jj)*cos(Dlt(ii)-Dlt(jj)));

end

dQ(ii)=Q(ii)-bb;

end

%求雅可比矩陣

fori=1:4

forj=1:4

if(j==i)

bb=0;

forjj=l:5

ifg)

bb=bb;

else

bb=bb+U(i)*UGj)*(G(ijj)*sin(Dlt(i)-Dlt(jj))-B(i,jj)*cos(Dlt(i)-Dlt(jj)));

end

end

H(iJ)=bb;

N(ij)=-U(i)*U(j)*G(ij)-P(i);

else

H(ij)=-U(i)*U(j)*(G(ij)*sin(Dlt(i)-Dlt(j))-B(ij)*cos(Dlt(i)-Dlt(j)))；

N(i,j)=-U(i)s)Ud)*(G(i,j)*cos(Dlt(i)-Dk(j))+B(ij)*sin(Dlt(i)-Dlt(j)));

end

J(i,j)=H(i,j);

ifj<=3

J(i,j+4尸N(i,j);

end

end

end

fori=l:3

forj=1:4

ifj==i

K(i,j)=U(i)*U(i)*G(ij)日(i);

L(iJ)=U(i)*U(i)*B(ij)-Q(i);

else

K(iJ)=-N(i,j);

L(ij)=H(i.j);

end

J(i+4j)=K(ij);

ifj<=3

J(i+4j+4)=L(i,j);

end

end

end

%解修正方程

XX=inv(J)*[dRdQ],;

fori=l:4

dDlt(i)=XX(i);

ifi<=3

dU(i)=XX(i+4)*U(i);

U(i)=U(i)-dU(i);

end

end

dP

dQ

Dl(

U

C=max(abs(dP));

CC=max(abs(dQ));

if(C<Ep&&CC<Ep)

break

end

end

%求平衡節(jié)點(diǎn)的輸入功率

Pn=0;

forjj=l:5

YY=IY11,Y12,Y13,Y14,Y15;Y21,Y22,Y23,Y24,Y25;Y31,Y32,Y33,Y34,Y35;Y41,Y42,Y43,Y4

4.Y45;Y5LY52,Y53,Y54,Y55]

%形成Bl矩陣(導(dǎo)納矩隴去掉平衡節(jié)點(diǎn)形成)

B1=[Y11,Y12,Y13,Y14;Y21,Y22,Y23,Y24;Y31,Y32,Y33,Y34;Y41,Y42,Y43,Y44];

Bll=imag(Bl)

%形成B2矩陣(導(dǎo)納矩I笆去抻平衡節(jié)點(diǎn)和PV節(jié)點(diǎn)形成)

B2=[Y11,Yl2,Y13;Y21,Y22,Y23;Y31,Y32,Y33];

B22=imag(B2)

G=reaI(YY);

B=imag(YY);

%賦電壓初值

U=[1.0,l.oJ.oJ.05J.05]';

Dk=[0,0,0,0,0]';

P=[-1.6,-2,-3.7,5T;

,

Q=[-0.8rl,-1.3J；

Ep=le-4%迭代精度為0.0001

fork=1:10()

forii=l:4

aa=0;

forjj=l:5

aa=aa+U(ii)*U(jj)*(G(ii,jj)*cos(Dlt(ii)-Dlt(jj))+B(ii,jj)*sin(Dlt(ii)-Dlt(jj)));

end

dP(ii)=P(ii)-aa;

end

forii=l:4

UP(ii)=dP(ii)/U(ii)；

end

dDlt=(inv(BU))*UP';

forii=l:4

dDlt(ii)=dDlt(ii)/U(ii);%求出相位角修正值

Dlt(ii)=Dlt(ii)-dDlt(ii);

end

forii=l:3

bb=0;

forjj=l:5

bb=bb+U(ii)*U(jj)*(G(iijj)*sin(Dlt(ii)-Dlt(jj))-B(ii,jj)*cos(Dlt(ii)-Dlt(jj)));

end

dQ(ii)=Q(ii)-bb;