1、山東交通學(xué)院電力分析課程設(shè)計(jì)報(bào)告書院(部)別 信息科學(xué)與電氣工程學(xué)院 班 級 xxx 學(xué) 號 xxx 姓 名 xxx 指導(dǎo)教師 xxx 時(shí) 間 2013.11.18-2013.11.29 山東交通學(xué)院 電力系統(tǒng)分析課程設(shè)計(jì)課 程 設(shè) 計(jì) 任 務(wù) 書題 目 電力系統(tǒng)分析課程設(shè)計(jì) 學(xué) 院 信息科學(xué)與電氣工程學(xué)院 專 業(yè) 電氣工程及其自動化 班 級 xxx 學(xué)生姓名 xxx 學(xué) 號 xxx 11 月 19 日至 12 月 2 日 共 2 周指導(dǎo)教師(簽字) 院 長(簽字) 年 月 日目錄第一節(jié) 設(shè)計(jì)任務(wù)2一、設(shè)計(jì)內(nèi)容及要求2二、設(shè)計(jì)原始資料2三、設(shè)計(jì)完成后提交的文件和圖表2四、進(jìn)程安排3五、主要參考

2、資料3第二節(jié) 復(fù)雜網(wǎng)絡(luò)N-R法潮流分析與計(jì)算的設(shè)計(jì)4一、設(shè)計(jì)題目4二、N-R法的求解過程6第三節(jié) 手算潮流計(jì)算8第四節(jié) 計(jì)算機(jī)算法潮流計(jì)算10一、計(jì)算機(jī)程序的實(shí)現(xiàn)10二、計(jì)算機(jī)程序的結(jié)果11總結(jié)18參考文獻(xiàn)19附錄20第一節(jié) 設(shè)計(jì)任務(wù)一、設(shè)計(jì)內(nèi)容及要求復(fù)雜網(wǎng)絡(luò)牛頓拉夫遜法潮流分析與計(jì)算的設(shè)計(jì)電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)中一項(xiàng)最基本的計(jì)算，設(shè)計(jì)內(nèi)容為復(fù)雜網(wǎng)絡(luò)潮流計(jì)算的計(jì)算機(jī)算法牛頓-拉夫遜法。首先，根據(jù)給定的電力系統(tǒng)簡圖，通過手算完成計(jì)算機(jī)算法的兩次迭代過程，從而加深對牛頓-拉夫遜法的理解，有助于計(jì)算機(jī)編程的應(yīng)用。其次，利用計(jì)算機(jī)編程對電力系統(tǒng)穩(wěn)態(tài)運(yùn)行的各參數(shù)進(jìn)行解析和計(jì)算；編程完成復(fù)雜網(wǎng)絡(luò)的節(jié)點(diǎn)

3、導(dǎo)納矩陣的形成；電力系統(tǒng)支路改變、節(jié)點(diǎn)增減的程序變化；編程完成各元件的功率損耗、各段網(wǎng)絡(luò)的電壓損耗、各點(diǎn)電壓、功率大小和方向的計(jì)算。二、設(shè)計(jì)原始資料1、給出一個(gè)六節(jié)點(diǎn)、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)；2、給出一個(gè)五節(jié)點(diǎn)、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)；參數(shù)給定，可以選用直角坐標(biāo)表示的牛拉公式計(jì)算，也可以選用極坐標(biāo)表示的牛拉公式計(jì)算。三、設(shè)計(jì)完成后提交的文件和圖表1計(jì)算說明書部分設(shè)計(jì)報(bào)告和手算潮流的步驟及結(jié)果2圖紙部分：電氣接線圖及等值電路；潮流計(jì)算的計(jì)算機(jī)算法，即程序；運(yùn)算結(jié)果等以圖片的形式附在設(shè)計(jì)報(bào)告中。四、進(jìn)程安排第一周：第一天上午：選題，查資料，制定設(shè)計(jì)方案；第一天下午第五天：復(fù)習(xí)潮流計(jì)

4、算的計(jì)算機(jī)算法，完成給定網(wǎng)絡(luò)的潮流計(jì)算；第二周：第一四天：利用matlab編程完成潮流計(jì)算，并對照手算結(jié)果，分析誤差第五天下午：答辯，交設(shè)計(jì)報(bào)告。五、主要參考資料電力系統(tǒng)分析（第三版） 于永源主編，中國電力出版社，2007年電力系統(tǒng)分析，何仰贊 溫增銀編著，華中科技大學(xué)出版社，2002年版；電力系統(tǒng)分析，韓楨祥主編，浙江大學(xué)出版社，2001年版；電力系統(tǒng)穩(wěn)態(tài)分析，陳珩 編，水利電力出版社；第二節(jié) 復(fù)雜網(wǎng)絡(luò)N-R法潮流分析與計(jì)算的設(shè)計(jì)電力系統(tǒng)的潮流計(jì)算是電力系統(tǒng)分析課程基本計(jì)算的核心部分之一。它既有自身的獨(dú)立意義，又有電力系統(tǒng)規(guī)劃設(shè)計(jì)、運(yùn)行和研究的理論基礎(chǔ)，因此課程設(shè)計(jì)的重要性自不待言。一、設(shè)計(jì)

