【正文】 part *amp。} else if(vallget_val()){insertline_part(tlchild, vall,u_end,u_begin,r,x,b,Pe,Xe,K,P0,Q0)。ofstream outfile1。的有功損耗P_cost=tget_Pcost()。tget_val()get_lchild()get_val()0amp。}else{tUend()。}else{} }PreOrder(tlchild)。amp。pass_w2(t,tlchild)。void main(){ofstream outfile。(“e:”)。(“e:”,ios::ate)。()。bb =aap_cost_all。電力系統(tǒng)潮流計算的結(jié)果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎(chǔ)。因此，潮流計算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運算。阻抗法改善了電力系統(tǒng)潮流計算問題的收斂性，解決了導(dǎo)納法無法解決的一些系統(tǒng)的潮流計算，在當(dāng)時獲得了廣泛的應(yīng)用，曾為我國電力系統(tǒng)設(shè)計、運行和研究作出了很大的貢獻(xiàn)。自從20世紀(jì)60年代中期采用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、計算速度方面都超過了阻抗法，成為直到目前仍被廣泛采用的方法。但是，到目前為止這些新的模型和算法還不能取代牛頓法和PQ分解法的地位。選取一典型模型進(jìn)行驗證，試驗程序是否可靠。第8周: 提出一種合理的程序設(shè)計方法。潮流計算的計算機(jī)算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計算原理較簡單、計算過程也不復(fù)雜，而且由于人們引入泰勒級數(shù)和非線性代數(shù)方程等在算法里從而進(jìn)一步提高了算法的收斂性和計算速度。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導(dǎo)，得出對應(yīng)的雅可比矩陣，給未知節(jié)點賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點電壓不平衡量，節(jié)點電壓加上節(jié)點電壓不平衡量構(gòu)成節(jié)點電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。y(1,5)=1/(+)。endY% Y為導(dǎo)納矩陣G=real(Y)。Q=imag(S)。N(i,j)=0。oQ(i)=oQ(i)+Q(i)。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)))。U(i)=U(i)+oU(i)。牛頓拉夫遜法是求解非線性方程的迭代過程，其計算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。S% 求節(jié)點注入電流I=Y*U39。oU(i)=x1(i+4)*U(i)。endendend% 當(dāng)i=j時，求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)))。oQ(i)=oQ(i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。k=0。S(5)=0。% 求節(jié)點導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。y(1,3)=1/(+)。牛頓迭代法其最大優(yōu)點是在方程f(x)=0的單根附近時誤差將呈平方減少，而且該法還可以用來求方程的重根、復(fù)根。北京：中國電力出版社，2001。第4周: 熟悉電力系統(tǒng)潮流計算的相關(guān)理論及計算機(jī)語言。四、本課題研究方案確定一種計算方法，如牛頓拉夫遜法。近20多年來，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進(jìn)牛頓法和PQ分解法進(jìn)行的。牛頓法是數(shù)學(xué)中求解非線性方程式的典型方法，有較好的收斂性。這就需要較大的內(nèi)存量。在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。**************************************************************************************************************************************************************習(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,第四篇：潮流計算畢業(yè)論文科學(xué)技術(shù)學(xué)院畢業(yè)設(shè)計（論文）開題報告題目：電力系統(tǒng)潮流分析計算機(jī)輔助設(shè)計學(xué) 科 部：信息學(xué)科部專業(yè)：電氣工程及其自動化班級：電氣082班學(xué)號：7022808070姓名：黃義軍指導(dǎo)教師：劉愛國填表日期：2011 年 月 日一、選題的依據(jù)及意義：電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本電氣計算。(”e:“,ios::ate)。ioutfile()。l++。()。tSbegin()。tSbegin()。PostOrder(trchild)。pass_U2(t,trchild)。tlchildUend()。t){//先序遍歷二叉樹修改電壓0)get_lchild()get_val()get_rchild()get_val()get_rchild()get_val()if(tget_val()0amp。//計算并存放各個支路的電壓損耗U_cost=tget_Ubegin()tget_Uend()。tget_val()display(tlchild)。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()。}。t,double vall,double u_end,double u_begin,double r, Pe,double Xe, Q0)。(*a).S_end[1]=(*b).S_begin[1]+(*a).S_org[1]。(*b).U_begin=(*a).U_end。double get_Qcost(){//返回支路無功損耗return S_begin[1]S_end[1]。double get_Pend(){//返回支路末端有功return S_end[0]。}else{ if(val0amp。//支路首段功率計算函數(shù) double Uend,I2。sqr(U_begin)*X[2]/2sqr(U_end)*X[2]/2。lchild = rchild= NULL。S_end[1]=S_org[1]=Xe。lchild = rchild= NULL。S_end[1]=S_org[1]=0。//支路電壓降落：U[0]電壓降落橫分量，U[1]電壓降落縱分量double k。//支路末端電壓double U_begin。主函數(shù)中定義的數(shù)據(jù)ofstream outfile。//支路的后繼兩個支路 }class BinTree{//定義樹類void PreOrder(){PreOrder(root)。//支路排序double U_end。例如習(xí)題34：++104SDST=++j30MVA20+j15MVA程序二叉樹結(jié)構(gòu)示意圖：本二叉樹中節(jié)點1即支路1為題目中節(jié)點1與節(jié)點2之間的部分； 本二叉樹中節(jié)點2即支路2為題目中節(jié)點2與節(jié)點3之間的部分； 對于習(xí)題33：程序二叉樹結(jié)構(gòu)示意圖：本二叉樹中節(jié)點2即支路2為題目中節(jié)點1與節(jié)點2之間的部分； 本二叉樹中節(jié)點1即支路1為題目中節(jié)點2與節(jié)點3之間的部分； ：如下多支路網(wǎng)絡(luò)：對于三節(jié)點網(wǎng)絡(luò)需先進(jìn)行以下網(wǎng)絡(luò)處理轉(zhuǎn)化為標(biāo)準(zhǔn)二叉樹，而后進(jìn)行計算；圖中三角形表明該項阻抗為零，為純導(dǎo)線，并進(jìn)行相應(yīng)參數(shù)補(bǔ)充進(jìn)行計算。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。EF=E+j*F。disp(E+j*F)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。%△Qend end for i=m+1:n1 %PV節(jié)點 FX(i)=PP(im)。%總的節(jié)點數(shù)FX=ones(2*n2,1)。B=[ 00。%PP=PV節(jié)點的P值 V=[]。這些特點將在電力系統(tǒng)的在線計算、靜態(tài)安全分析等領(lǐng)域發(fā)揮明顯優(yōu)勢。和過去的模型相比，該模型非常的簡潔、直觀，易于編程。它有一個特別的優(yōu)點是容許根據(jù)兩個子優(yōu)化問題各自的特性而采用不同的求解方法。0(20)(=x)253。0239。除上述之外，還有其它一些用途不同的特殊潮流問題，如諧波潮流、動態(tài)潮流等，在此不一一列舉。采用對稱分量坐標(biāo)，將系統(tǒng)各量分為正、負(fù)、零序分量，并對系統(tǒng)中不對稱元件的序分量之間的耦合，通過加電流源補(bǔ)償?shù)姆椒ㄊ怪怦?。Uj(GijsinqijBijcosqij)(9)j206。ij206。潮流算法多種多樣，但一般要滿足四個基本要求：(i)可靠收斂；(ii)計算速2n個變量作為已知量而預(yù)先給以指定?；趯?dǎo)納矩陣的高斯—塞德爾法是電力系統(tǒng)中最早得到應(yīng)用的潮流計算方法。首先簡單回顧了潮流計算的發(fā)展歷史，對當(dāng)前基于計算機(jī)的各種潮流算法的原理及其優(yōu)缺點，作了簡要介紹和比較，并介紹了它們采用的一些特別技術(shù)及程序設(shè)計技巧；接著簡要分析了三種新型的潮流計算方法的計算原理及優(yōu)缺點，它們分別是基于人工智能的潮流計算方法、基于L1范數(shù)和現(xiàn)代內(nèi)點理論的電力系統(tǒng)潮流計算方法、基于符號分析的潮流計算方法等。它根據(jù)給定系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)及運行條件來確定整個系統(tǒng)的運行狀態(tài)：主要是各節(jié)點電壓(幅值和相角)，網(wǎng)絡(luò)中功率分布及功率損耗等。同時，60年代初開始出現(xiàn)運用非線性規(guī)劃的最優(yōu)潮流算法。對應(yīng)于這些節(jié)點，分別對其注入有功、無功功率，有功功率及電壓模值以及電壓模值和相角加以指定；并且對平衡節(jié)點來說，其電壓相角一般作為系統(tǒng)電壓相角的基準(zhǔn)(即θ=0)。(Gijfj+Bijej)(7)j206。這兩種形式的潮流方程統(tǒng)稱為節(jié)點功率方程，是牛頓拉夫遜等潮流算法所采用的主要數(shù)學(xué)模型。后者則將交流系統(tǒng)潮流方程組和直流系統(tǒng)的方程組分開來求解，求解直流系統(tǒng)方程組時各換流站的交流母線電壓由交流系統(tǒng)潮流的解算