【正文】 f(L)=f(L)J(k1,N1)。 e(L)=e(L)J(k,N1)。disp(JJN1)。 %end %disp(JZ1)。 J(k1,k)=0。disp(J)。]。,num2str(a),39。JZ1=[39。)次39。Jacobi矩陣前代第(39。%消去運算 end %用當(dāng)前行K3列元素減去當(dāng)前行K列元素乘以第k行K3列元素 J(k2,k)=0。% 用K行K列對角元素去除K行K列后的非對角元素進(jìn)行規(guī)格化 end J(k,k)=1。disp(J)。]。,num2str(a),39。 end end end %(if B2(i,5)~=3 else) end %(if i~=isb) end %(for i=1:n)n個節(jié)點2n行(每節(jié)點兩個方程P和Q或U) JZ0=[39。 % PV節(jié)點有功誤差J(m,N)=DP。q=q+1。 m=p+1。J(p,q)=X5。 p=2*i1。% dP/df X5=2*e(i)。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 elseif j1==iamp。 % PV節(jié)點有功誤差J(m,N)=DP。q=q+1。 m=p+1。J(p,q)=X5。 p=2*i1。 % dP/df X5=0。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 end end else %if B2(i,5)~=3 否則（即為PV節(jié)點） %=============== 下面是針對PV節(jié)點來求取Jacobi矩陣的元素 =========== for j1=1:n if j1~=isbamp。%擴展列△P J(m,N)=DP。q=q+1。 m=p+1。J(p,q)=X3。% dQ/df p=2*i1。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 % X4=dQ/df X2=dp/df elseif j1==iamp。 J(p,q)=X4。q=q+1。 % X3=dQ/de J(p,N)=DQ節(jié)點無功功率差 J(p,N)=DQ。 J(p,q)=X3。 % X1=dP/de X4=dQ/df p=2*i1。 % X2=dP/df=dQ/de=X3 X3=X2。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。%Σ(Gij*fj+Bij*ej) end for j1=2:n %第i行共n列(2n個Jacobi矩陣元素dP/de及dP/df或dQ/de及dQ/df) if j1~=isbamp。 for j1=1:n %第i行共n列(n個節(jié)點間互導(dǎo)納及節(jié)點電壓相乘即電流) C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。 %當(dāng)前不滿足要求的節(jié)點數(shù)為零 break %退出迭代運算 end %以上為求取各個節(jié)點的功率及功率偏差及PV節(jié)點的電壓偏差 %================= 求取Jacobi矩陣形成修正方程 =================== for i=2:n %n個節(jié)點2n行(每節(jié)點兩個方程P和Q或U) if i~=isb %非平衡節(jié)點(PQ或PV節(jié)點) if B2(i,5)~=3 %下面是針對PQ節(jié)點來求取Jacobi矩陣的元素 =========== C(i)=0。 %不滿足要求的節(jié)點數(shù)；a為迭代次數(shù) ICT1=ICT1+1。 if DET=pr %PQ節(jié)點的功率偏差量及PV節(jié)點的電壓偏差量是否滿足要求 IT2=IT2+1。 %PV節(jié)點有功功率差J(m,N1)擴展列△P J(p,N1)=V(i)^2V2。 %PQ節(jié)點有功功率差J(m,N1)擴展列△P J(p,N1)=Q(i)Q1。%節(jié)點功率Q計算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej) V2=e(i)^2+f(i)^2。的計算值 P1=C(i)*e(i)+f(i)*D(i)。%Σ(Gij*fj+Bij*ej) end %求i節(jié)點有功和無功功率P39。 for j1=1:n %第i行共n列(n個節(jié)點間互導(dǎo)納及節(jié)點電壓相乘即電流) C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。C(i)=0。 %求取各個節(jié)點的功率及功率偏差及PV節(jié)點的電壓偏差 for i=1:n %n個節(jié)點2n行(每節(jié)點兩個方程P和Q或U) p=2*i1。 %迭代次數(shù)ICTa；不滿足收斂要求的節(jié)點數(shù)IT2while IT2~=0 % N0=2*n 雅可比矩陣的階數(shù)；N1=N0+1擴展列 IT2=0。N1=N0+1。IT2=1。Q=imag(S)。 %i節(jié)點注入功率SGSL B(i,i)=B(i,i)+B2(i,4)。 V(i)=abs(B2(i,3))。 %分解出導(dǎo)納陣的實部和虛部 for i=1:n %給定各節(jié)點初始電壓的實部和虛部 e(i)=real(B2(i,3))。disp(Y)%給定各節(jié)點初始電壓及給定各節(jié)點注入功率G=real(Y)。導(dǎo)納矩陣 Y=39。 %對角元j側(cè)+線路電納的一半 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./。 %非對角元 Y(q,p)=Y(p,q)。q=B1(i,2)。 %對角元K側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)。 %非對角元 Y(q,p)=Y(p,q)。q=B1(i,1)。q=B1(i,2)。S1=zeros(nl)。V=zeros(1,n)。e=zeros(1,n)。)。%0 10+6i 0 3] %input(39。%B2=[0 0 115 0 1。請輸入由支路參數(shù)形成的矩陣： B1=39。%1 3 4+16i 0 1 0 0。)。%input(39。)。%input(39。)。%input(39。)。%input(39。)。0 70+ 10 0 2] %input(39。0 40+ 35 0 2。0 50+ 10 0 2。0 50+ 35 0 2。150 0 230 0 3。)。 ] %input(39。 7 9 + 1 0 0。 5 6 + 0 220/11 1 1。 2 9 + 1 0 0。 1 5 + 1 0 0?！禡ATLAB實用指南》，蘇金明，王永利，電子工業(yè)出版社。同時感謝在課設(shè)過程中幫助我的老師們和同學(xué)，讓我在課設(shè)中能夠有很多收獲!七．參考文獻(xiàn)《電力系統(tǒng)分析課程設(shè)計綜合實驗》，祝淑萍 等，中國電力出版社。課設(shè)雖然結(jié)束了，但是我學(xué)習(xí)電力系統(tǒng)的過程仍在繼續(xù)，相信以后還會接觸到更多關(guān)于電力系統(tǒng)的知識，更深入的了解潮流計算以及其他計算分析。在潮流計算的過程中也是坎坷不斷，做了三次潮流計算的更改才形成最終這篇論文，但也是頗有收獲！第一次所有數(shù)值都采用的標(biāo)幺值，后來覺得不合理又重新做了一次；第二次由于先調(diào)的電廠電壓再調(diào)的變比，但是考慮到電廠電壓是不可輕易調(diào)的，要考慮實際情況，所以最后采用先調(diào)變比再調(diào)電廠電壓的方法，符合了要求。潮流計算是實現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)展的必要手段和重要環(huán)節(jié)。課設(shè)的過程在老師的指導(dǎo)下我知道了潮流計算是電力系統(tǒng)的最基本、最常用的計算。所以我們花了兩周的時間進(jìn)行潮流計算的課程設(shè)計，課程設(shè)計的過程我也是受益匪淺。了解了電力系統(tǒng)方面的基礎(chǔ)知識，也知道電力系統(tǒng)在日后的學(xué)習(xí)工作中占據(jù)很重要的地位。具體的潮流分布如下表：各條支路的首端功