【正文】 .......................................................... 錯誤 !未定義書簽。 《電力系統(tǒng)分析》課程設計 基于 Matlab 的電力系統(tǒng)潮流計算 專 業(yè) 電氣工程及其自動化 班 級 1402 班 姓 名 郭娟 學 號 14039208 同 組 成 員 孫帆 郭娟 陳朝輝 蘇新波 呂應 發(fā) 指 導 教 師 楚冰清 完 成 時 間 2021 .12 .30 1 目 錄 1 概述 .......................................................................................................... 錯誤 !未定義書簽。 2 電力系統(tǒng)潮流計算概述 .......................................................................... 錯誤 !未定義書簽。電力系統(tǒng)的出現(xiàn)，使高效、無污染、使用方便、易于調控的電能得到廣泛應用，推動了社會生產各個領域的變化，開創(chuàng)了電 力時代，發(fā)生了第二次技術革命。其目的是校驗系統(tǒng)是否能安全運行，即是否有過負荷的元件或電壓過低的母線等。解決電力系統(tǒng)潮流計算問題是以導納距陣為基礎的，因此，只要在迭代過程中盡可能保持方程式系數(shù)距陣的稀疏性，就可以大大提高牛頓法潮流程序的放率。它的任務是根據(jù)給定的運行條件和網(wǎng)路結構確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡中的功率分布以及功率損耗等。手工方式可用來計算一些接線較簡單的電力網(wǎng)，但若將其用于接線復雜的電力網(wǎng)則計算量過大，難于保證計算準確性；計算機方式從數(shù)學上看可歸結為用數(shù)值方法解非線性代數(shù)方程，數(shù)學邏輯簡單完整，可快速精確地完成計算，但其缺點是物理概念不明顯，物理規(guī)律被埋沒在循環(huán)往復的數(shù)值求解過程中，基本原理不太明顯。待求的運行狀態(tài)參量包括電網(wǎng)各母線節(jié)點的電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡的功率損耗等。對現(xiàn)有的電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計以及對電力系統(tǒng)進行靜態(tài)和穩(wěn)態(tài)分析都是以潮流計算為基礎。在電力系統(tǒng)調度運行的多個領域問題是研究電力系統(tǒng)穩(wěn)態(tài)問題的基礎和前提。阻抗法改善了系統(tǒng)潮流計算問題的收斂性，解決了導納無法求解的一些系統(tǒng)的潮流計算，在 60 年代獲得了廣泛的應用，阻抗法德主要缺點是占用計算機內存大，每次迭代的計算量大。這是數(shù)學中解決非線性方程式的典型方法，有較好的收斂性。速度方面都超過了阻抗法，成為了 60 年代末期以后廣泛采用的優(yōu)秀方法。它在數(shù)學類科技應用軟件中在 數(shù)值計算 方面首屈一指?？梢灾苯诱{用 ,用戶也可以將自己編寫的實用程序導入到 MATLAB 函數(shù)庫中方便自己以后調用，此外許多的 MATLAB 愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進行下載就可以用。2﹪） /， ΔPs=110kW， ΔP0=，U0﹪ =， I0﹪ =； 變壓器 T3： SFL18000/110，（ 110177。 潮流計算過程 各元件參數(shù)計算 ① 120Km 線路 7 SSlbBlxXlrR46101101101?? ?????????????????? （ 1） ② 100Km 線路 SSlbBlxXlrR46200202202?? ?????????????????? （ 2） ③ 70Km 線路 SSlbBlxXlrR46303303303?? ?????????????????? （ 3） ④ 變壓器 T1， T2 M V AjK V AjSIjPSSUUXXSUPRNNNkTTNNkT)()1001 6 0 0 (100%1 6 0 0 010010100%1 6 0 0 0121110100013232213223223???????????????????????????? （ 4） ⑤ 變壓器 T3 M V AjK V AjSIjPSSUUXSUPRNNNKTNKKT)()100(100%800010010100%80001105210003323233223223??????????????????????????? （ 5） ⑥ 變壓器 T4 K V AjK V AjSIjPSSUUXSUPRNNNKTNNKT)()1001 6 0 0 (100%1 6 0 0 010010100%1 6 0 0 01106210004323243223224??????????????????????????? （ 6） 8 錯誤 !未找到引用源。39。39。jjUBjXRUSSSM V AjUBjSSM V AjjSSNNLDLDNLDLDLD????????????????? M V AjM V Aj )8 6 7 8 (1 8 5 ??? （ 13） ⑧ 間 100Km 線路損耗 )()(2)()(2)()(221222122333222323)2(223232342339。1422)1(139。)1(111141 ??????????? U XQRPUUUU ??? （ 18） 由于 的輸入電壓大于 110 額定值， 所以調壓關系不滿足。)1(1239。 n=8。 B1=[1 2 + 1 0。 4 5 + 0 1。 8 9 + 0 1。 0 50+ 220 0 0 2。 0 50+ 220 0 0 2。 Y=zeros(n)。sida=zeros(1,n)。 if ym~=0 %SB=input(39。請輸入電壓基準值 :UB=39。 BB2=B2。B1 矩陣 B1=39。 B2(i,3)=B2(i,3)./UB。B2 矩陣 B2=39。 else p=B1(i,2)。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。)。 f(i)=imag(B2(i,3))。 %i 節(jié)點無功補償量 end %=================================================================== P=real(S)。N0=2*n。a=a+1。%Σ(Gij*ejBij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。,Q39。 %節(jié)點無功功率差 %=============== 以 上 為 除 平 衡 節(jié) 點 外 其 它 節(jié) 點 的 功 率 計 算 ================= %================= 求取 Jacobi 矩陣 =================== for j1=1:n if j1~=isbamp。 % dP/df=dQ/de X3=X2。 J(p,q)=X3。J(m,N)=DP。 elseif j1==iamp。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。J(p,q)=X3。q=q+1。 end end 15 else %=============== 下面是針對 PV 節(jié)點來求取 Jacobi 矩陣的元素 =========== DP=P(i)P1。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 p=2*i1。