【正文】 修正量 i Δ V (7)修正各節(jié)點(diǎn)電壓幅值 i V (i ) (i 1) (i 1) i i i V =V ? ? Δ V ? (18) (8)返回 (2)進(jìn)行迭代，直到各節(jié)點(diǎn)功率誤差及電壓誤差都滿足收斂條件。首先， 因?yàn)樾拚匠淌降南禂?shù)矩陣就是導(dǎo) 納矩陣的虛部，因此在迭代過程中不必象牛頓法那樣進(jìn)行形成雅可比矩陣的 計(jì)算，這樣不僅是僅減少了運(yùn)算量，而且也大大簡(jiǎn)化了程序。 isb=1。 pr=。 % n=6。 nl=6。 B1=[1 2 0+ 0 1 1。 4 5 + 0 1 0 0。 0 + 0 2。)。4 0。請(qǐng)輸入節(jié)點(diǎn)數(shù)： n=39。請(qǐng)輸入支路數(shù)： nl=39。2 3 2+8i 0 1 0 0。 %B2=[0 0 115 0 1。請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=39。f=zeros(1,n)。 % % %求導(dǎo)納矩陣 %for i=1:n % if X(i,2)~=0。q=B1(i,2)。 %非對(duì)角元 Y(q,p)=Y(p,q)。q=B1(i,2)。 %對(duì)角元 j側(cè) +線路電納的一半 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./。disp(Y)。 f(i)=imag(B2(i,3))。 %i節(jié)點(diǎn)無功補(bǔ)償量 (電納值 ) end %================== 用牛頓 拉 夫 遜 法 迭 代 求 解 非 線 性 代 數(shù) 方 程 （ 功 率方程）======================= P=real(S)。N0=2*n。a=a+1。) 次消去 運(yùn)算 39。,num2str(a),39。功率方程第 (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。 for j1=1:n %第 i 行共 n列 (n個(gè)節(jié)點(diǎn)間互導(dǎo)納及節(jié)點(diǎn)電壓相乘即電流 ) C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。的計(jì)算值 P1=C(i)*e(i)+f(i)*D(i)。 %PQ節(jié)點(diǎn)有功功率差 J(m,N1)擴(kuò)展列△ P J(p,N1)=Q(i)Q1。 end disp(JZ0)。 %PQ 節(jié)點(diǎn)的功率偏差量及 PV 節(jié)點(diǎn)的電壓偏差量是否滿足要求 IT2=IT2+1。 %當(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。%Σ(Gij*fj+Bij*ej) end for j1=2:n %第 i行共 n列 (2n個(gè) Jacobi矩陣元素 dP/de及 dP/df或 dQ/de及 dQ/df) if j1~=isbamp。 % X2=dP/df=dQ/de=X3 X3=X2。 J(p,q)=X3。q=q+1。 % X4=dQ/df X2=dp/df elseif j1==iamp。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(p,q)=X3。q=q+1。 end end else %if B2(i,5)~=3 % 否則（即為 PV節(jié)點(diǎn)） %=============== 下面是針對(duì) PV節(jié)點(diǎn)來求取 Jacobi矩陣的元素 =========== for j1=1:n if j1~=isbamp。 % dP/df X5=0。J(p,q)=X5。q=q+1。 elseif j1==iamp。% dP/df X5=2*e(i)。J(p,q)=X5。q=q+1。 end end end %(if B2(i,5)~=3 else) end %(if i~=isb) end %(for i=1:n)n 個(gè)節(jié)點(diǎn) 2n行 (每節(jié)點(diǎn)兩個(gè)方程 P和 Q或 U) JZ0=[39。]。% 用 K行 K列對(duì)角元素去除 K行 K列后的非對(duì)角元素進(jìn)行規(guī)格化 end J(k,k)=1。Jacobi 矩陣第 (39。JZ1=[39。]。 J(k1,k)=0。%disp(J)。 f(L)=f(L)J(k1,N1)。)。)。)。 %計(jì)算各節(jié)點(diǎn)電壓的角度 E(k)=e(k)+f(k)*j。 disp(E)。 disp(39。 %顯示各節(jié)點(diǎn)的電壓大小 V的模值 disp(39。各節(jié)點(diǎn)的電壓相角 sida為 (節(jié)點(diǎn)號(hào)從小到大排列 )： 39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。)。)。 for i=1:nl p=B1(i,1)。 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)))))。 end end ZF=[39。,num2str(q),39。 disp(39。各條支路的末端功率 Sj為 (順序同您輸入 B1時(shí)一致 )： 39。 if B1(i,7)==0 Sj(q,p)=E(q)*conj(E(q)*B1(i,4)./2+(E(q)E(p))./B1(i,3))。 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(q),39。,num2str(Sj(q,p))]。)。 for i=1:nl XVII p=B1(i,1)。DS(39。)=39。39。JDP=[]。 JDDYJD=strcat(JDDYJD,num2str(i),39。)。),39。,num2str(real(S(i))),39。(39。 end JDDYJD=strcat(39。節(jié)點(diǎn)電壓模值： 39。,JDP)。 ZF subplot(4,1,2)。節(jié)點(diǎn)電壓角 度 39。 P=real(S)。 ylabel(39。 subplot(4,1,4)。ylabel(39。%*********************************** figure(2)。 for i=1:nl JDH=strcat(JDH,num2str(i),39。,num2str(B1(i,2)),39。(39。), 39。 bar(P1)。)。 grid on。 ylabel(39。 subplot(3,2,3)。 title(JDH1)。 grid on。 ylabel(39。 subplot(3,2,5)。 xlabel(JDH)。 grid on。ylabel(39。 四、 結(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