MATLAB極坐標(biāo)求解潮流計(jì)算

1、電力系統(tǒng)分析大作業(yè)潮流計(jì)算班級(jí)：XXX姓名：XXX學(xué)號(hào)：XXX例題在圖所示的簡單電力系統(tǒng)中，網(wǎng)絡(luò)各元件參數(shù)的標(biāo)% %= =o.ioo.io十八.4.4。= =J7uoJ7uo= =yO.O1yO.O1528528q q學(xué)= =_/O-3_/O-39 9k k= =1.11.1= =O-O-1 12 2H-H-JO.SQJO.SQ，i4o3410i4o3410= =_ _/ /.0192.0192z z2424= =KO8KO8+ +yO4yO4= =JO.014.3JO.014.3系統(tǒng)中節(jié)點(diǎn) 1、 2 為 P。 節(jié)點(diǎn)， 節(jié)點(diǎn) 3 為叫節(jié)點(diǎn)，節(jié)點(diǎn) 4 為平衡節(jié)點(diǎn)，已給定P Pu uI I雞，=V

2、.3O=V.3O一用.18.18P P”必=。55-55-刊1313心=0.5=0.5,%,%=1.10=1.10乙=1.05X0*=1.05X0*容許誤差二1010一；試用牛頓法。rPQrPQ分解法計(jì)算潮流分布程序清單：帆角坐標(biāo)法求解潮流計(jì)算clear;close;clc;Branch=120.1+0.4j0.01528j10310.3jinf1.111 40.12+0.50j0.01920j10240.08+0.40j0.01413j10%Branch矩陣：1、支路首端號(hào)；2、支路末端號(hào)；3、支路阻抗；4、支路對(duì)地導(dǎo)納；5、支路的變化；6、支路首端處于K側(cè)為1,1側(cè)為0Y=zeros(4)

