【正文】 [J].電網(wǎng)技術 ,2021(6):10331038. [4] 陸敏政 .電力系統(tǒng)分析 [M].北京 。在進行課題設計的過程中，加深了我對潮流計算的認識，尤其是對牛頓拉夫遜潮流計算的求解思路有了比較透徹的理解。 DDS(i)=DS(i)*SB。,39。39。 if B1(i,6)==0 Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。各節(jié)點的功率 S 為 (節(jié)點號從小到大排列 )： 39。各節(jié)點的電 壓大小 V 為 (節(jié)點號從小到大排列 )： 39。沒有達到精度要求的個數(shù)： 39。%消去運算 end J(k3,k)=0。J(p,N)=DV。J(p,N)=DV。 J(m,q)=X1。 J(m,q)=X1。% 節(jié)點功率 Q 計算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej) %求 P39。 %i 節(jié)點注入功率 SGSL B(i,i)=B(i,i)+B2(i,5)。q=B1(i,2)。 BB1=B1。 200 0 229 229 0 3]。 2 3 + 0 1。)1(1239。0444 ）????????????M V Aj ) ??（ （ 9）④ 120Km 與 100Km 線路交點 3 末端功率損耗 M V AjSSSS LDT )] 8 6 ()560 0 9 7 4 1 [(0333 ????????????M VAj )( ?? （ 10） 1 234T 1T 2T 3T 4110/S 0 = 0 .0 1 8 5 + j0 .1 4 4S 0 = 0 .0 1 8 5 + j0 .1 4 4G 1 G 2S 0 = 0 .0 1 1 6 + j0 .1 7 6 S 0 = 0 .0 1 1 6 + j0 .1 7 6LD 4 = 1 0 .2 + j6 .2LD 3 = 5 .6 + j4 .4S 0 = 0 .0 0 9 7 6 + j0 .0 8 8R T1 +jX T1 R T2 +jX T2R T3 +jX T3R T4 +jX T4R T4 +jX T4R 1 +jX 1 R 3 +jX 3R 2 +jX 2 R4 +jX 4jB 1 /2jB 1 /2 jB 3 /2jB 3 /2jB 2 /2jB 2 /2 jB 4 /2 jB 4 /2 9 ⑤ 間 100Km 線路損耗 M V AjjUBjXRUSSSM V AjUBjSSM V AjjSSNNLDLDNLDLDLD)()(2)()(2)()(1771772222221414)1(122141441439。 MATLAB 和 Mathematica、 Maple 并稱為三大數(shù)學軟件。 運行方式管理中，潮流是確定電網(wǎng)運行方式的基本出發(fā)點；在規(guī)劃領域，需要進行潮流分析驗證規(guī)劃方案的合理性；在實時運行環(huán)境，調(diào)度員潮流提供了多個在預想操作情況下電網(wǎng)的潮流 分布以及校驗運行可靠性。 2 電力系統(tǒng)潮流計算概述 電力系統(tǒng)簡介 電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本電氣計算。 設計內(nèi)容 ........................................................................................... 錯誤 !未定義書簽。其目的是校驗系統(tǒng)是否能安全運行，即是否有過負荷的元件或電壓過低的母線等。待求的運行狀態(tài)參量包括電網(wǎng)各母線節(jié)點的電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡的功率損耗等。這是數(shù)學中解決非線性方程式的典型方法，有較好的收斂性。2﹪） /， ΔPs=110kW， ΔP0=，U0﹪ =， I0﹪ =； 變壓器 T3： SFL18000/110，（ 110177。jjUBjXRUSSSM V AjUBjSSM V AjjSSNNLDLDNLDLDLD????????????????? M V AjM V Aj )8 6 7 8 (1 8 5 ??? （ 13） ⑧ 間 100Km 線路損耗 )()(2)()(2)()(221222122333222323)2(223232342339。 n=8。 0 50+ 220 0 0 2。 if ym~=0 %SB=input(39。 B2(i,3)=B2(i,3)./UB。)。a=a+1。 % dP/df=dQ/de X3=X2。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 % N=N0+1 即 N=2*n+1 擴展列 △ P、△ Q for k2=k1:N1 % 擴展列 △ P、 △ Q J(k,k2)=J(k,k2)./J(k,k)。 %不滿足要求的節(jié)點數(shù)加 1 end end ICT2(a)=IT2。 disp(E)。各節(jié)點的電壓相角 sida 為 (節(jié)點號從小到大排列 )： 39。39。,39。 else Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 end disp(39。 %disp(DDS(i))。利用牛頓 拉夫遜法、采用 Matlab 語言進行潮流計算，提高了準確性，使用方便，已得到成功應用，而它在我們以后的專業(yè)課程的學習中應用也是很廣泛的，掌握好了 MATLAB，這為以后的專業(yè)課程的學習將有很大的幫助。 這次課程設計，雖然沒有想象中的那么簡單，花費了很多的時間來不斷地實驗，但是最后終于完成了預定的目的，收獲頗豐，也很好的鍛煉了自己運用知識的能力。39。)。 end disp(Sj(q,p))。)=39。 disp(39。 disp(sida)。 disp(EE)。 end %用高斯消去法解 w=J*V disp(39。 % 對角元規(guī)格化 if k~=3 % 不是第三行 %============================================================ k4=k1。% dP/df X5=2*e(i)。 % dP/df X5=0。% dQ/df p=2*i1。 % X1=dP/de X4=dQ/df p=2*i1。D(i)=0。B=imag(Y)。 B2(i,5)=B2(i,5)./SB。)。 0 40+ 220 0 0 2。 isb=1。39。2﹪） /， ΔPs=62kW， ΔP0=，Us﹪ =， I0﹪ =； 導線型號均為 LGJ150，參數(shù) r0=， x0=， b0=106S/km。自從 60 年代中期，牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算，即節(jié)點電壓和功率 分布，用以檢查系統(tǒng)各元件是否過負荷。但是一般靜態(tài)安全分析需要校驗的狀態(tài)數(shù)非常多，用嚴格的潮流計算來分析這些狀態(tài)往往計算量過大，因此不得不尋求一些特殊的算法以滿足要求。 設計目的 ........................................................................................... 錯誤 !未定義書簽。自從 20 世紀 60 年代中期利用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、速度方面都超過了阻抗法，成為直到目前仍在廣泛采用的優(yōu)秀方法。潮流計算結果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安 4 全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。 3 Matlab 概述 Matlab 簡介 MATLAB 是由美國 mathworks 公司發(fā)布的主要面對科學計算、可視化以及交互式程序設計的高科技計算環(huán)境。 繪制等效電路 功率分布計算 圖 2 環(huán)形網(wǎng)絡的等效圖 （ 1） 各元件功率損耗 ① 兩臺 T4 變壓器并聯(lián)損耗 M V