【正文】 J(p,q)=X2。j=jph m=p+1。J(m,q)=X4。De(k)=DX(i+1)。Df(p)=P(i)P1。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。Q=imag(S)。}第五篇：電力系統(tǒng)潮流計(jì)算程序設(shè)計(jì)電力系統(tǒng)潮流計(jì)算程序設(shè)計(jì)姓名：韋應(yīng)順學(xué)號(hào)：2011021052 電力工程學(xué)院牛頓—拉夫遜潮流計(jì)算方法具有能夠?qū)⒎蔷€性方程線性化的特點(diǎn)，而使用MATLAB語言是由于MATLAB語言的數(shù)學(xué)邏輯強(qiáng)，易編譯。}//全網(wǎng)數(shù)據(jù)// int high=1,low=1。QZL[i][j]=(x1*y2x2*y1)。}}//平衡節(jié)點(diǎn)//for(i=0。//調(diào)用Jacobi函數(shù)和高斯函數(shù)//for(int t=1。K[i][j]=1。ZLs++。B[i].Pd=m3[i]。Generator*G。Ja[2*i+1][2*j+1]=2*f[i]。jif(i!=j){if(type[i]==1){Ja[2*i][2*j]=(YG[i][j]*e[i]+YB[i][j]*f[i])。//deltaP deltaQ deltaV^2// void fun1(double YG[][50],double YB[][50],double e[],double f[],int type[],int N,double W[],double P[],double Q[],double V[]){ double dP=0,dQ=0,dV=0。//節(jié)點(diǎn)號(hào)類型負(fù)荷有功負(fù)荷無功母線數(shù)據(jù)（類型 1＝PV節(jié)點(diǎn)，2＝PQ節(jié)點(diǎn)，3＝平衡節(jié)點(diǎn)）struct BUS { int busno。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。oQ(i)=0。S(3)=。% 輸入原始數(shù)據(jù)，求節(jié)點(diǎn)導(dǎo)納矩陣y(1,2)=1/(+)。data[ii][7]=。data[ii][6]=。//參數(shù)分別為前節(jié)點(diǎn)的有功、無功、電壓data[ii][3]=Tm[j].GetPl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][5]=。x!=1amp。data[ii][3]=Pd。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。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。amp。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。statemp=line[i].start。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。int count=1。Complex aa,bb,cc,dd,B。//計(jì)算節(jié)點(diǎn)功率for(i=0。class Complex//定義復(fù)數(shù)類 { public: double real。//取極坐標(biāo)Complex(){RecPolar=0。=**。=1。}class gauss { public: static void gauss_slove(double **a,double *x,int NN)。pivrow[k]=i。}}} } x[n1]=a[n1][n]/a[n1][n1]。0節(jié)點(diǎn)25 。data[ii][0]=。coutcinx。data[ii][2]=。data[ii][6]=。Tf[k].Kk=f。data[ii][4]=data[ii1][4]+data[ii1][0]+data[ii][6]。2 基本步驟和設(shè)計(jì)流程圖形成了雅克比矩陣并建立了修正方程式，運(yùn)用牛頓拉夫遜法計(jì)算潮流的核心問題已經(jīng)解決，已有可能列出基本計(jì)算步驟并編制流程圖。B=imag(Y)。M(i,j)=0。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)。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。float B。} } //Jacobi矩陣// void Jacobi(double YG[][50],double YB[][50],double e[50],double f[50],int type[50],int N ,double Ja[100][101]){int i,j。if(type[i]==1){Ja[2*i+1][2*j]=b[i]+YB[i][i]*e[i]YG[i][i]*f[i]。jvoid main(){ ifstream fin。for(i=0。i++){finm1[i]。iYB[i][j]=0。iC[2*i]=1。for(i=0。y1=f[i]/L[k].kf[j]。pr[j][i]=I2*L[k].R。ki=L[k].busi1。V(i)=B2(i3)。V=abs(v)。Df(p)=P(i)P1。[x,y]=size(DX)。X2=G(i,j)*e(i)+B(i,j)*f(i)。J(m,q)=X4。end J(m,q)=X3。J(:,no)=[]。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。p=2*i1。%該子程序是用來求取jacci矩陣for i=1:n X3=0。end end end Df=Df’。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。[De,Di]=hxf(J,D,F,ph,n,no)。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個(gè)端點(diǎn) %第二列存貯支路的另一個(gè)端點(diǎn) %第三列存貯支路阻抗%第四列存貯支路的對(duì)地導(dǎo)納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號(hào)B2=input(‘n請(qǐng)輸入節(jié)點(diǎn)信息：B2=’)； %第一列為電源側(cè)的功率 %第二列為負(fù)荷側(cè)的功率 %第三列為該點(diǎn)的電壓值%第四列為該點(diǎn)的類型：1為PQ，2為PV節(jié)點(diǎn)，3為平衡節(jié)點(diǎn) A=input（‘n請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗：A=’）； ip=input（‘n請(qǐng)輸入修正值：ip=’）； %ip為修正值)。kVh){Vh=sqrt(e[k]*e[k]+f[k]*f[k])。pr[i][j]=I2*L[k].R。ja+=(YG[i][j]*e[j]YB[i][j]*f[j])。f[i]=C[2*i+1]。ifor(j=i。iL[i].busi=m1[i]。i++){finm1[i]。if(!fin){coutgetchar()。for(k=0。Ja[2*i+1][2*j+1]=Ja[2*i][2*j]。idouble A=0,B=0。float Qd。endfor i=1:4a(i)=a(i)+oa(i)。endoP(i)=oP(i)+P(i)。P=real(S)。y(1,4)=1/(+)。if(fabs(data[ii][2]data[ii1][2])第三篇：電力系統(tǒng)潮流計(jì)算南 京 理 工 大 學(xué)《電力系統(tǒng)穩(wěn)態(tài)分析》課程報(bào)告姓名XX學(xué) 號(hào)： 5*** 自動(dòng)化學(xué)院 電氣工程基于牛頓拉夫遜法的潮流計(jì)算例題編程報(bào)學(xué)院(系): 專業(yè): 題目: 任課教師 碩士導(dǎo)師 告楊偉 XX2015年6月10號(hào)基于牛頓拉夫遜法的潮流計(jì)算例題編程報(bào)告摘要：電力系統(tǒng)潮流計(jì)算的目的在于：確定電力系統(tǒng)的運(yùn)行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護(hù)的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計(jì)算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟(jì)運(yùn)行提供分析的基礎(chǔ)。//參數(shù)要改data[ii][3]=Tf[k].GetPt(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。j++。i++。amp。data[ii][6]=。0節(jié)點(diǎn)28 。 }x[i]=(a[i][n]sum)/a[i][i]。a[pivrow[k]][j]=a[k][j]。pivrow=new int[n]。}voidComplex::PrintfComplex(Complex c1)//按直角坐標(biāo)輸出 { if(==0){ (6)。return Node。return Node。//求兩個(gè)復(fù)數(shù)積static Complex divideComplex(Complex c1,Complex c2)。coutt(bus[i].)。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。//輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N)。X=new Complex *[N]。int _tmain(int argc, _TCHAR* argv[]){int i。//計(jì)算節(jié)點(diǎn)功率for(i=0。Complex aa,bb,cc,dd,B。int count=1。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。//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].)。//動(dòng)態(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。double data[300][8]。