【正文】 J(p,q)=X6。J(m,N)=DP。 m=p+1。J(p,q)=X5。 p=2*i1。 % dP/df X5=0。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 % PV 節(jié)點電壓誤差 for j1=1:n if j1~=isbamp。 end end else %=============== 下面是針對 PV 節(jié)點來求取 Jacobi 矩陣的元素 =========== DP=P(i)P1。J(m,N)=DP。q=q+1。%擴展列△ Q m=p+1。J(p,q)=X3。% dQ/df p=2*i1。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 elseif j1==iamp。 J(p,q)=X4。J(m,N)=DP。m=p+1。 J(p,q)=X3。 % X1=dP/de X4=dQ/df p=2*i1。 % dP/df=dQ/de X3=X2。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 %節(jié)點無功功率差 %=============== 以 上 為 除 平 衡 節(jié) 點 外 其 它 節(jié) 點 的 功 率 計 算 ================= %================= 求取 Jacobi 矩陣 =================== for j1=1:n if j1~=isbamp。 %電壓模平方 %========= 以下針對非 PV 節(jié)點來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV 節(jié)點 DP=P(i)P1。,Q39。% 節(jié)點功率 P 計算 ei Σ(Gij*ejBij*fj)+fiΣ (Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。% Σ(Gij*ejBij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。D(i)=0。a=a+1。a=0。N0=2*n。 ICT1=0。 %i 節(jié)點無功補償量 end %=================================================================== P=real(S)。 %PV 節(jié)點電壓給定模值 end for i=1:n %給定各節(jié)點注入功率 S(i)=B2(i,1)B2(i,2)。 f(i)=imag(B2(i,3))。B=imag(Y)。)。 %對角元 1側(cè) end %求導納矩陣 disp(39。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 else p=B1(i,2)。 disp(B2) end % % % for i=1:nl %支路數(shù) if B1(i,6)==0 %左節(jié)點處于低壓側(cè) p=B1(i,1)。B2 矩陣 B2=39。 B2(i,5)=B2(i,5)./SB。 B2(i,3)=B2(i,3)./UB。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。B1 矩陣 B1=39。 B1(i,4)=B1(i,4)./YB。 BB2=B2。 YB=SB./UB./UB。請輸入電壓基準值 :UB=39。)。 if ym~=0 %SB=input(39。您輸入的參數(shù)是標么值？（若不是則輸入一個不為零的數(shù)值）39。UB=220。 %－－－－－－－修改部分－－－－－－－－－－－－ ym=0。sida=zeros(1,n)。f=zeros(1,n)。 Y=zeros(n)。請輸入各節(jié)點參數(shù)形成的矩陣： B2=39。 0 0 220 0 2。 0 0 220 0 2。 0 0 220 0 2。 0 0 220 0 2。 B2=[0 0 235 0 1。 7 9 + 1 0 0。 5 6 + 0 231/11 1 1。 2 9 + 1 0 0。 1 5 + 1 0 0。)。%input(39。)。%input(39。)。%input(39。)。%input(39。 附錄： %本程序的功能是用牛頓 —— 拉夫遜法進行潮流計算 % B1 矩陣： 支路首端號； 末端號； 支路阻抗； 支路對地電納 % 支路的變比； 支路首端處于 K 側(cè)為 1， 1 側(cè)為 0 % B2 矩陣： 該節(jié)點發(fā)電機功率； 該節(jié)點負荷功率； 節(jié)點電壓初始值 % PV 節(jié)點電壓 V 的給定值； 節(jié)點所接的無功補償設備的容量 % 節(jié)點分類標號 clear。 《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩，中國電力出版社 。我們應該努力學好專業(yè)課，打好專業(yè)基礎，在以后的工作崗位上有所成就。 此次課程設計，讓我充分了解了解到精準的潮流計算不僅可以節(jié)省能 源并充分發(fā)揮能源的作用，還可以使電網(wǎng)處于穩(wěn)定且平穩(wěn)的狀態(tài)，這對電網(wǎng)的完全運行起到關鍵性的作用。在 DDTRS 仿真時，對參數(shù)的設置需要細心，有一點疏忽就會導致運行結(jié)果不合理，經(jīng)過調(diào)整使得兩種方法的計算結(jié)果相差甚小。 在潮流計算的過程中，為了滿足調(diào)壓要求，需要一次一次的嘗試，首選的調(diào)壓措施是改變機端電壓，一般不能滿足所有的節(jié)點的電壓要求，這就需要調(diào)節(jié)個別的變壓器分接頭，這是兩種較為經(jīng)濟的調(diào)壓措施。 七、心得體會 從本次電力系統(tǒng)潮流計算課程設計中，通過計算在不同的負荷下各支路節(jié)點潮流分布的變化，分析負荷變動時，對電壓，功率的影響，我們掌握了常用的調(diào)整電壓的措施 ， 鍛煉了自己實際操作分析能力 。雖然通過無功補償可以滿足要求，但其成本投入比較大，造價高。 無功補償設備的添加常見的有：串聯(lián)電抗器、并聯(lián)電容器、調(diào)相機等。用變壓器調(diào)壓的前提是必須保證系統(tǒng)的無功功率充足。適用于出廠母線附近沒有其他的用電設備網(wǎng)絡，調(diào)壓范圍在— 之間。 電力系統(tǒng)調(diào)壓措施中，常用的調(diào)壓手段有：借改變發(fā)電機機端電壓、改變變壓器的變比、無功補償設備的添加等。保證用戶的電壓接近額定值是電力系統(tǒng)調(diào)度的基本任務之一。 電壓是衡量電力系統(tǒng)電能質(zhì)量的標準之一。用以研究系統(tǒng)規(guī)劃和運行中提出的各種問題。 電壓 1 2 3 4 5 6 7 8 9 10 DDTRS MATLAB 由上述比較可觀察到 : 應用 DDTRS 仿真和 MATLAB 計算的到的結(jié)果大致相同，可以認為此潮流計算是正確的 。具體的潮流分布如下： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)= S(9,2)=+ S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)=+ S(9,7)= S(9,10)=+ S(10,9)= 各節(jié)點的功率 S為 (節(jié)點號從小到大排列 )： + + + 仿真比較 線路每次發(fā)生變化時相應的修改 DDTRS 中電路的線路參數(shù)，仿真并進行調(diào)節(jié)，可以得到相應的潮流分布圖，現(xiàn)與 由 MATLAB 得到的潮流分布數(shù)據(jù)圖比較（以斷開 9 節(jié)點間的支路為例）： 由 MATLAB 得到潮流數(shù)據(jù)圖： 由 DDTRS 得到潮流分布圖： 在用 DDTRS 繪制的系統(tǒng)潮流分布圖中，各個節(jié)點電壓在調(diào)試之后均能滿足要求，均達到系統(tǒng)電壓等級標準。具體的潮流分布如下： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)= S(9,2)=+ S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)=+ S(9,7)= S(9,10)=+ S(10,9)= 各節(jié)點的功率 S為 (節(jié)點號從小到大排列 )： + + + 斷開 9 支路的一條線 當斷開 9節(jié)點間的雙回線中的一條時 時， B1矩陣 的改變， B1 矩陣變?yōu)? 1 3 + 1 0 0 1 5 + 1 0 0 1 7 + 1 0 0 2 9 17+ 1 0 0 3 4 + 0 231/ 1 1 5 6 + 0 231/11 1 1 7 8 + 0 231/ 1 1 7 9 + 1 0 0 9 10 + 0 231/11 1 1 用修改后的 B2 矩陣帶入程序，用 MATLAB 進行潮流計算并進行調(diào)整， 對應的調(diào)節(jié)方法和結(jié)果見下表： 調(diào)節(jié)方法： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭四 未調(diào)整 235 235 231 231 231 231 調(diào)整后 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 調(diào)整后 由上表觀察到， 調(diào)節(jié)后 的電壓均在題目允許的范圍內(nèi)，對線路損耗進行分析，統(tǒng)計每次調(diào)整后各個支路的有功損耗，記錄于下表： 功率損耗 未調(diào)整 調(diào)整后 B1= (1,3)支路 (1,5)支路 (1,7)支路 (2,9)支路 + + (3,4)支路 + + (5,6)支路 + + (7,8)支路 + + (7,9)支路 (9,10)支路 + + 總損耗 由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認為是合理的。具體的潮流分布如下： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)= S(9,2)=+ S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)=+ S(9,7)= S(9,10)=+ S(10,9)= 各節(jié)點的功率 S為 (節(jié)點號從小到大排列 )： + + + 斷開 5 支路的一條線 當斷開 5節(jié)點間的雙回線中的一條時 時， B1矩陣 的改變， B1 矩陣變?yōu)?： 1 3 + 1 0