用matlab電力系統(tǒng)潮流計算(課程設計)_第1頁
已閱讀1頁,還剩21頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領

文檔簡介

1、<p><b>  信息工程學院</b></p><p><b>  課程設計報告書</b></p><p>  題目: 電力系統(tǒng)潮流計算 </p><p>  專 業(yè): 電氣工程及其自動化 </p><p>  班 級:

2、 </p><p>  學 號: </p><p>  學生姓名: </p><p>  指導教師: </p><p>  2011年 6 月 10 日</p><p>  201

3、1年6月12日 </p><p><b>  摘 要</b></p><p>  電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計算和靜態(tài)安全分析。本文主要運用的事潮流計算,潮流計算是電力網(wǎng)絡設計與運行中最基本的運算,對電力網(wǎng)絡的各種設計方案及各種運行方式進行潮流計算,可以得到各種電網(wǎng)各節(jié)點的電壓,并求得網(wǎng)絡的潮流及網(wǎng)絡中的各元件的電力損耗,進而求得

4、電能損耗。本位就是運用潮流計算具體分析,并有MATLAB仿真。</p><p>  關鍵詞: 電力系統(tǒng) 潮流計算 MATLAB仿真</p><p><b>  Abstract</b></p><p>  Electric power system steady flow calculation and analysis of the

5、static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and t

6、he operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The</p

7、><p>  Key words: Power system; Flow calculation; MATLAB simulation</p><p><b>  目 錄</b></p><p>  1 任務提出與方案論證2</p><p><b>  2 總體設計3</b></p>

8、<p>  2.1潮流計算等值電路3</p><p>  2.2建立電力系統(tǒng)模型3</p><p>  2.3模型的調(diào)試與運行3</p><p><b>  3 詳細設計4</b></p><p>  3.1 計算前提4</p><p><b>  3.2手工計算

9、7</b></p><p>  4設計圖及源程序11</p><p>  4.1MATLAB仿真11</p><p>  4.2潮流計算源程序11</p><p><b>  5 總結19</b></p><p><b>  參考文獻20</b><

10、/p><p>  1 任務提出與方案論證</p><p>  潮流計算是在給定電力系統(tǒng)網(wǎng)絡結構、參數(shù)和決定系統(tǒng)運行狀態(tài)的邊界條件的情況下確定系統(tǒng)穩(wěn)態(tài)運行狀態(tài)的一種基本方法,是電力系統(tǒng)規(guī)劃和運營中不可缺少的一個重要組成部分??梢哉f,它是電力系統(tǒng)分析中最基本、 最重要的計算,是系統(tǒng)安全、 經(jīng)濟分析和實時控制與調(diào)度的基礎。常規(guī)潮流計算的任務是根據(jù)給定的運行條件和網(wǎng)路結構確定整個系統(tǒng)的運行狀態(tài),如各母

11、線上的電壓(幅值及相角)、網(wǎng)絡中的功率分布以及功率損耗等。潮流計算的結果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎。在電力系統(tǒng)運行方式和規(guī)劃方案的研究中,都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟性。同時,為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài),也需要進行大量而快速的潮流計算。因此,潮流計算是電力系統(tǒng)中應用最廣泛、最基本和最重要的一種電氣運算。在系統(tǒng)規(guī)劃設計和安排系統(tǒng)的運行方式時,采用離線潮流計算;在電力系統(tǒng)運行狀態(tài)的實時監(jiān)控

12、中,則采用在線潮流計算。是電力系統(tǒng)研究人員長期研究的一個課題。它既是對電力系統(tǒng)規(guī)劃設計和運行方式的合理性、 可靠性及經(jīng)濟性進行定量分析的依據(jù) , 又是電力系統(tǒng)靜態(tài)和暫態(tài)穩(wěn)定計算的基礎。</p><p>  潮流計算經(jīng)歷了一個由手工到應用數(shù)字電子計算機的發(fā)展過程,現(xiàn)在的潮流算法都以計算機的應用為前提用計算機進行潮流計算主要步驟在于編制計算機程序,這是一項非常復雜的工作。對系統(tǒng)進行潮流分析,本文利用 MATLAB中的

