電力系統(tǒng)分析課程設(shè)計--電力系統(tǒng)潮流計算_第1頁
已閱讀1頁,還剩22頁未讀, 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

1、<p>  信息工程學院課程設(shè)計任務(wù)書</p><p>  年 月 日 </p><p><b>  摘 要</b></p><p>  潮流計算是電力網(wǎng)絡(luò)設(shè)計及運行中最基本的計算,對電力網(wǎng)絡(luò)的各種設(shè)計方案及各種運行方式進行潮流計算,可以得到各種電網(wǎng)各節(jié)點的電壓,并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中各元

2、件的電力損耗,進而求得電能損耗。本文以一個具體例子分析潮流計算的具體方法,并給出MATLAB仿真。</p><p>  關(guān)鍵詞: 電力系統(tǒng) 潮流計算 Mmatlab仿真</p><p><b>  目 錄</b></p><p><b>  1 需求分析6</b></p><p>  1

3、.1 所給電路圖及設(shè)計要求6</p><p>  1.2 潮流計算的意義6</p><p><b>  2 總體設(shè)計8</b></p><p>  2.1 等值電路圖8</p><p>  2.2 設(shè)計思路及方法8</p><p><b>  3 詳細設(shè)計9</b>

4、;</p><p><b>  3.1參數(shù)計算9</b></p><p>  3.2潮流計算12</p><p>  3.2 MATLAB仿真 14 </p><p><b>  4 總結(jié)23</

5、b></p><p><b>  參考文獻24</b></p><p><b>  1 需求分析</b></p><p><b>  1.1 二級標題</b></p><p>  1.2 潮流計算的意義</p><p>  潮流計算是研究電力系統(tǒng)

6、穩(wěn)態(tài)運行情況的一種基本電氣計算,常規(guī)潮流計算的任務(wù)是根據(jù)給定的運行條件和網(wǎng)路結(jié)構(gòu)確定整個系統(tǒng)的運行狀態(tài),如各母線上的電壓(幅值及相角)、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計算的結(jié)果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎(chǔ)。通過潮流計算可以判斷電網(wǎng)母線電壓、支路電流和功率是否越限,如果有越限,就應采取措施,調(diào)整運行方式。對于正在規(guī)劃的電力系統(tǒng),通過潮流計算,可以為選擇電網(wǎng)供電方案和電氣設(shè)備提供依據(jù)。潮流計算還可以為繼電保護和自動裝置定整計算、

7、電力系統(tǒng)故障計算和穩(wěn)定計算等提供原始數(shù)據(jù)。</p><p>  具體表現(xiàn)在以下方面:</p><p>  (1)在電網(wǎng)規(guī)劃階段,通過潮流計算,合理規(guī)劃電源容量及接入點,合理規(guī)劃網(wǎng)架,選擇無功補償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。</p><p>  (2)在編制年運行方式時,在預計負荷增長及新設(shè)備投運基礎(chǔ)上,選擇典型方式進行潮流計

8、算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對規(guī)劃、基建部門提出改進網(wǎng)架結(jié)構(gòu),加快基建進度的建議。</p><p>  (3)正常檢修及特殊運行方式下的潮流計算,用于日運行方式的編制,指導發(fā)電廠開機方式,有功、無功調(diào)整方案及負荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。</p><p>  (4)預想事故、設(shè)備退出運行對靜態(tài)安全的影響分析及作出預想的運行方式調(diào)整方案。&l

9、t;/p><p>  總結(jié)為在電力系統(tǒng)運行方式和規(guī)劃方案的研究中,都需要進行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟性。同時,為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài),也需要進行大量而快速的潮流計算。因此,潮流計算是電力系統(tǒng)中應用最廣泛、最基本和最重要的一種電氣運算。在系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式時,采用離線潮流計算;在電力系統(tǒng)運行狀態(tài)的實時監(jiān)控中,則采用在線潮流計算。 </p><p

10、><b>  2 總體設(shè)計</b></p><p><b>  2.1 等值電路圖</b></p><p><b>  設(shè)計思路及方法</b></p><p>  潮流計算(load flow calculation)根據(jù)電力系統(tǒng)接線方式、參數(shù)和運行條件計算電力系統(tǒng)穩(wěn)態(tài)運行狀態(tài)下的電氣量。通常給

11、定的運行條件有電源和負荷節(jié)點的功率、樞紐點電壓、平衡節(jié)點的電壓和相位角。待求的運行狀態(tài)量包括各節(jié)點電壓及其相位角和各支路(元件)通過的電流(功率)、網(wǎng)絡(luò)的功率損耗等。</p><p>  此題先按電網(wǎng)平均額定電壓VN計算功率損耗,再利用上面求得的功率分布,從電源點開始,順著功率的傳送方向,依次計算各段線路的電壓降落,求出各節(jié)點電壓。對于首端電壓和功率已知的線路,只需順著電路的方向算各個變壓器及線路的電壓降落和功率

