【正文】 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Jacobi矩陣第(3)次消去運(yùn)算 Columns 1 through 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 11 through 13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Jacobi矩陣第(3)次回代運(yùn)算 Columns 1 through 10 0 0 Columns 11 through 12 各個(gè)節(jié)點(diǎn)電壓模 0 功率方程第(4)次差值： Columns 1 through 10 0 0 Columns 11 through 12 形成的第(4)次Jacobi矩陣： Columns 1 through 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 11 through 13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Jacobi矩陣第(4)次消去運(yùn)算 Columns 1 through 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 。)。bar(Q3)。支路有功損耗39。Q3=imag(DS)。)。bar(Q2)。支路末端注入有功39。Q2=imag(Sjy)。)。bar(Q1)。支路號(hào)39。ylabel(39。endP1=real(Siz)。,39。)。,num2str(B1(i,1)),39。JDH=[]。)。title(JDQ)。)。bar(P)。grid on。title(JDDYJD)。節(jié)點(diǎn)無功：39。JDP=strcat(39。,JDDYJD)。),39。)。 JDP=strcat(JDP,num2str(i),39。(39。,num2str(sida(i)),39。JDDYJD=[]。 end zws=0。 disp(ZF)。,39。 DS(i)=Si(p,q)+Sj(q,p)。各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：39。 disp(39。,num2str(p),39。 end end ZF=[39。 else if B1(i,6)==0 Sj(q,p)=E(q)*conj((E(q)E(p)*B1(i,5))*(1./(B1(i,3)*B1(i,5)^2)))。for i=1:nl p=B1(i,1)。)。,num2str(Si(p,q))]。,num2str(p),39。 else Si(p,q)=E(p)*conj((E(p)E(q)*B1(i,5))*(1./(B1(i,3)*B1(i,5)^2)))。 if B1(i,7)==0 Si(p,q)=E(p)*conj(E(p)*B1(i,4)./2+(E(p)E(q))./B1(i,3))。各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：39。 %顯示各節(jié)點(diǎn)的注入功率disp(39。 %計(jì)算各節(jié)點(diǎn)的功率 S = 電壓 X 注入電流的共軛值enddisp(39。disp(sida)。)。)。39。各節(jié)點(diǎn)的電壓復(fù)數(shù)值E為(節(jié)點(diǎn)號(hào)從小到大排列)：39。for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。disp(39。 %============================== 結(jié)束一次迭代 ==============================end%********************** 下面為迭代計(jì)算結(jié)束后的有關(guān)輸出過程 *****************disp(39。 end disp(39。 e(L)=e(L)J(k,N1)。 end disp(JZ1)。disp(J)。,num2str(a),39。)次消去運(yùn)算39。%消去運(yùn)算 end %用當(dāng)前行K3列元素減去當(dāng)前行K列元素乘以第k行K3列元素 J(k2,k)=0。disp(J)。,num2str(a),39。 % PV節(jié)點(diǎn)有功誤差J(m,N)=DP。 m=p+1。 p=2*i1。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 % PV節(jié)點(diǎn)有功誤差J(m,N)=DP。 m=p+1。 p=2*i1。非對(duì)角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。%擴(kuò)展列△P J(m,N)=DP。 m=p+1。% dQ/df p=2*i1。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 J(p,q)=X4。 % X3=dQ/de J(p,N)=