【正文】 _end[0]=(*b).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*)。//電壓傳遞函數pass_w(line_part*,line_part*,line_part*)。line_part(0,0,0,0,0,0,0,0,0,0,0)。root = aa。}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)。private:*amp。t,double vall,double u_end,double u_begin,double r, Pe,double Xe, Q0)。}。friend void //功率傳遞函數 BinTree(){line_part *aa=new line_part *Getroot(){return root。} void insertline_part(doubledouble x=0,double b=0,doubledouble K=0,double insertline_part(root,vall,u_end,u_b} void PreOrder(){PreOrder(root)。}。void PostOrder(){PostOrder(root)。}。void display(){display(root)。}。line_part *root。void insertline_part(line_partdouble x,double b,doubledouble K,double P0,double void PreOrder(line_part *amp。t)。void PostOrder(line_part *t)。void display(line_part *amp。t)。void BinTree::insertline_part(line_part *amp。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()。if(t==0||tget_val()t=new line_part(vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。} else if(vallget_val()){insertline_part(tlchild, vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。} else {insertline_part(trchild, vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。} }。/////////////////////////////////////////////////////////////////////////////////////////////////////////////// void BinTree::display(line_part *amp。t){if(tget_val()0amp。amp。tget_val()display(tlchild)。display(trchild)。ofstream outfile1。(“e:”,ios::ate)。p_cost_all +=tget_Pcost()。double U。// 計算并存放各個節(jié)點的電壓相角（始端為零）U=atan2(tU[1],(tget_Ubegin()tU[0]))/*180。outfile1get_val()”outfile1””get_Ubegin()outfile1””get_Uend()//計算并存放各個支double U_cost。//計算并存放各個支路的電壓損耗U_cost=tget_Ubegin()tget_Uend()。outfile1”outfile1“outfile1”outfile1“outfile1”double P_cost。的有功損耗P_cost=tget_Pcost()。outfile1“double Q_cost。的無功損耗Q_cost=tget_Qcost()。outfile1”()。”get_Pbegin()get_Qbegin()//計算并存放各個支“get_Pend()get_Qend()//計算并存放各個支//計算并存放各個支路 “//計算并存放各個支路 ”}。/////////////////////////////////////////////////////////////////////////////////////////////////////////////////void BinTree::PreOrder(line_part *amp。t){//先序遍歷二叉樹修改電壓0)get_lchild()get_val()get_rchild()get_val()get_rchild()get_val()if(tget_val()0amp。amp。tget_val()get_lchild()get_val()0amp。amp。t { if(tget_rchild()get_val()0amp。amp。t{tUend()。pass_U(t,tlchild,trchild)。tlchildUend()。trchildUend()。}else{tUend()。pass_U2(t,tlchild)。tlchildUend()。} } else { if(tget_rchild()get_val()0amp。amp。t{tUend()。pass_U2(t,trchild)。trchildUend()。}else{} }PreOrder(tlchild)。PreOrder(trchild)。} }。void BinTree::PostOrder(line_part *t){//后序遍歷二叉樹修改潮流if(tget_val()0amp。amp。tget_val()get_lchild()get_val()get_rchild()get_val()get_rchild()get_val(){ PostOrder(tlchild)。PostOrder(trchild)。if(tget_lchild()get_val()0amp。amp。t { if(tget_rchild()get_val()0amp。amp。t{tlchildSbegin()。trchildSbegin()。pass_w(t,tlchild,trchild)。tSbegin()。}else{tlchildSbegin()。pass_w2(t,tlchild)。tSbegin()。} } else { if(tget_rchild()get_val()0amp。amp。t{trchildSbegin()。pass_w2(t,trchild)。tSbegin()。}else{}}。void main(){ofstream outfile。//輸入數據流定義infilea。cout”const int M(a)。量double *A=new double[10*M]。支路數據數組while(l}}//輸出數據流定義 (“e:”)。()。()。ifstream infile。(“e:”)。int a。int l(0)。char b。//支路節(jié)點數常//輸入流輸入各infileA[l]。(b)。l++。}。(“e:”,ios::ate)。outfile”i=0。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。for(int//添加節(jié)點(A[11*i+0],A[11 A[11*i+5],A[11*i+6],A[11*i+7],A[11} for(i = 0。ioutfile()。()。()。(“e:”,ios::ate)?！啊薄啊眔utfileline_part *ee=new line_part()。ee=()。double aa =eeget_Pbegin()。outfile()。(”e:“,ios::ate)。double bb。bb =aap_cost_all。outfileoutfiledouble cc。cc = p_cost_all/aa*100。outfile()。p_cost_all=。}}。**************************************************************************************************************************************************************習題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習題34 Input: 最大負荷 2 1,110,118,0,0,1,0,0 2,11,110,0,40,30,10, 最小負荷 2 1,110,113,0,0,1,0,0 2,11,110,0,20,15,10,第二篇：電力系統(tǒng)潮流計算南 京 理 工 大 學《電力系統(tǒng)穩(wěn)態(tài)分析》課程報告姓名XX學 號： 5*** 自動化學院 電氣工程基于牛頓拉夫遜法的潮流計算例題編程報學院(系): 專業(yè): 題目: 任課教師 碩士導師 告楊偉 XX2015年6月10號基于牛頓拉夫遜法的潮流計算例題編程報告摘要：電力系統(tǒng)潮流計算的目的在于：確定電力系統(tǒng)的運行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護的整定提供依據、為電力系統(tǒng)的穩(wěn)定計算提供初值、為電力系統(tǒng)規(guī)劃和經濟運行提供分析的基礎。潮流計算的計算機算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計算原理較簡單、計算過程也不復雜，而且由于人們引入泰勒級數和非線性代數方程等在算法里從而進一步提高了算法的收斂性和計算速度。同時基于MATLAB的計算機算法以雙精度類型進行數據的存儲和運算, 數據精確度高，能進行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。一 研究內容通過一道例題來認真分析牛頓拉夫遜法的原理和方法（采用極坐標形式的牛拉法），同時掌握潮流計算計算機算法的相關知識，能看懂并初步使用MATLAB軟件進行編程，培養(yǎng)自己電力系統(tǒng)潮流計算機算法編程能力。例題如下：用牛頓拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。計算精度要求各節(jié)點電壓修正量不大于106。二 牛頓拉夫遜法潮流計算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。牛頓迭代法其最大優(yōu)點是在方程f(x)=0的單根附近時誤差將呈平方減少，而且該法還可以用來求方程的重根、復根。電力系統(tǒng)潮流計算，一般來說，各個母線所供負荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據網絡結構形成節(jié)點導納矩陣，然后由節(jié)點導納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉化為求解非線性方程組的問題了。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導，得出對應的雅可比矩陣，給未知節(jié)點賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點電壓不平衡量（未知的）構成了誤差方程，解誤差方程，得到節(jié)點電壓不平衡量，節(jié)點電壓加上節(jié)點電壓不平衡量構成節(jié)點電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。2 基本步驟和設計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經解決，已有可能列出基本計算步驟并編制流程圖。由課本總結基本步驟如下：1)形成節(jié)點導納矩陣Y；2)設各節(jié)點電壓的初值，如果是直角坐標的話設電壓的實部e和虛部f；如果是極坐標的話則設電壓的幅值U和相角a；3)將各個節(jié)點電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點電壓的變化量，即修正量； 5)計算各個節(jié)點電壓的新值，即修正后的值；6)利用新值從第（3）步開始進入下一