3、;%節(jié)點(diǎn)導(dǎo)納矩陣fori=1:4ifBranch(i,6)=0%不含變壓器的支路j=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-1/Branch(i,3);Y(k,j)=Y(j,k);Y(j,j)=Y(j,j)+1/Branch(i,3)+Branch(i,4);Y(k,k)=Y(k,k)+1/Branch(i,3)+Branch(i,4);elsej=Branch(i,1);k=Branch(i,2);Y(j,k)=Y(j,k)-Branch(i,5)/Branch(i,3);Y(k,j)=Y(j,k);Y(j，j尸Y(j，j)+1/Branch(i,3)

4、;Y(k,k)=Y(k,k)+Branch(i,5)A2/Branch(i,3);endenddisp(節(jié)點(diǎn)導(dǎo)納矩陣,);YG=real(Y);B=imag(Y);V=1;1;1.1;1.05;%給定V的初始計(jì)算值e=real(V);f=imag(V);disp(節(jié)點(diǎn)電壓的實(shí)部：,)edisp(節(jié)點(diǎn)電壓的虛部：,)fdisp(節(jié)點(diǎn)注入有功功率：,)Ps=-0.3;-0.55;0.5;0disp(節(jié)點(diǎn)注入無功功率：,)Qs=-0.18;-0.13;0;0岫各節(jié)點(diǎn)電壓向量(狀態(tài)變量)可得各節(jié)點(diǎn)注入功率:P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G

5、*f+B*e);del_W=-0.30-P(1);-0.18-Q(1);-0.55-P(2);-0.13-Q(2);0.5用(3);1.1八2-(3)人2片(3)人2;n=0;%輔助循環(huán)計(jì)數(shù)變量while(any(del_W1e-5)&(n5)n=n+1;disp(迭代次數(shù)：);disp(n)%修正方程式:%-0.30-P(1)%-0.18-Q(1)%-0.55-P(2)%-0.13-Q(2)%0.5-P(3)del_e(1)del_f(1)del_e(2)+J*del_f(2)=0del_e(3)%1.1A2-e(3)A2-f(3)A2del_f(3)%雅克比矩陣為：J=P1_e1P

6、1_f1P1_e2P1_f2P1_e3P1_e3;%Q1_e1Q1_f1Q1_e2Q1_f2Q1_e3Q1_f3%P2_e1P2_f1P2_e2P2_f2P2_e3P2_f3%Q2_e1Q2_f1Q2_e2Q2_f2Q2_e3Q2_f3%P3_e1P3_f1P3_e2P3_f2P3_e3P3_f3%V3_e1V3_f1V3_e2V3_f2V3_e3V3_f3%求雅可比矩陣的參數(shù)GeBf=G*e-B*f;%輔助計(jì)算函數(shù)GfBe=G*f+B*e;%輔助計(jì)算函數(shù)fori=1:2forj=1:3ifi=jP_e(i,j)=-GeBf(i)-G(i,j)*e(i)-B(i,j)*f(i);P_f(i,j

7、)=-GfBe(i)-G(i,j)*f(i)+B(i,j)*e(i);Q_e(i,j)=GfBe(i)+B(i,j)*e(i)-G(i,j)*f(i);Q_f(i,j)=-GeBf(i)+G(i,j)*e(i)+B(i,j)*f(i);elseP_e(i,j)=-G(i,j)*e(i)-B(i,j)*f(i);Q_f(i,j尸-P_e(i,j);P_f(i,j)=B(i,j)*e(i)-G(i,j)*f(i);Q_e(i,j)=P_f(i,j);endendendforj=1:2P_e(3,j)=-G(3,j)*e(3)-B(3,j)*f(3);P_f(3,j)=B(3,j)*e(3)-G(3

8、,j)*f(3);V_e(3,j)=0;V_f(3,j)=0;endP_e(3,3)=-GeBf(3)-G(3,3)*e(3)-B(3,3)*f(3);P_f(3,3)=-GfBe(3)-G(3,3)*f(3)+B(3,3)*e(3);V_e(3,3)=-2*e(3);V_f(3,3)=-2*f(3);%求出雅克比矩陣%J=zeros(6,6);fori=1:3forj=1:3J(2*i-1,2*j-1)=P_e(i,j);J(2*i-1,2*j)=P_f(i,j);endendfori=1:2forj=1:3J(2*i,2*j-1)=Q_e(i,j);J(2*i,2*j)=Q_f(i,j);

9、endendforj=1:3J(6,2*j-1)=V_e(3,j);J(6,2*j)=V_f(3,j);enddisp(,雅可比矩陣為：,)J%解修正方程并得出修正后的電壓向量%del_V=-Jdel_W;cleari;fork=1:3V(k)=V(k)+del_V(2*k-1)+i*del_V(2*k);enddisp(節(jié)點(diǎn)電壓為：,)disp(V)%計(jì)算節(jié)點(diǎn)不平衡量%e=real(V);f=imag(V);P=e.*(G*e-B*f)+f.*(G*f+B*e);Q=f.*(G*e-B*f)-e.*(G*f+B*e);del_W=-0.30-P(1);-0.18-Q(1);-0.55-P(2

10、);-0.13-Q(2);0.5-P(3);1.1A2-e(3)A2-f(3)A2;disp(節(jié)點(diǎn)不平衡量為：)disp(del_W)end%最終Z論%disp(最終各節(jié)點(diǎn)電壓幅值為：,)disp(abs(V)disp(最終各節(jié)點(diǎn)電壓相角(度)為：)disp(180*angle(V)/pi)disp(最終各節(jié)點(diǎn)注入功率：,)S=P+i*Qdisp(最終各節(jié)點(diǎn)注入有功功率為：,)Pdisp(最終各節(jié)點(diǎn)注入無功功率為：,)Q運(yùn)行結(jié)果:Branch=1.00002.00000.1000+0.4000i0+0.0153i1.000003.00001.00000+0.3000iInf1.10001.00

11、001.00004.00000.1200+0.5000i0+0.0192i1.000002.00004.00000.0800+0.4000i0+0.0141i1.00000節(jié)點(diǎn)導(dǎo)納矩陣Y=1.0421-8.2429i-0.5882+2.3529i0+3.6667i-0.4539+1.8911i-0.5882+2.3529i1.0690-4.7274i0-0.4808+2.4038i0+3.6667i00-3.3333i0-0.4539+1.8911i-0.4808+2.4038i0節(jié)點(diǎn)電壓的實(shí)部：e=1.00001.00001.10001.0500節(jié)點(diǎn)電壓的虛部：f=0000節(jié)點(diǎn)注入有功功率:

12、Ps=-0.3000-0.55000.500000.9346-4.2616i節(jié)點(diǎn)注入無功功率:Qs=-0.1800-0.130000迭代次數(shù)：1雅可比矩陣為：-1.0194-8.37190.58822.352903.6667-8.11381.06482.3529-0.58823.666700.58822.3529-1.0450-4.8770002.3529-0.5882-4.57781.09300004.0333000-3.66670000-2.20000節(jié)點(diǎn)電壓為：0.9935-0.0088i0.9763-0.1078i1.1000+0.1267i1.0500節(jié)點(diǎn)不平衡量為：-0.0013-

13、0.0028-0.0135-0.05470.0030-0.0160迭代次數(shù)：2雅可比矩陣為：-0.8091-8.36130.60522.33250.03243.6429-7.99921.40712.3325-0.60523.6429-0.03240.82802.2338-1.0190-4.6364002.2338-0.8280-4.36412.087800-0.46444.033300-0.0324-3.64290000-2.2000-0.2533節(jié)點(diǎn)電壓為：0.9847-0.0086i0.9590-0.1084i1.0924+0.1289i1.0500節(jié)點(diǎn)不平衡量為：-0.0000-0.0000-0.0003-0.00110.