13、SimpowerSystems工具箱設計電力系統(tǒng),在simulink 環(huán)境下,不僅可以仿真系統(tǒng)的動態(tài)過程,還可以對系統(tǒng)進行穩(wěn)態(tài)潮流分析。</p><p><b>  2 總體設計</b></p><p>  SimpowerSystems使用Simulink環(huán)境,可以將該系統(tǒng)中的發(fā)電機、變壓器,線路等模型聯(lián)結起來,形成電力系統(tǒng)仿真模擬圖。在加人測量模塊,并對各元件的參

14、數(shù)進行設置后,用measurement和sink中的儀器可以觀察各元件的電壓、電流、功率的大小。</p><p>  2.1潮流計算等值電路</p><p>  2.2建立電力系統(tǒng)模型</p><p>  在Simulink中按照電力系統(tǒng)原型選擇元件進行建模。所建立的模型和建立的方法在詳細設計中詳述。</p><p>  在電力系統(tǒng)模型的建立

15、工程中主要涉及到的是:元器件的選擇及其參數(shù)的設置;發(fā)電機選型;變壓器選擇;線路的選擇;負荷模型的選擇;母線選擇。</p><p>  2.3模型的調(diào)試與運行</p><p>  建立系統(tǒng)模型,并設置好參數(shù)以后,就可以在Simulink環(huán)境下進行仿真運行。運行的具體結果和分析也在詳細設計中詳述。</p><p><b>  3 詳細設計</b>&

16、lt;/p><p><b>  3.1 計算前提</b></p><p>  首先是發(fā)電機的參數(shù)計算,先對5個發(fā)電廠簡化為5臺發(fā)電機來計算。</p><p><b>  發(fā)電機G1:</b></p><p><b>  發(fā)電機G2:</b></p><p>

17、<b>  發(fā)電機G3:</b></p><p><b>  發(fā)電機G4:</b></p><p><b>  發(fā)電機G5:</b></p><p>  其次是變電站的參數(shù)計算,我們還是對7個變電站簡化為7臺變壓器來計算。</p><p><b>  變壓器T1:&l

18、t;/b></p><p>  變壓器T2:(雙并聯(lián))</p><p>  變壓器T3:(四并聯(lián))</p><p>  變壓器T4:(雙并聯(lián))</p><p><b>  變壓器T5:</b></p><p>  變壓器T6:(兩個三繞組變壓器并聯(lián))</p><p>

19、  變壓器T7:(雙并聯(lián))</p><p>  再次是傳輸線參數(shù)計算,5條傳輸線的具體計算如下。</p><p><b>  根據(jù)教材查得</b></p><p><b>  線路L1:</b></p><p><b>  線路L2:</b></p><

20、p>  線路L3:(雙回路)</p><p><b>  線路L4:</b></p><p>  線路L5:(雙回路)</p><p><b>  3.2手工計算</b></p><p><b>  FLR1:</b></p><p><b&

21、gt;  FLR2:</b></p><p><b>  FLR3:</b></p><p><b>  FLR4:</b></p><p><b>  FLR5:</b></p><p>  計算每一個FLR的功率分布和電壓分布計算如下:</p>&

22、lt;p><b>  FLR1:</b></p><p><b>  FLR2:</b></p><p><b>  功率分布:</b></p><p><b>  電壓分布:</b></p><p><b>  FLR3:</b&g

23、t;</p><p><b>  功率分布:</b></p><p><b>  電壓分布:</b></p><p><b>  FLR4:</b></p><p><b>  功率分布:</b></p><p><b>

24、  電壓分布:</b></p><p><b>  FLR5:</b></p><p>  這里我們先將f點和發(fā)電機G5當做電源,經(jīng)過和構成兩端供電網(wǎng)絡以g點作為運算負荷進行計算。</p><p><b>  電壓分布:</b></p><p><b>  4設計圖及源程序&l

25、t;/b></p><p>  4.1MATLAB仿真</p><p>  相關的原始數(shù)據(jù)輸入格式如下:</p><p>  1、B1是支路參數(shù)矩陣,第一列和第二列是節(jié)點編號。節(jié)點編號由小到大編寫。</p><p>  2、對于含有變壓器的支路,第一列為低壓側節(jié)點編號,第二列為高壓側節(jié)點編號,將變壓器的串聯(lián)阻抗置于低壓側處理,第三列為支

