【正文】 .......................................................... 錯(cuò)誤 !未定義書簽。 《電力系統(tǒng)分析》課程設(shè)計(jì) 基于 Matlab 的電力系統(tǒng)潮流計(jì)算 專 業(yè) 電氣工程及其自動(dòng)化 班 級(jí) 1402 班 姓 名 郭娟 學(xué) 號(hào) 14039208 同 組 成 員 孫帆 郭娟 陳朝輝 蘇新波 呂應(yīng) 發(fā) 指 導(dǎo) 教 師 楚冰清 完 成 時(shí) 間 2021 .12 .30 1 目 錄 1 概述 .......................................................................................................... 錯(cuò)誤 !未定義書簽。 2 電力系統(tǒng)潮流計(jì)算概述 .......................................................................... 錯(cuò)誤 !未定義書簽。電力系統(tǒng)的出現(xiàn)，使高效、無污染、使用方便、易于調(diào)控的電能得到廣泛應(yīng)用，推動(dòng)了社會(huì)生產(chǎn)各個(gè)領(lǐng)域的變化，開創(chuàng)了電 力時(shí)代，發(fā)生了第二次技術(shù)革命。其目的是校驗(yàn)系統(tǒng)是否能安全運(yùn)行，即是否有過負(fù)荷的元件或電壓過低的母線等。解決電力系統(tǒng)潮流計(jì)算問題是以導(dǎo)納距陣為基礎(chǔ)的，因此，只要在迭代過程中盡可能保持方程式系數(shù)距陣的稀疏性，就可以大大提高牛頓法潮流程序的放率。它的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。手工方式可用來計(jì)算一些接線較簡單的電力網(wǎng)，但若將其用于接線復(fù)雜的電力網(wǎng)則計(jì)算量過大，難于保證計(jì)算準(zhǔn)確性；計(jì)算機(jī)方式從數(shù)學(xué)上看可歸結(jié)為用數(shù)值方法解非線性代數(shù)方程，數(shù)學(xué)邏輯簡單完整，可快速精確地完成計(jì)算，但其缺點(diǎn)是物理概念不明顯，物理規(guī)律被埋沒在循環(huán)往復(fù)的數(shù)值求解過程中，基本原理不太明顯。待求的運(yùn)行狀態(tài)參量包括電網(wǎng)各母線節(jié)點(diǎn)的電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡(luò)的功率損耗等。對現(xiàn)有的電力系統(tǒng)的運(yùn)行和擴(kuò)建，對新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)以及對電力系統(tǒng)進(jìn)行靜態(tài)和穩(wěn)態(tài)分析都是以潮流計(jì)算為基礎(chǔ)。在電力系統(tǒng)調(diào)度運(yùn)行的多個(gè)領(lǐng)域問題是研究電力系統(tǒng)穩(wěn)態(tài)問題的基礎(chǔ)和前提。阻抗法改善了系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納無法求解的一些系統(tǒng)的潮流計(jì)算，在 60 年代獲得了廣泛的應(yīng)用，阻抗法德主要缺點(diǎn)是占用計(jì)算機(jī)內(nèi)存大，每次迭代的計(jì)算量大。這是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。速度方面都超過了阻抗法，成為了 60 年代末期以后廣泛采用的優(yōu)秀方法。它在數(shù)學(xué)類科技應(yīng)用軟件中在 數(shù)值計(jì)算 方面首屈一指?？梢灾苯诱{(diào)用 ,用戶也可以將自己編寫的實(shí)用程序?qū)氲?MATLAB 函數(shù)庫中方便自己以后調(diào)用，此外許多的 MATLAB 愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進(jìn)行下載就可以用。2﹪） /， ΔPs=110kW， ΔP0=，U0﹪ =， I0﹪ =； 變壓器 T3： SFL18000/110，（ 110177。 潮流計(jì)算過程 各元件參數(shù)計(jì)算 ① 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 錯(cuò)誤 !未找到引用源。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 額定值， 所以調(diào)壓關(guān)系不滿足。)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。請輸入電壓基準(zhǔn)值 :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é)點(diǎn)無功補(bǔ)償量 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é)點(diǎn)無功功率差 %=============== 以 上 為 除 平 衡 節(jié) 點(diǎn) 外 其 它 節(jié) 點(diǎn) 的 功 率 計(jì) 算 ================= %================= 求取 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é)點(diǎn)來求取 Jacobi 矩陣的元素 =========== DP=P(i)P1。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 p=2*i1。