【正文】 ) ............ ... ... ...( 0) ( 0) ( 0)( , , ... , )...| | || | || | |nnnnn nn n nnf f fx x xfx x xf f ff x x xx x xf x x xf f fx x x? ? ? ??? ? ???????? ? ??????? ? ? ?????? ? ? ?? ? ??12( 0)( 0)...( 0)nxxx??? ???? ??? ??????? ???? ?????????? (412) 方程式 (318)是對(duì)于修正量1(0)x?,2(0)x?,…… , (0)nx? 的線性方程組 ,稱為牛頓法的修正方程式 .利用高斯消去法或三角分解法可以解出修正量1(0)x?,2(0)x?,…… , (0)nx?。取收斂判據(jù)為 |?Pi|?105 和 |?Qi(?Vi2)|?105。 潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，它的任務(wù)是對(duì)給定的運(yùn)行條件確定系 統(tǒng)的運(yùn)行狀態(tài)，如各母線 上的電壓、網(wǎng)絡(luò)中的功率分布以及功率損耗等 。在首先討論了電力系統(tǒng)節(jié)點(diǎn)的分類(lèi)以及潮流計(jì)算結(jié)果的約束條件后，具體介紹了常用于潮流計(jì)算的主要方法－牛頓－拉夫遜法。 在論文寫(xiě)作過(guò)程中，得到了吳茜瓊老師的親切關(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)+... xx 理 工學(xué)院畢業(yè)設(shè)計(jì)論文 29 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。 end A(n)=sum。 if NodeData(i,1)==1 xx 理 工學(xué)院畢業(yè)設(shè)計(jì)論文 35 NodeData(i,4)=NodeData(i,4)+X(m)。 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)))。 end end end m=m+1。 if NodeData(i,1)==1 YMatrix(m,n)=YMatrix(m,n)+NodeData(i,4)*A(m+1)。 end A(n)=sum。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。Select Node datafile39。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStruxx 理 工學(xué)院畢業(yè)設(shè)計(jì)論文 28 ctureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。 Temporary1=max(TopoStructureAndBranchPara(:,1))。謝謝 ! xx 理 工學(xué)院畢業(yè)設(shè)計(jì)論文 23 參考文獻(xiàn) [1] 楊以涵 .電力系統(tǒng)基礎(chǔ) .水利電力出版社， 1986. [2] 于永源 . 電力系統(tǒng)分析 .中國(guó)電力出版社， 1996 年 . [3] 諸俊偉等 .電力系統(tǒng)分析 [M].北京：中國(guó)電力出版社， 1995. [4] 楊少兵 ,駱平 . 電力系統(tǒng)分析的教學(xué)軟件開(kāi)發(fā)，電力系統(tǒng)潮流分析 .華北電力技術(shù) ， 2021(10). [5] 韓禎 祥，吳國(guó)炎等，電力系統(tǒng)分析，浙江大學(xué)出版社， 1993 年 . [6] 何仰贊等 .電力系統(tǒng)分析上冊(cè) [M]．武漢：華中理工大學(xué)出版社 . [7] 李久勝等編 .電氣自動(dòng)化英語(yǔ) .哈爾濱：哈爾濱工業(yè)大學(xué)出版社，1999. [8] 尹克寧 .電力工程 .中國(guó)水利電力出版社， 1989. [9] 韋剛等 .電力系統(tǒng)分析要點(diǎn)與習(xí)題 .中國(guó)電力出版社 ， 2021. [10] 陳珩 .電力系統(tǒng)穩(wěn)態(tài)分析 . 水利電力出版社 ， 1994. [11] 陸敏政 .電力工程 .北京：中國(guó)電力出版社， 1996. [12] Steven Holzner(美 )著 . 詳實(shí) 翻譯組譯 . Visual Basic 6 技術(shù)內(nèi)幕 . 機(jī)械工業(yè)出版社 ， 2021. [13] 諸俊偉等 .電力系統(tǒng)分析 [M].北京：中國(guó)電力出版社， 1995. [14] 張伯明 ,陳壽孫 .高等電力網(wǎng)絡(luò)分析 [M].北京：清華大學(xué)出版社1996. [15] 周全仁等 .電網(wǎng)計(jì)算與程序設(shè)計(jì) [M].長(zhǎng)沙：湖南科學(xué)技術(shù)出版社 ，1983. [16] 紀(jì)建偉等 .電力系統(tǒng)分析 .中國(guó)水利電力出版社， 1989. xx 理 工學(xué)院畢業(yè)設(shè)計(jì)論文 24 附 錄 輸 入 原 始 數(shù) 據(jù)形 成 節(jié) 點(diǎn) 導(dǎo) 納 矩 陣設(shè) 電 壓 初 值設(shè) 迭 代 次 數(shù) k ＝ 0( 0 )e( 0 )f計(jì) 算 誤 差 向 量()kp? ()kQ?2( )kU?收 斂 否否求 雅 可 比 矩 陣 元 素K = k + 1修 正 節(jié) 點(diǎn) 電 壓解 修 正 方 程 ，求 解()ke?()kf?( 1 ) ( ) ( )k k kf f f?? ? ?( 1 ) ( ) ( )k k ke e e? ? ? ?K m a xk不 收 斂 停 機(jī)是是求 p u 節(jié) 點(diǎn) 無(wú) 功 功 率 ， 求 平衡 節(jié) 點(diǎn) 功 率停 機(jī)求 支 路 功 率 分 布 和 損 耗停 機(jī)否 牛頓－拉 夫遜法計(jì)算潮流程序框圖 xx 理 工學(xué)院畢業(yè)設(shè)計(jì)論文 25 計(jì)算程序 %function OutPut=PowerFlowCalculation(handles) %is a subroutine of PowerSystemCalculation function OutPut=PowerFlow(handles) %the following program is open a data file and get the Number of % Node and Branch [fname,pname] = uigetfile(39。因此選擇開(kāi)發(fā)了潮流計(jì)算程序界面 ,VB 作為一種可視化的編程工具 ,其編程簡(jiǎn)單、界面友好 ,但它的計(jì)算能力相對(duì)較弱 ,因此將 MATLAB 與 VB 的優(yōu)勢(shì)結(jié)合起來(lái) ,讓它們?nèi)￠L(zhǎng)補(bǔ)短 ,進(jìn)行無(wú)縫鏈接 ,可使編程更專業(yè)、更靈活。這給以后的分析計(jì)算帶來(lái)了很大的方便，它 有利于節(jié)省內(nèi)存、提高 計(jì)算速度以及改善收斂等。 (5)為電力系統(tǒng)的規(guī)劃和擴(kuò)建提供依據(jù)。 xx 理工學(xué)院畢業(yè)設(shè)計(jì)（論文） 16 圖 (41)牛頓－拉夫遜方法的幾何意義 xx 理工學(xué)院畢業(yè)設(shè)計(jì)（論文） 17 第 5 章 計(jì)算實(shí)例 算例 圖 1 為一五結(jié)點(diǎn)系統(tǒng)，各支路參數(shù)均為標(biāo)么值。如果第 k 次迭代中得到 ()kx ，則過(guò) ()( ) ( ), ( )kkkfyxx???????點(diǎn)作一切線，此切線同 x 軸的交點(diǎn)便確定了下一個(gè)近似值 ( 1)kx? 。 ( 0 )()()f xx f x? ? ? (45) 用所求的 (0)x? 去修正近似解，變得 ( 0 )( 1 ) ( 0 ) ( 0 ) ( 0 )39。 39。 （ 3）在網(wǎng)絡(luò)的原有節(jié)點(diǎn) i、 j 之間切除一條導(dǎo)納為ijy的支路， （見(jiàn)圖 3－ 1，（ c） )） ，其相當(dāng)在 i、 j 之間增加一條導(dǎo)納為－ijy的支路，因此與 i、 j 有關(guān)的元素應(yīng)作以下修改： 1）節(jié)點(diǎn) i、 j 的自導(dǎo)納增量iiY?＝j(luò)jY?＝－ijy； 2）節(jié)點(diǎn) i、 j 之間的互導(dǎo) 納增量ijY?＝j(luò)iY?＝ijy； （ 4）原有網(wǎng)絡(luò)節(jié)點(diǎn) i、 j 之間的導(dǎo)納由 ijy 變成 39。在復(fù)雜電力網(wǎng)中，這中情況較多，從而使矩陣中出現(xiàn)了大量的零元素、節(jié)點(diǎn)導(dǎo)納矩陣稱為稀疏矩陣。 在電工原理課中，已講過(guò)用節(jié)點(diǎn)導(dǎo)納矩陣表示的節(jié)點(diǎn)電壓方程為 nn nY UI?? ???? ??? ???? ???? ?? （ 31） 對(duì)于 n 個(gè)節(jié)點(diǎn)的網(wǎng)絡(luò)，它可以展開(kāi)為 111 1 1 2 12 1 2 2 221 2 3............2. . .. .. . ...... .. .. . .. .......nnn n nn nUIY Y YY Y YI UY Y YI U????????????? ?????????????????? ?? ?????????? ?????????? ?????????????? ?? (32) 式（ 31）中的是節(jié)點(diǎn)注入電流的列向量。因此，這一約束條件主要是對(duì) PQ 節(jié)點(diǎn)而言。由于平衡節(jié)點(diǎn)的電壓已給定，只需要計(jì)算其余（ n－ 1）個(gè)節(jié)點(diǎn)的電壓。 3．平衡節(jié)點(diǎn) 在潮流分布算出以前，網(wǎng)絡(luò)中的功率損失是未知的，因此，網(wǎng)絡(luò)中至少有一個(gè)節(jié)點(diǎn)的有功功律 P 是