5、題目1.系統(tǒng)圖的確定選擇六節(jié)點(diǎn)、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)，簡化電力系統(tǒng)圖如圖1所示，等值阻抗圖如圖2所示。運(yùn)用以直角坐標(biāo)表示的牛頓-拉夫遜計(jì)算如圖1所示系統(tǒng)中的潮流分布。計(jì)算精度要求各節(jié)點(diǎn)電壓的誤差或修正量不大于。 圖1 電力系統(tǒng)圖圖2 電力系統(tǒng)等值阻抗圖2.各節(jié)點(diǎn)的初值及阻抗參數(shù)該系統(tǒng)中，節(jié)點(diǎn)為平衡節(jié)點(diǎn)，保持=1.05+j0為定值，節(jié)點(diǎn)為PV節(jié)點(diǎn)，其他四個(gè)節(jié)點(diǎn)都是PQ節(jié)點(diǎn)。給定的注入電壓標(biāo)幺值、線路阻抗標(biāo)幺值、線路阻抗標(biāo)幺值、輸出功率標(biāo)幺值和變壓器變比標(biāo)幺值如圖2所示的注釋。表1 各節(jié)點(diǎn)電壓標(biāo)幺值參數(shù)UUUUUU1.051.001.001.001.001.05表2 線路、變壓器阻抗標(biāo)幺

6、值 線路T1L2L3L4L5T2阻抗j0.030.06+j0.250.04+j0.250.08+j0.300.1+j0.35j0.015表3 節(jié)點(diǎn)輸出功率節(jié)點(diǎn)功率2+j11.8+j0.401.6+j0.83.7+j1.35注：各PQ節(jié)點(diǎn)的電壓取1是為了方便計(jì)算和最后驗(yàn)證程序的正確性。二、N-R法的求解過程1、 給定個(gè)節(jié)點(diǎn)電壓初始值2、 將以上電壓初始值代入下式（1）式，求出修正方程式常數(shù)項(xiàng)向量、 3、 將電壓初始值代入下式（2）式，求出修正方程式中系數(shù)矩陣（雅可比矩陣）的個(gè)元素（為2（n-1）階方陣）。4、 解修正方程式，求出修正量方法：將雅可比矩陣化簡為上三角陣或下三角陣求雅可比矩陣的逆矩陣

7、5、修正各節(jié)點(diǎn)電壓6、將再代入（1）式，求出7、校驗(yàn)是否收斂，其收斂條件為式中，為向量中大分量的絕對值。這個(gè)收斂條件比較直觀，它可以直接顯示出最終結(jié)果的功率誤差8、如果收斂，就進(jìn)一步計(jì)算各段電力線路潮流和平衡節(jié)點(diǎn)功率，并打印出計(jì)算結(jié)果；如果不收斂，轉(zhuǎn)回3步進(jìn)行下一次迭代計(jì)算，直到收斂為止。用圖1和圖2的數(shù)據(jù)和等值網(wǎng)絡(luò)形成節(jié)點(diǎn)導(dǎo)納矩陣第三節(jié) 手算潮流計(jì)算 1. 節(jié)點(diǎn)導(dǎo)納矩陣 由圖1可知，該系統(tǒng)以串聯(lián)支路的阻抗標(biāo)幺值和對地并聯(lián)導(dǎo)納標(biāo)幺值得等值電路如圖2所示。以圖2可得相應(yīng)的節(jié)點(diǎn)導(dǎo)納矩陣。對角線上的元素為：非對角線上的元素為：所以導(dǎo)納矩陣為2.計(jì)算各節(jié)點(diǎn)功率的修正方程的初始值（不平衡量）根據(jù)上述N-