12、損耗,以此求得各節(jié)點的電壓和功率。</p><p><b>  3 詳細設(shè)計</b></p><p><b>  3.1參數(shù)計算</b></p><p><b>  變壓器1:</b></p><p><b>  =</b></p><

13、p><b>  勵磁功率 :</b></p><p><b>  線路一</b></p><p><b>  設(shè),,</b></p><p><b>  產(chǎn)生的充電功率:</b></p><p><b>  同理可得</b>&l

14、t;/p><p>  變壓器2(2臺變壓器并聯(lián)運行):</p><p><b>  ;</b></p><p><b>  線路2:</b></p><p><b>  產(chǎn)生的充電功率:</b></p><p><b>  變壓器3:</b&

15、gt;</p><p><b>  線路3(雙回路):</b></p><p><b>  產(chǎn)生的充電功率:</b></p><p>  變壓器4(兩臺變壓器并聯(lián)):</p><p><b>  ;</b></p><p><b>  線路4:

16、</b></p><p><b>  產(chǎn)生的充電功率:</b></p><p>  變壓器5(兩臺三繞組變壓器并聯(lián)):</p><p><b>  線路5(雙回路):</b></p><p><b>  產(chǎn)生的充電功率:</b></p><p&g

17、t;  變壓器6(兩臺變壓器并聯(lián)):</p><p><b>  ;</b></p><p><b>  發(fā)電機A:</b></p><p><b>  發(fā)電機B:</b></p><p><b>  發(fā)電機C:</b></p><p&

18、gt;<b>  發(fā)電機D:</b></p><p><b>  發(fā)電機E:</b></p><p><b>  3.2潮流計算</b></p><p>  變壓器1的功率損耗:</p><p><b>  a點功率</b></p><

19、p><b>  線路1的功率損耗:</b></p><p>  變壓器T2的功率損耗:</p><p><b>  所以b點功率為</b></p><p><b>  c點功率:</b></p><p>  計算線路I的電壓分布:</p><p>

20、<b>  VA=115V;</b></p><p>  變壓器2中的電壓降落:</p><p>  線路I中的電壓降落:</p><p><b>  線路II中:</b></p><p>  線路III的功率損耗和電壓降落:</p><p><b>  線路IV的

21、損耗:</b></p><p><b>  點f的功率和電壓:</b></p><p>  變壓器6的功率和損耗:</p><p><b>  點g的功率及電壓:</b></p><p>  線路5 的功率損耗及電壓降落:</p><p><b>  點

22、h的功率:</b></p><p>  變壓器5-2的功率損耗和電壓降落:</p><p><b>  點m的功率和電壓:</b></p><p>  變壓器5-3的功率損耗和電壓降落:</p><p><b>  點n的功率和電壓:</b></p><p>  

23、3.3MATLAB仿真(由陳志彪完成)</p><p>  相關(guān)的原始數(shù)據(jù)輸入格式如下:</p><p>  1、B1是支路參數(shù)矩陣,第一列和第二列是節(jié)點編號。節(jié)點編號由小到大編寫。</p><p>  2、對于含有變壓器的支路,第一列為低壓側(cè)節(jié)點編號,第二列為高壓側(cè)節(jié)點編號,將變壓器的串聯(lián)阻抗置于低壓側(cè)處理,第三列為支路的串列阻抗參數(shù),第四列為支路的對地導納參數(shù),

24、第五烈為含變壓器支路的變壓器的變比,第六列為變壓器是否是否含有變壓器的參數(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>  4、X為節(jié)點號和對地參數(shù)矩陣。其中第

25、一列為節(jié)點編號,第二列為節(jié)點對地參數(shù)。</p><p><b>  程序如下:</b></p><p><b>  close all</b></p><p><b>  clear all</b></p><p>  n=input('請輸入節(jié)點數(shù):n=');

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

27、);</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><p&

28、gt;<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,q)=Y(p,q)

29、-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 %含有變壓器的支路&l

30、t;/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><

31、;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(2*n-2,1);

32、</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><p>  fo

33、r 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,1)=OrgS(2

34、*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))*(real(Y(i,j

35、))*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><

36、;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><b&g

37、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

38、><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<

39、/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=zeros(n-

40、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>  PVU(t,1)=

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

42、lt;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)=real(B2(i,2

43、))-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></

44、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><

45、;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><p>&

46、lt;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><b&g

47、t;  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><p

48、><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><p><

49、;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><p>

50、;<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-1)=-imag

51、(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)=-Jacbi(2*h

52、-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))+real(Y(i,j)

53、)*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(2*h-1,2*

54、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><p&

55、gt;<b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  k=0;</b></p>

56、<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~=isb</b&g

57、t;</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><p

58、>  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><p>

59、;  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,3));</

60、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><b>  k

61、=0;</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b&

62、gt;  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><p>&l

63、t;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><p&

64、gt;<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>  f

