【正文】 Q=imag(S)。V(i)=B2(i,4)。disp(Y)%G=real(Y)。%對角元K側(cè)Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。q=B1(i,1)。S1=zeros(nl)。e=zeros(1,n)。% B1矩陣：支路首端號；末端號；支路阻抗；支路對地電納 %支路的變比；支路首端處于K側(cè)為1，1側(cè)為0 % B2矩陣：該節(jié)點發(fā)電機(jī)功率；該節(jié)點負(fù)荷功率；節(jié)點電壓初始值 %PV節(jié)點電壓V的給定值；節(jié)點所接的無功補(bǔ)償設(shè)備的容量 %節(jié)點分類標(biāo)號：1為平衡節(jié)點（應(yīng)為1號節(jié)點）；2為PQ節(jié)點； %3為PV節(jié)點；%input(39。00002。00002。+01。+0。+01。%input(39。%input(39。%input(39。%input(39。X為節(jié)點號和對地參數(shù)矩陣。4==ViDVL5==4設(shè)計圖及源程序相關(guān)的原始數(shù)據(jù)輸入格式如下：B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點編號。==37DVT63==180。35=(+)()180。+180。VN=()180。VN=()180。180。VN=()180。+45180。(+j4)=(+)MVA352Si=Sh+DSL5+DS06+jDQL5+5=(+)MVA 152+=180。Se=SG3+Sd3025DST4DS04jQL3DSL3=(+)MVAFLR4： DST5P2+Q2502+312=(RT5+jXT5)=(+)=(+)MVA22VN110P2++=(RL4+jXL4)=(+j24)=(+)MVA22VN110Se39。104S1DQL5=BL3VN2=FLR1：P2102DST1=2(RT1+jXT1)=(+)=(+)MVA2VN110Sa=10MW+DST1+DS01+jDQL1=(+)MVAP2++=(RL1+jXL1)=(+j16)=(+)MVAVN21102DST2P2+Q2402+452=(RT2+jXT2)=(+)=(+)MVAVN21102FLR2Sb=SG120DST2=60+=(+)MVASc=SbSa25jQL1DSL1=(+)MVA：DST3P2+Q22522+1562=(RT3+jXT3)=(+)=(+)MVA22VN110P++(R+jX)=(+j52)=(+)MVAL2L222VN1102222Sc39。l5=2180。l5=180。l5=180。106180。60==x0180。106180。70=14W BL3=2180。70==2180。104S DQ1L2=2BL2V2N==12180。130=52WB6L2=b0180。104S DQ1L1=2BL1V2N==r0180。40=16WB6L1=b0180。10S/km 6線路L1：線路L2：線路L3：（雙回路）線路L4：線路L5：（雙回路）RL1=r0180。103DS07=2180。180。103=(10180。3523=180。10= 2SNDS06=2180。180。VNXT61=180。10)變壓器T6：（兩個三繞組變壓器并聯(lián)）RT61=RT62=RT631180。(Dp0+j變壓器T4：（雙并聯(lián)）1RT1= XT4=XT1==2DS01=(+)MVART4=變壓器T5：RT5=4RT3==4XT3==1DS03=(+)MVA4163180。180。180。10=180。SN)=(+)MVA100 1Dps180。180。180。10=180。SN=(+)MVA 100RT2XT221Dps180。1102XT1=180。10=180。tan()=其次是變電站的參數(shù)計算，我們還是對7個變電站簡化為7臺變壓器來計算。tan()=27MVar發(fā)電機(jī)G4：P4=1180。tan()=45MVar發(fā)電機(jī)G2：P2=4180。I0%=%=(12)80MW239。30km35kV0km31236。所建立的模型和建立的方法在詳細(xì)設(shè)計中詳述。50MW“236。Dps=121kW10kV237。I%=%===73kW237。10MWVA35kV32km25MW110kV80km25MW110kV70km110kVYN,d112180。63MW“236。10MWYN,d114180。潮流計算經(jīng)歷了一個由手工到應(yīng)用數(shù)字電子計算機(jī)的發(fā)展過程，現(xiàn)在的潮流算法都以計算機(jī)的應(yīng)用為前提用計算機(jī)進(jìn)行潮流計算主要步驟在于編制計算機(jī)程序，這是一項非常復(fù)雜的工作。因此，潮流計算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。常規(guī)潮流計算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。關(guān)鍵詞： 電力系統(tǒng) 潮流計算 MATLABAbstract Electric power system steady flow calculation and analysis of the static safety paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide putation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the ponents of the power loss, and getting electric standard is to use the power flow calculation and analysis, the specific have MATLAB words: Power system。S %求節(jié)點注入電流 I=Y*U39。end ox,U,e,count %求節(jié)點注入的凈功率 i=5。oU(i)=ox(i+4)*U(i)。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。N(i,i)=N(i,i)U(i)*U(j)*(G(i, j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))。N(i,j)=U(i)*U(j)*(G(i,j)*cos(e(i)e(j))+B(i,j)*sin(e(i)e(j)))。end oP(i)=oP(i)+P(i)。oP(i)=0。count=0 %計算迭代次數(shù) while max(fx)1e5 for i=1:4for j=1:4H(i,j)=0。U(5)=。S(5)=0。37電力系統(tǒng)潮流分析計算的MATLAB仿真%原始節(jié)點功率S(1)=+。%求互導(dǎo)納 for i=1:5 for j=1:5 if i~=j Y(i,j)=y(i,j)。y(3,4)=1/(+)。y(1,4)=1/(+)。感謝老師們百忙中審閱我的論文！在論文即將完成之際，我的心情非常激動，從開始進(jìn)入課題到論文的順利完成，有多少可敬的師長、同學(xué)、朋友給了我無言的幫助，在這里請接受我誠摯的謝意！35電力系統(tǒng)潮流分析計算的MATLAB仿真參考文獻(xiàn)[1] 于永源,(第二版).中國電力出版社, 2004.[2] , 1997.[3] 《電氣工程師手冊》, 2000.[4] 何仰贊, 2002.[5] 吳際舜, 2000.[6] (第三版).中國電力出版社, 2007.[7] , 2006.[8] , 2007.[9] [J].電力系統(tǒng)及其自動化學(xué)報, 1999.[10] , 2003.[11] 張志涌, 附件 MATLAB計算實例程序clear。四年的時光，仿佛就在昨天。她的教誨，令我終身難忘。34電力系統(tǒng)潮流分析計算的MATLAB仿真致謝大學(xué)四年的學(xué)習(xí)生活即將結(jié)束了，本次畢業(yè)設(shè)計也已經(jīng)接近尾聲，作為一個本科生的畢業(yè)設(shè)計，由于經(jīng)驗的匱乏，難免有許多考慮不周全的地方，但在各位老師的督促指導(dǎo)下，以及同學(xué)們的關(guān)心和幫助下，我的各方面都有所提高，整體素質(zhì)也有了很大的進(jìn)步。通過實例計算分析，取得了比較滿意的效果。27電力系統(tǒng)潮流分析計算的MATLAB仿真Y值：迭代過程：28電力系統(tǒng)潮流分析計算的MATLAB仿真電力系統(tǒng)潮流分析計算的MATLAB仿真30電力系統(tǒng)潮流分析計算的MATLAB仿真31電力系統(tǒng)潮流分析計算的MATLAB仿真電壓值：32電力系統(tǒng)潮流分析計算的MATLAB仿真平衡節(jié)點注入功率及電流：圖電力系統(tǒng)潮流分析計算的MATLAB仿真結(jié)論在電力系統(tǒng)的潮流計算算法中牛頓一拉夫遜法是得到電力系統(tǒng)研究人員認(rèn)可的算法之一，在本文中我們采用牛頓一拉夫遜法，主要是同時考慮到計算的準(zhǔn)確和程序的運(yùn)行速度。有了“對應(yīng)元素”的規(guī)定，數(shù)組的運(yùn)算實質(zhì)上就是針對數(shù)組內(nèi)部的每個元素進(jìn)行的。而乘除法運(yùn)算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對應(yīng)元素之間的乘除法，它們的運(yùn)算符為“.*”和“./”或“.”。但需注意進(jìn)行數(shù)除時，常數(shù)通常只能做除數(shù)。在傳統(tǒng)的MATLAB算法中，右除是先計算矩陣的逆再相乘，而左除則不需要計算逆矩陣直接進(jìn)行除運(yùn)算。原始數(shù)據(jù)輸入格式的設(shè)計，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運(yùn)算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強(qiáng)大的功能。j238。DU2S=iij182。DU182。(GBijejijfj)+Giiei+Biifi(j=i)j206。f=237。Gijei+Bijfi)L182。(Gijfj+Bijej)+BiieiGiifi(j=i)238。DQi239。j206。 Nij=(318)(Gijfj+Bijej)+BiieiGiifi(j=i)182。(j185。(GijejBijfj)GiieiBiifi(j=i)239。(j185。i185。254。182。=2ei239。182。n182。182。n239。DPi=229。182。DPi182。ej182。22182。(315)182。182。182。182。DUn1182。DUm+1182。fn1182。182。fm+1L182。f1L182。en1M182。fm+1M182。em+1M182。fmM182。emM182。f1M182。en1182。en1M182。DU2m+1L182。DQm182。DQ1182。em+1182。em+1182。em+1182。Pm+1182。DPm182。DP1182。em182。emM182。DU2m+1L182。DQm182。DQ1182。e1235。2234。182。e1234。2234。182。DQm234。182。M234。DQ1234。182。其中。234。234。m+1234。234。234。M234。Df1234。233。2234。234。m+1234。DP234。234。M234。DQ1234。nni=m+1,m+2,L,n1(3)對于平衡節(jié)點平衡節(jié)點只設(shè)一個,電壓為已知,不參見迭代,其電壓為:Un=en+jfn(312)2：修正方程兩組迭代式中包括2(n1),并將其按泰勒級數(shù)展開,略去Dei,Dfi二次方程及以后各項,得到修正方程如下:DW=JDU(313)20電力系統(tǒng)潮流分析計算的MATLAB仿真233。(Gijfj+Bijej)239。j=1j=1254。(310)nnDQi=Qifi229。(GijejBijfj)fi229。以直角坐標(biāo)系形式表示：1：迭代推算式采用直角坐標(biāo)時,節(jié)點電壓相量及復(fù)數(shù)導(dǎo)納可表示為:amp。253。i2ijejBfij)f229。)=0239。(Gfijj206。(Gfjijj206。事實上，除了平衡節(jié)點的功率方程式在迭代過程中沒有約束作用以外，其余每個節(jié)點都可以列出兩個方程式。解決這個問題的辦法可以用高斯法迭代1~2次，以此迭代結(jié)果作為牛頓法的初值。對于正常運(yùn)行的系統(tǒng)，各節(jié)點電壓一般均在額定值附近，偏移不會太大，并且各節(jié)點間的相位角差也不大，所以對各節(jié)點可以采用統(tǒng)一的電壓初值(也稱為平直電壓)，如假定：Ui(0)(0)=1 qi(0)=0或ei(0)=1 fi(0)=0(i=q1,2,L,n。牛頓法也具有良好的收斂可靠性，對于對以節(jié)點導(dǎo)納矩陣為基礎(chǔ)的高斯法呈病態(tài)的系統(tǒng)，牛頓法也能可靠收斂。有上式可見，牛頓法的核心便是反復(fù)形式并求解修正方程式。接著就從x(0)出發(fā)，重復(fù)上述計算過程。對于非線性代數(shù)方程組：f(x)=0 即fi(x1,x2,L,xn)=0(i=1,2L,n,(31)在待求量x的某一個初始估計值x(0)附近，將上式展開成泰勒級數(shù)并略去二階及以上的高階項，得到如下的經(jīng)線性化的方程組：39。實際電力系統(tǒng)的潮流技術(shù)那主要采用牛頓拉夫遜法。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負(fù)荷。常用的方法是迭代法和牛頓法，在計算過程15電力系統(tǒng)潮流分析計算的MATLAB仿真中，或得出結(jié)果之后用約束條件進(jìn)行檢驗。axGmiaxPQ節(jié)點的有功功率和無功功率，以及PU節(jié)點的有功功率，在給定是就必須滿足上述條件，因此，對平衡節(jié)點的P和Q以及PU節(jié)點的Q應(yīng)按上述條件進(jìn)行檢驗。Q163。因此，這一約束條件對PQ節(jié)點而言。這些要求夠成了潮流問題中某些變量的約束條件，常用的約束條件如下：① 節(jié)點電壓應(yīng)滿足Uimin163。關(guān)于平衡節(jié)點的選擇，一般選擇系統(tǒng)中擔(dān)任調(diào)頻調(diào)壓的某一發(fā)電廠(或發(fā)電機(jī))，有時也可能按其他原則選擇，例如，為提高計算的收斂性。PU節(jié)點上的發(fā)電機(jī)稱為PU機(jī)(或PU給定型發(fā)電機(jī))平衡節(jié)點在潮流計算中，這類節(jié)點一般只設(shè)一個。PU節(jié)點這類節(jié)點給出的參數(shù)是該節(jié)點的有功功率P及電壓幅值U，待求量為該節(jié)點的無功功率Q及電壓向量的相角q。通常變電所母線都是PQ節(jié)點，當(dāng)某些發(fā)電機(jī)的輸出功率P。 電力系統(tǒng)節(jié)點分類用一般的電路理論求解網(wǎng)絡(luò)方程，目的是給出電壓源(或電流源)研究網(wǎng)絡(luò)內(nèi)的電流(或電壓)分布，作為基礎(chǔ)的方程式，一般用線性代數(shù)方程式表示。j=1254。(217)namp。Pi=Ui229。nnamp。(Gi