【正文】 PQ分解法在計(jì)算速度方面有顯著的提高，迅速得到了推廣。在用數(shù)字計(jì)算機(jī)求解電力系統(tǒng)潮流問題的開始階段，人們普遍采用以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的高斯賽德爾迭代法（一下簡稱導(dǎo)納法）。0節(jié)點(diǎn)25 。}}} } x[n1]=a[n1][n]/a[n1][n1]。pivrow[k]=i。}class gauss { public: static void gauss_slove(double **a,double *x,int NN)。=1。=**。//取極坐標(biāo)Complex(){RecPolar=0。class Complex//定義復(fù)數(shù)類 { public: double real。//計(jì)算節(jié)點(diǎn)功率for(i=0。Complex aa,bb,cc,dd,B。int count=1。//動態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。statemp=line[i].start。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。amp。if ph=n q=q+1。J(p,q)=X1。J(p,q)=X2。For i=1:2:x K=k+1。p=p+1。jd=angle(v)*180/p。p=B1(i2)。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。L(i,i)=L(i,i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))endendN(i,i)=N(i,i)2*(U(i))^2*G(i,i)。M(i,j)=0。B=imag(Y)。2 基本步驟和設(shè)計(jì)流程圖形成了雅克比矩陣并建立了修正方程式，運(yùn)用牛頓拉夫遜法計(jì)算潮流的核心問題已經(jīng)解決，已有可能列出基本計(jì)算步驟并編制流程圖。采用牛頓算法，比較容易收斂，迭代次數(shù)設(shè)置為100次。自從60年代中期，牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。此外，電力系統(tǒng)的潮流計(jì)算也是計(jì)算機(jī)系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)，所以潮流計(jì)算是研究電力系統(tǒng)的一種和重要和基礎(chǔ)的計(jì)算。本文計(jì)算采用直角坐標(biāo)形式下的牛頓拉夫遜法，牛頓拉夫遜法有很好的收斂性，但要求有合適的初值。除它自身的重要作用之外，潮流計(jì)算還是網(wǎng)損計(jì)算、靜態(tài)安全分析、暫態(tài)穩(wěn)定計(jì)算、小干擾靜態(tài)穩(wěn)定計(jì)算、短路計(jì)算、靜態(tài)和動態(tài)等值計(jì)算的基礎(chǔ)。各元件中流過的功率，系統(tǒng)的功率損耗等等。克服阻抗法缺點(diǎn)是另一個(gè)途徑是采用牛頓拉夫遜法。圖1電力系統(tǒng)模型 Figure 1 Power System Model單線圖完成后需要進(jìn)行方案定義，通常定義成常規(guī)型。牛頓迭代法其最大優(yōu)點(diǎn)是在方程f(x)=0的單根附近時(shí)誤差將呈平方減少，而且該法還可以用來求方程的重根、復(fù)根。% 求節(jié)點(diǎn)導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。k=0。endendend% 當(dāng)i=j時(shí)，求H,N,M,Lfor i=1:4for j=1:5if i~=j H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。S% 求節(jié)點(diǎn)注入電流I=Y*U39。For i=1:m S(i)=B2(i1)B2(i2)。end S(ph)=sum(hh)*v(ph)。V1=e(i)^2+f(i)^2。End ElseDX=[DX1,0,0]。X1=G(i)*f(i)B(i,j)*e(i)。J(p,q)=X2。end J(p,q)=X1。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。for(i=0。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。coutt(bus[i].)。//動態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。for(i=0。goto LOOP。dd=Complex::CaddC(aa,cc)。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。//顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。=。=。cout”。//列n*2矩陣for(i=k。i{aik=a[i][k]。====節(jié)點(diǎn)電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點(diǎn)10。非線性代數(shù)方程組的解法離不開迭代，因此，潮流計(jì)算方法首先要求它是能可靠的收斂，并給出正確答案。解決電力系統(tǒng)潮流計(jì)算問題是以導(dǎo)納矩陣為基礎(chǔ)的，因此，只要在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計(jì)算效率。另外，為了解決病態(tài)潮流計(jì)算，出現(xiàn)了將潮流計(jì)算表示為一個(gè)無約束非線性規(guī)劃問題的模型，即非線性規(guī)劃潮流算法。阻抗矩陣是滿矩陣，阻抗法要求計(jì)算機(jī)儲存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。0 00 節(jié)點(diǎn)5。i){sum=0。exit(0)。void gauss::gauss_slove(double **a,double *x,int NN){int n,i,j,k,*pivrow,**pivcol。return Node。} Complex Complex::divideComplex(Complex c1,Complex c2)//復(fù)數(shù)除法 {Complex Node。=sqrt(*+*)。//0表示直角坐標(biāo)，1表示極坐標(biāo)static Complex CaddC(Complex c1,Complex c2)。bus[i].V=Complex::Rec2Polar(bus[i].V)。statemp=line[i].start。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。}Bus::ScanfBusData(bus)。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。dd=Complex::CaddC(aa,cc)。goto LOOP。for(i=0。//動態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。J(no:)=[]。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。X4=X1。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。p=p+1。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。【】【】Function tisco %這是一個(gè)電力系統(tǒng)潮流計(jì)算的程序 n=input（‘n請輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請輸入支路數(shù)：m=’)。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。oQ(i)=0。S(3)=。% 輸入原始數(shù)據(jù)，求節(jié)點(diǎn)導(dǎo)納矩陣y(1,2)=1/(+)。Figure 3 flow calculation results輸出報(bào)表如下圖。PSASP簡介及潮流結(jié)仿真研究 PSASP簡介基于電網(wǎng)基礎(chǔ)數(shù)據(jù)庫、固定模型庫以及用戶自定義模型庫的支持，PSASP可進(jìn)行電力系統(tǒng)(輸電、供電和配電系統(tǒng))的各種計(jì)算分析。各點(diǎn)電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等?！峨娏ο到y(tǒng)分析綜合程序》（Power System Analysis SoftwarePackage）PSASP)是由中國電力科學(xué)研究院研發(fā)的電力系統(tǒng)分析程序。關(guān)鍵詞： 電力系統(tǒng)；潮流計(jì)算；psasp 中圖分類號：N031 文獻(xiàn)標(biāo)識碼：APower System Flow Calculation Based On PSASPJiShujie(College of Science, Liaoning Technical University, HuLudao 125105, Liaoning, China)Abstract: Power flow calculation is a calculation of the steadystate operation of the power system, which determines the status of each part of the entire power system according to the given operating conditions and system system flow calculation is very important to the analysis and calculation to the various issues for system planning and operation of research paper conducts simple simulation of power flow calculation by psasp power system analysis software, flowing calculation by Newton words: power system。在電力系統(tǒng)調(diào)度運(yùn)行的多個(gè)領(lǐng)域都涉及到電網(wǎng)潮流計(jì)算。阻抗法改善了系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納無法求解的一些系統(tǒng)的潮流計(jì)算，在60年代獲得了廣泛的應(yīng)用，阻抗法德主要缺點(diǎn)是占用計(jì)算機(jī)內(nèi)存大，每次迭代的計(jì)算量大。靈活多樣的計(jì)算結(jié)果報(bào)表、圖形和圖示化輸出。例題如下：用牛頓拉夫遜法計(jì)算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點(diǎn)，節(jié)點(diǎn)5電壓保持U=，其他四個(gè)節(jié)點(diǎn)分別為PQ節(jié)點(diǎn)，給定的注入功率如圖所示。y(3,4)=1/(+)。a=zeros(1,5)。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))+P(i)。p=A(i1)；Y(p p)=1./A(i2)。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)；J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)；[De,Df]=hxf(J,Df,ph,n,no)。D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j)。if ph=n DX(no)=0。for j=1:n m=p+1。X4=0。J(m,q)=X4。75．[14]IEEE Power Engmeering Socie哆．IEEE std 421．5．2005 IEEE玎ccOmmended practice for excitation syst