26、路的串列阻抗參數(shù),第四列為支路的對地導納參數(shù),第五烈為含變壓器支路的變壓器的變比,第六列為變壓器是否是否含有變壓器的參數(shù),其中“1”為含有變壓器,“0”為不含有變壓器。</p><p>  3、B2為節(jié)點參數(shù)矩陣,其中第一列為節(jié)點注入發(fā)電功率參數(shù);第二列為節(jié)點負荷功率參數(shù);第三列為節(jié)點電壓參數(shù);第六列為節(jié)點類型參數(shù),其中“1”為平衡節(jié)點,“2”為PQ節(jié)點,“3”為PV節(jié)點參數(shù)。</p><p&

27、gt;  4、X為節(jié)點號和對地參數(shù)矩陣。其中第一列為節(jié)點編號,第二列為節(jié)點對地參數(shù)。</p><p>  4.2潮流計算源程序</p><p><b>  close all</b></p><p><b>  clear all</b></p><p>  n=input('請輸入節(jié)點數(shù):

28、n=');</p><p>  n1=input('請輸入支路數(shù):n1=');</p><p>  isb=input('請輸入平衡節(jié)點號:isb=');</p><p>  pr=input('請輸入誤差精度:pr=');</p><p>  B1=input('請輸入支路參數(shù)

29、:B1=');</p><p>  B2=input('請輸入節(jié)點參數(shù):B2=');</p><p>  X=input('節(jié)點號和對地參數(shù):X=');</p><p>  Y=zeros(n);</p><p>  Times=1; %置迭代次數(shù)為初始值</p>

30、;<p><b>  %創(chuàng)建節(jié)點導納矩陣</b></p><p>  for i=1:n1</p><p>  if B1(i,6)==0 %不含變壓器的支路</p><p>  p=B1(i,1);</p><p>  q=B1(i,2);</p><p>  Y(p,

31、q)=Y(p,q)-1/B1(i,3);</p><p>  Y(q,p)=Y(p,q);</p><p>  Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);</p><p>  Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);</p><p>  else %含

32、有變壓器的支路</p><p>  p=B1(i,1);</p><p>  q=B1(i,2);</p><p>  Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));</p><p>  Y(q,p)=Y(p,q);</p><p>  Y(p,p)=Y(p,p)+1/B1(i,3);</p

33、><p>  Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));</p><p><b>  end</b></p><p><b>  end</b></p><p><b>  Y</b></p><p>  OrgS=zeros(

34、2*n-2,1);</p><p>  DetaS=zeros(2*n-2,1); %將OrgS、DetaS初始化</p><p>  %創(chuàng)建OrgS,用于存儲初始功率參數(shù)</p><p><b>  h=0;</b></p><p><b>  j=0;</b></p><

35、p>  for i=1:n %對PQ節(jié)點的處理</p><p>  if i~=isb&B2(i,6)==2 </p><p><b>  h=h+1;</b></p><p><b>  for j=1:n</b></p><p>  OrgS(2*h-1,

36、1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p>  OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(r

37、eal(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>  end</b></p><p><b>  end</b></p

38、><p><b>  end</b></p><p>  for i=1:n %對PV節(jié)點的處理,注意這時不可再將h初始化為0</p><p>  if i~=isb&B2(i,6)==3</p><p><b>  h=h+1;</b></p><p&g

39、t;<b>  for j=1:n</b></p><p>  OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)

40、));</p><p>  OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>

41、  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  OrgS</b></p><p>  %創(chuàng)建PVU 用于存儲PV節(jié)點的初始電壓</p><p>  PVU

42、=zeros(n-h-1,1);</p><p><b>  t=0;</b></p><p><b>  for i=1:n</b></p><p>  if B2(i,6)==3</p><p><b>  t=t+1;</b></p><p>  

