【正文】 整的Jacobi矩陣 %====下面用高斯消去法對由Jacobi矩陣形成的修正方程進行求解(按列消去、回代) for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） for k1=k+1:N1 % 從k+1列的Jacobi元素到擴展列的△P、△Q 或 △U J(k,k1)=J(k,k1)./J(k,k)。 end J(k,k)=1。 % 對角元規(guī)格化K行K列對角元素賦1 %================== 按列消去運算 for k2=k+1:N0 % 從k+1行到2*n最后一行 for k3=k+1:N1 %從k2+1列到擴展列消去k+1行后各行下三角元素 J(k2,k3)=J(k2,k3)J(k2,k)*J(k,k3)。%消去運算 end %用當前行K3列元素減去當前行K列元素乘以第k行K3列元素 J(k2,k)=0。 %當前行第k列元素已消為0 end end %JZ=[39。Jacobi矩陣前代第(39。,num2str(a),39。)次39。]。JZ1=[39。Jacobi矩陣回代第(39。,num2str(a),39。)次39。]。 %disp(JZ)。disp(J)。 %==================== 按列回代運算 for k=N0:1:3 for k1=k1:1:3 J(k1,N1)=J(k1,N1)J(k1,k)*J(k,N1)。 J(k1,k)=0。 end end %for m=1:N0 % JJN1(m)=J(m,N1)。 %end %disp(JZ1)。disp(J)。disp(JJN1)。 %修改節(jié)點電壓 for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N1)。 %修改節(jié)點電壓實部 k1=k+1。 f(L)=f(L)J(k1,N1)。 %修改節(jié)點電壓虛部 end %===== 結束一次迭代 ===end%********************** 下面為迭代計算結束后的有關輸出過程 *****************disp(39。迭代次數(shù)：39。)。disp(ICT11)。disp(39。沒有達到精度要求的個數(shù)：39。)。disp(ICT2)。for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 %計算各節(jié)點電壓的模值 sida(k)=atan(f(k)./e(k))*180./pi。 %計算各節(jié)點電壓的角度 E(k)=e(k)+f(k)*j。 %將各節(jié)點電壓用復數(shù)表示end%=============== 計算各輸出量 ========disp(39。各節(jié)點的電壓復數(shù)值E為(節(jié)點號從小到大排列)：39。)。disp(E)。 %顯示各節(jié)點的實際電壓值E用復數(shù)表示disp(39。39。)。disp(39。各節(jié)點的電壓模值大小V為(節(jié)點號從小到大排列)：39。)。disp(V)。 %顯示各節(jié)點的電壓大小V的模值disp(39。39。)。disp(39。各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)：39。)。disp(sida)。 %顯示各節(jié)點的電壓相角for p=1:n C(p)=0。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 %計算各節(jié)點注入電流的共軛值 end S(p)=E(p)*C(p)。 %計算各節(jié)點的功率 S = 電壓 X 注入電流的共軛值enddisp(39。各節(jié)點的功率S為(節(jié)點號從小到大排列)：39。)。disp(S)。 %顯示各節(jié)點的注入功率disp(39。39。)。disp(39。各條支路的首端功率Si為(順序同您輸入B1時一致)：39。)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 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))))。 Siz(i)=Si(p,q)。 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))))。 Siz(i)=Si(p,q)。 end ZF=[39。S(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(Si(p,q))]。 disp(ZF)。 disp(39。39。)。enddisp(39。各條支路的末端功率Sj為(順序同您輸入B1時一致)：39。)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 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))))。 Sjy(i)=Sj(q,p)。 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))))。 Sjy(i)=Sj(q,p)。 end ZF=[39。S(39。,num2str(q),39。,39。,num2str(p),39。)=39。,num2str(Sj(q,p))]。 disp(ZF)。 disp(39。39。)。enddisp(39。各條支路的功率損耗DS為(順序同您輸入B1時一致)：39。)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 DS(i)=Si(p,q)+Sj(q,p)。 ZF=[39。DS(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(DS(i))]。 disp(ZF)。 disp(39。39。)。 end zws=0。JDDY=[]。JDP=[]。JDQ=[]。JDDYJD=[]。 for i=1:n %總網(wǎng)損為所有節(jié)點注入功率的代數(shù)和 zws=zws+S(i)。 JDDYJD=strcat(JDDYJD,num2str(i),39。(39。,num2str(sida(i)),39。),39。)。 JDDY=strcat(JDDY,num2str(i),39。(39。,num2str(V(i)),39。),39。)。 JDP=strcat(JDP,num2str(i),39。(39。,num2str(real(S(i))),39。),39。)。 JDQ=strcat(JDQ,num2str(i),39。(39。,num2str(imag(S(i))),39。),39。)。endJDDYJD=strcat(39。節(jié)點電壓角度：39。,JDDYJD)。JDDY=strcat(39。節(jié)點電壓模值：39。,JDDY)。JDP=strcat(39。節(jié)點有功：39。,JDP)。JDQ=strcat(39。節(jié)點無功：39。,JDQ)。 ZF=[39?？偩W(wǎng)損S=39。,num2str(zws)]。 ZFS=[39。有功總損耗= 39。,num2str(real(zws)),39。 無功總損耗= 39。,num2str(imag(zws))]。 disp(ZF)。figure(1)。subplot(4,1,1)。bar(V)。title(JDDY)。ylabel(39。節(jié)點電壓模值39。)。%xlabel(39。節(jié)點號39。)。grid on。subplot(4,1,2)。plot(sida)。title(JDDYJD)。ylabel(39。節(jié)點電壓角度39。)。grid on。subplot(4,1,3)。P=real(S)。Q=imag(S)。bar(P)。title(JDP)。ylabel(39。節(jié)點注入有功39。)。grid on。subplot(4,1,4)。bar(Q)。title(JDQ)。xlabel(ZFS)。ylabel(39。節(jié)點注入無功39。)。grid on。%***********************************figure(2)。subplot(3,2,1)。JDH=[]。JDH1=[]。for i=1:nl JDH=strcat(JDH,num2str(i),39。(39。,num2str(B1(i,1)),39。,39。,num2str(B1(i,2)),39。), 39。)。 JDH1=strcat(JDH1,num2str(i),39。(39。,num2str(B1(i,2)),39。,39。,num2str(B1(i,1)),39。), 39。)。endP1=real(Siz)。Q1=imag(Siz)。bar(P1)。title(JDH)。ylabel(39。支路首端注入有功39。)。%xlabel(39。支路號39。)。grid on。subplot(3,2,2)。bar(Q1)。title(JDH)。ylabel(39。支路首端注入無功39。)。grid on。subplot(3,2,3)。P2=real(Sjy)。Q2=imag(Sjy)。bar(P2)。title(JDH1)。ylabel(39。支路末端注入有功39。)。grid on。subplot(3,2,4)。bar(Q2)。title(JDH1)。ylabel(39。支路末端注入無功39。)。grid on。subplot(3,2,5)。P3=real(DS)。Q3=imag(DS)。bar(P3)。xlabel(JDH)。ylabel(39。支路有功損耗39。)。grid on。subplot(3,2,6)。bar(Q3)。xlabel(JDH)。ylabel(39。支路無功損耗39。)。grid on。