【正文】 方法，越來越廣泛的應用到電力系統(tǒng)潮流計算中。國內(nèi)專家學者對解決此問題也進行了許多有益的探討。與現(xiàn)代內(nèi)點算法相結合的求解過程表現(xiàn)出了良好的收斂性和快速性，計算結果準確、可靠，計算各種病態(tài)系統(tǒng)均可良好的收斂，基于L1范數(shù)的數(shù)學規(guī)劃模型將傳統(tǒng)電力系統(tǒng)潮流的直接迭代求解轉化為對一簡單規(guī)劃問題的求解后，對系統(tǒng)運行中各部分的控制可更加簡便?；诜柗治龇椒ǖ某绷饔嬎惴椒ㄍㄟ^建立電力網(wǎng)絡的拓撲模型生成拓撲網(wǎng)絡的全部樹和2樹，應用網(wǎng)絡的k樹樹支導納乘積對電力網(wǎng)絡的節(jié)點電壓方程進行拓撲求解，進而得出所求變量(即各節(jié)點電壓)的符號表達式(即關于元件參數(shù)符號的顯式表達式)。另外，傳統(tǒng)的潮流計算方法都是純“數(shù)值計算”，利用這些方法計算出來的結果是數(shù)字而不是函數(shù)，它們的特點是逐點進行完整的數(shù)值計算，因此不可避免地存在收斂性問題、冗余項對消問題、計算機有效字長效應問題和相近數(shù)值求差時發(fā)生的浮點運算誤差問題。%v=PV節(jié)點個數(shù)P=[ 0]。%V=PV節(jié)點的U值E=[1 1 1 ]39。 0 0。 15 0。Y=G+j*B。%F(x)矩陣F1=zeros(n1,n1)。%△P FX(n+i1)=Q(i)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,w)=0。F1(n+i1,w)=2*E(i)。endendend end X=inv(F1)*(FX)。fprintf(39。fprintf(39。enddisp(39。for i=1:n1 fprintf(39。%d號節(jié)點電壓的相角度數(shù)為39。end fprintf(39。運行結果：運行結果復制如下：第1次迭代后的節(jié)點電壓分別為:第1次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:00第2次迭代后的節(jié)點電壓分別為:第2次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第3次迭代后的節(jié)點電壓分別為:第3次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第4次迭代后的節(jié)點電壓分別為:第4次迭代后功率偏差△P △Q電壓偏差△V的平方分別為: *收斂后的節(jié)點電壓用極坐標表示為: 1號節(jié)點電壓的幅值為:2號節(jié)點電壓的幅值為:3號節(jié)點電壓的幅值為:4號節(jié)點電壓的幅值為:平衡節(jié)點的功率 第三篇：潮流計算電力系統(tǒng)大作業(yè)(C++編寫)程序設計所涉及二叉樹解釋任意一棵樹或一個森林都能唯一地對應一棵二叉樹，由此而編寫本程序。程序說明文檔******************************************************************************************************************************************************* 本程序測試使用方法：在E盤根目錄下建立輸入文件：輸入文件名：。********************************************************************************************************************************************************* 程序中變量定義說明：類中定義的變量class line_part{//定義支路類double U[2]。//支路末端電壓double U_begin。//支路首段功率：S_begin[0]有功，S_begin[1]無功double S0[2]。}。}。//輸出數(shù)據(jù)流定ifstream infile。}。//支路首段電壓double X[3]。//變壓器勵磁損耗：S0[0]有功，S0[1]無功double S_org[2]。//變壓器變比public: line_part(){ val=0。X[1]=0。S_begin[0]=0。S0[1]=0。}line_part(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){val=vall。X[1]=x。S_begin[0]=0。S0[1]=Q0。}。//功率傳遞函數(shù)friend void pass_U2(line_part*a,line_part*b)。}。//U_heng即是△u2 0])/U_begin。Uend=k*U_end。U[1]=U_zong。amp。}。}。}。}。}。}。void pass_U2(line_part *a,line_part *b){(*b).U_begin=(*a).U_end。}。root = aa。}。void PostOrder(){PostOrder(root)。line_part *root。void display(line_part *amp。if(t==0||tget_val()t=new line_part(vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。/////////////////////////////////////////////////////////////////////////////////////////////////////////////// void BinTree::display(line_part *amp。display(trchild)。double U。outfile1”outfile1“outfile1”outfile1“outfile1”double P_cost。outfile1”()。amp。amp。trchildUend()。} } else { if(tget_rchild()get_val()0amp。trchildUend()。void BinTree::PostOrder(line_part *t){//后序遍歷二叉樹修改潮流if(tget_val()0amp。if(tget_lchild()get_val()0amp。t{tlchildSbegin()。}else{tlchildSbegin()。amp。}else{}}。量double *A=new double[10*M]。ifstream infile。char b。}。如下“()。()。ee=()。double bb。outfile()。它的任務是根據(jù)給定的運行條件和網(wǎng)路結構確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡中的功率分布以及功率損耗等。利用電子計算機進行潮流計算從20世紀50年代中期就已經(jīng)開始。同時，為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài)，也需要進行大量而快速的潮流計算。這個方法的原理比較簡單，要求的數(shù)字計算機的內(nèi)存量也比較小，適應當時的電子數(shù)字計算機制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計算人員轉向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個元素進行計算，因此，每次迭代的計算量很大。為了克服阻抗法在內(nèi)存和速度方面的缺點，后來發(fā)展了以阻抗矩陣為基礎的分塊阻抗法。解決電力系統(tǒng)潮流計算問題是以導納矩陣為基礎的，因此，只要在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計算效率。牛頓法的特點是將非線性方程線性化。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡、模糊算法也逐漸被引入潮流計算。，運用所學專業(yè)知識，提出一種合理高效的潮流計算算法。結合C語言，編寫一套適用的程序完成潮流計算。進度安排：第1周: 收集相關參考資料和相關文獻。第5周: 實習第6周: 寫實習報告第7周: 確定一種計算方法。六、參考文獻：[1] Tankut Yalcinoz, Onur Ko168。第五篇：電力系統(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)濟運行提供分析的基礎。例題如下：用牛頓拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。電力系統(tǒng)潮流計算，一般來說，各個母線所供負荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據(jù)網(wǎng)絡結構形成節(jié)點導納矩陣，然后由節(jié)點導納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉化為求解非線性方程組的問題了。% 如圖所示1，2，3，4為PQ節(jié)點,5為平衡節(jié)點y=0。y(1,4)=1/(+)。y(3,4)=1/(+)。endendend% 求節(jié)點導納矩陣中自導納for i=1:5Y(i,i)=sum(y(i,:))。S(2)=。P=real(S)。a=zeros(1,5)。while max(x2)1e6for i=1:4for j=1:4H(i,j)=0。oP(i)=0。endoP(i)=oP(i)+P(i)。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(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)))。endJ=[H,N。endfor i=1:4a(i)=a(i)+oa(i)。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)。四運行結果節(jié)點導納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結果分析在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設計流程。