【正文】 dP/df X5=2*e(i)。 elseif j1==iamp。J(m,N)=DP。J(p,q)=X5。 % dP/df X5=0。 % PV 節(jié)點電壓誤差 for j1=1:n if j1~=isbamp。J(m,N)=DP。%擴展列△ Q m=p+1。% dQ/df p=2*i1。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 J(p,q)=X4。m=p+1。 % X1=dP/de X4=dQ/df p=2*i1。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 %電壓模平方 %========= 以下針對非 PV節(jié)點來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV節(jié)點 DP=P(i)P1。%節(jié)點功率 P 計算 eiΣ (Gij*ejBij*fj)+fiΣ (Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。D(i)=0。a=0。 ICT1=0。 %PV 節(jié)點電壓給定模值 end for i=1:n %給定各節(jié)點注入功率 S(i)=B2(i,1)B2(i,2)。B=imag(Y)。 %對角元 1 側 end %求導納矩陣 disp(39。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 disp(B2) end % % % for i=1:nl %支路數 if B1(i,6)==0 %左 節(jié)點處于低壓側 p=B1(i,1)。 B2(i,5)=B2(i,5)./SB。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。 B1(i,4)=B1(i,4)./YB。 YB=SB./UB./UB。)。您輸入的參數是標么值？（若不是則輸入一個不為零的數值） 39。 %－－－－－－－修改部分－－－－－－－－－－－－ ym=0。f=zeros(1,n)。請輸入各節(jié)點參數形成的矩陣： B2=39。 0 + 1 0 0 2。 0 0 1 0 0 2。)。 5 9 + 0 1。 3 7 + 0 1。 1 4 + 1 0。%input(39。%input(39。%input(39。%input(39。在潮流計算的大量數值計算過程中更顯示出其優(yōu)勢，是我們應該掌握的一門基本技術。潮流計算是實現電力系統(tǒng)安全經濟發(fā)展的必要 手段和重要環(huán)節(jié)。 0 + 1 0 0 2。 0 0 1 0 0 2。 6 10 + 0 1]。 3 7 + 0 1。 1 4 + 1 0。 6 10 + 0 1]。 3 7 + 0 1。 1 4 + 1 0。 5 9 + 0 1。 3 7 + 0 1。 調節(jié)方 法 發(fā)電廠1 機端電壓 發(fā)電廠2 機端電壓 變電所1 變比 變電所2 變比 變電所3 變比 變電所4 變比 未調整 1 1 1 1 調整后 調節(jié)后的輸出結果： B1=[1 3 + 1 0。 5 9 + 0 1。 3 7 + 0 1。 調整前后各節(jié)點的 電壓 標幺值如下： 各支路的功率損耗如下表： 功率損耗 未調整 調整后 （ 1,3）支路 （ 1,4）支路 (2,6) 支路 (3,4) 支路 (3,7) 支路 + + (4,5) 支路 (4,8) 支路 + + (5,6) 支路 (5,9) 支路 + + (6,10) 支路 + + 各節(jié)點的功率 S 為 (節(jié)點號從小到大排列 )： + + + 調整后的具體潮流分布 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,4)=+ S(4,1)= S(2,6)=110+ S(6,2)= S(3,4)= S(4,3)= S(3,7)=+ S(7,3)=5031i 節(jié)點電壓 1 2 3 4 5 6 7 8 9 10 調整前 0,9973 0,9998 調整后 S(4,5)=+ S(5,4)= S(4,8)=+ S(8,4)= S(5,6)=+ S(6,5)= S(5,9)=+ S(9,5)=5031i S(6,10)=+ S(10,6)= 圖像如下： 斷開 6 節(jié)點之間一回線時的潮流計算 B1 陣變?yōu)椋? B1=[1 3 + 1 0。 5 6 + 1 0。 3 4 + 1 0。 電壓調節(jié)方法：通過改變發(fā)電機極端電壓、變壓器變比或增加無功功率補償來調整電壓。 5 6 + 1 0。 3 4 + 1 0。 0 + 1 0 0 2]。 0 0 1 0 0 2。 B2=[0 0 0 1。 4 8 + 0 1。 2 5 + 1 0。 B2矩陣不變。 4 8 + 0 1。 2 5 + 1 0。 0 + 1 0 0 2。 0 0 1 0 0 2。 5 9 + 0 1。 3 7 + 0 1。 調節(jié)方 法 發(fā)電廠1 機端電壓 發(fā)電廠2 機端電壓 變電所1 變比 變電所2 變 比 變電所3 變比 變電所4 變比 未調整 1 1 1 1 調整后 調整后的輸出結果： B1=[1 3 + 1 0。 5 9 + 0 1。 3 7 + 0 1。 調整前后各節(jié)點 電壓 的標幺值如下： 各支路的功率損耗如下表： 功率損耗 未調整 調整后 （ 1,4）支路 （ 2,5）支路 (2,6) 支路 (3,4) 支路 (3,7) 支路 + + (4,5) 支路 (4,8) 支路 + + (5,6) 支路 (5,9) 支路 + + (6,10) 支路 + + 各節(jié)點的功率 S 為 (節(jié)點號從小到大排列 )： + + + + 0 + 調整后的具體潮流分布 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,4)=+ S(4,1)= 節(jié)點電壓 1 2 3 4 5 6 7 8 9 10 調整前 0,9716 調整后 S(2,5)=+ S(5,2)= S(2,6)=+ S(6,2)= S(3,4)= S(4,3)=+ S(3,7)=+ S(7,3)=5031i S(4,5)= S(5,4)=+ S(4,8)=+ S(8,4)= S(5,6)= S(6,5)= S(5,9)=+ S(9,5)=5031i S(6,10)=+ S(10,6)= 圖像如下： 斷開 4節(jié)點之間一回線時的潮流計算 B1矩陣變?yōu)椋? B1=[1 3 + 1 0。 5 6 + 1 0。 3 4 + 1 0。 0 + 1 0 0 2] 電壓調節(jié)方法：通過改變發(fā)電機極端電壓、變壓器變比或增加無功功率補償來調整電壓。 0 0 1 0 0 2。 B2=[0 0 0 1。 4 8 + 0 1。 2 6 + 1 0。 0 + 1 0 0 2。 0 0 1 0 0 2。 5 9 + 0 1。 3 7 + 0 1。 1 4 + 1 0。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 0 1 0 0 2。 5 9 + 0 1。 3 7 + 0 1。 1 4 + 1 0。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 0 1 0 0 2。 5 9 + 0 1。 3 7 + 0 1。 1 4 + 1 0。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 0 1 0 0 2。 6 10 + 0 1]。 4 5 + 1 0。 節(jié)點 1 2 3 4 5 6 7 8 9 10 電壓 2 5 + 1 0。 關鍵詞 牛頓 拉夫遜法（ NewtonRaphson） 潮流計算與 DDRT S 仿真 潮流計算結果分析與調控 利用 MATLAB 程序運行的結果 未調整時各節(jié)點電壓的標么值如下： 圖像如下： 由調壓范圍可知，變電所 4 的節(jié)點電壓均不滿足電壓范圍。 牛頓 — 拉夫遜 NewtonRaphson 法是數學上解非線性方程組的有效方法，有較好的收斂性。 4. 對不同的負荷變化，進行潮流的調節(jié)控制，并說明調節(jié)控制的方法，并列表表示調節(jié)控制的參數變化。（變電所低壓母線電壓 10KV 要求調整范圍在 之間；電壓 35KV要求調整范圍在 3536 之間） 5. 輪流斷開環(huán)網一回線，分析潮流的分布。 變電所資料： (一 ) 變電所 4 低壓母線的電壓等級分別為： 35KV 35KV 35KV 10KV (二 ) 變電所的負荷分別為： （ 1） 50MW 50MW 70MW 60MW （ 2） 60MW 70MW 50MW 45MW （ 3） 50MW 40MW 50MW 70MW 電廠一 電廠二 （ 4） 70MW 50MW 60MW 60MW （ 5） 70MW 60MW 50MW 50MW （ 6） 45MW 55MW 60MW 70MW （ 7） 70MW 45MW 50MW 60MW （ 8） 50MW 40MW 50MW 60MW （ 9） 40MW 40MW 50MW 60MW （ 10） 50MW 40MW 40MW 60MW (三 ) 每個變電所的功率因數均為 cosφ =； (四 ) 變電所 1 和變電所 2 分別配有兩臺容量為 75MVA 的變壓器，短路損耗