電力系統(tǒng)潮流計算的matlab輔助程序設(shè)計-潮流計算程序

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

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

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

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

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

6、程序lfybuslfybus，busoutbusout和lineflowlineflow。decoupledecouple是快速解耦法對實際電力系統(tǒng)潮流計算開發(fā)的程序，同高斯法和牛頓法一樣需要用到三個程序：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