【正文】 if ICT3~=0 L=0。 for i=1:na DQ(i)=A(i,i)*DQ(i)。 if i==na for LZ=2:i L=i+2LZ。 IC4=L1。 for MZ=1:IC4 I=IC4+1MZ。 DQ(I)=DQ(I)A(I,L)*DQ(L)。 end end else IC1=i+1。 for k=IC1:na DQ(k)=DQ(k)A(k,i)*DQ(i)。 end end end L=0。 for i=1:n if B2(i,6)==2 L=L+1。 V(i)=V(i)DQ(L)。 end end kp=1。 K=K+1。 else kq=0。 if kp~=0 K=K+1。 end end for i=1:n Dy(K1,i)=V(i)。 end end disp(39。迭代次數(shù) 39。) disp(K)。 disp(39。每次沒有達(dá)到精度要求的有功功率個(gè)數(shù)為 39。)。 disp(Np)。 disp(39。每次沒有達(dá)到精度要求的無(wú)功功率個(gè)數(shù)為 39。)。 disp(Nq)。 for k=1:n E(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j。 O(k)=O(k)*180./pi。 end disp(39。各節(jié)點(diǎn)的電壓標(biāo)么值 E 為 39。)。 disp(E)。 disp(39。各節(jié)點(diǎn)的電壓 V 大小 39。)。 disp(V)。 disp(39。各節(jié)點(diǎn)的電壓相角 O39。)。 disp(O)。 for p=1:n C(p)=0。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 end S(p)=E(p)*C(p)。 end disp(39。各節(jié)點(diǎn)的功率為 39。)。 disp(S)。 disp(39。各條支路的首端功率為 39。)。 for i=1:nl if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。 disp(Si(p,q))。 end disp(39。各條支路的末端功率為 39。)。 for i=1:nl if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 disp(Sj(q,p))。 end disp(39。各條支路的功率損耗為 39。)。 for i=1:nl if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。 else p=B1(i,2)。q=B1(i,1)。 end DS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。 end for i=1:K Cs(i)=i。 for j=1:n Dy(K,j)=Dy(K1,j)。 end end disp(39。每次迭代后各節(jié)點(diǎn)的電壓值如圖所示 39。)。 plot(Cs,Dy) xlabel(39。迭代次數(shù) 39。) ylabel(39。電壓 39。) title(39。電壓迭代次數(shù)曲線 39。)。 題號(hào)： 2 節(jié)點(diǎn)數(shù)： 4 支路數(shù)： 4 計(jì)算精度： 支路 1： + 1┠ ————— □————— ┨ 3 支路 2： + 1┠ ————— □————— ┨ 4 支路 3： + 2┠ ————— □————— ┨ 4 支路 4： + 3┠ ————— □————— ┨ 4 節(jié)點(diǎn) 1： PQ 節(jié)點(diǎn)， S(1)= 節(jié)點(diǎn) 2： PQ 節(jié)點(diǎn)， S(2)= 節(jié)點(diǎn) 3： PV 節(jié)點(diǎn)， P(3)= V(3)= 節(jié)點(diǎn) 4：平衡節(jié)點(diǎn)， U(4)=∠ include include define M 20 /*節(jié)點(diǎn)數(shù)、支路數(shù)極限值 */ define N 5 /*迭代次數(shù)極限值 */ int n,m,dd=0,pq=0,pv=0,tt,qy。 float eps。 /*節(jié)點(diǎn)數(shù)、支路數(shù)、對(duì)地支路數(shù)、 pq 節(jié)點(diǎn)數(shù)、 pv 節(jié)點(diǎn)數(shù)、控制打印、互聯(lián)網(wǎng)數(shù)、精度 */ struct jiedian /*節(jié)點(diǎn)數(shù)據(jù)結(jié)構(gòu)體 */ {int s。 /*節(jié)點(diǎn)類型： 1PQ 節(jié)點(diǎn)； 2PV 節(jié)點(diǎn)； 3平衡節(jié)點(diǎn) */ float p,q,e,f,v。 /*節(jié)點(diǎn)的有功、無(wú)功、電壓實(shí)部、虛部、電壓輻值 */ }jiedian[M]。 /*如引用節(jié)點(diǎn) 1 的無(wú)功，則為 jiedian[1].q */ struct zhilu /*支路數(shù)據(jù)結(jié)構(gòu)體 */ {int p1,p2,s。 /*支路兩端節(jié)點(diǎn)號(hào)，支路類型： 1普通支路； 2變壓器支路； 3對(duì)地支路 */ float r,x,b,kt。 /*支路的電阻、電抗、導(dǎo)納、變壓器非標(biāo)準(zhǔn)變比 */ }zhilu[M]。 struct hulianwang /*互聯(lián)網(wǎng)結(jié)構(gòu)體 */ {int num,pv。 /*互聯(lián)區(qū)域號(hào) ,指定的 pv 節(jié)點(diǎn)號(hào) */ float p,eps2。 /*規(guī)定的有功功率及其允許誤差 */ int count。 /*每個(gè)互聯(lián)區(qū)域包括的節(jié)點(diǎn)數(shù) */ int a[M]。 /*包括的節(jié)點(diǎn) */ } hulianwang[M]。 static float G[M][M],B[M][M],G1[M][M],B1[M][M]。 /*節(jié)點(diǎn)導(dǎo)納陣矩陣 */ static float ykb[2*M][2*M]。 /*節(jié)點(diǎn)導(dǎo)納陣矩陣 */ static float yinzi[2*M][2*M]。 /*因子表 */ static float P[M][M],Q[M][M]。 /*潮流計(jì)算結(jié)果 */ FILE *fp1,*fp2。 /*文件輸入、輸出指針 */ void input() /*從文件 輸入線路基本參數(shù)、節(jié)點(diǎn)數(shù)據(jù)、支路數(shù)據(jù) */ { int i,j,h。 /*循環(huán)變量、節(jié)點(diǎn)（支路）類型 */ /**********************************打開電網(wǎng)參數(shù)文件 **********************/ fp1=fopen(,r)。 if (fp1==NULL) { printf(Can not open file !\n)。