【正文】 K側(cè)標識為1時歸算至末端側(cè)，0歸算至首端側(cè)B2矩陣：該節(jié)點發(fā)電機功率； 該節(jié)點負荷功率；PQ節(jié)點電壓初始值，或平衡節(jié)點及PV節(jié)點電壓的給定值節(jié)點所接無功補償并聯(lián)電容（感）的電納節(jié)點分類標號：1為平衡節(jié)點（應(yīng)為1號節(jié)點）；2為PQ節(jié)點；3為PV節(jié)點；B1=[ 1 3 + 1 01 5 + 1 01 7 + 1 02 9 + 1 0 3 4 + 0 220/ 15 6 + 0 220/11 17 8 + 0 220/ 17 9 + 1 09 10 + 0 220/11 1]B2=[ 0 0 231 0 1200 0 231 0 30 0 220 0 20 50+ 35 0 20 0 220 0 20 50+ 10 0 20 0 220 0 20 40+ 35 0 20 0 220 0 20 70+ 10 0 2]由于各節(jié)點電壓不應(yīng)超過各自的允許范圍，由題目要求知，變壓器低壓側(cè)電壓有名值變化范圍如下表：節(jié)點46108 電壓下限（KV）3535電壓上限(KV)3636同時保證發(fā)電廠發(fā)出的有功功率要大于變電所有功功率，并且兩者差值越小表明線路有功功率損耗越小，越滿足實際要求。第二次調(diào)變比的具體的潮流分布如下：各條支路的首端功率Si各條支路的末端功率SjS(1,3)= +S(3,1)= S(1,5)= +S(5,1)= S(1,7)= + S(7,1)= S(2,9)= S(9,2)= +S(3,4)= + S(4,3)= S(5,6)= + S(6,5)= S(7,8)= +S(8,7)= S(7,9)= + S(9,7)= S(9,10)= +S(10,9)= 各節(jié)點的功率S為(節(jié)點號從小到大排列)：Columns 1 through 5 + + + Columns 6 through 10 + 得到的圖像如下：由于調(diào)節(jié)變壓器分接頭達到極限仍未達到題目要求的范圍，所以接下來調(diào)節(jié)發(fā)電機機端電壓。了解了電力系統(tǒng)方面的基礎(chǔ)知識，也知道電力系統(tǒng)在日后的學(xué)習(xí)工作中占據(jù)很重要的地位。課設(shè)的過程在老師的指導(dǎo)下我知道了潮流計算是電力系統(tǒng)的最基本、最常用的計算。在潮流計算的過程中也是坎坷不斷，做了三次潮流計算的更改才形成最終這篇論文，但也是頗有收獲！第一次所有數(shù)值都采用的標幺值，后來覺得不合理又重新做了一次；第二次由于先調(diào)的電廠電壓再調(diào)的變比，但是考慮到電廠電壓是不可輕易調(diào)的，要考慮實際情況，所以最后采用先調(diào)變比再調(diào)電廠電壓的方法，符合了要求?！峨娏ο到y(tǒng)穩(wěn)態(tài)分析》，陳珩，中國電力出版社。 %input(39。 %input(39。%input(39。%input(39。 1 5 + 1 0 0。 5 6 + 0 220/11 1 1。 ] %input(39。200 0 230 0 3。0 50+ 10 0 2。0 70+ 10 0 2] %input(39。%input(39。%input(39。%input(39。%input(39。%1 3 4+16i 0 1 0 0。%B2=[0 0 115 0 1。)。V=zeros(1,n)。q=B1(i,2)。 %非對角元 Y(q,p)=Y(p,q)。q=B1(i,2)。 %對角元j側(cè)+線路電納的一半 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./。disp(Y)%給定各節(jié)點初始電壓及給定各節(jié)點注入功率G=real(Y)。 V(i)=abs(B2(i,3))。Q=imag(S)。N1=N0+1。 %求取各個節(jié)點的功率及功率偏差及PV節(jié)點的電壓偏差 for i=1:n %n個節(jié)點2n行(每節(jié)點兩個方程P和Q或U) p=2*i1。 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)。的計算值 P1=C(i)*e(i)+f(i)*D(i)。 %PQ節(jié)點有功功率差J(m,N1)擴展列△P J(p,N1)=Q(i)Q1。 if DET=pr %PQ節(jié)點的功率偏差量及PV節(jié)點的電壓偏差量是否滿足要求 IT2=IT2+1。 %當前不滿足要求的節(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。%Σ(Gij*fj+Bij*ej) end for j1=2:n %第i行共n列(2n個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é)點） %=============== 下面是針對PV節(jié)點來求取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個節(jié)點2n行(每節(jié)點兩個方程P和Q或U) JZ0=[39。]。% 用K行K列對角元素去除K行K列后的非對角元素進行規(guī)格化 end J(k,k)=1。Jacobi矩陣前代第(39。JZ1=[39。]。 J(k1,k)=0。disp(JJN1)。 f(L)=f(L)J(k1,N1)。disp(ICT11)。disp(ICT2)。 %將各節(jié)點電壓用復(fù)數(shù)表示end%=============== 計算各輸出量 ===========================disp(39。 %顯示各節(jié)點的實際電壓值E用復(fù)數(shù)表示disp(39。各節(jié)點的電壓模值大小V為(節(jié)點號從小到大排列)：39。39。)。 %計算各節(jié)點注入電流的共軛值 end S(p)=E(p)*C(p)。disp(S)。disp(39。q=B1(i,2)。 Siz(i)=Si(p,q)。,39。 disp(ZF)。enddisp(39。q=B1(i,2)。 Sjy(i)=Sj(q,p)。,39。 disp(ZF)。enddisp(39。q=B1(i,2)。,num2str(p),39。,num2str(DS(i))]。)。JDQ=[]。(39。 JDDY=strcat(JDDY,num2str(i),39。)。),39。,num2str(imag(S(i))),39。節(jié)點電壓角度：39。,JDDY)。JDQ=strcat(39?？偩W(wǎng)損S=39。,num2str(real(zws)),39。figure(1)。ylabel(39。節(jié)點號39。plot(sida)。)。Q=imag(S)。節(jié)點注入有功39。bar(Q)。節(jié)點注入無功39。subplot(3,2,1)。(39。), 39。,num2str(B1(i,2)),39。)。title(JDH)。%xlabel(39。subplot(3,2,2)。支路首端注入無功39。P2=real(Sjy)。ylabel(39。subplot(3,2,4)。支路末端注入無功39。P3=real(DS)。ylabel(39。subplot(3,2,6)。支路無功損耗3