43、PVU(t,1)=B2(i,3);</p><p><b>  end</b></p><p><b>  end</b></p><p><b>  PVU</b></p><p>  %創(chuàng)建DetaS,用于存儲有功功率、無功功率和電壓幅值的不平衡量</p>&

44、lt;p><b>  h=0;</b></p><p>  for i=1:n %對PQ節(jié)點的處理</p><p>  if i~=isb&B2(i,6)==2</p><p><b>  h=h+1;</b></p><p>  DetaS(2*h-1,1)=re

45、al(B2(i,2))-OrgS(2*h-1,1);</p><p>  DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);</p><p><b>  end</b></p><p><b>  end</b></p><p><b>  t=0;</b

46、></p><p>  for i=1:n %對PV節(jié)點的處理,注意這時不可再將h初始化為0</p><p>  if i~=isb&B2(i,6)==3</p><p><b>  h=h+1;</b></p><p><b>  t=t+1;</b></p

47、><p>  DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);</p><p>  DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;</p><p><b>  end</b></p>&

48、lt;p><b>  end</b></p><p><b>  DetaS</b></p><p>  %創(chuàng)建I,用于存儲節(jié)點電流參數(shù)</p><p>  i=zeros(n-1,1);</p><p><b>  h=0;</b></p><p&g

49、t;<b>  for i=1:n</b></p><p><b>  if i~=isb</b></p><p><b>  h=h+1;</b></p><p>  I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));</p&g

50、t;<p><b>  end</b></p><p><b>  end</b></p><p><b>  I</b></p><p>  %創(chuàng)建Jacbi(雅可比矩陣)</p><p>  Jacbi=zeros(2*n-2);</p><

51、;p><b>  h=0;</b></p><p><b>  k=0;</b></p><p>  for i=1:n %對PQ節(jié)點的處理</p><p>  if B2(i,6)==2</p><p><b>  h=h+1;</b></p>

52、<p><b>  for j=1:n</b></p><p><b>  if j~=isb</b></p><p><b>  k=k+1;</b></p><p>  if i==j %對角元素的處理</p><p>  Jacbi(2*h-1,2*k

53、-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>  Jacbi(2*h,2*k-1)=-

54、Jacbi(2*h-1,2*k)+2*real(I(h,1));</p><p>  Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));</p><p>  else %非對角元素的處理</p><p>  Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+re

55、al(Y(i,j))*imag(B2(i,3));</p><p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>  Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);</p><p>  Jacbi(2*h,2*k)=Jacbi

56、(2*h-1,2*k-1);</p><p><b>  end</b></p><p>  if k==(n-1) %將用于內(nèi)循環(huán)的指針置于初始值,以確保雅可比矩陣換行</p><p><b>  k=0;</b></p><p><b>  end</b></p>

57、;<p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  k=0;</b><

58、;/p><p>  for i=1:n %對PV節(jié)點的處理</p><p>  if B2(i,6)==3</p><p><b>  h=h+1;</b></p><p><b>  for j=1:n</b></p><p><b>  if j~=is

59、b</b></p><p><b>  k=k+1;</b></p><p>  if i==j %對角元素的處理</p><p>  Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p&g

60、t;<p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>  Jacbi(2*h,2*k-1)=2*imag(B2(i,3));</p><p>  Jacbi(2*h,2*k)=2*real(B2(i,3));</p>

61、<p>  else %非對角元素的處理</p><p>  Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,

62、3));</p><p>  Jacbi(2*h,2*k-1)=0;</p><p>  Jacbi(2*h,2*k)=0;</p><p><b>  end</b></p><p>  if k==(n-1) %將用于內(nèi)循環(huán)的指針置于初始值,以確保雅可比矩陣換行</p><p><

63、;b>  k=0;</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p&

64、gt;<b>  end</b></p><p><b>  Jacbi</b></p><p>  %求解修正方程,獲取節(jié)點電壓的不平衡量</p><p>  DetaU=zeros(2*n-2,1);</p><p>  DetaU=inv(Jacbi)*DetaS;</p>&l

65、t;p><b>  DetaU</b></p><p><b>  %修正節(jié)點電壓</b></p><p><b>  j=0;</b></p><p>  for i=1:n %對PQ節(jié)點處理</p><p>  if B2(i,6)==2</p>

