【正文】 V =V ? ? ΔV ? (18)(8)返回(2)進(jìn)行迭代，直到各節(jié)點(diǎn)功率誤差及電壓誤差都滿足收斂條件。isb=1。%n=6。B1=[1 2 0+ 0 1 1。 0 + 0 2。4 0。請(qǐng)輸入支路數(shù)：nl=39。%B2=[0 0 115 0 1。f=zeros(1,n)。q=B1(i,2)。q=B1(i,2)。disp(Y)。 %i節(jié)點(diǎn)無(wú)功補(bǔ)償量(電納值)end%==================用牛頓拉夫遜法迭代求解非線性代數(shù)方程（功率方程）=======================P=real(S)。a=a+1。,num2str(a),39。 %求取各個(gè)節(jié)點(diǎn)的功率及功率偏差及PV節(jié)點(diǎn)的電壓偏差 for i=1:n %n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U) p=2*i1。的計(jì)算值 P1=C(i)*e(i)+f(i)*D(i)。 end disp(JZ0)。 %當(dāng)前不滿足要求的節(jié)點(diǎn)數(shù)為零 break %退出迭代運(yùn)算 end %以上為求取各個(gè)節(jié)點(diǎn)的功率及功率偏差及PV節(jié)點(diǎn)的電壓偏差 %================= 求取Jacobi矩陣形成修正方程 =================== for i=2:n %n個(gè)節(jié)點(diǎn)2n行(每節(jié)點(diǎn)兩個(gè)方程P和Q或U) if i~=isb %非平衡節(jié)點(diǎn)(PQ或PV節(jié)點(diǎn)) if B2(i,5)~=3 %下面是針對(duì)PQ節(jié)點(diǎn)來求取Jacobi矩陣的元素 =========== C(i)=0。 % X2=dP/df=dQ/de=X3 X3=X2。q=q+1。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。q=q+1。 % dP/df X5=0。q=q+1。% dP/df X5=2*e(i)。q=q+1。]。Jacobi矩陣第(39。]。%disp(J)。)。)。disp(E)。 %顯示各節(jié)點(diǎn)的電壓大小V的模值disp(39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。)。 else if B1(i,6)==0 Si(p,q)=E(p)*(conj(E(p)*B1(i,4)... +(E(p)*B1(i,5)E(q))*(1./(B1(i,3)*B1(i,5)))))。,num2str(q),39。各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：39。 else Sj(q,p)=E(q)*(conj(E(q)*B1(i,4)... +(E(q)*B1(i,5)E(p))*(1./(B1(i,3)*B1(i,5)))))。,num2str(Sj(q,p))]。for i=1:nl p=B1(i,1)。)=39。JDP=[]。)。,num2str(real(S(i))),39。endJDDYJD=strcat(39。,JDP)。節(jié)點(diǎn)電壓角度39。ylabel(39。ylabel(39。for i=1:nl JDH=strcat(JDH,num2str(i),39。(39。bar(P1)。grid on。subplot(3,2,3)。grid on。subplot(3,2,5)。grid on。四、結(jié)果分析及結(jié)論 B1 = Columns 1 through 5 0 + 0 + 0 + + 0 + + 0 + + 0 0 + 0 Columns 6 through 7 0 0 0 0 0 0 0 0 B2 = 0 0 0 0 + 0 0 + 0 0 + 0 0 + 0 0 0 導(dǎo)納矩陣 Y= Columns 1 through 5 0 0 + 0 0 0 0 + + 0 + 0 + + 0 0 0 + + 0 + 0 + 0 0 0 0 0 + Column 6 0 0 0 0 0 + 0 功率方程第(1)次差值： Columns 1 through 10 0 0 Columns 11 through 12 0 形成的第(1)次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 0 Columns 11 through 13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Jacobi矩陣第(1)次消去運(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 thr