【正文】 name] = uigetfile(39。 endend for CircleNumber=1:NumberOfBranch if TopoStructureAndBranchPara(CircleNumber,5) NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1))=... NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1))+... 1/... ((TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))) 。Branchname=strcat(pname,fname)。多少個日日夜夜，吳老師不僅在學(xué)業(yè)上給我們以精心指導(dǎo)，除了敬佩劉孝書老師的專業(yè)水平外，他的治學(xué)嚴謹和科學(xué)研究的精神也是我永遠學(xué)習(xí)的榜樣，并將積極影響我今后的學(xué)習(xí)和工作。從計算的速度來看，速度比較快。在復(fù)雜電力網(wǎng)中，在各個節(jié)點中沒有直接相連的節(jié)點很多，從而使矩陣中有很多零元素、節(jié)點導(dǎo)納矩陣成為稀疏矩陣。(2)檢查電力系統(tǒng)各節(jié)點的電壓是否滿足電壓質(zhì)量的要求。由于節(jié)點電壓可以采用不同的坐標(biāo)系表示，牛頓－拉夫遜法潮流計算也將相應(yīng)的采用不同的計算公式。f(x)＝0的解相當(dāng)于曲線與x軸的交點。先給出解的近似值它與真解的誤差為，則將滿足方程，即 (42)將(38)式左邊的函數(shù)在附近展成泰勒級數(shù)，于是便得 (43)式中，,……分別為函數(shù)在處的一階導(dǎo)數(shù)，….，n階導(dǎo)數(shù)。因此節(jié)點導(dǎo)納矩陣是一個對稱、稀疏且具有對角線優(yōu)勢的方陣。網(wǎng)絡(luò)中有接地支路時，通常以大地為參考點，節(jié)點電壓就是各節(jié)點的對地電壓。因此對平衡節(jié)點的P和Q以及PU節(jié)點Q應(yīng)按上述條件進行檢驗。 通過求解方程得到了全部節(jié)點電壓以后，就可以進一步計算各類節(jié)點的功率以及網(wǎng)絡(luò)中功率的分布?；鶞?zhǔn)節(jié)點的電壓幅值是給定的。若系統(tǒng)中某些發(fā)電廠送出的功率在一定時間內(nèi)為固定時，則該發(fā)電廠母線可作為PQ節(jié)點。因此，對于每個節(jié)點，必須給定這6個變量中的4個，使待求量的數(shù)目同方程的數(shù)目相等，才能對方程求解。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡(luò)、模糊算法也逐漸被引入潮流計算。解決電力系統(tǒng)潮流計算問題是以導(dǎo)納矩陣為基礎(chǔ)的，因此，只要在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計算效率。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個元素進行計算，因此，每次迭代的計算量很大。在計算機技術(shù)還未發(fā)展以前，電力系統(tǒng)的潮流分布計算多采用“手工”近似計算，即按照電路的基本關(guān)系，用手工來推算各節(jié)點的功率和電壓。而MATLAB以復(fù)數(shù)矩陣為基本運算單元，且內(nèi)置眾多高精度、高可靠性矩陣、數(shù)組運算函數(shù)、數(shù)值計算方法，可大大提高編程的效率。關(guān)鍵詞：牛頓－拉夫遜法，復(fù)雜電力系統(tǒng)，潮流計算The method of Newton Raphson based on polarABSTRACTPower system load flow calculation is the most basic puting functions, the basic idea is based on some of the electricity network nodes to solve the unknown quantity of known volume,In power system, power flow, which can ensure that electrical net can work well and give the high quality power, but also later provide the necessary datas in the enlargement of the power system. has special function.There are lots of methods about power flow. We mainly use the method of NewtonRaphson based on polar in my design. According to the basic knowledge of the electrical network, we established the mathematics model which can presents the power system ,then puted again and again unknown members of the each bus with the method of NewtonRaphSon until the unknown numbers meet the demand of the power system. We can write down the block diagram and write the order with the Matlab in place of the traditional methods. Complex power system is a large system which involves lots of bus bars and branches. We also chose a fivebus power system for testing the advantages in the relity.KEY WORDS: NewtonRaphson，power system，power flow4目　錄前　言 1第一章 電力系統(tǒng)潮流計算的基本知識 2 2 2第二章 潮流計算的節(jié)點 5 節(jié)點的分類 5 7第三章 電力網(wǎng)絡(luò)的數(shù)學(xué)模型 8 節(jié)點導(dǎo)納矩陣的形成 9 節(jié)點導(dǎo)納矩陣的修改 9第四章 潮流計算的原理 12 牛頓－拉夫遜法 12第五章 計算實例 17 17 節(jié)點導(dǎo)納的形成 17 計算結(jié)果 18結(jié)　論 20謝 辭 22參考文獻 23附　錄 24計算程序 25外文資料翻譯 41前　言 潮流計算是電力系統(tǒng)中應(yīng)用最廣泛和最重要的一種電氣計算。它不僅能確保電力網(wǎng)絡(luò)能夠正常的運行工作、提供較高質(zhì)量的電能，還能在以后的電力系統(tǒng)擴建中各種計算提供必要的依據(jù)。對這樣的系統(tǒng)進行潮流分析時，采用人工計算的方法已經(jīng)不再適用。在復(fù)雜的電力系統(tǒng)潮流計算中需要對電力系統(tǒng)網(wǎng)絡(luò)進行必要的計算，用來獲得必要的數(shù)據(jù)。(4) 根據(jù)功率分布，可以選折電力系統(tǒng)的電氣設(shè)備和導(dǎo)線截面積，可以為電力系統(tǒng)繼電保護整定計算提供必要的數(shù)據(jù)等。20世紀60年代初，數(shù)字計算機已經(jīng)發(fā)展到第二代，計算機的內(nèi)存和計算速度發(fā)生了很大的飛躍，從而為高斯賽德爾迭代法的采用創(chuàng)造了條件。這個方法把一個大系統(tǒng)分割為幾個小的地區(qū)系統(tǒng)，在計算機內(nèi)只需存儲各個地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡(luò)線的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時也提高了節(jié)省速度。20世紀70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。圖22表示某個三節(jié)點的簡單電力系統(tǒng)及其等值電路，其網(wǎng)絡(luò)方程為即 （i＝1，2，3） (21) 因為，所以節(jié)點電流用功率和電壓可以表示為 (22)將式(22)帶入式(21)可得這是一組復(fù)數(shù)方程式，如果把實部和虛部分開，便得到6個實數(shù)方程。這類節(jié)點的有功功率P和無功功率Q是給定的。在電力系統(tǒng)只能中，這一類的數(shù)目很少。根據(jù)以上所述可以看到，盡管網(wǎng)絡(luò)方程是線性方程但是由于在定解條件中不能給定節(jié)點電流，只能給出節(jié)點功率，這就使潮流方程變?yōu)榉蔷€性方程了。對于平衡節(jié)點的PU節(jié)點，其電壓幅值必須按上述條件給定。也可以說是對電力系統(tǒng)的運行狀態(tài)、變量和網(wǎng)絡(luò)參數(shù)在電力系統(tǒng)潮流分布的計算中，廣泛采用的是節(jié)點電壓方程。亦等于節(jié)點i，j之間所連支路元件導(dǎo)納的負值，其表示式為 (34)依據(jù)互導(dǎo)納的物理意義可知＝，即＝；特別地，當(dāng)節(jié)點i、j之間無直接支路相連時，＝＝0。（2）在原有網(wǎng)絡(luò)的節(jié)點i與j之間增加一條導(dǎo)納為的支路(見圖3－1，（b）)，則與i、j有關(guān)的元素應(yīng)作如下修改：1）節(jié)點i、j的自導(dǎo)納增量＝＝；2）節(jié)點i、j的互導(dǎo)納增量＝＝－。但是，這樣的迭代計算可以反復(fù)進行下去，迭代計算的通式是 (47)迭代過程的收斂判據(jù)為 (48)或 (49)式中，為預(yù)先給定的小正數(shù)。 迭代過程一直到滿足收斂判據(jù) (418)或 (419)為止。105。潮流計算是電力系統(tǒng)分析中的一種最基本的計算，它的任務(wù)是對給定的運行條件確定系統(tǒng)的運行狀態(tài)，如各母線上的電壓、網(wǎng)絡(luò)中的功率分布以及功率損耗等。在首先討論了電力系統(tǒng)節(jié)點的分類以及潮流計算結(jié)果的約束條件后，具體介紹了常用于潮流計算的主要方法－牛頓－拉夫遜法。在論文寫作過程中，得到了吳茜瓊老師的親切關(guān)懷和耐心的指導(dǎo)。,39。if Temporary1 Temporary2 NumberOfNode=Temporary1。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))=... NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))+... +1/(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))+j*TopoStructureAndBranchPara(CircleNumber,5)。Nodename=strcat(pname,fname)。 for w=1:NumberOfNode sum=sum+NodeData(w,4)*(real(NAM(i,w))*cosd(NodeData(i,5)NodeData(w,5))+imag(NAM(i,w))*sind(NodeData(i,5)NodeData(w,5)))。 m=m+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 else YMatrix(m,n)=NodeData(i,4)*NodeData(k,4)*(real(NAM(i,k))*sind(NodeData(i,5)NodeData(k,5))imag(NAM(i,k))*cosd(NodeData(i,5)NodeData(k,5)))。 for l=1:NumberOfNode1 if i==l YMatrix(m,n)=NodeData(i,4)^2*real(NAM(i,i))+NodeData(i,4)*A(m1)。 end end end m=m+1。 end end m=1。 A(n)=sum1。 else sum=0。 end sum=0。 S(TopoStructureAndBranchPara(i,2),TopoStructureAndBranchPara(i,1))=(abs(V(TopoStructureAndBranchPara(i,2))))^2*(j)*TopoStructureAndBranchPara(i,5)+V(TopoStructureAndBranchPara(i,2))*... conj((V(TopoStructureAndBranchPara(i,2))V(TopoStructureAndBranchPara(i,1)))*1/(TopoStructureAndBranchPara(i,3)+j*TopoStructureAndBranchPara(