【正文】 。(3)解耦最優(yōu)潮流把最優(yōu)潮流的整體最優(yōu)化問題分解為有功優(yōu)化和無功優(yōu)化兩個子優(yōu)化問題。)239。(19)h(u,x)163。其突出優(yōu)點是通過一次計算就提供了系統(tǒng)運行和規(guī)劃的全面信息。(ii)序分量法(或?qū)ΨQ分量法)。iBqni j(8)isji(i=1,2, ,nQi=Ui229。(Gijfj+Bijej)(6)j206。至于用于特殊用途的潮流算法如諧波潮流、適于低壓配電網(wǎng)的潮流算法也得到了較快的發(fā)展。1956年ward等人編制成實用的計算機(jī)潮流計算程序，標(biāo)志著電子計算機(jī)開始在電力系統(tǒng)潮流計算中應(yīng)用。本文對電力系統(tǒng)潮流計算進(jìn)行了綜述。它既是對電力系統(tǒng)規(guī)劃設(shè)計和運行方式的合理性、可靠性及經(jīng)濟(jì)性進(jìn)行定量分析的依據(jù)，又是電力系統(tǒng)靜態(tài)和暫態(tài)穩(wěn)定計算的基礎(chǔ)。60年代末Dom8mel和Tinney提出最優(yōu)潮流的簡化梯度法，70年代有人提出海森矩陣法，80年代SunDl提出最優(yōu)潮流牛頓算法，還可把解耦技術(shù)應(yīng)用于最優(yōu)潮流,從而形成解耦型最優(yōu)潮流牛頓算法，還可把解禍技術(shù)應(yīng)用于最優(yōu)潮流，從而形成解耦型最優(yōu)潮流牛頓算法。交流電力系統(tǒng)中的復(fù)數(shù)電壓變量可以用兩種坐標(biāo)形式來表示Ui=Uiejqi(3)或Ui=ei+jfi(4)而復(fù)數(shù)導(dǎo)納為Yij=Gij+jBij(5)將(3)、式(4)以及式(5)代入以導(dǎo)納矩陣為基礎(chǔ)的式(1)，并將實部與虛部分開，可得到以下兩種形式的潮流方程。ij206。對于以上潮流方程中的有關(guān)運行變量，還可以按其性質(zhì)的不同加以分類，這對于進(jìn)行例如靈敏度分析以及最優(yōu)潮流的研究等都是比較方便的。這樣交替迭代計算，直到收斂。主要方法有最小二乘估計法、支路潮流狀態(tài)估計法、遞推狀態(tài)估計法等。優(yōu)點是原理簡單，設(shè)計簡便。先不考慮不等式約束，構(gòu)造拉格朗日函數(shù)：L(x,l)=f(x)+lTg(x)，定義向量Z=x,l，則應(yīng)用海森矩陣法求最優(yōu)解點Z*的迭代方程為：W△z=d，式中：W，d分別為L對于Z的海森矩陣及梯度向量。因此，即使是在病態(tài)系統(tǒng)的情況下，計算過程不會發(fā)散。根據(jù)電力網(wǎng)絡(luò)在實際運行中的特點，結(jié)合網(wǎng)絡(luò)圖論理論提出了運用符號分析方法求解電力網(wǎng)絡(luò)潮流的新思路，有望克服傳統(tǒng)數(shù)值計算方法在收斂性、冗余項對消、計算機(jī)有效字長效應(yīng)等方面的不足。%m=PQ節(jié)點個數(shù) v=1。 0 0。000]。%最后的節(jié)點電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點FX(i)=P(i)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。end a=a+1。disp(FX)。fprintf(39。disp(PPH)。******************************************************************************************************************************************************* 程序功能說明：本程序可以計算任意長度線型開始網(wǎng)絡(luò)潮流；支持多電壓等級下的計算；可在除供電節(jié)點外任意節(jié)點引出負(fù)載；*******************************************************************************************************************************************************輸入格式說明：以支路為基本單位，按潮流方向輸入數(shù)據(jù)：以下例示意：//支路個數(shù)1,110,118,0,0,1,0,0//線路等效的支路2,11,110,0,40,30,10,//變壓器等效的支路.//按此方式知道輸入所有的支路.//輸入從上到下的順序為潮流在.//線型開式網(wǎng)絡(luò)中的流動方向/*第一行輸入支路個數(shù)，回車第二行至后輸入各個支路參數(shù)，回車分隔不同支路；各行輸入的支路參數(shù)順序是：支路編號，末端電壓，始端電壓，線路等效電阻，線路等效感抗，線路等效容納，末端輸入有功，末端輸入無功，變比，變壓器有功勵磁損耗，變壓器無功勵磁損耗對于線路等效電路：變壓器有功/無功損耗輸入零，變比輸入1；對于變壓器等效電路：所有的參數(shù)均歸算至高壓側(cè)，Rt,Xt對應(yīng)輸入線路等效電阻/感抗的位置，線路等效容抗為零；對于個節(jié)點的引出負(fù)荷：輸入至以此節(jié)點為末節(jié)點的支路的末端輸入有功/無功部分；*/ ******************************************************************************************************************************************************** 輸出文件格式說明：輸出的內(nèi)容包括（1）支路信息：每個支路元件的始端有功、無功和末端有功、無功；有功損耗無功損耗；電壓損耗；（2）全網(wǎng)信息：全網(wǎng)的總電源有功、總負(fù)荷有功、有功損耗、網(wǎng)損率；（3）迭代信息：每次完整迭代后的所有內(nèi)容；具體在輸出文件中都明確標(biāo)出。//支路末端功率：S_end[0]有功，S_end[1]無功double S_begin[2]。//樹的后序遍歷修改潮流void display(){display(root)。double sqr(double x){//平方計算函數(shù) return x*x。//支路首段功率：S_begin[0]有功，S_begin[1]無功double S0[2]。X[0]=0。S0[0]=0。X[0]=r。S0[0]=P0。//電壓傳遞函數(shù)friend void pass_w(line_part*a,line_part*b,line_part*c)。1])/U_begin。q_begin=S_begin[1]+sqr(U_begin)*U_heng=(p_begin*X[0]+q_begin*X[U_zong=(p_begin*X[1]q_begin*X[U_end=sqrt(sqr(U_beginU_heng)+U[0]=U_heng。double get_Ubegin(){//返回支路首段電壓return U_begin。double get_Qend(){//返回支路末端無功return S_end[1]。line_part *get_rchild(){//返回支路無功損耗return rchild。}。line_part(0,0,0,0,0,0,0,0,0,0,0)。}。void PostOrder(line_part *t)。} }。p_cost_all +=tget_Pcost()。的無功損耗Q_cost=tget_Qcost()。t { if(tget_rchild()get_val()0amp。tlchildUend()。} }。amp。} } else { if(tget_rchild()get_val()0amp。cout”const int M(a)。int l(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])。“”“”outfileline_part *ee=new line_part()。cc = p_cost_all/aa*100。現(xiàn)在的潮流算法都以計算機(jī)的應(yīng)用為前提。二、國內(nèi)外研究現(xiàn)狀及發(fā)展趨勢（含文獻(xiàn)綜述）：在用數(shù)字計算機(jī)求解電力系統(tǒng)潮流問題的開始階段，人們普遍采用以節(jié)點導(dǎo)納矩陣為基礎(chǔ)的高斯賽德爾迭代法（一下簡稱導(dǎo)納法）。當(dāng)系統(tǒng)不斷擴(kuò)大時，這些缺點就更加突出。PQ分解法在計算速度方面有顯著的提高，迅速得到了推廣。三、本課題研究內(nèi)容。從而提高電力系統(tǒng)運行的經(jīng)濟(jì)性。第10周: 將整體程序模塊化，并定義出每個模塊的功能。一 研究內(nèi)容通過一道例題來認(rèn)真分析牛頓拉夫遜法的原理和方法（采用極坐標(biāo)形式的牛拉法），同時掌握潮流計算計算機(jī)算法的相關(guān)知識，能看懂并初步使用MATLAB軟件進(jìn)行編程，培養(yǎng)自己電力系統(tǒng)潮流計算機(jī)算法編程能力。由課本總結(jié)基本步驟如下：1)形成節(jié)點導(dǎo)納矩陣Y；2)設(shè)各節(jié)點電壓的初值，如果是直角坐標(biāo)的話設(shè)電壓的實部e和虛部f；如果是極坐標(biāo)的話則設(shè)電壓的幅值U和相角a；3)將各個節(jié)點電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點電壓的變化量，即修正量； 5)計算各個節(jié)點電壓的新值，即修正后的值；6)利用新值從第（3）步開始進(jìn)入下一次迭代，直至達(dá)到精度退出循環(huán)； 7)計算平衡節(jié)點的功率和線路功率，輸出最后計算結(jié)果； ① 公式推導(dǎo)② 流程圖三matlab編程代碼clear。y(2,5)=1/(+)。% 輸入原始節(jié)點的給定注入功率S(1)=+。U(5)=。L(i,j)=0。% x2為不平衡量列向量% 求雅克比矩陣% 當(dāng)i~=j時，求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)))。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。endk,U,a% 求節(jié)點注入功率i=5。六 結(jié)論通過這個任務(wù)，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，達(dá)到一個新的層次。Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))+Q(i)。M,L]% J為雅克比矩陣x1=((inv(J))*x2)。L(i,j)=H(i,j)。oQ(i)=0。x1=ones(8,1)。S(3)=。for i=1:5for j=i:5y(j,i)=y(i,j)。% 輸入原始數(shù)據(jù)，求節(jié)點導(dǎo)納矩陣y(1,2)=1/(+)。計算精度要求各節(jié)點電壓修正量不大于106。 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] 許主平，周少武，鄒軍安。第2周: 總結(jié)整理資料，熟習(xí)課題。20世紀(jì)70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進(jìn)來