【正文】 (417) 式中的 X 和 X? 分別是由 n 個(gè)變量和修正量組成的 n 維列向量； F(X)是由 n個(gè)多元函數(shù)組成的 n 維列項(xiàng)量； J 是 n 階方陣，稱為雅可比矩陣，它的第 i、j 個(gè)元素 iij ifJ x??? 是第 n 個(gè)函數(shù) 12( , ,..., , )nif x x x對(duì)第 j 個(gè)變量jx的偏導(dǎo)數(shù)；上角標(biāo) (k)表示 J 陣的每一個(gè)元素都在點(diǎn), ,()( ) ( )( .. ., )12i kkk nf xxx處取值。 1? 和 2? 為預(yù)先給定的小正數(shù)。由于節(jié)點(diǎn)電壓可以采用不同的坐標(biāo)系表示，牛頓－拉夫遜法潮流計(jì)算也將相應(yīng)的采用不同的計(jì)算公式。假定結(jié)點(diǎn) 3 為 PQ節(jié)點(diǎn)，結(jié)點(diǎn) 4 為 PV 節(jié)點(diǎn)、結(jié)點(diǎn) 5 為平衡結(jié)點(diǎn)，試分別用 直角坐標(biāo)和極坐標(biāo)牛頓－拉夫遜法 計(jì)算其潮流。 給定： S1=?? S2=?? S3=?? P4= |V1（ 0） |=|V2（ 0） |=|V3（ 0） |=|V4（ 0） |= |V4|=|V5|= (0)3(0)2(0)1 ??? 5(0)4 ?? 1 . 0 5 : 14j 0 . 0 1 51 : 1 . 0 52J 0 . 2 50 . 0 8 + j 0 . 3 0j 0 . 2 5J 0 . 2 50 . 0 4 + j 0 . 2 50 . 1 + j 0 . 3 5J 0 . 2 51j 0 . 0 335 圖 1 5 節(jié)點(diǎn)系統(tǒng) 節(jié)點(diǎn)導(dǎo)納的形成 根據(jù)節(jié)點(diǎn)導(dǎo)納矩陣的定義，可求的節(jié)點(diǎn)導(dǎo)納矩陣各元素，即 05( 0 )4( 0 )3( 0 )2( 0 )1 ????? fffffxx 理工學(xué)院畢業(yè)設(shè)計(jì)（論文） 18 11 1 0 1 2 1 3 110 . 2 5 0 . 0 4 0 . 2 5 0 . 1 0 . 3 5j jjyyyY ? ? ? ? ? ??? =++ = 與節(jié)點(diǎn)１有關(guān)的互導(dǎo)納為 12 21 12yYY? ? ?=+ 31 13 13yYY? ? ?=0754717+ 支路 24 為變壓器支路，可以求出節(jié)點(diǎn) 2 的自導(dǎo)納為 222 2 0 1 2 2 3 4 5 /y y y y kY ? ? ? ? ＝ ＋ ++ ＝ 與節(jié)點(diǎn) 2 有關(guān)的互導(dǎo)納為 23 3224 42 4242 29 87 6 30 11 20 33/ 63 .49 20 64jjYY y kYY? ? ? ?? ? ? ? ? 用類似的方法可以求出導(dǎo)納矩陣的其他元素，最后可得到節(jié)點(diǎn)導(dǎo)納矩陣為 742 025 717 0 0 665 156 509 025 901 876 206 3 0 156 082 1 033 717 876 593 0 603 2 509 033 785 80j j jjj j jYjj j j??? ? ?? ? ?? ? ??? ? ?? ? ?? 206 3 0 666 7 00 0 603 2 0 333jjj??????????? 計(jì)算結(jié)果 xx 理工學(xué)院畢業(yè)設(shè)計(jì)（論文） 19 節(jié)點(diǎn)原始數(shù)據(jù)： scanf the 1th numbers: 2 4 0 scanf the 2th numbers: 4 3 scanf the 3th numbers: 4 5 scanf the 4th numbers: 5 3 0 scanf the 5th numbers: 3 1 0 計(jì)算結(jié)果 v=+ v=+ v=+ v=+ the balance node:+ s=+ s=+ s=+xx 理工學(xué)院畢業(yè)設(shè)計(jì)（論文） 20 結(jié) 論 電力系統(tǒng)潮流計(jì)算分布計(jì)算，是指電力系統(tǒng)在某一穩(wěn)定狀態(tài)的正常運(yùn)行方式下，電力網(wǎng)絡(luò)各節(jié)點(diǎn)的電壓和功率分布的計(jì)算。 (2)檢查電力系統(tǒng)各節(jié)點(diǎn)的電壓是否滿足電壓質(zhì)量的要求。 (4) 檢查電力系統(tǒng)各元件是否過(guò)負(fù)荷。 (6)為調(diào)整計(jì)算、經(jīng)濟(jì)運(yùn)行計(jì)算、短路計(jì)算和穩(wěn)定計(jì)算提供必要的數(shù)據(jù)。潮流計(jì)算的數(shù)學(xué)模型是以 節(jié)點(diǎn)的方程為基礎(chǔ)，推導(dǎo)出相應(yīng)的功率方程。在復(fù)雜電力網(wǎng)中，在各個(gè)節(jié)點(diǎn)中沒(méi)有直接相連的節(jié)點(diǎn)很多，從而使矩陣中有很多零元素、節(jié)點(diǎn)導(dǎo)納矩陣成為稀疏矩陣。因此，節(jié)點(diǎn)導(dǎo)納矩陣是一個(gè)對(duì)稱、稀疏且具有對(duì)角線優(yōu)勢(shì)的方陣。功率方程是非線性代數(shù) 方程組，必須采用數(shù)值求解的方法進(jìn)行計(jì)算。在實(shí)例計(jì)算這一章節(jié)中的驗(yàn)證，從計(jì)算量來(lái)看，計(jì)算量不是特別的龐大，計(jì)算結(jié)果可以迅速的收斂，可以快速準(zhǔn)確的計(jì)算出各個(gè)節(jié)點(diǎn)的待求量。從計(jì)算的速度來(lái)看，速度比較快。 另外， MATLAB 的現(xiàn)有界面對(duì)用戶而言 不直接、不方便。 xx 理工學(xué)院畢業(yè)設(shè)計(jì)論文 22 謝 辭 經(jīng)過(guò)幾個(gè)月的忙碌和工作，本次畢業(yè)論文設(shè)計(jì)已經(jīng)接近尾聲，在這次的畢業(yè)設(shè)計(jì)中，我們做的是牛頓－拉夫遜法潮流計(jì)算，課題類型為工程設(shè)計(jì)類的題目，由于經(jīng)驗(yàn)的匱乏，難免有許多考慮不周全的地方，如果沒(méi)有導(dǎo)師的督促指導(dǎo)，以及一起工作的同學(xué)們的支持，想要完成這個(gè)設(shè)計(jì)是難以想象的。她嚴(yán)肅的科學(xué)態(tài)度，嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神，精益求精的工作作風(fēng)，深深地感染和激勵(lì)著我。多少個(gè)日日夜夜，吳老師不僅在學(xué)業(yè)上給我們以精心指導(dǎo)，除了敬佩劉孝書(shū)老師的專業(yè)水平外，他的治學(xué)嚴(yán)謹(jǐn)和科學(xué)研究的精神也是我永遠(yuǎn)學(xué)習(xí)的榜樣，并將積極影響我今后的學(xué)習(xí)和工作。 在論文即將完成之際，我的心情無(wú)法平靜，從開(kāi)始進(jìn)入課題到論文的順利完成，有多少可敬的師長(zhǎng)、同學(xué)、朋友給了 我無(wú)言的幫助，在這里請(qǐng)接受我誠(chéng)摯的謝意 ! 最后我還要感謝自動(dòng)化系和我的母校 — 洛陽(yáng)理工學(xué)院這三年來(lái)對(duì)我的栽培。*.dat39。Select Topo datafile39。 Branchname=strcat(pname,fname)。 [NumberOfBranch,NumberOfPara]=size(TopoStructureAndBranchPara)。 Temporary2=max(TopoStructureAndBranchPara(:,2))。 else NumberOfNode=Temporary2。 end end 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)+... xx 理工學(xué)院畢業(yè)設(shè)計(jì)論文 27 j*TopoStructureAndBranchPara(CircleNumber,4))) 。 NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))=... 1/... (TopoStructureAndBranchPara(CircleNumber,5)*(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))) 。 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)。 NAM(TopoStructureAndBranchPara( CircleNumber,1),TopoStructureAndBranchPara( CircleNumber,2))=... 1/(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))。 end end % [fname,pname] = uigetfile(39。,39。)。 NodeData = csvread(Nodename)。 datamax=0。 xx 理工學(xué)院畢業(yè)設(shè)計(jì)論文 30 n=1。 sum1=0。 sum1=sum1+NodeData(w,4)*(real(NAM(i,w))*sind(NodeData(i,5)NodeData(w,5))imag(NAM(i,w))*cosd(NodeData(i,5)NodeData(w,5)))。 n=n+1。 n=n+1。 datamax=max(datamax,B(m))。 xx 理工學(xué)院畢業(yè)設(shè)計(jì)論文 31 B(m)=NodeData(i,3)NodeData(i,4)*sum1。 m=m+1。 for r=1:NumberOfNode sum=sum+NodeData(r,4)*(real(NAM(i,r))*cosd(NodeData(i,5)NodeData(r,5))+imag(NAM(i,r))*sind(NodeData(i,5)NodeData(r,5)))。 n=n+1。 datamax=max(datamax,B(m))。 end end while datamaxfloatmin m=1。 for k=1:NumberOfNode1 if i==k YMatrix(m,n)=NodeData(i,4)^2*imag(NAM(i,i))。 end n=n+1。 n=n+1。 n=n+1。 n=n+1。 if NodeData(i,1)==1 n=1。 n=n+1。 n=n+1。 n=n+1。 n=n+1。 end end X=inv(YMatrix)*B39。 m=1。 m=m+1。 m=m+1。 n=1。 sum1=0。 sum1=sum1+NodeData(v,4)*(real(NAM(i,v))*sind(NodeData(i,5)NodeData(v,5))imag(NAM(i,v))*cosd(NodeData(i,5)NodeData(v,5)))。 n=n+1。