【正文】 ranchPara(i,5)*(TopoStructureAndBranchPara(i,3)+... j*TopoStructureAndBranchPara(i,4)))))。 for u=1:NumberOfNode sum=sum+real(NAM(NumberOfNode,u))*real(V(u))imag(NAM(NumberOfNode,u))*imag(V(u))。 m=m+1。 end A(n)=sum。 datamax=max(datamax,B(m))。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 end A(n)=sum。 for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。 if NodeData(i,1)==1 NodeData(i,4)=NodeData(i,4)+X(m)。 datamax=0。 else YMatrix(m,n)=NodeData(i,4)*(real(NAM(i,l))*sind(NodeData(i,5)NodeData(l,5))imag(NAM(i,l))*cosd(NodeData(i,5)NodeData(l,5)))。 else YMatrix(m,n)=NodeData(i,4)*NodeData(l,4)*(real(NAM(i,l))*cosd(NodeData(i,5)NodeData(l,5))+imag(NAM(i,l))*sind(NodeData(i,5)NodeData(l,5)))。 end end end m=m+1。 end if NodeData(k,1)==1 if i==k YMatrix(m,n)=NodeData(i,4)*real(NAM(i,i))A(m)。 if NodeData(i,1)==1 YMatrix(m,n)=YMatrix(m,n)+NodeData(i,4)*A(m+1)。 m=m+1。 end A(n)=sum。 datamax=max(datamax,B(m))。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 end A(n)=sum。for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。 floatmin=。Select Node datafile39。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。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。Temporary1=max(TopoStructureAndBranchPara(:,1))。)。謝謝!22參考文獻(xiàn)[1] ，1986.[2] 于永源. ，1996年.[3] [M].北京：中國(guó)電力出版社，1995.[4] 楊少兵,駱平. 電力系統(tǒng)分析的教學(xué)軟件開發(fā)，2000(10).[5] 韓禎祥，吳國(guó)炎等，電力系統(tǒng)分析，浙江大學(xué)出版社，1993年.[6] [M]．武漢：華中理工大學(xué)出版社.[7] ：哈爾濱工業(yè)大學(xué)出版社，1999.[8] ，1989.[9] ，2004.[10] . 水利電力出版社，1994.[11] ：中國(guó)電力出版社，1996.[12] Steven Holzner(美)著. 詳實(shí)翻譯組譯. Visual Basic 6 技術(shù)內(nèi)幕. 機(jī)械工業(yè)出版社，2000.[13] [M].北京：中國(guó)電力出版社，1995.[14] 張伯明,[M].北京：清華大學(xué)出版社1996.[15] [M].長(zhǎng)沙：湖南科學(xué)技術(shù)出版社，1983.[16] ，1989.47附　錄 牛頓－拉夫遜法計(jì)算潮流程序框圖 計(jì)算程序%function OutPut=PowerFlowCalculation(handles)%is a subroutine of PowerSystemCalculationfunction OutPut=PowerFlow(handles)%the following program is open a data file and get the Number of % Node and Branch[fname,pname] = uigetfile(39。從課題的選擇到項(xiàng)目的最終完成，吳老師都始終給予我們細(xì)心的指導(dǎo)和不懈的支持。因此選擇開發(fā)了潮流計(jì)算程序界面,VB 作為一種可視化的編程工具,其編程簡(jiǎn)單、界面友好,但它的計(jì)算能力相對(duì)較弱,因此將MATLAB 與VB 的優(yōu)勢(shì)結(jié)合起來,讓它們?nèi)￠L(zhǎng)補(bǔ)短,進(jìn)行無縫鏈接,可使編程更專業(yè)、更靈活。從效果來看，牛頓－拉夫遜法的迭代次數(shù)較少。這給以后的分析計(jì)算帶來了很大的方便，它有利于節(jié)省內(nèi)存、提高計(jì)算速度以及改善收斂等。當(dāng)電力系統(tǒng)中必需的已知條件給定后潮流分布，取決于網(wǎng)絡(luò)的結(jié)構(gòu)，而網(wǎng)絡(luò)結(jié)構(gòu)在功率方程中的反映是節(jié)點(diǎn)導(dǎo)納矩陣或節(jié)點(diǎn)阻抗矩陣。(5)為電力系統(tǒng)的規(guī)劃和擴(kuò)建提供依據(jù)。它的主要目的：(1) 根據(jù)功率分布，可以選折電力系統(tǒng)的電氣設(shè)備和導(dǎo)線截面積，可以為電力系統(tǒng)繼電保護(hù)整定計(jì)算提供必要的數(shù)據(jù)等。取收斂判據(jù)為|DPi|163。 將牛頓－拉夫遜法用于潮流計(jì)算，要求將潮流方程寫成形如方程式(31)的形式。 應(yīng)用牛頓法求解多變量非線性方程組(31)時(shí)，假定已給出各變量的初值，…. ，令，….. 分別為各變量的修正量，使其滿足方程(31)即 (410)將上式中的n個(gè)多元函數(shù)在初始值附近分別展成泰勒級(jí)數(shù),并略去含有,，……，二次及以上階次的各項(xiàng)，便得 (411)方程式(317)也可以寫成矩陣形式 (412)方程式(318)是對(duì)于修正量,……, 的線性方程組,……, 。函數(shù)y＝f(x)為圖中的曲線。解此方程可得修正量 (45)用所求的去修正近似解，變得 (46)由于(310)是略去高次項(xiàng)的簡(jiǎn)化式，因此所解出的修正量也只是近似值。圖(31)11洛陽(yáng)理工學(xué)院畢業(yè)設(shè)計(jì)（論文）第4章 潮流計(jì)算的原理 牛頓－拉夫遜法設(shè)有單變量非線性方程 (41)求解此方程時(shí)?，F(xiàn)在幾種典型的接線變化說明具體的修改方法。一般來說，即對(duì)角線元素的絕對(duì)值大于非對(duì)角線元素的絕對(duì)值，使節(jié)點(diǎn)導(dǎo)納矩陣稱為具有對(duì)角線優(yōu)勢(shì)的矩陣。其表示式為 (33)節(jié)點(diǎn)導(dǎo)納矩陣的非對(duì)角線元素 (i＝…n ,j=…n但ij)稱為互導(dǎo)納。是節(jié)點(diǎn)電壓的列向量。如果計(jì)算出來的結(jié)果不滿足這些約束條件，必須修改某些變量的給定值，甚至修改系統(tǒng)的運(yùn)行方式。2. 所有電源節(jié)點(diǎn)的有功功率和無功功率必須滿足的條件 和 PQ節(jié)點(diǎn)的有功功率和無功功率以及PU節(jié)點(diǎn)的有功功率，在給定時(shí)就必須滿足上式條件。這些要求構(gòu)成了潮流問題中某些變量的約束條件，通常的約束條件有：1. 所有節(jié)點(diǎn)電壓必須滿足 (i＝1，2，……n) 這個(gè)條件是說各節(jié)點(diǎn)電壓的幅值應(yīng)限制在一定的范圍之內(nèi)。所以方程式的數(shù)目實(shí)際上只有2（n－1）個(gè)。一般選折主調(diào)頻發(fā)電廠為平衡節(jié)點(diǎn)比較合適。另外，必須選定一個(gè)節(jié)點(diǎn)，指定其電壓相位為零，作為計(jì)算各節(jié)點(diǎn)電壓相位的參考，這個(gè)節(jié)點(diǎn)稱為基準(zhǔn)節(jié)點(diǎn)。這類節(jié)點(diǎn)必須有足夠的可調(diào)無功容量，用以維持給定的電壓幅值，因而又稱之為電壓控制節(jié)點(diǎn)。通常變電所都是在這一類型的節(jié)點(diǎn)，由于沒有發(fā)電機(jī)設(shè)備，故發(fā)電機(jī)功率為零。還要給定其中的2個(gè)，將剩下的2個(gè)作為待求變量，方程組才可以求解。對(duì)于n個(gè)節(jié)點(diǎn)的網(wǎng)絡(luò)，可以寫2n個(gè)方程，但是確有6n個(gè)變量。在電氣網(wǎng)絡(luò)理論中，一般是給出電壓源或電流源，為求得網(wǎng)絡(luò)內(nèi)電流和電壓的分布，只要直接求解網(wǎng)絡(luò)方程就可以了。近20多年來，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進(jìn)牛頓法和PQ分解法進(jìn)行的。PQ分解法在計(jì)算速度方面有顯著的提高，迅速得到了推廣。牛頓拉夫遜法是數(shù)學(xué)中求解非線性方程式的典型方法，有較好的收斂性。當(dāng)系統(tǒng)不斷擴(kuò)大時(shí)，這些缺點(diǎn)就更加突出。這就需要較大的內(nèi)存量。在用數(shù)字計(jì)算機(jī)求解電力系統(tǒng)潮流問題的開始階段，人們普遍采用以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的高斯賽德爾迭代法。(6) 為調(diào)整計(jì)算、經(jīng)濟(jì)運(yùn)行計(jì)算、短路計(jì)算和穩(wěn)定計(jì)算提供必要的數(shù)據(jù)。(2) 檢查電力系統(tǒng)各節(jié)點(diǎn)的電壓是否滿足電壓質(zhì)量的要求。在這個(gè)設(shè)計(jì)中，我們選折了MATLAB開發(fā)潮流計(jì)算程序，是因?yàn)槌绷饔?jì)算在數(shù)學(xué)上一般屬于多元非線性代數(shù)方程組的求解，必須采用迭代計(jì)算其中涉及大量的向量和矩陣運(yùn)算，使用傳統(tǒng)的編程語言將十分麻煩。簡(jiǎn)單網(wǎng)絡(luò)的潮流計(jì)算，比如：輻射型網(wǎng)絡(luò)的潮流計(jì)算和閉式網(wǎng)絡(luò)的潮流計(jì)算。本設(shè)計(jì)中還用了一個(gè)五節(jié)點(diǎn)的電力系統(tǒng)網(wǎng)絡(luò)來驗(yàn)證本設(shè)計(jì)在實(shí)際運(yùn)行中的優(yōu)越性。我們可以畫出計(jì)算框圖，用MATLAB編寫出程序，來代替?zhèn)鹘y(tǒng)的手算算法?；跇O坐標(biāo)的牛頓拉夫遜法潮流計(jì)算 洛陽(yáng)理工學(xué)院畢業(yè)設(shè)計(jì)（論文）畢業(yè)設(shè)計(jì)基于極坐標(biāo)的牛頓－拉夫遜法潮流計(jì)算摘 要潮流計(jì)算是電力系統(tǒng)最基本的計(jì)算功能，其基本思想是根據(jù)電力網(wǎng)絡(luò)上某些節(jié)點(diǎn)的已知量求解未知量，潮流計(jì)算在電力系統(tǒng)中有著獨(dú)特的作用。根據(jù)電力系統(tǒng)網(wǎng)絡(luò)的基本知識(shí)，構(gòu)建出能代表電力系統(tǒng)系統(tǒng)網(wǎng)絡(luò)的數(shù)學(xué)模型，然后用牛頓—拉夫遜法反復(fù)計(jì)算出各個(gè)接點(diǎn)的待求量，直到各個(gè)節(jié)點(diǎn)的待求量滿足電力系統(tǒng)的要求。計(jì)算機(jī)計(jì)算已逐漸成為分析復(fù)雜系統(tǒng)潮流分布的主要方法。潮流計(jì)算可以分為簡(jiǎn)單網(wǎng)絡(luò)的潮流計(jì)算和復(fù)雜系統(tǒng)的潮流計(jì)算。潮流計(jì)算在電力系統(tǒng)規(guī)劃設(shè)計(jì)及運(yùn)行方式分析的離線及在線計(jì)算中都發(fā)揮著重要的作用。它的主要目的：(1) 檢查電力系統(tǒng)各元件是否過負(fù)荷。(5) 為電力系統(tǒng)的規(guī)劃和擴(kuò)建提供依據(jù)。潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，它的任務(wù)是對(duì)給定的運(yùn)行條件確定系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓、網(wǎng)絡(luò)中的功率分布以及功率損耗等。阻抗矩陣是滿矩陣，高斯賽德爾迭代法要求計(jì)算機(jī)儲(chǔ)存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。但是，高斯賽德爾迭代法的主要缺點(diǎn)就是占用計(jì)算機(jī)的內(nèi)存很大，每次迭代的計(jì)算量很大?？朔咚官惖聽柕ㄈ秉c(diǎn)的另一途徑是采用牛頓拉夫遜法。