【正文】 bi(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)W[2*i+1]=dQ。W[2*i]=dP。}dV=V[i]*V[i]e[i]*e[i]f[i]*f[i]。jA+=YG[i][j]*e[j]YB[i][j]*f[j]。idouble A=0,B=0。int i,j。}。float B。float R。//支路信息 節(jié)點I 節(jié)點J R X B/2 k struct Line { int busi。float Vg。//發(fā)電機數(shù)據(jù) 節(jié)點號有功發(fā)電電壓幅值 struct Generator { int busno。float Qd。int type。第四篇：電力系統(tǒng)通用潮流計算C語言程序include include include include using namespace std。利用x(*)=x(k+1)+DX(k+1)進行多次迭代，通過迭代判據(jù)得到所需要的精度值即準確值x(*)。四運行結(jié)果節(jié)點導納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結(jié)果分析在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設(shè)計流程。endS(5)=P(5)+Q(5)*sqrt(1)。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)。endk=k+1。endfor i=1:4a(i)=a(i)+oa(i)。% x1為所求△x的列向量% 求節(jié)點電壓新值，準備下一次迭代for i=1:4oa(i)=x1(i)。endJ=[H,N。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)。N(i,i)=N(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。M(i,j)=N(i,j)。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。endx2=[oP,oQ]39。endoP(i)=oP(i)+P(i)。endend% 求有功、無功功率不平衡量for i=1:4for j=1:5oP(i)=oP(i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。oP(i)=0。M(i,j)=0。while max(x2)1e6for i=1:4for j=1:4H(i,j)=0。x2=ones(8,1)。a=zeros(1,5)。% 賦初值，U為節(jié)點電壓的幅值，a為節(jié)點電壓的相位角U=ones(1,5)。P=real(S)。S(4)=。S(2)=。B=imag(Y)。endendend% 求節(jié)點導納矩陣中自導納for i=1:5Y(i,i)=sum(y(i,:))。endendY=0。y(3,4)=1/(+)。y(2,3)=1/(+)。y(1,4)=1/(+)。y(4,5)=0。% 如圖所示1，2，3，4為PQ節(jié)點,5為平衡節(jié)點y=0。2 基本步驟和設(shè)計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經(jīng)解決，已有可能列出基本計算步驟并編制流程圖。電力系統(tǒng)潮流計算，一般來說，各個母線所供負荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)形成節(jié)點導納矩陣，然后由節(jié)點導納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉(zhuǎn)化為求解非線性方程組的問題了。二 牛頓拉夫遜法潮流計算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。例題如下：用牛頓拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。同時基于MATLAB的計算機算法以雙精度類型進行數(shù)據(jù)的存儲和運算, 數(shù)據(jù)精確度高，能進行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。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)濟運行提供分析的基礎(chǔ)。int jj。data[ii][6]=。data[ii][4]=data[ii1][4]+data[ii1][0]+data[ii][6]。} data[ii][2]=U1。data[ii][0]=。ii++。data[ii][5]=Tf[k].Kk。//參數(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])。data[ii][7]=。data[ii][1]=Tf[k].GetQ0()。Tf[k].Kk=f。Tf[k].P0=d。Tf[k].Pk=b。int k=0。}if(x==3){coutcinabcdef。j++。data[ii][4]=Tm[j].GetQl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][2]=Tm[j].GetUl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][6]=。data[ii][0]=。Tm[j].Xl=b。int j=0。}if(x==2){coutcinabc。i++。data[ii][6]=。data[ii][4]=N[i].Qn。data[ii][2]=。data[ii][0]=。N[i].Pn=a。Node N[100]。}if(x==0)break。amp。amp。amp。coutcinx。double f。double d。double b。ii++。data[ii][6]=。data[ii][4]=Qd。data[ii][2]=Ud。data[ii][0]=。//分別為0前面有功，1前面無功，2壓降，3中間有功，4中間無功,5變比,6后面有功，7后面無功int ii。coutUd。0 20節(jié) 點30 。。0節(jié)點29 。0節(jié)點28 。節(jié) 點27 。0節(jié)點26 。0節(jié)點25 。節(jié)點2。}coutcout}====節(jié)點功率計算值==== icon=1進行第2次迭代 節(jié)點電壓修正量================= ====節(jié)點功率計算值==== icon=1進行第3次迭代 節(jié)點電壓修正量================= ====節(jié)點功率計算值==== icon=0,迭代結(jié)束。x[pivcol[k][1]]=x[pivcol[k][0]]。k=0。 }x[i]=(a[i][n]sum)/a[i][i]。for(j=(i+1)。i=0。}}} } x[n1]=a[n1][n]/a[n1][n1]。for(j=k。}}if(k!=(n1))//將矩陣化為上三角形式{for(i=(k+1)。a[i][pivcol[k][1]]=a[i][k]。}}if(pivcol[k][1]!=k)//列變換{for(i=0。a[pivrow[k]][j]=a[k][j]。}if(pivrow[k]!=k)//行變換{for(j=k。//列}}}if(pivot{coutgetchar()。pivrow[k]=i。i{for(j=k。pivcol[k][1]=k。pivrow[k]=k。for(i=0。pivrow=new int[n]。double eps,pivot,sum,aik,al。}。}class gauss { public: static void gauss_slove(double **a,double *x,int NN)。} else {coutComplex::zeroComplex(Complex c1)//清零 { =0。(ios::left)。(6)。(ios::right)。}voidComplex::PrintfComplex(Complex c1)//按直角坐標輸出 { if(==0){ (6)。}double Complex::getComplexReal(Complex c1)//取實部 {return 。Node=(Complex::divideComplex(Node,c1))。=1。return Node。=。return Node。=(r1)。return Node。=(*+*)/(pow(,2)+pow(,2))。return Node。=**。return Node。=。return Node。=+。return Node。=atan2(,)*180/。Complex Complex::Rec2Polar(Complex c1)//極坐標表示 { Complex Node。//取極坐標Complex(){RecPolar=0。//顯示多個復數(shù)static void zeroComplex(Complex c1)。//求一個復數(shù)虛部static void PrintfComplex(Complex c1)。//取倒數(shù)static double getComplexReal(Complex c1)。//除數(shù)static Complex getconj(Complex c1)。//求兩個復數(shù)積static Complex divideComplex(Complex c1,Complex c2)。//求兩個復數(shù)和static Complex subComplex(Co