【正文】 s method in optimization. As the algorithm based on the Newton’s method in optimization has considerable putational cost and power control cannot be realized in all nodes, the algorithm based on the bination of the NewtonRaphson methods and the Newton’s method in optimization is offered to be utilized for calculating speed, enhancing the power flow calculation. II. THEORETICAL BACKGROUND A． Steadystate equations 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 29 The system of steadystate equations, in general, can be expressed as follows: ( , ) = 0W W X Y where W is the vector of parameters given for power flow calculation. In power flow calculation, real and reactive powers are set in each bus except for the slack bus. Ingeneration buses, the modulus of voltage can be fixed. W(X,Y) is the nonlinear vector function of steadystate equations. Variables Y define the quasiconstant parameters associated with an equivalent circuit of an electrical work. X is a required state vector, it defines steady state of EPS. The dimension of the state vector coincides with the number of nonlinear equations of the system (1). There are various known forms of notation of the steadystate equations. Normally, they are nodalvoltage equations in the form of power balance or in the form of current balance. Complex quantities in these equations can be presented in polar or rectangular coordinates, which leads to a sufficiently large variety forms of the steadystate equations notation. There are variable methods of a nonlin。 最后，衷心地感謝在百忙之中評(píng)閱論文和參加答辯的各位專家、教授！ 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 27 參考文獻(xiàn) 1 陳燕萍，王啟，趙彩虹等 ． 2021． 基于改進(jìn)前推回代法的輻射配電網(wǎng)潮流計(jì)算方法 [J]． 南京師范大學(xué)學(xué)報(bào)， 2021． 8(1)： 24~29． 2 劉健 ,畢鵬翔 ,董海鵬 ． 復(fù)雜配電網(wǎng)簡(jiǎn)化分析與優(yōu)化 [M]． 北京 ： 中國(guó)電力出版社 ， 2021． 3 王守相，王成山．現(xiàn)代配電系統(tǒng)分析 [M]．北京：高等教育出 版社 ． 2021． 4 王丹，常寶立．一種用于配網(wǎng)潮流計(jì)算的節(jié)點(diǎn)編號(hào)新方法 [J]．電力系統(tǒng)及其自動(dòng)化學(xué)報(bào) ． 2021， 15(1)： 2226． 5 何麗 ． 地區(qū)配電網(wǎng)潮流計(jì)算方法的研究 [D] ． 哈爾濱 ： 東北農(nóng)業(yè)大學(xué) 2021． 6 李光琦 .電力系統(tǒng)暫態(tài)分析 [M]． 北京 ： 中國(guó)電力出版社 1998． 7 宋文南，李樹鴻，張堯 ． 1990． 電力系統(tǒng)潮流計(jì)算 ． 天津 ： 天津大學(xué)出版社 ． 8 劉耀年，豈小梅，李國(guó)鵬，等．基于回路阻抗法的配電網(wǎng)潮流計(jì)算[J]．繼電器， 2021， 32(8)： 810． 9 汪宇霆，張焰，張益波．基于改進(jìn)回路電流法的配電網(wǎng)潮流 通用算法[J]．電力系統(tǒng)保護(hù)與控制， 2021， 38(20)： 5761， 68． 10 索南加樂，李懷強(qiáng)，羅云照，等．一種新的配網(wǎng)潮流常 Jacobian牛頓算法 [J]．西安交通大學(xué)學(xué)報(bào) ， 2021， 36(12)： 12221226． 11 汪芳宗，葉婧，李燕山．一種新的少環(huán)配電網(wǎng)潮流計(jì)算方法 [J]．電網(wǎng)技術(shù) ， 2021， 32(S1)： 4750． 12 孫健，江道灼．基于牛頓法的配電網(wǎng)絡(luò) Zbus潮流計(jì)算方法 [J]．電網(wǎng)技術(shù)， 2021， 28(15)： 4044． 13 顧晨，樂秀璠，張曉明．基于改進(jìn)前推回代法的弱環(huán)配電網(wǎng)三相潮流計(jì)算 [J]．電力系統(tǒng)保護(hù)與控制 ， 2021， 38(19)： 160164． 14 王明俊，于爾鏗，劉廣 .配電系統(tǒng)自動(dòng)化及其發(fā)展 [M]． 北京：中國(guó)電力出版社 1998． 15 何仰贊，溫增銀 ． 電力系統(tǒng)分析 (上冊(cè) )(第三版 )[M]． 武漢 ： 華中科技大學(xué)出版社 ， 2021． 1． 16 王明崗，畢睿華 ． 基于前推后代法解決環(huán)網(wǎng)的配電網(wǎng)潮流計(jì)算 [J]． 南京工程學(xué)院學(xué)報(bào)， 2021， 3(1)： 36~43． 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 28 附錄 A 英文文獻(xiàn) Power Flow Calculation by Combination of NewtonRaphson Method and Newton’s Method in Optimization. Andrey Pazderin, Sergey Yuferev URAL STATE TECHNICAL UNIVERSITY – UPI Email: , AbstractIn this paper, the application of the Newton’s method in optimization for power flow calculation is considered. Convergence conditions of the suggested method using an example of a threemachine system are investigated. It is shown, that the method allows to calculate nonexistent state points and automatically pulls them onto the boundary of power flow existence domain. A bined method which is posed of NewtonRaphson method and Newton’s method in optimization is presented in the paper. Index Terms—Newton method, Hessian matrix, convergence of numerical methods, steady state stability Ⅰ .INTRODUCTION The solution of the power flow problem is the basis on which other problems of managing the operation and development of electrical power systems (EPS) are solved. The plexity of the problem of power flow calculation is attributed to nonlinearity of steadystate equations system and its high dimensionality, which involves iterative methods. The basic problem of the power flow calculation is that of the solution feasibility and iterative process convergence [1]. The desire to find a solution which would be on the boundary of the existence domain when the given nodal capacities are outside the existence domain of the solution, and it is required to pull the state point back onto the feasibility boundary, motivates to develop methods and algorithms for power flow calculation, providing reliable convergence to the solution. The algorithm for the power flow calculation based on the Newton39。 同時(shí)我要感謝研究生 師兄的指點(diǎn)和幫助，在與師兄的交流中，我學(xué)到了很多知識(shí)，獲得了不少啟發(fā)。她 淵博的學(xué)識(shí)水平、嚴(yán)謹(jǐn)踏實(shí)的工作作風(fēng)、兢兢業(yè)業(yè)的科研精神和處處為人師表的作風(fēng)使我獲益匪淺。 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 26 致謝 本課題的選題、課題研究及論文撰寫工作是在 付 老師細(xì)心指導(dǎo)下完成的。 充分利用配電網(wǎng)絡(luò)的結(jié)構(gòu)特點(diǎn) ， 對(duì) 節(jié)點(diǎn)進(jìn)行編號(hào) ,直接采用支路參數(shù)進(jìn)行計(jì)算 ,不需要形成導(dǎo)納矩陣和進(jìn)行三角矩陣 ， 簡(jiǎn)單處理代數(shù)方程就可以計(jì)算節(jié)點(diǎn)電壓和功率分布 。 a 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 21 b 圖 41 程序編輯 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 22 a 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 23 b 圖 42 程序運(yùn)行 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 24 運(yùn)算結(jié)果如表 43 所示，各節(jié)點(diǎn)電壓幅值逐漸降低，有功功率，無功功率也逐漸降低，符合邏輯。 圖 41 8 節(jié)點(diǎn)主干饋線配電網(wǎng) 表 41 支路參數(shù) 支路序號(hào) 電阻 R/Ω 電抗 X/Ω 1 2 3 4 5 6 7 表 42 節(jié)點(diǎn)負(fù)荷功率 節(jié)點(diǎn)序號(hào) 有功負(fù)荷 /kW 無功負(fù)荷 /Kvar 1 2 3 16 4 5 16 6 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 20 7 8 程序 運(yùn)行 下面為所編程序在軟件 Visual C++ 中運(yùn)行截圖。 哈爾濱理工大學(xué)學(xué)士學(xué)位論文 19 第 4章 配電網(wǎng)潮流計(jì)算程序仿真 算例分析 利用上面提出的潮流計(jì)算方法取 8 個(gè)節(jié)點(diǎn)時(shí)，如 圖 41 所示配電網(wǎng)絡(luò)盡行潮流計(jì)算網(wǎng)絡(luò) 。 cout節(jié)點(diǎn) i+1的電壓相角偏差 ai+1=a[i]\n。 cout節(jié)點(diǎn) i+1的電壓幅值 Vi+1=V[i+1]\n。 deltaV[i+1]=(P[i+1]*X[i]+Q[i+1]*R[i])/V[i+1]。 cout節(jié)點(diǎn) i的注入無功功率 Qi=Q[i] *1000\n。 Q[i]=Q[i+1]+LQ[i]+(P[i+1]*P[i+1]+Q[i+1]*Q[i+1])/(Vn*Vn)*X[i]。i=1。 cout節(jié)點(diǎn) n的注入無功功率 Qn=Q[n] *1000\n。 cinX[i]。 cinR[i]。i=n1。 LQ[1]=0。 Q[n]=LQ[n]。 LQ[i]=LQ[i]/1000。 cout請(qǐng)輸入第 i個(gè)節(jié)點(diǎn)的無功負(fù)荷功率 LQi=。 cinLP[i]。i=n。 cinVn。 cinV[1]。 a=(double *)malloc((n+1)*sizeof(double))。 V=(double *)malloc((n+1)*sizeof(double))。 LP=(double *)malloc((n+1)*sizeof(double))。 R=(double *)malloc