8、R法的求解過程，將各節(jié)點(diǎn)電壓的初始值代入（1）式和（2）式，進(jìn)行潮流計(jì)算。設(shè)各節(jié)點(diǎn)電壓初始值為：e1=1.05 f1=0 e2=1 f2=0 e3=1 f3=0 e4=1 f4=0 e5=1 f5=0 e6=1.05 f6=0求出修正方程式常數(shù)項(xiàng)向量=3、將電壓初始值代入（3）式 求得雅克比矩陣各元素。雅科比矩陣JJ= 4解修正方程，求出修正量如下：=5、修正各節(jié)點(diǎn)電壓 結(jié)果如下： =1.0500 =1.14267 =0.93619 =0.86686 =1.09717 =1.0500 =0 =-0.11248 =-0.57420 =-0.58412 =-0.07762 =0.00275第四節(jié)

9、計(jì)算機(jī)算法潮流計(jì)算一、計(jì)算機(jī)程序的實(shí)現(xiàn)圖2.2 程序流程圖二、計(jì)算機(jī)程序的結(jié)果總結(jié)通過這次課程設(shè)計(jì)，我對牛頓-拉夫遜法算法流程有了深刻的認(rèn)識，復(fù)雜網(wǎng)絡(luò)潮流計(jì)算的計(jì)算機(jī)算法非常實(shí)用，我們應(yīng)該更好的掌握計(jì)算機(jī)算法。在課程設(shè)計(jì)的過程中，我們遇到了很多困難。同學(xué)們相互交流，共同合作，有利于解決問題。同時(shí)讓我感受到了團(tuán)隊(duì)合作的力量。學(xué)會如何跟團(tuán)隊(duì)成員溝通是至關(guān)重要，我們要虛心聽取他人意見，發(fā)現(xiàn)問題及時(shí)向團(tuán)隊(duì)成員反饋，這樣有利于發(fā)現(xiàn)問題并且避免進(jìn)一步的錯(cuò)誤。在設(shè)計(jì)過程中，我通過查閱大量的有關(guān)資料，得出一個(gè)結(jié)論，資料非常多，各個(gè)版本的資料闡述牛頓-拉夫遜法方式也各不相同。所以如何學(xué)會篩選資料對于提高課程設(shè)計(jì)

10、的效率和成功率至關(guān)重要。團(tuán)隊(duì)合作固然重要，但是個(gè)人的獨(dú)立思考也不能忽視。善于獨(dú)立思考有利于樹立自己在團(tuán)隊(duì)合作中的自信心，也為后續(xù)課程設(shè)計(jì)增添活力。這次潮流計(jì)算的手算部分讓我深深的感覺到了，粗心大意是我們最大的敵人。這也為我今后的學(xué)習(xí)和工作敲響了警鐘，我們做事情要保持嚴(yán)謹(jǐn)踏實(shí)的心態(tài)。踏實(shí)認(rèn)真，端正態(tài)度，積極上進(jìn)，是我們的方向。在設(shè)計(jì)的過程中所學(xué)到的東西是這次畢業(yè)設(shè)計(jì)最大的收獲和財(cái)富，使我終身受益！對我們的團(tuán)隊(duì)我要表示深深的感謝!同時(shí)，對給過我?guī)椭乃型瑢W(xué)和老師再次表示忠心的感謝! 參考文獻(xiàn)電力系統(tǒng)分析（第三版） 于永源主編，中國電力出版社，2007年電力系統(tǒng)分析，何仰贊 溫增銀編著，華中科技大

11、學(xué)出版社，2002年版；電力系統(tǒng)分析，韓楨祥主編，浙江大學(xué)出版社，2001年版；電力系統(tǒng)穩(wěn)態(tài)分析，陳珩 編，水利電力出版社。附錄y11=1/0.03j;y12=-1/(1.05*0.03j);y13=0;y14=0;y15=0;y16=0;y21=y12;y22=0.5j+1/(0.06+0.025j)+1/(0.04+0.25j)+1/(1.052)*0.03j);y23=-1/(0.06+0.025j);y24=0;y25=-1/(0.04+0.25j);y26=0;y31=y13;y32=y23;y33=0.5j+1/(0.06+0.025j)+1/(0.08+0.3j);y34=-1/

12、(0.08+0.3j);y35=0; y36=0;y41=y14;y42=y24;y43=y34;y44=0.5j+1/(0.08+0.3j)+1/(0.1+0.35j);y45=-1/(0.1+0.35j);y46=0;y51=y15;y52=y25;y53=y35;y54=y45;y55=0.5j+1/(0.1+0.35j)+1/(0.04+0.25j)+1/(1.052)*0.015j);y56=-1/(1.05*0.015j);y61=0;y62=0;y63=0;y64=0;y65=y56;y66=1/0.015j;%輸入導(dǎo)納參數(shù)YB=y11 y12 y13 y14 y15 y16;

