【正文】 。)。ylabel(39。bar(Q3)。grid on。支路有功損耗39。xlabel(JDH)。Q3=imag(DS)。subplot(3,2,5)。)。ylabel(39。bar(Q2)。grid on。支路末端注入有功39。title(JDH1)。Q2=imag(Sjy)。subplot(3,2,3)。)。ylabel(39。bar(Q1)。grid on。支路號(hào)39。)。ylabel(39。bar(P1)。endP1=real(Siz)。), 39。,39。(39。)。,num2str(B1(i,2)),39。,num2str(B1(i,1)),39。for i=1:nl JDH=strcat(JDH,num2str(i),39。JDH=[]。%***********************************figure(2)。)。ylabel(39。title(JDQ)。subplot(4,1,4)。)。ylabel(39。bar(P)。P=real(S)。grid on。節(jié)點(diǎn)電壓角度39。title(JDDYJD)。subplot(4,1,2)。)。%xlabel(39。節(jié)點(diǎn)電壓模值39。title(JDDY)。subplot(4,1,1)。 disp(ZF)。 無功總損耗= 39。有功總損耗= 39。,num2str(zws)]。 ZF=[39。節(jié)點(diǎn)無功：39。,JDP)。JDP=strcat(39。節(jié)點(diǎn)電壓模值：39。,JDDYJD)。endJDDYJD=strcat(39。),39。(39。)。,num2str(real(S(i))),39。 JDP=strcat(JDP,num2str(i),39。),39。(39。)。,num2str(sida(i)),39。 JDDYJD=strcat(JDDYJD,num2str(i),39。JDDYJD=[]。JDP=[]。 end zws=0。39。 disp(ZF)。)=39。,39。DS(39。 DS(i)=Si(p,q)+Sj(q,p)。for i=1:nl p=B1(i,1)。各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：39。)。 disp(39。,num2str(Sj(q,p))]。,num2str(p),39。,num2str(q),39。 end ZF=[39。 else 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))))。 if B1(i,6)==0 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))))。for i=1:nl p=B1(i,1)。各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：39。)。 disp(39。,num2str(Si(p,q))]。,num2str(q),39。,num2str(p),39。 end ZF=[39。 else 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))))。 if B1(i,6)==0 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))))。for i=1:nl p=B1(i,1)。各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：39。)。 %顯示各節(jié)點(diǎn)的注入功率disp(39。)。 %計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值enddisp(39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。disp(sida)。各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。 %顯示各節(jié)點(diǎn)的電壓大小V的模值disp(39。)。disp(39。39。disp(E)。各節(jié)點(diǎn)的電壓復(fù)數(shù)值E為(節(jié)點(diǎn)號(hào)從小到大排列)：39。 %計(jì)算各節(jié)點(diǎn)電壓的角度 E(k)=e(k)+f(k)*j。for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。)。disp(39。)。 %修改節(jié)點(diǎn)電壓虛部 end %===== 結(jié)束一次迭代 ===end%********************** 下面為迭代計(jì)算結(jié)束后的有關(guān)輸出過程 *****************disp(39。 %修改節(jié)點(diǎn)電壓實(shí)部 k1=k+1。 %修改節(jié)點(diǎn)電壓 for k=3:2:N01 L=(k+1)./2。disp(J)。 end end %for m=1:N0 % JJN1(m)=J(m,N1)。 %==================== 按列回代運(yùn)算 for k=N0:1:3 for k1=k1:1:3 J(k1,N1)=J(k1,N1)J(k1,k)*J(k,N1)。 %disp(JZ)。)次39。Jacobi矩陣回代第(39。]。,num2str(a),39。 %當(dāng)前行第k列元素已消為0 end end %JZ=[39。 % 對(duì)角元規(guī)格化K行K列對(duì)角元素賦1 %================== 按列消去運(yùn)算 for k2=k+1:N0 % 從k+1行到2*n最后一行 for k3=k+1:N1 %從k2+1列到擴(kuò)展列消去k+1行后各行下三角元素 J(k2,k3)=J(k2,k3)J(k2,k)*J(k,k3)。 %=============================== 以上為形成完整的Jacobi矩陣 %====下面用高斯消去法對(duì)由Jacobi矩陣形成的修正方程進(jìn)行求解(按列消去、回代) for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點(diǎn)） for k1=k+1:N1 % 從k+1列的Jacobi元素到擴(kuò)展列的△P、△Q 或 △U J(k,k1)=J(k,k1)./J(k,k)。 disp(JZ0)。)次Jacobi矩陣：39。形成的第(39。 J(m,q)=X2。J(p,q)=X6。 J(m,q)=X1。 % PV節(jié)點(diǎn)電壓誤差J(p,N)=DV。q=2*j11。 X6=2*f(i)。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。j1~=isb %非平衡節(jié)點(diǎn)amp。 J(m,q)=X2。J(p,q)=X6。 J(m,q)=X1。 % PV節(jié)點(diǎn)電壓誤差J(p,N)=DV。q=2*j11。X6=0。 % dP/de X2=B(i,j1)*e(i)G(i,j1)*f(i)。j1~=i %非平衡節(jié)點(diǎn)amp。 J(m,q)=X2。J(p,q)=X4。 J(m,q)=X1。%擴(kuò)展列△Q J(p,N)=DQ。q=2*j11。 % dQ/de X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。j1~=isb %非平衡節(jié)點(diǎn)amp。J(m,q)=X2。 % X1=dP/de J(m,N)=DP節(jié)點(diǎn)有功功率差 J(m,N)=DP。 J(m,q)=X1。m=p+1。q=2*j11。 % X2=dp/df X3=dQ/de X4=X1。 % X1=dP/de=dQ/df=X4 X2=B(i,j1)*e(i)G(i,j1)*f(i)。j1~=i %非平衡節(jié)點(diǎn)amp。%Σ(Gij*ejBij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。D(i)=0。 %迭代次數(shù) if ICT2(a)==0。 %不滿足要求的節(jié)點(diǎn)數(shù)加1 end end ICT2(a)=IT2。 %PV節(jié)點(diǎn)電壓模平方差J(p,N1)擴(kuò)展列△U end end %(if i~=isb) 非平衡節(jié)點(diǎn)(PQ或PV節(jié)點(diǎn)) end %(for i=1:n) n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U) %判斷功率偏差量及PV節(jié)點(diǎn)的電壓偏差量是否滿足要求 for k=3:N0 %除去平衡節(jié)點(diǎn)2號(hào)以外的所有節(jié)點(diǎn) DET=abs(J(k,N1))。 %PQ節(jié)點(diǎn)無功功率差J(p,N1)擴(kuò)展列△Q else %PV節(jié)點(diǎn)================== J(m,N1)=P(i)P1。 %電壓模平方 %===求取功率差及PV節(jié)點(diǎn)電壓模平方差 ========= if i~=isb %非平衡節(jié)點(diǎn)(PQ或PV節(jié)點(diǎn)) if B2(i,5)~=3 %非PV節(jié)點(diǎn)(只能是PQ節(jié)點(diǎn)) J(m,N1)=P(i)P1。%節(jié)點(diǎn)功率P計(jì)算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。%Σ(Gij*ejBij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。D(i)=0。m=p+1。a=a+1。a=0。N0=2*n。 ICT1=0。 %i節(jié)點(diǎn)無功補(bǔ)償量(電納值)end%====用牛頓拉夫遜法迭代求解非線性代數(shù)方程（功率方程）P=real(S)。 %PV、平衡節(jié)點(diǎn)及PQ節(jié)點(diǎn)電壓模值 endfor i=1:n %給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)B2(