【正文】 ， 在求節(jié)點(diǎn)導(dǎo)納矩陣時(shí)的要用到反饋的指令，所以在 MATLAB 中有下 面幾種命令可以解決兩個(gè)系統(tǒng)間的連接問題。示意圖如下： u1 y1 u + y u2 y2 并聯(lián)后的系統(tǒng)傳遞函數(shù)表示式為： 1 2 2 112 12( ) ( ) ( ) n d n dg s g s g s dd? ? ?? ? ? ?其中 n d1和 n d2分別為 g1（ s）、 g2（ s）的傳遞函數(shù)分子、分母系數(shù)行矢量。 y2 u2 連接后系統(tǒng)的傳遞函數(shù)表示為： 1 1 22 1 2 1 2()() 1 ( )g s n dgs g s d d n n??? ?? 命令格式： [ n， d ] = feedback（ n1， d1， n2， d2） 或： [ n， d ] = feedback（ n1， d1， n2， d2， sign） [ a， b， c， d ] = feedback（ a1， b1， c1， d1， a2， b2， c2， d2， sign） 其中 sign 是指示 y2到 u1連接的符號，缺省時(shí)默認(rèn)為負(fù)（即 sign = 1）。 正、負(fù)反饋后閉環(huán)系統(tǒng)為： 111 1 1()() ()g s ngs m g s n d?? ? 命令格式： [ n， d ] = cloop （ n1， d1， sign） [ ac， bc， cc， dc] = cloop（ a， b， c， d， sign） 通過以上對 MATLAB 基本指令的了解，我們就可以對所求的電力系統(tǒng)網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣進(jìn)行畫編程框架圖。其牛頓 拉夫遜潮流計(jì)算程序框圖如下所示。電力系統(tǒng)極坐標(biāo)下的牛頓 拉夫遜法潮流計(jì)算 :39。 clear n=input(39。)。請輸入 PV 結(jié)點(diǎn)數(shù)： n1=39。 n2=input(39。)。請輸入平衡結(jié)點(diǎn)： isb=39。 輸入原始數(shù)據(jù) 形成節(jié)點(diǎn)導(dǎo)納矩陣 給定節(jié)點(diǎn)電壓初值 ei(0),fi(0) 用公式計(jì)算不平衡功率 △ Pi（ k） i △ Qi（ k） i△ Vi2（ k） i max（ |△ Pi（ K） i△ Qi（ i） i△ Vi2（ k） i|ε ） ε 解修正方程求 △δ （ k） △ V（ k） δ （ k+1） =δ （ k） +△δ （ k） △ V（ k+1） =V（ k） +△ V（ k） K 0 計(jì)算平衡節(jié)點(diǎn)功率及全部路線功率 輸出 K+1=k 是 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 20 pr=input(39。)。請輸入變比矩陣看 :K=39。 C=input(39。)。請輸入支路導(dǎo)納矩陣： y=39。 U=input(39。)。請輸入各結(jié)點(diǎn)的功率： S=39。 Z=zeros(1,n)。L=zeros(n2,n2)。 for m=1:n for R=1:n C(m,m)=C(m,m)+y(m,R)。 C(R,R)=C(R,R)+1/(C(m,R)/(1K(m,R)))。 C(R,m)=C(m,R)。 end 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 21 end end for m=1:n for R=1:n if m==R Y(m,m)=C(m,m)+Z(m)。 end end end disp(39。)。 disp(39。)。 B=imag(Y)。 U1=abs(U)。 PR=1。 Q=imag(S)。 end for m=1:n1+n2 for R=1:n PT(R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)O(R))+B(m,R)*sin(O(m)O(R)))。 PP(m)=P(m)PT1(m)。 end for m=1:n2 for R=1:n QT(R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)O(R))B(m,R)*cos(O(m)O(R)))。 (m)=Q(m)QT1(m)。 end PR1=max(abs(PP))。 PR=max(PR1,PR2)。 else H(m,R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)O(R))B(m,R)*cos(O(m)O(R)))。 else N(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)O(R))+B(m,R)*sin(O(m)O(R)))。 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 24 else J(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)O(R))+B(m,R)*sin(O(m)O(R)))。 else L(m,R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)O(R))B(m,R)*cos(O(m)O(R)))。J L]。 PQ=[PP39。]。 DA1=DA39。 end 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 25 for m=n:n1+n2+n2 UU1(mn1n2)=DA1(m)。 UU=UU1*UD2。 end for m=1:n2 U1(m)=U1(m)+UU(m)。 end for m=1:n2 u(k+1,m)=U1(m)。 end for m=1:n b(m)=U1(m)*cos(O(m))。 end 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 26 U=b+i*c。 end PH=sum(PH1)。 else C1(m,m)=C(m,m)。(m~=R) SS(m,R)=U1(m)^2*conj(C1(m,m))+U(m)*(conj(U(m))conj(U(R)))*conj(C1(m,R))。迭代中的△ P： 39。disp(PP1)。迭代中的△ Q： 39。disp(1)。迭代中相角 :39。disp(o)。迭代中電壓的模 :39。disp(u)。平衡結(jié)點(diǎn)的功率 :39。disp(PH)。全部線路 功率分布 :39。disp(SS)。0 0 0 。0 0 0 0] 請輸入支路阻抗矩陣： C=[0 + + inf。+ + 0 +。 0 0 0。0 0 0 0] 請輸入結(jié)點(diǎn)電壓矩陣： U=[1+0i 1+0i +0i +0i] 請輸入各結(jié)點(diǎn)的功率： S=[ 0] 結(jié)點(diǎn)導(dǎo)納矩陣 : + + 0 + + 0 + + + + 0 0 + + 迭代中的雅克比矩陣 : 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 28 迭代中的△ P： 迭代中的△ Q： 迭代中相角 : 迭代中電壓的模 : 平衡結(jié)點(diǎn)的功率 : 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 29 + 全部線路功率分布 : 0 + 0 0 0 0 + + 0 第七章 總 結(jié) 此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過精密的設(shè)計(jì)和計(jì)算。同時(shí)由于求解過程中用到求節(jié)點(diǎn)導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對以前所學(xué)的知識有了一次很好的溫習(xí)。 第八章 參考文獻(xiàn) [1] 何仰贊等 .電力系統(tǒng)分析 [M]. 武漢：華中理工大學(xué)出版社， [2] 西安交通大學(xué)等 .電力系統(tǒng)計(jì)算 [M].北京：水利電力出版社， [3] 陳 衍 .電力系統(tǒng)穩(wěn)態(tài)分析 [M].北京：水利電力出版社， [4] 李光琦 .電力系統(tǒng)暫態(tài)分析 [M].北京： 水利電力出版社， [5] 于永源，楊綺雯 . 電力系統(tǒng)分析（第二版） [M]. 北京：中國電力出版社，