13、y21 y22 y23 y24 y25 y26;y31 y32 y33 y34 y35 y36; y41 y42 y43 y44 y45 y46;y51 y52 y53 y54 y55 y56;y61 y62 y63 y64 y65 y66%形成導(dǎo)納矩陣%*牛頓拉夫遜法潮流計(jì)算* E(2)=1.00;E(3)=1.00;E(4)=1.00;E(5)=1.00;E(6)=1.1;F(2)=0;F(3)=0;F(4)=0;F(5)=0; F(6)=0; %設(shè)定電壓初始值G=real(YB);B=imag(YB); %提取導(dǎo)納實(shí)部,虛部分別為G,BS(2)=-2-1j;S(3)=-1.8-0.4j;

14、S(4)=-1.6-0.8j;S(5)=-3.7-1.3j; S(6)=5;%設(shè)定注入功率初始值P=real(S);Q=imag(S); %提取注入功率實(shí)部,虛部分別為P,Qk=0;precision=0.1; %設(shè)定參數(shù)k為迭代次數(shù)初值為0及精度大于0.00001N1=5;%設(shè)定除平衡節(jié)點(diǎn)外節(jié)點(diǎn)數(shù)之和while precision > 0.00001 E(1)=1.1;F(1)=0;%輸入平衡節(jié)點(diǎn)電壓實(shí)部,虛部e,f for m=2:N1+1 for n=1:N1+1 Pt(n)=(E(m)*(G(m,n)*E(n)-B(m,n)*F(n)+F(m)*(G(m,n)*F(n)+B(m,

15、n)*E(n); Qt(n)=(F(m)*(G(m,n)*E(n)-B(m,n)*F(n)-E(m)*(G(m,n)*F(n)+B(m,n)*E(n);%以設(shè)定值計(jì)算的注入功率 end dP(m)=P(m)-sum(Pt); dQ(m)=Q(m)-sum(Qt); %計(jì)算功率誤差endfor m=2:N1+1 for n=1:N1+1 Ai(n)=G(m,n)*E(n)-B(m,n)*F(n); Bi(n)=G(m,n)*F(n)+B(m,n)*E(n); end N(m,m)=sum(Ai)+G(m,m)*E(m)+B(m,m)*F(m);%計(jì)算dP/de H(m,m)=sum(Bi)-B(

16、m,m)*E(m)+G(m,m)*F(m);%計(jì)算dP/df L(m,m)= -sum(Bi)-B(m,m)*E(m)+G(m,m)*F(m);%計(jì)算dQ/de J(m,m)= sum(Ai)-G(m,m)*E(m)-B(m,m)*F(m);%計(jì)算dQ/df endfor m=2:N1+1 JJ(2*m-3,2*m-3)=N (m,m); JJ(2*m-3,2*m-2)=H(m,m); JJ(2*m-2,2*m-3)=L(m,m); JJ(2*m-2,2*m-2)=J(m,m); end %給雅可比矩陣的對角元素賦值for m=2:N1+1 for n=2:N1+1 if m=n else N

17、(m,n)=G(m,n)*E(m)+B(m,n)*F(m); %計(jì)算dP/de H(m,n)=-B(m,n)*E(m)+G(m,n)*F(m); %計(jì)算dP/df L(m,n)=-B(m,n)*E(m)+G(m,n)*F(m); %計(jì)算dQ/de J(m,n)=-B(m,n)*F(m)-G(m,n)*E(m);%計(jì)算dQ/df%*給雅可比矩陣的非對角元素賦值* JJ(2*m-3,2*n-3)=N(m,n); JJ(2*m-3,2*n-2)=H(m,n); JJ(2*m-2,2*n-3)=L(m,n); JJ(2*m-2,2*n-2)=J(m,n); end end endfor i=1:8JJ

18、(10,i)=0;%修正雅可比矩陣第十行1,到8元素為0endJJ(10,9)=2*E(6); %修正雅可比矩陣第十行9元素為2eJJ(10,10)=2*F(6);%修正雅可比矩陣第十行10元素為2f for m=1:N1 PQ(2*m-1)=dP(m+1); PQ(2*m)=dQ(m+1);end%寫PQ矩陣，奇數(shù)行為dP，偶數(shù)行為dQPQ(10)=1.12-(E(6)2+F(6)2);%修正PQ矩陣第十個(gè)元素為du2 dU=inv(JJ)*PQ'%求dU為JJ的逆矩陣乘以PQ' precision=max(abs(PQ); %取精度為PQ的絕對值的最大值 for n=1:N1 E(n+1)=E(n+1)+dU(2*n-1); F(n+1)