【正文】 _begin,double r, Pe,double Xe, Q0)。}。}。void PostOrder(line_part *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é)點(diǎn)double ass=tget_val()。} }。tget_val()display(tlchild)。p_cost_all +=tget_Pcost()。//計(jì)算并存放各個(gè)支路的電壓損耗U_cost=tget_Ubegin()tget_Uend()。的無功損耗Q_cost=tget_Qcost()。t){//先序遍歷二叉樹修改電壓0)get_lchild()get_val()get_rchild()get_val()get_rchild()get_val()if(tget_val()0amp。t { if(tget_rchild()get_val()0amp。tlchildUend()。tlchildUend()。pass_U2(t,trchild)。} }。PostOrder(trchild)。amp。tSbegin()。} } else { if(tget_rchild()get_val()0amp。tSbegin()。cout”const int M(a)。()。int l(0)。l++。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])。ioutfile()?！啊薄啊眔utfileline_part *ee=new line_part()。(”e:“,ios::ate)。cc = p_cost_all/aa*100。**************************************************************************************************************************************************************習(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,第四篇：潮流計(jì)算畢業(yè)論文科學(xué)技術(shù)學(xué)院畢業(yè)設(shè)計(jì)（論文）開題報(bào)告題目：電力系統(tǒng)潮流分析計(jì)算機(jī)輔助設(shè)計(jì)學(xué) 科 部：信息學(xué)科部專業(yè)：電氣工程及其自動(dòng)化班級(jí)：電氣082班學(xué)號(hào)：7022808070姓名：黃義軍指導(dǎo)教師：劉愛國填表日期：2011 年 月 日一、選題的依據(jù)及意義：電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種基本電氣計(jì)算?，F(xiàn)在的潮流算法都以計(jì)算機(jī)的應(yīng)用為前提。在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。二、國內(nèi)外研究現(xiàn)狀及發(fā)展趨勢（含文獻(xiàn)綜述）：在用數(shù)字計(jì)算機(jī)求解電力系統(tǒng)潮流問題的開始階段，人們普遍采用以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的高斯賽德爾迭代法（一下簡稱導(dǎo)納法）。這就需要較大的內(nèi)存量。當(dāng)系統(tǒng)不斷擴(kuò)大時(shí)，這些缺點(diǎn)就更加突出。牛頓法是數(shù)學(xué)中求解非線性方程式的典型方法，有較好的收斂性。PQ分解法在計(jì)算速度方面有顯著的提高，迅速得到了推廣。近20多年來，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進(jìn)牛頓法和PQ分解法進(jìn)行的。三、本課題研究內(nèi)容。四、本課題研究方案確定一種計(jì)算方法，如牛頓拉夫遜法。從而提高電力系統(tǒng)運(yùn)行的經(jīng)濟(jì)性。第4周: 熟悉電力系統(tǒng)潮流計(jì)算的相關(guān)理論及計(jì)算機(jī)語言。第10周: 將整體程序模塊化，并定義出每個(gè)模塊的功能。北京：中國電力出版社，2001。一 研究內(nèi)容通過一道例題來認(rèn)真分析牛頓拉夫遜法的原理和方法（采用極坐標(biāo)形式的牛拉法），同時(shí)掌握潮流計(jì)算計(jì)算機(jī)算法的相關(guān)知識(shí)，能看懂并初步使用MATLAB軟件進(jìn)行編程，培養(yǎng)自己電力系統(tǒng)潮流計(jì)算機(jī)算法編程能力。牛頓迭代法其最大優(yōu)點(diǎn)是在方程f(x)=0的單根附近時(shí)誤差將呈平方減少，而且該法還可以用來求方程的重根、復(fù)根。由課本總結(jié)基本步驟如下：1)形成節(jié)點(diǎn)導(dǎo)納矩陣Y；2)設(shè)各節(jié)點(diǎn)電壓的初值，如果是直角坐標(biāo)的話設(shè)電壓的實(shí)部e和虛部f；如果是極坐標(biāo)的話則設(shè)電壓的幅值U和相角a；3)將各個(gè)節(jié)點(diǎn)電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點(diǎn)電壓的變化量，即修正量； 5)計(jì)算各個(gè)節(jié)點(diǎn)電壓的新值，即修正后的值；6)利用新值從第（3）步開始進(jìn)入下一次迭代，直至達(dá)到精度退出循環(huán)； 7)計(jì)算平衡節(jié)點(diǎn)的功率和線路功率，輸出最后計(jì)算結(jié)果； ① 公式推導(dǎo)② 流程圖三matlab編程代碼clear。y(1,3)=1/(+)。y(2,5)=1/(+)。% 求節(jié)點(diǎn)導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。% 輸入原始節(jié)點(diǎn)的給定注入功率S(1)=+。S(5)=0。U(5)=。k=0。L(i,j)=0。oQ(i)=oQ(i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。% x2為不平衡量列向量% 求雅克比矩陣% 當(dāng)i~=j時(shí)，求H,N,M,Lfor i=1:4for j=1:4if i~=jH(i,j)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。endendend% 當(dāng)i=j時(shí)，求H,N,M,Lfor i=1:4for j=1:5if i~=j H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。oU(i)=x1(i+4)*U(i)。endk,U,a% 求節(jié)點(diǎn)注入功率i=5。S% 求節(jié)點(diǎn)注入電流I=Y*U39。六 結(jié)論通過這個(gè)任務(wù)，自己在matlab編程，潮流計(jì)算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識(shí)儲(chǔ)備不足，對(duì)matlab的性能和特點(diǎn)還不能有一個(gè)全面的把握，對(duì)word軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，達(dá)到一個(gè)新的層次。牛頓拉夫遜法是求解非線性方程的迭代過程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))+Q(i)。U(i)=U(i)+oU(i)。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。M(i,i)=M(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。L(i,j)=H(i,j)。oQ(i)=oQ(i)+Q(i)。oQ(i)=0。N(i,j)=0。x1=ones(8,1)。Q=imag(S)。S(3)=。endY% Y為導(dǎo)納矩陣G=real(Y)。for i=1:5for j=i:5y(j,i)=y(i,j)。y(1,5)=1/(+)。% 輸入原始數(shù)據(jù)，求節(jié)點(diǎn)導(dǎo)納矩陣y(1,2)=1/(+)。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對(duì)功率平衡方程求偏導(dǎo)，得出對(duì)應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點(diǎn)電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點(diǎn)電壓不平衡量，節(jié)點(diǎn)電壓加上節(jié)點(diǎn)電壓不平衡量構(gòu)成節(jié)點(diǎn)電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計(jì)算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于106。潮流計(jì)算的計(jì)算機(jī)算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計(jì)算原理較簡單、計(jì)算過程也不復(fù)雜，而且由于人們引入泰勒級(jí)數(shù)和非線性代數(shù)方程等在算法里從而進(jìn)一步提高了算法的收斂性和計(jì)算速度。 multiobjective optimizationmethod to environmental economic ,29（1）:4250 [2] ， economic dispatch model and Power and Energy , 29(1):113120 [3] 何仰贊，：華中科技大學(xué)出版社，2002 [4] 王錫凡，方萬良，：科學(xué)出版社，2003 [5] 宋文南，李樹鴻，：天津大學(xué)出版社，1990 [6] 王晶，翁國慶，：西安電子科技大學(xué)出版社，2008.[7] :水利電力出版社，1987.[8] 周全仁，:湖南科學(xué)技術(shù)出版社，1983.[9] 許主平，周少武，鄒軍安。第8周: 提出一種合理的程序設(shè)計(jì)方法。第2周: 總結(jié)整理資料，熟習(xí)課題。選取一典型模型進(jìn)行驗(yàn)證，試驗(yàn)程序是否可靠。但是，到目前為止這些新的模型和算法還不能取代牛頓法和PQ分解法的地位。20世紀(jì)70年代后期，有人提出采用更精確的模型，即將泰勒級(jí)數(shù)的高階項(xiàng)也包括進(jìn)來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。自從20世紀(jì)60年代中期采用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、計(jì)算速度方面都超過了阻抗法，成為直到目前仍被廣泛采用的方法。這個(gè)方法把一個(gè)大系統(tǒng)分割為幾個(gè)小的地區(qū)系統(tǒng)，在計(jì)算機(jī)內(nèi)只需存儲(chǔ)各個(gè)地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡(luò)線的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時(shí)也提高了計(jì)算速度。阻抗法改善了電力系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納法無法解決的一些系統(tǒng)的潮流計(jì)算，在當(dāng)時(shí)獲得了廣泛的應(yīng)用，曾為我國電力系統(tǒng)設(shè)計(jì)、運(yùn)行和研究作出了很大的貢獻(xiàn)。20世紀(jì)60年代初，數(shù)字計(jì)算機(jī)已經(jīng)發(fā)展到第二代，計(jì)算機(jī)的內(nèi)存和計(jì)算速度發(fā)生了很大的飛躍，從而為阻抗法的采用創(chuàng)造了條件。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。此后，潮流計(jì)算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對(duì)潮流計(jì)算的一些基本要求進(jìn)行的。電力系統(tǒng)潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。p_cost_all=。bb =aap_cost_all。double aa =eeget_Pbegin()。()。BinTree elec。(“e:”,ios::ate)。//支路節(jié)點(diǎn)數(shù)常//輸入流輸入各infileA[l]。(