【正文】 eAndBranchPara(i,2))))^2*(conj((1TopoStructureAndBranchPara(i,5))/(TopoStructureAndBranchPara(i,5)^2*(TopoStructureAndBranchPara(i,3)+j*TopoStructureAndBranchPara(i,4)))))+... V(TopoStructureAndBranchPara(i,2))*conj((V(TopoStructureAndBranchPara(i,2))V(TopoStructureAndBranchPara(i,1)))*(1/(TopoStructureAndBranchPara(i,5)*(TopoStructureAndBranchPara(i,3)+... j*TopoStructureAndBranchPara(i,4)))))。 end Balance=(NodeData(NumberOfNode,4)*sumNodeData(NumberOfNode,5)*sum1)+j*(NodeData(NumberOfNode,5)*sum+NodeData(NumberOfNode,4)*sum1)。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 m=m+1。 for v=1:NumberOfNode sum=sum+NodeData(v,4)*(real(NAM(i,v))*cosd(NodeData(i,5)NodeData(v,5))+imag(NAM(i,v))*sind(NodeData(i,5)NodeData(v,5)))。 for i=1:NumberOfNode1 NodeData(i,5)=NodeData(i,5)+X(m)。 end if NodeData(l,1)==1 if i==l YMatrix(m,n)=NodeData(i,4)*imag(NAM(i,i))A(m)。 else YMatrix(m,n)=NodeData(i,4)*(real(NAM(i,k))*cosd(NodeData(i,5)NodeData(k,5))+imag(NAM(i,k))*sind(NodeData(i,5)NodeData(k,5)))。 for i=1:NumberOfNode1 n=1。 else sum=0。 A(n)=sum1。m=1。*.dat39。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))=... NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))+... 1/... (TopoStructureAndBranchPara(CircleNumber,5)^2*(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))) 。TopoStructureAndBranchPara= csvread(Branchname)。在此謹(jǐn)向吳老師致以誠摯的謝意和崇高的敬意。在計(jì)算初始時(shí)，要選定比較適合的初始值才能滿足計(jì)算結(jié)果的迅速收斂，如果選值不合適，計(jì)算結(jié)果不會(huì)收斂，可能成為發(fā)散型的算式。一般來說，對(duì)角元素的絕對(duì)值大于非對(duì)角元素的絕對(duì)值，使節(jié)點(diǎn)導(dǎo)納矩陣成為具有對(duì)角線優(yōu)勢(shì)的矩陣。(3)根據(jù)對(duì)各種運(yùn)行方式的潮流分布計(jì)算，可以幫助我們正確地選擇系統(tǒng)的接線方式，合理調(diào)整負(fù)荷，以保證電力系統(tǒng)安全、可靠地的運(yùn)行，向用戶供給高質(zhì)量的電能。圖(41)牛頓－拉夫遜方法的幾何意義16 第5章 計(jì)算實(shí)例圖1為一五結(jié)點(diǎn)系統(tǒng)，各支路參數(shù)均為標(biāo)么值。如果第k次迭代中得到，則過點(diǎn)作一切線，此切線同x軸的交點(diǎn)便確定了下一個(gè)近似值。 如果差值很小，(39)式右端的二次及以上階次的各項(xiàng)均可略去。 節(jié)點(diǎn)導(dǎo)納矩陣的修改在電力系統(tǒng)中，接線方式或運(yùn)行狀態(tài)等均會(huì)發(fā)生變化，從而使網(wǎng)絡(luò)接線改變。 節(jié)點(diǎn)導(dǎo)納矩陣的形成節(jié)點(diǎn)導(dǎo)納矩陣的對(duì)角線元素稱為自導(dǎo)納。3. 某些節(jié)點(diǎn)之間電壓的相位差應(yīng)滿足 為了保證系統(tǒng)運(yùn)行的穩(wěn)定性，要求某些輸電線路兩端電壓相位差不超過一定的數(shù)值。這些計(jì)算結(jié)果代表了潮流方程在數(shù)學(xué)上的一組解答。（亦稱為松弛節(jié)點(diǎn)、搖擺節(jié)點(diǎn)）?？梢婋娏ο到y(tǒng)的絕大多數(shù)節(jié)點(diǎn)屬于這一類型。通常把負(fù)荷功率作已知量，并把節(jié)點(diǎn)功率和引入網(wǎng)絡(luò)方程。但是，到目前為止這些新的模型和算法還不能取代牛頓拉夫遜法和PQ分解法的地位。自從20世紀(jì)60年代中期采用了最佳順序消去法以后，牛頓拉夫遜法在收斂性、內(nèi)存要求、計(jì)算速度方面都超過了阻抗法，成為直到目前仍被廣泛采用的方法。高斯賽德爾迭代法改善了電力系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納法無法解決的一些系統(tǒng)的潮流計(jì)算，在當(dāng)時(shí)獲得了廣泛的應(yīng)用，曾為我國電力系統(tǒng)設(shè)計(jì)、運(yùn)行和研究做出了很大的貢獻(xiàn)。隨著電子計(jì)算機(jī)技術(shù)的進(jìn)步，電力系統(tǒng)潮流分布的計(jì)算幾乎已普遍采用計(jì)算機(jī)來進(jìn)行，通過求解描述電力系統(tǒng)狀態(tài)的數(shù)學(xué)模型，而得到較精確的解。第一章 電力系統(tǒng)潮流計(jì)算的基本知識(shí)電力系統(tǒng)潮流計(jì)算分布計(jì)算，是指電力系統(tǒng)在某一穩(wěn)定狀態(tài)的正常運(yùn)行方式下，電力網(wǎng)絡(luò)各節(jié)點(diǎn)的電壓和功率分布的計(jì)算。其任務(wù)是根據(jù)給定的網(wǎng)絡(luò)結(jié)構(gòu)及運(yùn)行條件，求出整個(gè)網(wǎng)絡(luò)的運(yùn)行狀態(tài)，其中個(gè)母線的電壓、網(wǎng)絡(luò)中的功率分布以及整個(gè)系統(tǒng)的功率損耗等。計(jì)算潮流分布的方法很多，本設(shè)計(jì)主要用的是基于極坐標(biāo)的牛頓－拉夫遜法。復(fù)雜電力系統(tǒng)是一個(gè)包括大量母線、支路的龐大系統(tǒng)。它們是復(fù)雜電力系統(tǒng)潮流計(jì)算的基礎(chǔ)。(3) 根據(jù)對(duì)各種運(yùn)行方式的潮流分布計(jì)算，可以幫助我們正確地選擇系統(tǒng)的接線方式，合理調(diào)整負(fù)荷，以保證電力系統(tǒng)安全、可靠地的運(yùn)行，向用戶供給高質(zhì)量的電能。這個(gè)方法的原理比較簡(jiǎn)單，要求的數(shù)字計(jì)算機(jī)的內(nèi)存量也比較小，適應(yīng)當(dāng)時(shí)的電子數(shù)字計(jì)算機(jī)制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計(jì)算人員轉(zhuǎn)向以阻抗矩陣為主的逐次代入法。為了克服阻抗法在內(nèi)存和速度方面的缺點(diǎn)，后來發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。牛頓拉夫遜法的特點(diǎn)是將非線性方程線性化。但是，在潮流計(jì)算中，在網(wǎng)絡(luò)的運(yùn)行狀態(tài)求出以前，無論是電源的電勢(shì)值，還是節(jié)點(diǎn)的注入的電流，都是無法準(zhǔn)確給定的。根據(jù)電力系統(tǒng)的實(shí)際運(yùn)行條件，按給定變量的不同，一般將節(jié)點(diǎn)分為以下三種類型。一般是選折有一定無功儲(chǔ)備的發(fā)電廠和具有可調(diào)無功電源設(shè)備的變電所作為PU節(jié)點(diǎn)。但在進(jìn)行潮流計(jì)算時(shí)也可以按照慣例的原則來選折，例如，為了提高導(dǎo)納矩陣法潮流程序收斂性，也可以選折出線較多的發(fā)電廠母線做為平衡節(jié)點(diǎn)。從保證電能質(zhì)量和供電安全的要求來看，電力系統(tǒng)的所有電氣設(shè)備都必須運(yùn)行在額定電壓附近。第3章 電力網(wǎng)絡(luò)的數(shù)學(xué)模型電力網(wǎng)絡(luò)的數(shù)學(xué)模型指的是將網(wǎng)絡(luò)有關(guān)參數(shù)和變量及其相互關(guān)系歸納起來所組成的、可以反映網(wǎng)絡(luò)性能的數(shù)學(xué)方程式組。由式（3－3）可見，互導(dǎo)納在數(shù)值上就等于節(jié)點(diǎn)i施加單位電壓，其它節(jié)點(diǎn)全部接地時(shí)，經(jīng)節(jié)點(diǎn)j注入網(wǎng)絡(luò)的電流。（1）從原有網(wǎng)絡(luò)的節(jié)點(diǎn)i引出一條導(dǎo)納為的支路(見圖3－1，（a）)，j為新增加的節(jié)點(diǎn)，由于新增加了一個(gè)節(jié)點(diǎn)，所以節(jié)點(diǎn)導(dǎo)納矩陣增加一階，矩陣作如下修改：1）原有節(jié)點(diǎn)i的自導(dǎo)納的增量＝；2）新增節(jié)點(diǎn)j的自導(dǎo)納＝；3）新增的非對(duì)角元素＝＝－；其它新增的非對(duì)角元均為零。修正后的近似解同真解仍然有誤差。然后對(duì)初始近似值進(jìn)行修正 (i=1,2,….,n) (413)如此反復(fù)迭代，在進(jìn)行k＋1次迭代時(shí)，從求解修正方程式 (414)得到修正量，并對(duì)各變量進(jìn)行修正 (i=1,2,…,n) (415)式(320)和(321)也可以縮寫為 (416)和 (417) 式中的X和分別是由n個(gè)變量和修正量組成的n維列向量；F(X)是由n個(gè)多元函數(shù)組成的n維列項(xiàng)量；J是n階方陣，稱為雅可比矩陣，它的第i、j個(gè)元素是第n個(gè)函數(shù)對(duì)第j個(gè)變量的偏導(dǎo)數(shù)；上角標(biāo)(k)表示陣的每一個(gè)元素都在點(diǎn)處取值。105和|DQi(DVi2)|163。(6)為調(diào)整計(jì)算、經(jīng)濟(jì)運(yùn)行計(jì)算、短路計(jì)算和穩(wěn)定計(jì)算提供必要的數(shù)據(jù)。功率方程是非線性代數(shù)方程組，必須采用數(shù)值求解的方法進(jìn)行計(jì)算。洛陽理工學(xué)院畢業(yè)設(shè)計(jì)論文謝 辭經(jīng)過幾個(gè)月的忙碌和工作，本次畢業(yè)論文設(shè)計(jì)已經(jīng)接近尾聲，在這次的畢業(yè)設(shè)計(jì)中，我們做的是牛頓－拉夫遜法潮流計(jì)算，課題類型為工程設(shè)計(jì)類的題目，由于經(jīng)驗(yàn)的匱乏，難免有許多考慮不周全的地方，如果沒有導(dǎo)師的督促指導(dǎo)，以及一起工作的同學(xué)們的支持，想要完成這個(gè)設(shè)計(jì)是難以想象的。*.dat39。Temporary2=max(TopoStructureAndBranchPara(:,2))。 else NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1))=... NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1))+... +1/(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))+j*TopoStructureAndBranchPara(CircleNumber,5)。)。 sum1=0。 datamax=max(datamax,B(m))。 n=n+1。 end n=n+1。 if NodeData(i,1)==1 n=1。 n=n+1。 m=m+1。 n=n+1。 m=m+1。 end end end for i=1:NumberOfNode V(i)=NodeData(i,4)*cosd(NodeData(i,5))+j*NodeData(i,4)*sind(NodeData(i,5))。 else S(TopoStructureAndBranchPara(i,1),TopoStructureAndBranchPara(i,2))=(abs(V(TopoStructureAndBranchPara(i,1))))^2*(j)*TopoStructureAndBranchPara(i,5)+V(TopoStructureAndBranchPara(i,1))*... conj((V(TopoStructureAndBranchPara(i,1))V(TopoStructureAndBranchPara(i,2)))*1/(TopoStructureAndBranchPara(i,3)+j*TopoStructureAndBranchPara(i,4)))。 sum1=sum1real(NAM(NumberOfNode,u))*imag(V(u))imag(NAM(NumberOfNode,u))*real(V(u))。 n=n+1。 datamax=max(datamax,B(m))。 sum1=0。 m=1。 n=n+1。 n=n+1。 endend while datamaxfloatmin m=1。 m=m+1。 n=n+1。datamax=0。 endend%[fname,p