【正文】 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)潮流計算程序設計電力系統(tǒng)潮流計算程序設計姓名：韋應順學號：2011021052 電力工程學院牛頓—拉夫遜潮流計算方法具有能夠將非線性方程線性化的特點，而使用MATLAB語言是由于MATLAB語言的數(shù)學邏輯強，易編譯。}//全網(wǎng)數(shù)據(jù)// int high=1,low=1。QZL[i][j]=(x1*y2x2*y1)。}}//平衡節(jié)點//for(i=0。//調用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é)點號類型負荷有功負荷無功母線數(shù)據(jù)（類型 1＝PV節(jié)點，2＝PQ節(jié)點，3＝平衡節(jié)點）struct BUS { int busno。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。oQ(i)=0。S(3)=。% 輸入原始數(shù)據(jù)，求節(jié)點導納矩陣y(1,2)=1/(+)。data[ii][7]=。data[ii][6]=。//參數(shù)分別為前節(jié)點的有功、無功、電壓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é)點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)//復數(shù)除法 {Complex Node。=sqrt(*+*)。//0表示直角坐標，1表示極坐標static Complex CaddC(Complex c1,Complex c2)。bus[i].V=Complex::Rec2Polar(bus[i].V)。statemp=line[i].start。//計算節(jié)點功率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)分配支路結構體Line::ScanfLineData(line)。amp。//輸出節(jié)點導納矩陣int NN=(N1)*2。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。statemp=line[i].start。//動態(tài)分配結點結構體for(int i=0。int count=1。Complex aa,bb,cc,dd,B。//計算節(jié)點功率for(i=0。class Complex//定義復數(shù)類 { public: double real。//取極坐標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é)點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 基本步驟和設計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經(jīng)解決，已有可能列出基本計算步驟并編制流程圖。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ù)得到所需要的精度值即準確值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)。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個端點 %第二列存貯支路的另一個端點 %第三列存貯支路阻抗%第四列存貯支路的對地導納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號B2=input(‘n請輸入節(jié)點信息：B2=’)； %第一列為電源側的功率 %第二列為負荷側的功率 %第三列為該點的電壓值%第四列為該點的類型：1為PQ，2為PV節(jié)點，3為平衡節(jié)點 A=input（‘n請輸入節(jié)點號及對地阻抗：A=’）； ip=input（‘n請輸入修正值：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)潮流計算南 京 理 工 大 學《電力系統(tǒng)穩(wěn)態(tài)分析》課程報告姓名XX學 號： 5*** 自動化學院 電氣工程基于牛頓拉夫遜法的潮流計算例題編程報學院(系): 專業(yè): 題目: 任課教師 碩士導師 告楊偉 XX2015年6月10號基于牛頓拉夫遜法的潮流計算例題編程報告摘要：電力系統(tǒng)潮流計算的目的在于：確定電力系統(tǒng)的運行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟運行提供分析的基礎。//參數(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é)點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)//按直角坐標輸出 { if(==0){ (6)。return Node。return Node。//求兩個復數(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é)點功率差值int icon=wehcon1(bus,N)。X=new Complex *[N]。int _tmain(int argc, _TCHAR* argv[]){int i。//計算節(jié)點功率for(i=0。Complex aa,bb,cc,dd,B。int count=1。//動態(tài)分配結點結構體for(int i=0。//i作為整個程序的循環(huán)變量int N=Bus::ScanfBusNo()。for(i=0。//whether converbence看迭代是否結束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。coutt(bus[i].)。//動態(tài)分配支路結構體Line::ScanfLineData(line)。for(i=0。goto LOOP。dd=Complex::CaddC(aa,cc)。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。//顯示一個復數(shù)static void PrintfmultiComplex(Complex C,int N)。=。=。cout”。//列n*2矩陣for(i=k。i{aik=a[i][k]。====節(jié)點電壓===============發(fā)電機發(fā)出功率======節(jié)點10。double data[300][8]。