【正文】 obi(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é)點導納矩陣經過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結果分析在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設計流程。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 基本步驟和設計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經解決，已有可能列出基本計算步驟并編制流程圖。電力系統(tǒng)潮流計算，一般來說，各個母線所供負荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據(jù)網絡結構形成節(jié)點導納矩陣，然后由節(jié)點導納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉化為求解非線性方程組的問題了。二 牛頓拉夫遜法潮流計算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。例題如下：用牛頓拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。同時基于MATLAB的計算機算法以雙精度類型進行數(shù)據(jù)的存儲和運算, 數(shù)據(jù)精確度高，能進行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請按任意鍵繼續(xù)... 20to19 20to10 22to21 26to25 27to25 28to27 高等電力系統(tǒng)分析 IEEE30節(jié)點潮流程序班級：電研114班姓名：王大偉學號：2201100151第三篇：電力系統(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ī)劃和經濟運行提供分析的基礎。。00 節(jié)點8 節(jié)點9 0 節(jié)點10 0 節(jié)點11 節(jié)點12 0 節(jié)點13 節(jié)點14 0 節(jié)點15 0 節(jié)點16 0 節(jié)點17 0 節(jié)點18 0 節(jié)點19 0 節(jié)點20 0 節(jié)點21 0 節(jié)點22 0 節(jié)點23 0 節(jié)點24 0 。00 0 節(jié)點7。0 0 節(jié)點6。0 00 節(jié)點5。00 0 節(jié)點4。節(jié)點3。====節(jié)點電壓===============發(fā)電機發(fā)出功率======節(jié)點10。x[pivcol[k][0]]=al。k){al=x[pivcol[k][1]]。} for(k=(n2)。j{sum +=a[i][j]*x[j]。i){sum=0。//解方程for(i=(n2)。j{a[i][j]=aik*a[k][j]/a[k][k]。i{aik=a[i][k]。a[i][k]=al。i{al=a[i][pivcol[k][1]]。a[k][j]=al。j{al=a[pivrow[k]][j]。exit(0)。//行pivcol[k][1]=j。j{if(pivot{pivot=fabs(a[i][j])。//列n*2矩陣for(i=k。//行pivcol[k][0]=k。ipivot= fabs(a[k][k])。pivcol=new int *[n]。n=NN。void gauss::gauss_slove(double **a,double *x,int NN){int n,i,j,k,*pivrow,**pivcol。static void gauss_output()。=0。cout”。(8)。cout”。(8)。}doubleComplex::getCompleximage(Complex c1)//取虛部 {return 。return Node。=0。}Complex Complex::getinverse(Complex c1)//取倒數(shù) { Complex Node。=。} Complex Complex::getconj(Complex c1)//取共軛 {Complex Node。=(r1)。} Complex Complex::ComDivRea(Complex c1,double r1)//復數(shù)除數(shù) { Complex Node。=(**)/(pow(,2)+pow(,2))。} Complex Complex::divideComplex(Complex c1,Complex c2)//復數(shù)除法 {Complex Node。=*+*。}Complex Complex::productComplex(Complex c1,Complex c2)//復數(shù)乘法 {Complex Node。=。}Complex Complex::subComplex(Complex c1,Complex c2)//復數(shù)減法 {Complex Node。=+。}Complex Complex::CaddC(Complex c1,Complex c2)//復數(shù)加法 {Complex Node。=1。=sqrt(*+*)。} }。//將復數(shù)復零static Complex Rec2Polar(Complex c1)。//顯示一個復數(shù)static void PrintfmultiComplex(Complex C,int N)。//求一個復數(shù)實部static double getCompleximage(Complex c1)。//求一個復數(shù)共軛static Complex getinverse(Complex c1)。//求兩個復數(shù)商static Complex ComDivRea(Complex c1,double r2)。//求兩個復數(shù)差static Complex productComplex(Complex c1,Complex c2)。//0表示直角坐標，1表示極坐標static Complex CaddC(Complex c1,Complex c2)。double image。}include using namespace std。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。endtemp=line[i].end。i{int statemp,endtemp。coutt(bus[i].)。i{coutComplex::PrintfComplex(bus[i].V)。bus[i].V=Complex::Rec2Polar(bus[i].V)。i{bus[i]. = bus[i]. + bus[i].。}coutBus::JisuanNodeScal(X,bus,N)。dd=Complex::CaddC(aa,cc)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)