電力系統(tǒng)分析潮流計(jì)算課程序設(shè)計(jì)及其MATLAB程序設(shè)計(jì)

電力系統(tǒng)分析潮流計(jì)算課程序設(shè)計(jì)及其MATLAB程序設(shè)計(jì)目錄TOC\o"1-3"\h\u一、程序設(shè)計(jì)目的 2二、程序設(shè)計(jì)要求 4三、13節(jié)點(diǎn)配網(wǎng)潮流計(jì)算 43.1主要流程 4第一步的前推公式如下（1-1）-（1-5）： 4第二步的回代公式如下（1-6）—（1-9）： 43.2配網(wǎng)前推后代潮流計(jì)算的原理 53.3配網(wǎng)前推后代潮流計(jì)算迭代過程 53.3計(jì)算原理 6四、計(jì)算框圖流程 7五、確定前推回代支路次序 7六、前推回代計(jì)算輸入文件 8主程序： 8輸入文件清單： 9計(jì)算結(jié)果： 10數(shù)據(jù)分析： 10七、配電網(wǎng)潮流計(jì)算的要點(diǎn) 11八、自我總結(jié) 11九、參考文獻(xiàn) 12附錄一MATLAB的簡介 12華南理工大學(xué)廣州學(xué)院電氣工程學(xué)院課程設(shè)計(jì)報(bào)告PAGEPAGE13一、程序設(shè)計(jì)目的開式網(wǎng)絡(luò)潮流計(jì)算：配電網(wǎng)的結(jié)構(gòu)特點(diǎn)呈輻射狀，在正常運(yùn)行時(shí)是開環(huán)的；配電網(wǎng)的潮流計(jì)算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收斂性。(1)在電網(wǎng)規(guī)劃階段，通過潮流計(jì)算，合理規(guī)劃電源容量及接入點(diǎn)，合理規(guī)劃網(wǎng)架，選擇無功補(bǔ)償方案，滿足規(guī)劃水平年的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。

(2)在編制年運(yùn)行方式時(shí)，在預(yù)計(jì)負(fù)荷增長及新設(shè)備投運(yùn)基礎(chǔ)上，選擇典型方式進(jìn)行潮流計(jì)算，發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié)，供調(diào)度員日常調(diào)度控制參考，并對(duì)規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu)，加快基建進(jìn)度的建議。

