【正文】 J(m,q)=X4。J(m,q)=X3。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。X2=G(i,j)*e(i)+B(i,j)*f(i)。Df(p)=Q(i)Q1。c(i)=v(B1(i2))*a(i)j*B1(i4)*v(B1(i2))^2/2。J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)。ph=input(‘n請輸入平衡母線的節(jié)點(diǎn)號：ph=’)； B1=input（‘n請輸入支路信號：B1=’）。%% generator data% bus PgQgQmaxQmin Vg mbase status PmaxPmin =[ 4 500 0 999909999 100 1 600 0。牛頓拉夫遜法是求解非線性方程的迭代過程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。M(i,i)=M(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。N(i,j)=0。endY% Y為導(dǎo)納矩陣G=real(Y)。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導(dǎo)，得出對應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點(diǎn)電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點(diǎn)電壓不平衡量，節(jié)點(diǎn)電壓加上節(jié)點(diǎn)電壓不平衡量構(gòu)成節(jié)點(diǎn)電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計(jì)算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。不同領(lǐng)域、不同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進(jìn)行計(jì)算、分析及設(shè)計(jì)工作。③正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開機(jī)方式,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。（2）隨機(jī)抽取5條電壓調(diào)整后的支路功率進(jìn)行對比，如表9所示。（三）電網(wǎng)的節(jié)點(diǎn)設(shè)置與分類。最后把Matpower潮流計(jì)算的最終結(jié)果與Matlab編程和Simulink仿真兩種方法的結(jié)果做比較。（2）變電所資料：①變電所4低壓母線的電壓等級分別為：10KV 35KV 10KV 35KV。（5）系統(tǒng)節(jié)點(diǎn)電壓如表4所示，―，變電所低壓母線35KV在35―36KV之間，經(jīng)過折算發(fā)現(xiàn)節(jié)點(diǎn)9不在指定范圍內(nèi)，需要進(jìn)行電壓調(diào)整，電壓調(diào)整后折算發(fā)現(xiàn)9節(jié)點(diǎn)已滿足要求，且相角隨電壓幅值而變化。所以考留到經(jīng)及調(diào)度、電網(wǎng)規(guī)劃、電力系統(tǒng)可靠性分析。MATLAB與C語言和FORTRAN語言相比更容易被掌握。二 牛頓拉夫遜法潮流計(jì)算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。endendY=0。x2=ones(8,1)。M(i,j)=N(i,j)。endS(5)=P(5)+Q(5)*sqrt(1)。4 2 0 0 0 0 1 0 100 1 ?！尽俊尽縁unction tisco %這是一個電力系統(tǒng)潮流計(jì)算的程序 n=input（‘n請輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請輸入支路數(shù)：m=’)。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。p=p+1。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。X4=X1。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。J(no:)=[]。J(no:)=[]。Else if j=amp。X3=X2。Df(k)=DX(i)。else p=2*i1。resulte1=[A(:,1),real(v),imag(v),V,jd,real(S’)，imag(S’),real(B2(:1)),imag(B2(:1)),real(B2(:2)),imag(B2(:,2))]。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。)計(jì)算輸出結(jié)果如下：MATPOWER Version , 07Feb2011AC Power Flow(Newton)Newton39。2 1 200 100 0 0 1 1 0 100 1 。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))+P(i)。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。a=zeros(1,5)。y(3,4)=1/(+)。例題如下：用牛頓拉夫遜法計(jì)算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點(diǎn)，節(jié)點(diǎn)5電壓保持U=，其他四個節(jié)點(diǎn)分別為PQ節(jié)點(diǎn)，給定的注入功率如圖所示。MATLAB程序設(shè)計(jì)語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。為此，論文根據(jù)理想變壓器對功率和電壓的變換原理，提出了一種有效的電壓變換模型來處理變壓器支路，從而改善了潮流算法的收斂特性。（3），說明選擇1號節(jié)點(diǎn)為平衡節(jié)點(diǎn)是正確的。（二）系統(tǒng)圖。電力系統(tǒng)潮流計(jì)算的目的是計(jì)算系統(tǒng)在給定狀態(tài)下的節(jié)點(diǎn)電壓及功率分布，來檢查系統(tǒng)中各電壓是否滿足要求，系統(tǒng)中各元件是否過負(fù)荷以及功率分配的合理性等。并設(shè)平衡節(jié)點(diǎn)為節(jié)點(diǎn)1，PV節(jié)點(diǎn)為節(jié)點(diǎn)10，剩余節(jié)點(diǎn)為PQ節(jié)點(diǎn)。四、結(jié)語對給定系統(tǒng)通過Matpower進(jìn)行潮流計(jì)算后，其運(yùn)算結(jié)果從節(jié)點(diǎn)、支路等方面與Matlab編程與Simulink仿真的結(jié)果進(jìn)行了對比，得出系統(tǒng)運(yùn)行穩(wěn)定，且在符合系統(tǒng)要求的情況下，三種方法的潮流計(jì)算結(jié)果基本一致?？偨Y(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。原始數(shù)據(jù)輸入格式的設(shè)計(jì)，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。由課本總結(jié)基本步驟如下：1)形成節(jié)點(diǎn)導(dǎo)納矩陣Y；2)設(shè)各節(jié)點(diǎn)電壓的初值，如果是直角坐標(biāo)的話設(shè)電壓的實(shí)部e和虛部f；如果是極坐標(biāo)的話則設(shè)電壓的幅值U和相角a；3)將各個節(jié)點(diǎn)電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點(diǎn)電壓的變化量，即修正量； 5)計(jì)算各個節(jié)點(diǎn)電壓的新值，即修正后的值；6)利用新值從第（3）步開始進(jìn)入下一次迭代，直至達(dá)到精度退出循環(huán)； 7)計(jì)算平衡節(jié)點(diǎn)的功率和線路功率，輸出最后計(jì)算結(jié)果； ① 公式推導(dǎo)② 流程圖三matlab編程代碼clear。% 輸入原始節(jié)點(diǎn)的給定注入功率S(1)=+。L(i,j)=0。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。六 結(jié)論通過這個任務(wù)，自己在matlab編程，潮流計(jì)算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點(diǎn)還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，達(dá)到一個新的層次。]。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。t=0。printcut(result1,S,b,c,result2)。If ph=n Df(n