65、or 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>  end&

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

67、hile 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></p>&

68、lt;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(2*h-1,1)

69、=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))*(rea

70、l(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&g

71、t;<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>  for j

72、=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><

73、;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</b><

74、;/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>  h=0;&

75、lt;/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(2*h-1,1)

76、;</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><p>

77、<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)=real(B2(i

78、,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></p>

79、;<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</b>&

80、lt;/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>  end</b

81、></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>  h=0;

82、</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><b>

83、;  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)=-imag(Y(

84、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)=-Jacbi(2*h-1,

85、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))+real(Y(i,j))

86、*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(2*h-1,2*k

87、-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</b>&

88、lt;/p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  k=0;</b></p><p><b>  for i=1:

89、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><p>&l

90、t;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>  Jacbi(2*h-1,2

91、*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>  else<

92、;/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><p> 

93、 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>  end

94、</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>

95、  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><p

96、><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*j,1)+Deta

97、U(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><p>&l

98、t;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><b>

99、;  B2</b></p><p>  Times=Times+1; %迭代次數(shù)加1</p><p><b>  end</b></p><p><b>  Times</b></p><p><b>  4 總結(jié)</b></p><p

100、>  透過這次電力系統(tǒng)分析課程設(shè)計,對自己以往的學習有了很深的認識。電力系統(tǒng)分析專業(yè)課的系統(tǒng)學習已經(jīng)過去了三個多月,由于電力系統(tǒng)分析課程本身較難,在學習過程中沒有做到足夠的理解學習以及在以后的時間沒有做到充分的溫習,因此在這次設(shè)計作業(yè)的解題過程中感遇到了很大的阻力。</p><p>  首先拿到圖就感到無從下手,已經(jīng)超出我們平時作業(yè)的水平,也超出了我們的能力范圍,以至好長時間都沒有動筆,只是每次拿出來看幾眼

101、。慢慢看出了些眉頭后,又用過請教同學和老師把部分不懂的細節(jié)弄清了。完成這個艱難的第一步后,對整個線路圖有了宏觀的把握,對解題方法也有了大致思路。第二步就應該對電路圖做詳細分析,計算線路參數(shù)作出等值電路,這一步相對進行的比較順利。接下來就是理清思路認真做潮流計算,在做潮流計算的時候由于計算量比較大,線路參數(shù)比較多,所花的時間很多。完成整個潮流計算后,面臨這一個更艱巨的任務(wù),就是手寫稿轉(zhuǎn)化成電子稿。由于公式比較復雜,上下標多,對公式編輯及畫

102、圖軟件應用甚少,操作起來很費時??偟脕碚f,還是基本上完成了任務(wù),現(xiàn)將在此次設(shè)計的過程中遇到的問題及收獲總結(jié)如下:</p><p><b>  存在的問題:</b></p><p>  不能熟記電力系統(tǒng)課本中有關(guān)公式,應用時需要查看課本。</p><p>  對三繞組的變壓器應用不深。</p><p>  對潮流計算的幾種

103、模型及基本運算了解不透,不能熟練應用。</p><p><b>  計算能力不強。</b></p><p>  對公式編輯及畫圖軟件viso操作不熟。</p><p>  由于在計算過程中忽略了變壓器的勵磁損耗,而且只做了一次迭代,沒有保證精度,與MATLAB仿真存在很一定的差異,有的甚至出現(xiàn)錯誤的數(shù)據(jù)。</p><p>

104、;<b>  收獲:</b></p><p>  為了較好的完成此次作業(yè),對主要的相關(guān)章節(jié)又進行較全面的學習,尤其是加深了對有關(guān)公式的應用。</p><p>  提高了對有關(guān)軟件的操作能力。</p><p>  由于此模板類似于畢業(yè)論文的模板,讓我們提前熟悉了畢業(yè)論文的要求和格式。這次的課程設(shè)計為明年的畢業(yè)論文打下了很好的基礎(chǔ)。</p&g

105、t;<p>  同時,通過此次的課程設(shè)計,讓我認識到一個很嚴重的問題:我們的學習不扎實??纪曛蠡蚴墙?jīng)過一段時間我們就對前面學過的知識沒有映像了,這主要是因為沒有理解性學習。所以,在以后的學習中,一定要學的扎實,學的牢靠,并且能學以致用,活學活用。</p><p><b>  參考文獻</b></p><p>  1.何仰贊 溫増銀 《電力系統(tǒng)分析》第三

106、版</p><p>  3.王守相 劉玉田 電力系統(tǒng)潮流計算研究現(xiàn)狀--《山東電力技術(shù)》1996年05期 </p><p>  4.飛思科技產(chǎn)品研發(fā)中心.MATLAB基礎(chǔ)與提高.北京:電子工業(yè)出版社,2005.</p><p>  5.Stephen J. Chapman. MATLAB Programming for Engineers Second Edi

溫馨提示

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

評論

0/150

提交評論