【正文】 jjjSSS ) 6 0 ()() 5 0 )2(1)1(11 ???????? （ （ 15）⑩ 位置 2 點(diǎn)總損耗： M V AjjjSSS )()()()2(2)1(22 ???????? （ 16） 調(diào)壓計(jì)算 計(jì)算 線路上的電壓值 MV AjXRUSSUBjSSNLDLDN )()(2 222 21439。39。 電網(wǎng)潮流計(jì)算思路 （ 1） 計(jì)算各元件參數(shù)，畫出等值電路； （ 2） 進(jìn)行網(wǎng)絡(luò)潮流計(jì)算； （ 3） 不滿足供電要求，進(jìn)行調(diào)壓計(jì)算。在新的版本中也加入了對(duì) C， FORTRAN， C++ ， JAVA的支持。內(nèi)存要求。這就迫使電力系統(tǒng)的計(jì)算人員轉(zhuǎn)向以阻抗矩陣為基礎(chǔ)的逐次代入法。各點(diǎn)電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。兩種方法各有優(yōu)缺點(diǎn)。牛頓法是 數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。電力系統(tǒng)建設(shè)往往是國(guó)家及地區(qū)國(guó)民經(jīng)濟(jì)發(fā)展規(guī)劃的重要組成部分。 《電力系統(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ò)誤 !未定義書簽。電力系統(tǒng)的出現(xiàn)，使高效、無污染、使用方便、易于調(diào)控的電能得到廣泛應(yīng)用，推動(dòng)了社會(huì)生產(chǎn)各個(gè)領(lǐng)域的變化，開創(chuàng)了電 力時(shí)代，發(fā)生了第二次技術(shù)革命。解決電力系統(tǒng)潮流計(jì)算問題是以導(dǎo)納距陣為基礎(chǔ)的，因此，只要在迭代過程中盡可能保持方程式系數(shù)距陣的稀疏性，就可以大大提高牛頓法潮流程序的放率。手工方式可用來計(jì)算一些接線較簡(jiǎn)單的電力網(wǎng)，但若將其用于接線復(fù)雜的電力網(wǎng)則計(jì)算量過大，難于保證計(jì)算準(zhǔn)確性；計(jì)算機(jī)方式從數(shù)學(xué)上看可歸結(jié)為用數(shù)值方法解非線性代數(shù)方程，數(shù)學(xué)邏輯簡(jiǎn)單完整，可快速精確地完成計(jì)算，但其缺點(diǎn)是物理概念不明顯，物理規(guī)律被埋沒在循環(huán)往復(fù)的數(shù)值求解過程中，基本原理不太明顯。對(duì)現(xiàn)有的電力系統(tǒng)的運(yùn)行和擴(kuò)建，對(duì)新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)以及對(duì)電力系統(tǒng)進(jìn)行靜態(tài)和穩(wěn)態(tài)分析都是以潮流計(jì)算為基礎(chǔ)。阻抗法改善了系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納無法求解的一些系統(tǒng)的潮流計(jì)算，在 60 年代獲得了廣泛的應(yīng)用，阻抗法德主要缺點(diǎn)是占用計(jì)算機(jī)內(nèi)存大，每次迭代的計(jì)算量大。速度方面都超過了阻抗法，成為了 60 年代末期以后廣泛采用的優(yōu)秀方法。可以直接調(diào)用 ,用戶也可以將自己編寫的實(shí)用程序?qū)氲?MATLAB 函數(shù)庫中方便自己以后調(diào)用，此外許多的 MATLAB 愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進(jìn)行下載就可以用。 潮流計(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。1422)1(139。)1(1239。 B1=[1 2 + 1 0。 8 9 + 0 1。 0 50+ 220 0 0 2。sida=zeros(1,n)。請(qǐng)輸入電壓基準(zhǔn)值 :UB=39。B1 矩陣 B1=39。B2 矩陣 B2=39。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 f(i)=imag(B2(i,3))。N0=2*n。%Σ(Gij*ejBij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。 %節(jié)點(diǎn)無功功率差 %=============== 以 上 為 除 平 衡 節(jié) 點(diǎn) 外 其 它 節(jié) 點(diǎn) 的 功 率 計(jì) 算 ================= %================= 求取 Jacobi 矩陣 =================== for j1=1:n if j1~=isbamp。 J(p,q)=X3。 elseif j1==iamp。J(p,q)=X3。 end end 15 else %=============== 下面是針對(duì) PV 節(jié)點(diǎn)來求取 Jacobi 矩陣的元素 =========== DP=P(i)P1。 p=2*i1。J(p,q)=X6。 p=2*i1。J(p,q)=X6。%消去運(yùn)算 end J(k3,k)=0。 %修改節(jié)點(diǎn)電壓實(shí)部 k1=k+1。)。 E(k)=e(k)+f(k)*j。39。 disp(39。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。)。)。 SSi(p,q)=Si(p,q)*SB。 disp(ZF)。 for i=1:nl p=B1(i,1)。 ZF=[39。 %disp(SSj(q,p))。q=B1(i,2)。,num2str(q),39。 end 運(yùn)行結(jié)果 21 22 23 6 心得體會(huì) 這次課程設(shè)計(jì)，雖然這次課程設(shè)計(jì)沒有想象中的那么簡(jiǎn)單，花費(fèi)了很多的時(shí)間來不斷地實(shí)驗(yàn)，雖說花費(fèi)了大量的時(shí)間但是最后終于完成了預(yù)定的目 的，收獲頗豐，也很好的鍛煉了自己運(yùn)用知識(shí)的能力。同時(shí)在設(shè)計(jì)過程中要求我們有足夠的耐心。從中收獲知識(shí) 、 收獲閱歷 、 收獲成熟，在此過程中，通過查找大量資料，以及不懈的努力，不僅培養(yǎng)了獨(dú)立思考、動(dòng)手操作的能力，在各種其它能力上也都有了提高。只有不斷的測(cè)試自己，挑戰(zhàn)自己，才能擁有更多的成功和快樂！享受過程，而不是結(jié)果！認(rèn)真對(duì)待每一件事，珍惜每一分一秒，學(xué)到最多的知識(shí)和方法，鍛煉自己的能力，這個(gè)是我課程設(shè)計(jì)中學(xué)到的 最重要的東西，也是以后都將受益。傳統(tǒng)的手工計(jì)算 工作量大，易出錯(cuò)，潮流計(jì)算費(fèi)時(shí)費(fèi)力。 此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過精密的設(shè) 24 計(jì)和計(jì)算。,num2str(DDS(i))]。 disp(DS(i))。39。,num2str(q),39。 if B1(i,6)==0 Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q).