66、;<p><b>  j=j+1;</b></p><p>  B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>  end</b></p><p><b>  end</b></p><

67、;p>  for i=1:n %對PV節(jié)點的處理</p><p>  if B2(i,6)==3</p><p><b>  j=j+1;</b></p><p>  B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b&

68、gt;  end</b></p><p><b>  end</b></p><p><b>  B2</b></p><p>  %開始循環(huán)**********************************************************************</p><

69、;p>  while abs(max(DetaU))>pr</p><p>  OrgS=zeros(2*n-2,1); %!!!初始功率參數(shù)在迭代過程中是不累加的,所以在這里必須將其初始化為零矩陣</p><p><b>  h=0;</b></p><p><b>  j=0;</b><

70、/p><p><b>  for i=1:n</b></p><p>  if i~=isb&B2(i,6)==2</p><p><b>  h=h+1;</b></p><p><b>  for j=1:n</b></p><p>  OrgS

71、(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p>  OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i

72、,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>  end</b></p><p><b>  end</b>

73、;</p><p><b>  end</b></p><p><b>  for i=1:n</b></p><p>  if i~=isb&B2(i,6)==3</p><p><b>  h=h+1;</b></p><p><b&g

74、t;  for j=1:n</b></p><p>  OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p

75、><p>  OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>  end<

76、/b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  OrgS</b></p><p><b>  %創(chuàng)建DetaS</b></p><p><b&g

77、t;  h=0;</b></p><p><b>  for i=1:n</b></p><p>  if i~=isb&B2(i,6)==2</p><p><b>  h=h+1;</b></p><p>  DetaS(2*h-1,1)=real(B2(i,2))-OrgS

78、(2*h-1,1);</p><p>  DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);</p><p><b>  end</b></p><p><b>  end</b></p><p><b>  t=0;</b></p>

79、<p><b>  for i=1:n</b></p><p>  if i~=isb&B2(i,6)==3</p><p><b>  h=h+1;</b></p><p><b>  t=t+1;</b></p><p>  DetaS(2*h-1,1)=

80、real(B2(i,2))-OrgS(2*h-1,1);</p><p>  DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;</p><p><b>  end</b></p><p><b>  end</b>

81、</p><p><b>  DetaS</b></p><p><b>  %創(chuàng)建I</b></p><p>  i=zeros(n-1,1);</p><p><b>  h=0;</b></p><p><b>  for i=1:n&l

82、t;/b></p><p><b>  if i~=isb</b></p><p><b>  h=h+1;</b></p><p>  I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));</p><p><b>  

83、end</b></p><p><b>  end</b></p><p><b>  I</b></p><p><b>  %創(chuàng)建Jacbi</b></p><p>  Jacbi=zeros(2*n-2);</p><p><b&

84、gt;  h=0;</b></p><p><b>  k=0;</b></p><p><b>  for i=1:n</b></p><p>  if B2(i,6)==2</p><p><b>  h=h+1;</b></p><p>

85、;<b>  for j=1:n</b></p><p><b>  if j~=isb</b></p><p><b>  k=k+1;</b></p><p><b>  if i==j</b></p><p>  Jacbi(2*h-1,2*k-1)

86、=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>  Jacbi(2*h,2*k-1)=-Jac

87、bi(2*h-1,2*k)+2*real(I(h,1));</p><p>  Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));</p><p><b>  else</b></p><p>  Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+rea

88、l(Y(i,j))*imag(B2(i,3));</p><p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>  Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);</p><p>  Jacbi(2*h,2*k)=Jacbi(

89、2*h-1,2*k-1);</p><p><b>  end</b></p><p>  if k==(n-1)</p><p><b>  k=0;</b></p><p><b>  end</b></p><p><b>  end&l

90、t;/b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  k=0;</b></p><p><b> 

91、 for i=1:n</b></p><p>  if B2(i,6)==3</p><p><b>  h=h+1;</b></p><p><b>  for j=1:n</b></p><p><b>  if j~=isb</b></p>&l

92、t;p><b>  k=k+1;</b></p><p><b>  if i==j</b></p><p>  Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>  Jacb

