【正文】 sum1=0。 n=n+1。 end end X=inv(YMatrix)*B39。 n=n+1。 B(m)=NodeData(i,3)NodeData(i,4)*sum1。NodeData = csvread(Nodename)。else NumberOfNode=Temporary2。她嚴肅的科學態(tài)度，嚴謹?shù)闹螌W精神，精益求精的工作作風，深深地感染和激勵著我。潮流計算的數(shù)學模型是以節(jié)點的方程為基礎，推導出相應的功率方程。和為預先給定的小正數(shù)。 （3）在網(wǎng)絡的原有節(jié)點i、j之間切除一條導納為的支路，（見圖3－1，（c）)），其相當在i、j之間增加一條導納為－的支路，因此與i、j有關的元素應作以下修改：1）節(jié)點i、j的自導納增量＝＝－；2）節(jié)點i、j之間的互導納增量＝＝；（4）原有網(wǎng)絡節(jié)點i、j之間的導納由變成見圖3－1，（d）)，相當于在節(jié)點i、j之間切除一條導納為的支路，在增加一條導納為的支路，則與i、j有關的元素應作如下修改： 1）節(jié)點i、j的自導納增量＝＝－； 2）節(jié)點i、j的互導納增量=＝－。在電工原理課中，已講過用節(jié)點導納矩陣表示的節(jié)點電壓方程為 （31）對于n個節(jié)點的網(wǎng)絡，它可以展開為 (32)式（31）中的是節(jié)點注入電流的列向量。由于平衡節(jié)點的電壓已給定，只需要計算其余（n－1）個節(jié)點的電壓。節(jié)點電壓（U，）是待求量。另外，為了解決病態(tài)潮流計算，出現(xiàn)了將潮流計算表示為一個無約束非線性規(guī)劃問題的模型，即非線性規(guī)劃潮流算法。阻抗矩陣是滿矩陣，高斯賽德爾迭代法要求計算機儲存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。潮流計算在電力系統(tǒng)規(guī)劃設計及運行方式分析的離線及在線計算中都發(fā)揮著重要的作用?；跇O坐標的牛頓拉夫遜法潮流計算 洛陽理工學院畢業(yè)設計（論文）畢業(yè)設計基于極坐標的牛頓－拉夫遜法潮流計算摘 要潮流計算是電力系統(tǒng)最基本的計算功能，其基本思想是根據(jù)電力網(wǎng)絡上某些節(jié)點的已知量求解未知量，潮流計算在電力系統(tǒng)中有著獨特的作用。在這個設計中，我們選折了MATLAB開發(fā)潮流計算程序，是因為潮流計算在數(shù)學上一般屬于多元非線性代數(shù)方程組的求解，必須采用迭代計算其中涉及大量的向量和矩陣運算，使用傳統(tǒng)的編程語言將十分麻煩。這就需要較大的內(nèi)存量。近20多年來，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進牛頓法和PQ分解法進行的。通常變電所都是在這一類型的節(jié)點，由于沒有發(fā)電機設備，故發(fā)電機功率為零。所以方程式的數(shù)目實際上只有2（n－1）個。是節(jié)點電壓的列向量。圖(31)11洛陽理工學院畢業(yè)設計（論文）第4章 潮流計算的原理 牛頓－拉夫遜法設有單變量非線性方程 (41)求解此方程時。 將牛頓－拉夫遜法用于潮流計算，要求將潮流方程寫成形如方程式(31)的形式。當電力系統(tǒng)中必需的已知條件給定后潮流分布，取決于網(wǎng)絡的結構，而網(wǎng)絡結構在功率方程中的反映是節(jié)點導納矩陣或節(jié)點阻抗矩陣。從課題的選擇到項目的最終完成，吳老師都始終給予我們細心的指導和不懈的支持。end%The following program is to form the Nodal Admittance Matrix % and the Topologic structure and Branch Parametres are arranged% I,J,R,X,C/K, and pay attention to the inpedence is in the side % Node I and the ratio of transformer is in the side of Node Jfor CircleNumber1=1:NumberOfBranch for CircleNumber2=1:NumberOfBranch NAM(CircleNumber1,CircleNumber2)=0。 floatmin=。 datamax=max(datamax,B(m))。 end if NodeData(k,1)==1 if i==k YMatrix(m,n)=NodeData(i,4)*real(NAM(i,i))A(m)。 datamax=0。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 for u=1:NumberOfNode sum=sum+real(NAM(NumberOfNode,u))*real(V(u))imag(NAM(NumberOfNode,u))*imag(V(u))。 end sum=0。 A(n)=sum1。 end end end m=m+1。 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)))。 m=m+1。Nodename=strcat(pname,fname)。if Temporary1 Temporary2 NumberOfNode=Temporary1。在論文寫作過程中，得到了吳茜瓊老師的親切關懷和耐心的指導。潮流計算是電力系統(tǒng)分析中的一種最基本的計算，它的任務是對給定的運行條件確定系統(tǒng)的運行狀態(tài)，如各母線上的電壓、網(wǎng)絡中的功率分布以及功率損耗等。 迭代過程一直到滿足收斂判據(jù) (418)或 (419)為止。（2）在原有網(wǎng)絡的節(jié)點i與j之間增加一條導納為的支路(見圖3－1，（b）)，則與i、j有關的元素應作如下修改：1）節(jié)點i、j的自導納增量＝＝；2）節(jié)點i、j的互導納增量＝＝－。也可以說是對電力系統(tǒng)的運行狀態(tài)、變量和網(wǎng)絡參數(shù)在電力系統(tǒng)潮流分布的計算中，廣泛采用的是節(jié)點電壓方程。根據(jù)以上所述可以看到，盡管網(wǎng)絡方程是線性方程但是由于在定解條件中不能給定節(jié)點電流，只能給出節(jié)點功率，這就使潮流方程變?yōu)榉蔷€性方程了。這類節(jié)點的有功功率P和無功功率Q是給定的。20世紀70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。20世紀60年代初，數(shù)字計算機已經(jīng)發(fā)展到第二代，計算機的內(nèi)存和計算速度發(fā)生了很大的飛躍，從而為高斯賽德爾迭代法的采用創(chuàng)造了條件。在復雜的電力系統(tǒng)潮流計算中需要對電力系統(tǒng)網(wǎng)絡進行必要的計算，用來獲得必要的數(shù)據(jù)。它不僅能確保電力網(wǎng)絡能夠正常的運行工作、提供較高質量的電能，還能在以后的電力系統(tǒng)擴建中各種計算提供必要的依據(jù)。而MATLAB以復數(shù)矩陣為基本運算單元，且內(nèi)置眾多高精度、高可靠性矩陣、數(shù)組運算函數(shù)、數(shù)值計算方法，可大大提高編程的效率。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個元素進行計算，因此，每次迭代的計算量很大。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡、模糊算法也逐漸被引入潮流計算。若系統(tǒng)中某些發(fā)電廠送出的功率在一定時間內(nèi)為固定時，則該發(fā)電廠母線可作為PQ節(jié)點。 通過求解方程得到了全部節(jié)點電壓以后，就可以進一步計算各類節(jié)點的功率以及網(wǎng)絡中功率的分布。網(wǎng)絡中有接地支路時，通常以大地為參考點，節(jié)點電壓就是各節(jié)點的對地電壓。先給出解的近似值它與真解的誤差為，則將滿足方程，即 (42)將(38)式左邊的函數(shù)在附近展成泰勒級數(shù)，于是便得 (43)式中，,……分別為函數(shù)在處的一階導數(shù)，….，n階導數(shù)。由于節(jié)點電壓可以采用不同的坐標系表示，牛頓－拉夫遜法潮流計算也將相應的采用不同的計算公式。在復雜電力網(wǎng)中，在各個節(jié)點中沒有直接相連的節(jié)點很多，從而使矩陣中有很多零元素、節(jié)點導納矩陣成為稀疏矩陣。多少個日日夜夜，吳老師不僅在學業(yè)上給我們以精心指導，除了敬佩劉孝書老師的專業(yè)水平外，他的治學嚴謹和科學研究的精神也是我永遠學習的榜樣，并將積極影響我今后的學習和工作。 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))) 。datamax=0。 m=m+1。 n=n+1。 m=1。 datamax=max(datamax,B(m))。 sum1=sum1real(NAM(NumberOfNode,u))*imag(V(u))imag(NAM(NumberOfNode,u))*real(V(u))。 end end end for i=1:NumberOfNode V(i)=NodeData(i,4)*cosd(NodeData(i,5))+j*NodeData(i,4)*sind(NodeData(i,5))。 n=n+1。 n=n+1。 end n=n+1。 datamax=max(datamax,B(m))。)。Temporary2=max(TopoStructureAndBranchPara(:,2))。洛陽理工學院畢業(yè)設計論文謝 辭經(jīng)過幾個月的忙碌和工作，本次畢業(yè)論文設計已經(jīng)接近尾聲，在這次的畢業(yè)設計中，我們做的是牛頓－拉夫遜法潮流計算，課題類型為工程設計類的題目，由于經(jīng)驗的匱乏，難免有許多考慮不周全的地方，如果沒有導師的督促指導，以及一起工作的同學們的支持，想要完成這個設計是難以想象的。(6)為調(diào)整計算、經(jīng)濟運行計算、短路計算和穩(wěn)定計算提供必要的數(shù)據(jù)。然后對初始近似值進行修正 (i=1,2,….,n) (413)如此反復迭代，在進行k＋1次迭代時，從求解修正方程式 (414)得到修正量，并對各變量進行修正 (i=1,2,…,n) (415)式(320)和(321)也可以縮寫為 (416)和 (417) 式中的X和分別是由n個變量和修正量組成的n維列向量；F(X)是由n個多元函數(shù)組成的n維列項量；J是n階方陣，稱為雅可比矩陣，它的第i、j個元素是第n個函數(shù)對第j個變量的偏導數(shù)；上角標(k)表示陣的每一個元素都在點處取值。（1）從原有網(wǎng)絡的節(jié)點i引出一條導納為的支路(見圖3－1，（a）)，j為新增加的節(jié)點，由于新增加了一個節(jié)點，所以節(jié)點導納矩陣增加一階，矩陣作如下修改：1）原有節(jié)點i的自導納的增量＝；2）新增節(jié)點j的自導納＝；3）新增的非對角元素＝＝－；其它新增的非對角元均為零。第3章 電力網(wǎng)絡的數(shù)學模型電力網(wǎng)絡的數(shù)學模型指的是將網(wǎng)絡有關參數(shù)和變量及其相互關系歸納起來所組成的、可以反映網(wǎng)絡性能的數(shù)學方程式組。但在進行潮流計算時也可以按照慣例的原則來選折，例如，為了提高導納矩陣法潮流程序收斂性，也可以選折出線較多的發(fā)電廠母線做為平衡節(jié)點。根據(jù)電力系統(tǒng)的實際運行條件，按給定變量的不同，一般將節(jié)點分為以下三種類型。牛頓拉