【正文】 .............................. 23 參考文獻(xiàn) ................................................................................................................................... 25 2 1 概述 設(shè)計目的 在如今的社會，電力已經(jīng)成為人們必不可少的需求，而建立結(jié)構(gòu)合理的大型電力系統(tǒng)不僅便于電能生產(chǎn)與消費的集中管理、統(tǒng)一調(diào)度和分配，減少總裝機容量節(jié)省動力設(shè)施投資，且有利于地區(qū)能源資源的合理開發(fā)利用，更大限度地滿足地區(qū)國民經(jīng)濟日益增長的用電需要。但是一般靜態(tài)安全分析需要校驗的狀態(tài)數(shù)非常多，用嚴(yán)格的潮流計算來分析這些狀態(tài)往往計算量過大，因此不得不尋求一些特殊的算法以滿足要求。潮流計算可以用傳統(tǒng)的手工方式進(jìn)行，也可以計算機軟件完成。潮流計算的目標(biāo)是求取電力系統(tǒng)在給定運行狀態(tài)的計算，即節(jié)點電壓和功率 分布，用以檢查系統(tǒng)各元件是否過負(fù)荷。這個方法的原理比較簡單，要求的數(shù)字計算機內(nèi)存量比較差下，適應(yīng) 50 年代電子計算機制造水平和當(dāng)時電力系統(tǒng)理論水平，但它的收斂性較差，當(dāng)系統(tǒng)規(guī)模變大時，迭代次數(shù)急劇上升，在計算中往往出現(xiàn)迭代不收斂的情況。自從 60 年代中期，牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。 Matlab 的應(yīng)用 MATLAB 的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，故用 MATLAB 來解算問題要比用 C， FORTRAN 等語言完成相同的事情簡捷得多，并且 mathwork 也吸收了像 Maple 等軟件的優(yōu)點 ,使 MATLAB成為一個強大的數(shù)學(xué)軟件。2﹪） /， ΔPs=62kW， ΔP0=，Us﹪ =， I0﹪ =； 導(dǎo)線型號均為 LGJ150，參數(shù) r0=， x0=， b0=106S/km。?????????????????M V Aj )( ?? （ 11） ⑥ 間 120Km 線路損耗 )()(2)()(2)()(221022102111221313)2(121131331339。39。 ) 2 1(1 ??T 的變壓器取輸出 %)51(121 ??? 錯誤 !未找到引用源。 isb=1。 6 7 + 0 1。 0 40+ 220 0 0 2。f=zeros(1,n)。)。 B1(i,4)=B1(i,4)./YB。 B2(i,5)=B2(i,5)./SB。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。B=imag(Y)。 ICT1=0。D(i)=0。 %電壓模平方 %========= 以下針對非 PV 節(jié)點來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV 節(jié)點 DP=P(i)P1。 % X1=dP/de X4=dQ/df p=2*i1。 J(p,q)=X4。% dQ/df p=2*i1。J(m,N)=DP。 % dP/df X5=0。J(m,N)=DP。% dP/df X5=2*e(i)。J(m,N)=DP。 % 對角元規(guī)格化 if k~=3 % 不是第三行 %============================================================ k4=k1。 end end end %====上面是用 線性變換方式將 Jacobi 矩陣化成單位矩陣 ===== for k=3:2:N01 L=(k+1)./2。 end %用高斯消去法解 w=J*V disp(39。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 disp(EE)。 VV=V*UB。 disp(sida)。 disp(39。 disp(39。 Siz(i)=Si(p,q)。)=39。各條支路的末端功率 Sj 為 (順序同您輸入 B1 時一致 )： 39。 end disp(Sj(q,p))。,num2str(SSj(q,p))]。)。,num2str(p),39。39。尤其是在剛開始的時候，不知道該從哪里入手，更是意識到這一點的重要性。 這次課程設(shè)計，雖然沒有想象中的那么簡單，花費了很多的時間來不斷地實驗，但是最后終于完成了預(yù)定的目的，收獲頗豐，也很好的鍛煉了自己運用知識的能力。水利電力出版社 ,1990. [5]陶海英 .基于 MATLAB 的電力系統(tǒng)穩(wěn)態(tài)仿真分析 [J].江西電力職業(yè)技術(shù)學(xué)院學(xué)報 ,2021(3):8082. [6]徐敏 ,彭瑜 .MATLAB 在《電力系統(tǒng)分析》教學(xué)中的應(yīng)用 [J].電力系統(tǒng)及其自動化學(xué)報 ,2021,(7):152155. 。利用牛頓 拉夫遜法、采用 Matlab 語言進(jìn)行潮流計算，提高了準(zhǔn)確性，使用方便，已得到成功應(yīng)用，而它在我們以后的專業(yè)課程的學(xué)習(xí)中應(yīng)用也是很廣泛的，掌握好了 MATLAB，這為以后的專業(yè)課程的學(xué)習(xí)將有很大的幫助。同時由于求解過程中用到求節(jié)點導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對以前所學(xué)的知識有了一次很好的溫習(xí)。 %disp(DDS(i))。 ZF=[39。 end disp(39。,num2str(p),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))))。)。,39。 Siz(i)=Si(p,q)。39。)。各節(jié)點的電壓相角 sida 為 (節(jié)點號從小到大排列 )： 39。)。 disp(E)。)。 %不滿足要求的節(jié)點數(shù)加 1 end end ICT2(a)=IT2。 end 17 else for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 % N=N0+1 即 N=2*n+1 擴展列 △ P、△ Q for k2=k1:N1 % 擴展列 △ P、 △ Q J(k,k2)=J(k,k2)./J(k,k)。 m=p+1。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 m=p+1。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。q=q+1。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(m,N)=DP。 % dP/df=dQ/de X3=X2。,Q39。a=a+1。 %i 節(jié)點無功補償量 end %=================================================================== P=real(S)。)。 else p=B1(i,2)。 B2(i,3)=B2(i,3)./UB。 BB2=B2。 if ym~=0 %SB=input(39。 Y=zeros(n)。 0 50+ 220 0 0 2。 4 5 + 0 1。 n=8。)1(111141 ??????????? U XQRPUUUU ??? （ 18） 由于 的輸入電壓大于 110 額定值， 所以調(diào)壓關(guān)系不滿足。jjUBjXRUSSSM V AjUBjSSM V AjjSSNNLDLDNLDLDLD????????????????? M V AjM V Aj )8 6 7 8 (1 8 5 ??? （ 13） ⑧ 間 100Km 線路損耗 )()(2)()(2)()(221222122333222323)2(223232342339。39。2﹪） /， ΔPs=110kW， ΔP0=，U0﹪ =， I0﹪ =； 變壓器 T3：