93、i(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>  Jacbi(2*h,2*k-1)=2*imag(B2(i,3));</p><p>  Jacbi(2*h,2*k)=2*real(B2(i,3));</p><p><b>

94、  else</b></p><p>  Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>  Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p>&

95、lt;p>  Jacbi(2*h,2*k-1)=0;</p><p>  Jacbi(2*h,2*k)=0;</p><p><b>  end</b></p><p>  if k==(n-1)</p><p><b>  k=0;</b></p><p><b

96、>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p>

97、<b>  Jacbi</b></p><p>  DetaU=zeros(2*n-2,1);</p><p>  DetaU=inv(Jacbi)*DetaS;</p><p><b>  DetaU</b></p><p><b>  %修正節(jié)點電壓</b></p&g

98、t;<p><b>  j=0;</b></p><p><b>  for i=1:n</b></p><p>  if B2(i,6)==2</p><p><b>  j=j+1;</b></p><p>  B2(i,3)=B2(i,3)+DetaU(2*

99、j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>  end</b></p><p><b>  end</b></p><p><b>  for i=1:n</b></p><p>  if B2(i,6)==3</p>&l

100、t;p><b>  j=j+1;</b></p><p>  B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>  end</b></p><p><b>  end</b></p><p>

101、;<b>  B2</b></p><p>  Times=Times+1; %迭代次數(shù)加1</p><p><b>  end</b></p><p><b>  Times</b></p><p><b>  5 總結</b></p&g

102、t;<p>  通過本次課程設計讓我有復習了一次潮流計算的相關知識,跟家清晰了什么事潮流計算以及潮流計算的在電力系統(tǒng)的重要性。電力系統(tǒng)的穩(wěn)定運行狀況即是正常運行狀況,是指電力系統(tǒng)在穩(wěn)定運行條件下電壓、功率的分布,也稱為潮流分布。電力系統(tǒng)分析的潮流計算是電力系統(tǒng)分析的一個重要的部分。通過對電力系統(tǒng)潮流分布的分析和計算,可進一步對系統(tǒng)運行的安全性,經(jīng)濟性進行分析、評估,提出改進措施。同時潮流分布也是電力系統(tǒng)規(guī)劃設計的一項基礎工

103、作。</p><p>  整個計算過程的模型建立并不是十分復雜,但計算過程十分繁瑣、計算量相當?shù)拇螅矣捎谥?jié)太多很容易算錯。不過在計算潮流計算的過程中卻對以往學過的電力系統(tǒng)分析的相關知識進行了一次較為深入的復習。而且整個計算對計算量的要求很大,鍛煉了我們的計算能力。而且對細節(jié)的把握也得到了鍛煉,做題的精細程度得到了提高。</p><p><b>  參考文獻</b>

104、;</p><p>  [1]何仰贊, 溫增銀《電力系統(tǒng)分析》(第三版)[M]. 華中科技大學,2002</p><p>  [2]http://baike.baidu.com/view/627420.</p><p>  [3]王守相,劉玉田 電力系統(tǒng)潮流計算研究現(xiàn)狀--《山東電力技術》1996年05期</p><p>  [4]何仰贊,溫

105、增銀. 電力系統(tǒng)分析(上冊)第三版[M]. 湖北:華中科技大學出版社,2002</p><p>  [5] 劉同娟. MATLAB在電路分析中的應用.電氣電子教學學報.2002</p><p>  [6] 西安交通大學等.電力系統(tǒng)計算[M].北京:水利電力出版社,1993.12</p><p>  [7] 李光琦.電力系統(tǒng)暫態(tài)分析[M].北京: 水利電力出版社,20

106、02.5</p><p>  [8]何仰贊,溫增銀. 電力系統(tǒng)分析(下冊)第三版[M]. 湖北:華中科技大學出版社,2002</p><p>  [9]韋化,李濱,杭乃善,等.大規(guī)模水一火電力系統(tǒng)最優(yōu)潮流的現(xiàn)代內(nèi)點算法實現(xiàn)[J].中國電機工程學報,2003.23(6):13一l8.</p><p>  [10]Chen Luo—nan,Suzuki Hideki,K

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論