【正文】 返回支路無功損耗return rchild。}。(*b).U_begin=(*a).U_end。void pass_w(line_part *a,line_part *b,line_part *c){(*a).S_end[0]=(*b).S_begin[0]+(*c).S_begin[0]+(*a).S_org[0]。}。}。(*a).S_end[1]=(*b).S_begin[1]+(*a).S_org[1]。class BinTree{ public:friend void pass_U(line_part*,line_part*,line_part*)。line_part(0,0,0,0,0,0,0,0,0,0,0)。}vall,double u_end=0,double u_begin=0,double r=0, Pe=0,double Xe=0, P0=0,double Q0=0){ egin,r,x,b,Pe,Xe,K,P0,Q0)。t,double vall,double u_end,double u_begin,double r, Pe,double Xe, Q0)。friend void //功率傳遞函數(shù) BinTree(){line_part *aa=new line_part *Getroot(){return root。}。}。}。void insertline_part(line_partdouble x,double b,doubledouble K,double P0,double void PreOrder(line_part *amp。void PostOrder(line_part *t)。t)。t, double vall,double u_end=0,double u_begin=0,double r=0，double x=0,double b=0,double Pe=0,double Xe=0，double K=0,double P0=0,double Q0=0){//插入節(jié)點double ass=tget_val()。} else if(vallget_val()){insertline_part(tlchild, vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。} }。t){if(tget_val()0amp。tget_val()display(tlchild)。ofstream outfile1。p_cost_all +=tget_Pcost()。// 計算并存放各個節(jié)點的電壓相角（始端為零）U=atan2(tU[1],(tget_Ubegin()tU[0]))/*180。//計算并存放各個支路的電壓損耗U_cost=tget_Ubegin()tget_Uend()。的有功損耗P_cost=tget_Pcost()。的無功損耗Q_cost=tget_Qcost()?！眊et_Pbegin()get_Qbegin()//計算并存放各個支“get_Pend()get_Qend()//計算并存放各個支//計算并存放各個支路 “//計算并存放各個支路 ”}。t){//先序遍歷二叉樹修改電壓0)get_lchild()get_val()get_rchild()get_val()get_rchild()get_val()if(tget_val()0amp。tget_val()get_lchild()get_val()0amp。t { if(tget_rchild()get_val()0amp。t{tUend()。tlchildUend()。}else{tUend()。tlchildUend()。amp。pass_U2(t,trchild)。}else{} }PreOrder(tlchild)。} }。amp。PostOrder(trchild)。amp。amp。trchildSbegin()。tSbegin()。pass_w2(t,tlchild)。} } else { if(tget_rchild()get_val()0amp。t{trchildSbegin()。tSbegin()。void main(){ofstream outfile。cout”const int M(a)。支路數(shù)據(jù)數(shù)組while(l}}//輸出數(shù)據(jù)流定義 (“e:”)。()。(“e:”)。int l(0)。//支路節(jié)點數(shù)常//輸入流輸入各infileA[l]。l++。(“e:”,ios::ate)。i*i+1],A[11*i+2],A[11*i+3],A[11*i+4]，*i+8],A[11*i+9],A[11*i+10])。BinTree elec。ioutfile()。()?！啊薄啊眔utfileline_part *ee=new line_part()。double aa =eeget_Pbegin()。(”e:“,ios::ate)。bb =aap_cost_all。cc = p_cost_all/aa*100。p_cost_all=。**************************************************************************************************************************************************************習(xí)題33 input: 3 2,35,0,3,2,1,0,0 1,35,35,1,2,0,5,3,1,0,0 3,35,35,2,4,0,2,3,1,0,0習(xí)題34 Input: 最大負(fù)荷 2 1,110,118,0,0,1,0,0 2,11,110,0,40,30,10, 最小負(fù)荷 2 1,110,113,0,0,1,0,0 2,11,110,0,20,15,10,第二篇：電力系統(tǒng)潮流計算南 京 理 工 大 學(xué)《電力系統(tǒng)穩(wěn)態(tài)分析》課程報告姓名XX學(xué) 號： 5*** 自動化學(xué)院 電氣工程基于牛頓拉夫遜法的潮流計算例題編程報學(xué)院(系): 專業(yè): 題目: 任課教師 碩士導(dǎo)師 告楊偉 XX2015年6月10號基于牛頓拉夫遜法的潮流計算例題編程報告摘要：電力系統(tǒng)潮流計算的目的在于：確定電力系統(tǒng)的運行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護(hù)的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟(jì)運行提供分析的基礎(chǔ)。同時基于MATLAB的計算機(jī)算法以雙精度類型進(jìn)行數(shù)據(jù)的存儲和運算, 數(shù)據(jù)精確度高，能進(jìn)行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。例題如下：用牛頓拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。二 牛頓拉夫遜法潮流計算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。電力系統(tǒng)潮流計算，一般來說，各個母線所供負(fù)荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)形成節(jié)點導(dǎo)納矩陣，然后由節(jié)點導(dǎo)納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉(zhuǎn)化為求解非線性方程組的問題了。2 基本步驟和設(shè)計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經(jīng)解決，已有可能列出基本計算步驟并編制流程圖。% 如圖所示1，2，3，4為PQ節(jié)點,5為平衡節(jié)點y=0。y(4,5)=0。y(1,4)=1/(+)。y(2,3)=1/(+)。y(3,4)=1/(+)。endendY=0。endendend% 求節(jié)點導(dǎo)納矩陣中自導(dǎo)納for i=1:5Y(i,i)=sum(y(i,:))。B=imag(Y)。S(2)=。S(4)=。P=real(S)。% 賦初值，U為節(jié)點電壓的幅值，a為節(jié)點電壓的相位角U=ones(1,5)。a=zeros(1,5)。x2=ones(8,1)。while max(x2)1e6for i=1:4for j=1:4H(i,j)=0。M(i,j)=0。oP(i)=0。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)))。endoP(i)=oP(i)+P(i)。endx2=[oP,oQ]39。N(i,j)=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,i)=N(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。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)。endJ=[H,N。% x1為所求△x的列向量% 求節(jié)點電壓新值，準(zhǔn)備下一次迭代for i=1:4oa(i)=x1(i)。endfor i=1:4a(i)=a(i)+oa(i)。endk=k+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)。endS(5)=P(5)+Q(5)*sqrt(1)。四運行結(jié)果節(jié)點導(dǎo)納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結(jié)果分析在這次學(xué)習(xí)和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機(jī)算法中的牛頓拉夫遜法進(jìn)行深入的研讀，弄明白了其原理、計算過程、公式推導(dǎo)以及設(shè)計流程。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。第三篇：電力系統(tǒng)潮流計算程序設(shè)計電力系統(tǒng)潮流計算程序設(shè)計姓名：韋應(yīng)順學(xué)號：2011021052 電力工程學(xué)院牛頓—拉夫遜潮流計算方法具有能夠?qū)⒎蔷€性方程線性化的特點，而使用MATLAB語言是由于MATLAB語言的數(shù)學(xué)邏輯強(qiáng)，易編譯。ph=input(‘n請輸入平衡母線的節(jié)點號：ph=’)； B1=input（‘n請輸入支路信號：B1=’）。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。f=zeros（1，n）；no=2*ph=1； Y(q,q)=Y(q,q)+1./B1(i3)+B1(i4)/2；End for i=1：nG=real(Y)。p=A(i1)；Y(p p)=1./A(i2)。End f(i)=imag(B2(i3))。p=B1(i1)。p=B1(i2)。Q=imag(S)。J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)。t=0。max(abs(Dfi)ipt=t+1。[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)。for i=1:n hh(i)=conj(Y(ph,i)*v(i))。B2(ph,1)=S(ph)。jd=angle(v)*180/p。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。c(i)=v(B1(i2))*a(i)j*B1(i4)*v(B1(i2))^2/2。printcut(result1,S,b,c,result2)。D(i)=0。D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j)。Q1=C(i)*f(i)D(i)*e(i)。IfB2(i4)=2 p=2*i1。p=p+1。Df(p)=P(i)P1。Df(p)=Q(i)Q1。If ph=n Df(noJ=[]。DX1=DX。if ph=n DX(no)=0。For i=(no+2):(x1+2)DX(i)=DX1(i2)。End k=0。For i=1:2:x K=k+1。De(k)=DX(i+1)。X2=G(i,j)*e(i)+B(i,j)*f(i)。switch B2(i4)X4=0。continue q=2*j1。for j=1:n m=p+1。J=ph J(m,q)=X3。q=q+1。J(p,q)=X2。J(m,q)=X4。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。X3=0。X4=0。P=2*i1。q=2*j1。J(p,q)=X1。j=jph m=p+1。J(m,q)=X3。q=q+1。J(p,q)=X2。J(m,q)=X4。end q=2*j1。end m=p+1。if ph=n q=q+1。J(p,q)=X2。J(m,q)=X4。End J(:,no)=[]。cess dyn鋤ic Voltage s