(3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算，用于日運(yùn)行方式的編制，指導(dǎo)發(fā)電廠開機(jī)方式，有功、無功調(diào)整方案及負(fù)荷調(diào)整方案，滿足線路、變壓器熱穩(wěn)定要求。

及電壓質(zhì)量要求。(4)預(yù)想事故、設(shè)備退出運(yùn)行對(duì)靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案。圖113節(jié)點(diǎn)配電網(wǎng)結(jié)構(gòu)圖表1系統(tǒng)支路參數(shù)支路R（Ω）X（Ω）B/2(TK)（S）1~23.3673.6850.02~32.3562.5410.03~41.1451.280.04~54.5245.040.02~60.8561.140.06~72.7452.9650.02~83.7434.2510.08~92.2372.7560.03~104.144.6960.03~111.3281.7630.011~122.4362.8660.04~133.5213.9660.0表2系統(tǒng)負(fù)荷參數(shù)節(jié)點(diǎn)編號(hào)節(jié)點(diǎn)類型節(jié)點(diǎn)初始電（kV）Pi（MVA）Qi（MVA）1根節(jié)點(diǎn)10.4002中間節(jié)點(diǎn)10.00.03420.03013中間節(jié)點(diǎn)10.00.06930.06424中間節(jié)點(diǎn)10.00.08450.07635葉節(jié)點(diǎn)10.00.02950.02616中間節(jié)點(diǎn)10.00.04740.04097葉節(jié)點(diǎn)10.00.11760.09578中間節(jié)點(diǎn)10.00.09460.08579葉節(jié)點(diǎn)10.00.09160.085910葉節(jié)點(diǎn)10.00.02710.022911中間節(jié)點(diǎn)10.00.06960.064312葉節(jié)點(diǎn)10.00.06760.057913葉節(jié)點(diǎn)10.00.02980.0242二、程序設(shè)計(jì)要求看懂前推回代法計(jì)算程序；報(bào)告敘述計(jì)算原理及計(jì)算流程；繪制計(jì)算流程框圖；確定前推回代支路次序（廣度優(yōu)先，或深度優(yōu)先）并編寫前推回代計(jì)算輸入文件，然后進(jìn)行潮流計(jì)算；整理專利計(jì)算結(jié)果；總結(jié)配電網(wǎng)潮流計(jì)算的要點(diǎn)；三、設(shè)計(jì)內(nèi)容1.根據(jù)電力系統(tǒng)網(wǎng)絡(luò)推導(dǎo)電力網(wǎng)絡(luò)數(shù)學(xué)模型，寫出節(jié)點(diǎn)導(dǎo)納矩陣；2.賦予各節(jié)點(diǎn)電壓變量（直角坐標(biāo)系形式）初值后，求解不平衡量；3.形成雅可比矩陣；4.求解修正量后，重新修改初值，從2開始重新循環(huán)計(jì)算；5.求解的電壓變量達(dá)到所要求的精度時(shí)，再計(jì)算各支路功率分布、功率損耗和節(jié)點(diǎn)6.上機(jī)編程調(diào)試；連調(diào)；7.計(jì)算分析給定系統(tǒng)潮流分析并與手工計(jì)算結(jié)果作比較分析。8.準(zhǔn)備計(jì)算機(jī)演示答辯，書寫該課程設(shè)計(jì)說明書（必須計(jì)算機(jī)打印）四、13節(jié)點(diǎn)配網(wǎng)潮流計(jì)算4.1牛拉法的原理及其基本方程牛頓迭代法（Newton'smethod）又稱為牛頓-拉夫遜（拉弗森）方法（Newton-Raphsonmethod），它是牛頓在17世紀(jì)提出的一種在實(shí)數(shù)域和復(fù)數(shù)域上近似求解方程的方法。設(shè)r是f(x)=0的根，選取x0作為r初始近似值，過點(diǎn)（x0,f(x0)）做曲線y=f(x)的切線L，L的方程為y=f(x0)f'(x0)(x-x0)，求出L與x軸交點(diǎn)的橫坐標(biāo)x1=x0-f(x0)/f'(x0)，稱x1為r的一次近似值。過點(diǎn)（x1,f(x1)）做曲線y=f(x)的切線，并求該切線與x軸的橫坐標(biāo)x2=x1-f(x1)/f'(x1)，稱x2為r的二次近似值。重復(fù)以上過程，得r的近似值序列，其中x(n+1)=x(n)－f(x(n))/f'(x(n))，稱為r的n+1次近似值，上式稱為牛頓迭代公式。解非線性方程f(x)=0的牛頓法是把非線性方程線性化的一種近似方法。把f(x)在x0點(diǎn)附近展開成泰勒級(jí)數(shù)f(x)=f(x0)+(x－x0)f'(x0)+(x－x0)^2*f''(x0)/2!+…取其線性部分，作為非線性方程f(x)=0的近似方程，即泰勒展開的前兩項(xiàng)，則有f(x0)+f'(x0)(x－x0)=f(x)=0設(shè)f'(x0)≠0則其解為x1=x0－f(x0)/f'(x0)這樣，得到牛頓法的一個(gè)迭代序列：x(n+1)=x(n)－f(x(n))/f'(x(n))。4.2PQ分解法的原理及其基本方程PQ分解法的基本思想是根據(jù)電力系統(tǒng)實(shí)際運(yùn)行特點(diǎn)：通常網(wǎng)絡(luò)上的電抗遠(yuǎn)大于電阻值，則系統(tǒng)母線電壓幅值的微小變化對(duì)母線有功功率的改變影響很小。同樣，母線電壓相角的少許改變，也不會(huì)引起母線有功功率的明顯改變，因此，節(jié)點(diǎn)功率方程在用極坐標(biāo)形式表示時(shí)，它的修正方程式可簡化為（2-1）這就是把2（n-1）階的線性方程組變成了兩個(gè)n-1階的線性方程組，將P和Q分開來進(jìn)行迭代計(jì)算，因而大大地減少了計(jì)算工作量。但是H、L在迭代過程中仍然在不斷的變化而且又都是不對(duì)稱矩陣。對(duì)牛頓法的進(jìn)一步簡化，即把式（2-1）中的系數(shù)矩陣簡化為在迭代過程中不變的對(duì)稱矩陣。在一般情況下，線路兩端的電壓相角是不大的（不超過10°～20°），因此，可以認(rèn)為（2-2）此外，與系統(tǒng)各節(jié)點(diǎn)無功功率相應(yīng)的導(dǎo)納遠(yuǎn)小于該節(jié)點(diǎn)自導(dǎo)納的虛部，即因而（2-3）考慮到以上關(guān)系，式（2-1）的系數(shù)矩陣中的個(gè)元素可以表示為(i,j=1,2,…,n-1)(i,j=1,2,…,m)而系數(shù)矩陣H和L則可以分別寫成：(2-4)(2-5)將式（2-4）和（2-5）代人式中，得到用和分別左乘以上兩式，便得(2-6)(2-7)這就是簡化了的修正方程式，它們也可以展開寫成(2-8)(2-9)在這兩個(gè)修正方程式中系數(shù)矩陣元素就是系統(tǒng)導(dǎo)納矩陣的虛部，因而系數(shù)矩陣是對(duì)稱矩陣，且在迭代過程中保持不變。這就大大減少了計(jì)算工作量。用極坐標(biāo)表示的節(jié)點(diǎn)功率增量為(2-10)式（2-8）、（2-9）和（2-10）構(gòu)成了PQ分解法迭代過程的基本方程式。4.3配網(wǎng)前推后代潮流計(jì)算的原理前推回代法在配電網(wǎng)潮流計(jì)算中簡單實(shí)用，所有的數(shù)據(jù)都是以矢量形式存儲(chǔ)，因此節(jié)省了大量的計(jì)算機(jī)內(nèi)存，對(duì)于任何種類的配電網(wǎng)只要有合理的R/X值，此方法均可保證收斂。算法的穩(wěn)定性也是評(píng)價(jià)配電網(wǎng)潮流算法的重要指標(biāo)。一般情況下，算法的收斂階數(shù)越高，算法的穩(wěn)定性越差，前推回代法的收斂階數(shù)為一階，因此它也具有較好的穩(wěn)定性。比較而言，前推回代法充分利用了網(wǎng)絡(luò)呈輻射狀的結(jié)構(gòu)特點(diǎn)，數(shù)據(jù)處理簡單，計(jì)算效率高，具有較好的收斂性，被公認(rèn)是求解輻射狀配電網(wǎng)潮流問題的最佳算法之一。4.4配網(wǎng)前推后代潮流計(jì)算迭代過程配網(wǎng)潮流前推回代潮流算法第n+1步的迭代公式如下：節(jié)點(diǎn)i的前推計(jì)算公式為由圖1可知節(jié)點(diǎn)i+1的注入有功功率和無功功率分別為：NBNB-1P(i+1)=∑PL(j)+∑LP(j)Q(i+1)=∑QL(j)+∑LQ(j)其中：i=1，2，。。。NB-1，NB為節(jié)點(diǎn)數(shù)，PL(i)為第i節(jié)點(diǎn)所帶負(fù)荷有功功率；QL(i)為第i節(jié)點(diǎn)所帶負(fù)荷無功功率；LP(i)為第i條線段上的有功功率損耗；LQ(i)為第i條線段上的無功功率損耗。收斂條件以前后兩次迭代的電壓偏差作為迭代收斂條件，(8)式表明節(jié)點(diǎn)電壓幅值最大偏差小于設(shè)定閾值，即認(rèn)為迭代收斂，則迭代過程結(jié)束。max{|Vk(i)-Vk-1(i)|}<ε.(8)在開始計(jì)算的時(shí)候先把各節(jié)點(diǎn)上的LP(i)、LQ(i)分別設(shè)為零，利用式(1)計(jì)算出P(i+1)、Q(i+1)，分別代入式(6)計(jì)算各節(jié)點(diǎn)電壓幅值，作為初始化的節(jié)點(diǎn)電壓值.然后把此電壓值代入公式(2)來修正線路的損耗LP(i)、LQ(i)，用經(jīng)過修正的線路損耗求新的節(jié)點(diǎn)注入功率，再根據(jù)節(jié)點(diǎn)注入功率求新的節(jié)點(diǎn)電壓，根據(jù)前后兩次迭代的電壓偏差是否小于設(shè)定閾值判斷是否收斂，如果滿足，則輸出結(jié)果，如果不滿足,必須重復(fù)上述過程直到滿足收斂條件為止。4.5計(jì)算原理第一步，從與葉節(jié)點(diǎn)聯(lián)系的支路開始，該支路的末端功率即等于葉節(jié)點(diǎn)功率，利用這個(gè)功率和對(duì)應(yīng)的節(jié)點(diǎn)電壓計(jì)算支路功率損耗，求得支路的首端功率。當(dāng)以某節(jié)點(diǎn)為始節(jié)點(diǎn)的各支路都計(jì)算完畢后，便想象將這些支路都拆去，是的該節(jié)點(diǎn)成為新的葉節(jié)點(diǎn)，其節(jié)點(diǎn)功率等于原有的負(fù)荷功率與該節(jié)點(diǎn)為始節(jié)點(diǎn)的各支路首端功率之和。于是計(jì)算便可延續(xù)下去，直到全部支路計(jì)算完畢。第二步，利用第一步所得的支路首端功率和本步驟剛算出的本支路始節(jié)點(diǎn)的電壓（對(duì)電源點(diǎn)為已知電壓），從電源點(diǎn)開始逐條支路進(jìn)行計(jì)算，求得各支路終節(jié)點(diǎn)的電壓。在迭代計(jì)算開始之前，要先處理好支路的計(jì)算順序問題。介紹兩種確定之路計(jì)算順序的方法。第一種方法是，按與葉節(jié)點(diǎn)聯(lián)接的之路順序，并將已排序的之路拆除，在此過程中將不斷出現(xiàn)新的葉節(jié)點(diǎn)，而與其聯(lián)接的支路又加入排序行列。這樣就可以全部排列好從葉節(jié)點(diǎn)向電源點(diǎn)計(jì)算功率損耗的支路順序。其逆序就是進(jìn)行電壓計(jì)算的支路順序。第二種是逐條追加支路的方法。首先從根節(jié)點(diǎn)（電源點(diǎn)）開始接出第一條支路，引出一個(gè)新節(jié)點(diǎn)，以后每次追加的支路都必須從已出現(xiàn)的節(jié)點(diǎn)接出遵循這個(gè)原則逐條追加支路，直到全部支路追加完畢。所得到的支路追加順序即是進(jìn)行電壓計(jì)算的支路順序，其逆序便是功率損耗計(jì)算的支路順序。計(jì)算框圖流程六、前推回代計(jì)算輸入文件主程序：[PQ,FT,RX]=case115();%調(diào)用數(shù)據(jù)文件NN=size(PQ,1);%節(jié)點(diǎn)數(shù)NB=size(FT,1);%支路數(shù)數(shù)V=PQ(:,1);%V初始電壓相量maxd=1k=0;whilemaxd>0.0001PQ2=PQ;%每一次迭代各節(jié)點(diǎn)的注入有功和無功相同PL=0.0;fori=1:NBkf=FT(i,1);%前推始節(jié)點(diǎn)號(hào)kt=FT(i,2);%前推終節(jié)點(diǎn)號(hào)x=(PQ2(kf,2)^2+PQ2(kf,3)^2)/(V(kf)^2);%計(jì)算沿線電流平方APQ1(i,1)=RX(i,1)*x+PQ2(kf,2);%計(jì)算支路首端有功/MWRX(i,1)~RPQ1(i,2)=RX(i,2)*x+PQ2(kf,3);%計(jì)算沿支路的無功損耗/MvarRX(i,2)~XPQ2(kt,2)=PQ2(kt,2)+PQ1(i,1);%用PQ1去修正支路末端節(jié)點(diǎn)的有功P單位MWPQ2(kt,3)=PQ2(kt,3)+PQ1(i,2);%用PQ1去修正支路末端節(jié)點(diǎn)的有功Q單位MvarPL=PL+RX(i,1)*x;%累積計(jì)算整個(gè)網(wǎng)絡(luò)的功率損耗endangle(1)=0.0;fori=NB:-1:1kf=FT(i,2);%回代始節(jié)點(diǎn)號(hào)kt=FT(i,1);%回代終節(jié)點(diǎn)號(hào)dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2))/V(kf);%計(jì)算支路電壓損耗的縱分量dv1dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1))/V(kf);%計(jì)算支路電壓損耗的橫分量dv2V2(kt)=sqrt((V(kf)-dv1)^2+(dv2)^2);%計(jì)算支路末端電壓/kVangle(kt)=atand(dv2/(V(kf)-dv1))+angle(kf);%計(jì)算支路電壓相角（單位度）endmaxd=abs(V2(2)-V(2));V2(1)=V(1);fori=3:1:NNifabs(V2(i)-V(i))>maxd;maxd=abs(V2(i)-V(i));endendmaxdk=k+1PQ1%潮流分布即支路首端潮流MVAV=V2%節(jié)點(diǎn)電壓模計(jì)算結(jié)果kVangle%節(jié)點(diǎn)電壓角度計(jì)算結(jié)果單位度PL%網(wǎng)損單位MWendclear輸入文件清單：function[PQ,FT,RX]=case115()PQ=[%節(jié)點(diǎn)電壓有功無功10.4 0 010.0 0.0342 0.030110.0 0.0693 0.064210.0 0.0845 0.076310.0 0.0295 0.026110.0 0.0474 0.040910.0 0.1176 0.095710.0 0.0946 0.085710.00.09160.085910.00.02710.022910.00.06960.064310.00.06760.057910.00.02980.0242];FT=[%首端末端5 413 412 114 31031137 69832628221];RX=[%RX4.5245.043.5213.9662.4362.8661.1451.284.144.6961.3281.7632.7452.9652.2372.7562.3562.5410.8561.143.7434.2513.3673.685];計(jì)算結(jié)果：k=8PQ1=0.02960.02620.02990.02430.06780.05810.14440.12730.02720.02300.13790.12310.11830.09640.09200.08640.38520.34450.16610.13790.18910.17500.81160.7278V=10.40009.87959.69919.66529.63779.84929.78729.73269.68709.67649.65789.62359.6444angle=00.30110.40110.42500.44380.34310.39420.38990.42660.42090.44980.48230.4452PL=0.0488數(shù)據(jù)分析：經(jīng)過8次迭代，各節(jié)點(diǎn)的電壓誤差均在0.0001kv以內(nèi)，計(jì)算到此結(jié)束。七、配電網(wǎng)潮流計(jì)算的要點(diǎn)需要注意系統(tǒng)中所有母線的電壓是否在允許的范圍以內(nèi)。系統(tǒng)中各種元件(線路、變壓器等)是否會(huì)出現(xiàn)過負(fù)荷。在進(jìn)行電壓和功率分布計(jì)算以前，先要對(duì)網(wǎng)絡(luò)等值電路作些簡化處理。確立前推回代的計(jì)算次序。在開始編寫程序之前，要了解matlab仿真軟件，包括其程序的編寫，程序的運(yùn)行，我們都要清楚的知道為了提高計(jì)算精度，重復(fù)計(jì)算，給出一個(gè)容許誤差作為計(jì)算的收斂的判據(jù)。八、自我總結(jié)此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過精密的設(shè)計(jì)和計(jì)算。在進(jìn)行課題設(shè)計(jì)的過程中，加深了我對(duì)潮流計(jì)算的認(rèn)識(shí)，尤其是對(duì)牛頓拉夫遜潮流計(jì)算的求解思路有了比較透徹的理解。在此次程序設(shè)計(jì)中，我在MATLAB編程，13節(jié)點(diǎn)配電網(wǎng)潮流計(jì)算，WORD文檔的編輯方面均有所提高，之前學(xué)過的潮流計(jì)算僅限于書本，通過完成這個(gè)作業(yè)使我對(duì)這塊知識(shí)了解得更加全面，但也暴露出了一些問題：理論知識(shí)還是不足，基本知識(shí)不充足，做起這個(gè)大作業(yè)難度還是很大的，對(duì)MATLAB的性能和特點(diǎn)還不能有一個(gè)全面的把握，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，從而達(dá)到一個(gè)新的層次。在過程設(shè)計(jì)中我發(fā)現(xiàn)了一些要點(diǎn)，要求得支路的首端功率，從與葉節(jié)點(diǎn)聯(lián)接的支路開始，該支路的末端功率即等于葉節(jié)點(diǎn)功率，利用這個(gè)功率和對(duì)應(yīng)的節(jié)點(diǎn)電壓計(jì)算支路功率損耗，才能求首段功率。還有在迭代計(jì)算開始之前，先要確定支路的計(jì)算順序問題，即確定好是深度優(yōu)先還是廣度優(yōu)先。要求得各支路終結(jié)點(diǎn)的電壓的話利用上述所得的支路首端功率和本步驟剛算出的本支路始節(jié)點(diǎn)的電壓，從電源點(diǎn)開始逐條支路進(jìn)行計(jì)算。而且在此次課程設(shè)計(jì)中，我發(fā)現(xiàn)了自己的基礎(chǔ)知識(shí)有很多的不足。這些基礎(chǔ)的缺乏給我的設(shè)計(jì)計(jì)劃造成了不小的障礙。在這個(gè)過程中，我明白了，只要用心去做，認(rèn)直去做，持之以恒，就會(huì)有新的發(fā)現(xiàn)，有意外的收獲。雖然很多錯(cuò)誤被克服了，但是系統(tǒng)中難免還有很多不足之處，有些地方感覺做的不是很正確，希望在以后的學(xué)習(xí)中能夠在同學(xué)們的幫助和老師的指導(dǎo)下做得更好。九、參考文獻(xiàn)[1]何仰贊溫增銀.電力系統(tǒng)分析.華中科技大學(xué)出版社，2002[2]蔡旭暉、劉衛(wèi)國、蔡立燕.

MATLAB基