1、電力系統(tǒng)潮流計(jì)算的電力系統(tǒng)潮流計(jì)算的MATLABMATLAB輔助程序設(shè)計(jì)輔助程序設(shè)計(jì)潮流計(jì)算，通常指負(fù)荷潮流，是電力系統(tǒng)分析和設(shè)計(jì)的主要組成部分，對(duì)系統(tǒng)規(guī)劃、安全運(yùn)行、經(jīng)濟(jì)調(diào)度和電力公司的功率交換非常重要。此外，潮流計(jì)算還是其它電力系統(tǒng)分析的基礎(chǔ)，比如暫態(tài)穩(wěn)定，突發(fā)事件處理等。現(xiàn)代電力系統(tǒng)潮流計(jì)算的方法主要:高斯法、牛頓法、快速解耦法和MATLAB的M語(yǔ)言編寫的MATPOWER4.1MATPOWER4.1，這里主要介紹高斯法、牛頓法和快

2、速解耦法。高斯法的程序是lfgausslfgauss，其與lfybuslfybus、busoutbusout和lineflowlineflow程序聯(lián)合使用求解潮流功率。lfybuslfybus、busoutbusout和lineflowlineflow程序也可與牛頓法的lfnewtonlfnewton程序和快速解耦法的decoupledecouple程序聯(lián)合使用。(讀者可以到MATPOWER主頁(yè)下載MATPOWER4.1然后將其解壓到M

3、ATLAB目錄下，即可使用該軟件進(jìn)行潮流計(jì)算)一、高斯賽德?tīng)柗ǔ绷饔?jì)算使用的程序：高斯賽德法的具體使用方法讀者可參考后面的實(shí)例，這里僅介紹各程序的編寫格式:lfgausslfgauss：該程序是用高斯法對(duì)實(shí)際電力系統(tǒng)進(jìn)行潮流計(jì)算，需要用到busdatabusdata和linedatalinedata兩個(gè)文件。程序設(shè)計(jì)為輸入負(fù)荷和發(fā)電機(jī)的有功MW和無(wú)功Mvar，以及節(jié)點(diǎn)電壓標(biāo)幺值和相角的角度值。根據(jù)所選復(fù)功率為基準(zhǔn)值將負(fù)荷和發(fā)電機(jī)的功率轉(zhuǎn)

4、換為標(biāo)幺值。對(duì)于PV節(jié)點(diǎn)，如發(fā)電機(jī)節(jié)點(diǎn)，要提供一個(gè)無(wú)功功率限定值。當(dāng)給定電壓過(guò)高或過(guò)低時(shí)，無(wú)功功率可能超出功率限定值。在幾次迭代之后(高斯塞德?tīng)柕鸀?0次)，需要檢查一次發(fā)電機(jī)節(jié)點(diǎn)的無(wú)功出力，如果接近限定值，電壓幅值進(jìn)行上下5%的調(diào)整，使得無(wú)功保持在限定值內(nèi)。lfybuslfybus：這個(gè)程序需要輸入線路參數(shù)、變壓器參數(shù)以及變壓器分接頭參數(shù)。并將這些參數(shù)放在名為linedatalinedata的文件中。這個(gè)程序?qū)⒆杩罐D(zhuǎn)換為導(dǎo)納，并得到

5、節(jié)點(diǎn)導(dǎo)納矩陣。busoutbusout：該程序以表格形式輸出結(jié)果，節(jié)點(diǎn)輸出包括電壓幅值和相角，發(fā)電機(jī)和負(fù)荷的有功和無(wú)功功率，以及并聯(lián)電容器或電抗器的有功和無(wú)功功率。lineflowlineflow：該程序輸出線路的相關(guān)數(shù)據(jù)，程序設(shè)計(jì)輸出流入線路終端的有功和無(wú)功的功率、線損以及節(jié)點(diǎn)功率，還包含整個(gè)系統(tǒng)的有功和無(wú)功損耗。lfnewtonlfnewton是牛頓拉夫遜法對(duì)實(shí)際電力系統(tǒng)潮流計(jì)算開(kāi)發(fā)的程序，數(shù)據(jù)準(zhǔn)備和程序格式和高斯賽德?tīng)柗ㄒ粯?，包?/p>

6、程序lfybuslfybus，busoutbusout和lineflowlineflow。decoupledecouple是快速解耦法對(duì)實(shí)際電力系統(tǒng)潮流計(jì)算開(kāi)發(fā)的程序，同高斯法和牛頓法一樣需要用到三個(gè)程序：lfybuslfybus、busoutbusout、lineflowlineflow。Pk=[]P=[]Qk=[]Q=[]S=[]V=[]fk=1:nbusn=busdata(k1)kb(n)=busdata(k2)Vm(n)=bus

7、data(k3)delta(n)=busdata(k4)Pd(n)=busdata(k5)Qd(n)=busdata(k6)Pg(n)=busdata(k7)Qg(n)=busdata(k8)Qmin(n)=busdata(k9)Qmax(n)=busdata(k10)Qsh(n)=busdata(k11)ifVm(n)=0Vm(n)=1.0V(n)=1j0elsedelta(n)=pi180delta(n)V(n)=Vm(n)(cos

8、(delta(n))jsin(delta(n)))P(n)=(Pg(n)Pd(n))basemvaQ(n)=(Qg(n)Qd(n)Qsh(n))basemvaS(n)=P(n)jQ(n)endDV(n)=0endnum=0AcurBus=0converge=1Vc=zeros(nbus1)jzeros(nbus1)Sc=zeros(nbus1)jzeros(nbus1)whileexist(accel)~=1accel=1.3endwh