電力系統短路計算課程設計 (2)_第1頁
已閱讀1頁,還剩22頁未讀, 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

1、<p><b>  課程設計(論文)</b></p><p>  題 目 名 稱 電力系統短路計算課程設計 </p><p>  課 程 名 稱 電力系統暫態(tài)分析 </p><p>  學 生 姓 名

2、 </p><p>  學 號 </p><p>  系 、專 業(yè) 07級電力系統運行及其控制 </p><p>  指 導 教 師 </p><p>

3、  2010年 7 月 2 日</p><p><b>  目 錄</b></p><p>  1手工計算...............................................1</p><p>  2 程序設計......................................7</p>

4、<p>  2.1 程序流程圖..........................................7</p><p>  2.2MATLAB程序設計....................... 8</p><p>  2.3 運行結果..........................................19 結束語............

5、................................. 20</p><p>  課程設計(論文)任務書</p><p>  指導教師(簽字): 學生(簽字):</p><p>  注:1.此表由指導教師填寫,經系、教研室審批,指導教師、學生簽字后生效;</p><p>  2.此表1式

6、3份,學生、指導教師、教研室各1份。</p><p><b>  2程序設計</b></p><p>  2.1 程序設計框圖</p><p>  利用MATLAB程序編程計算故障點的短路故障電流及故障點相電壓。兩相短路接地故障計算的原理框圖如圖2.1所示。</p><p>  圖2.1兩相短路接地故障計算的原理框圖&l

7、t;/p><p>  2.2MATLAB程序設計</p><p>  clear all;</p><p>  format long e</p><p>  Y1=[-j*23.619417 0 0 j*4.444444 j*9.523810 j*4.651163</p>

8、<p>  0 -j*35.238095 0 j*28.571429 0 0</p><p>  0 0 -j*27.809524 0 j*23.80524 0</p><p>  j*4.444444 j*28

9、.571429 0 -j*41.015873 0 j*8.000000</p><p>  j*9.523810 0 j*23.809524 0 -j*39.047619 j*5.714286</p><p>  j*4.651163 0 0

10、 j*8.000000 j*5.714286 -j*18.365449] %輸入正序負序導納矩陣</p><p>  Y0=[-j*31.3736 0 0 j*2.5 j*5.000000 j*2.564103 </p><p>  0 -j*33.8346

11、 0 j*28.5714 0 0</p><p>  0 0 -j*9.2083 0 j*5.2083 0</p><p>  j*2.5 j*28.5714 0 -j*35.0714 0

12、 j*4 </p><p>  j*5 0 j*5.2083 0 -j*5.2083 j2.857182 </p><p>  0 0 0 j*4 j*2.8571 -j*3.2784] %輸入零序導納矩陣、

13、</p><p>  disp('正序導納矩陣');</p><p><b>  Y1</b></p><p>  disp('負序導納矩陣');</p><p><b>  Y0</b></p><p>  N1=input('ple

14、ase input 正負序網絡節(jié)點數 :');</p><p>  N2=input('please input 零序網絡節(jié)點數 :');</p><p>  N3=input('please input 正負序短路編號(正負序導納矩陣中短路點對應的編號) :');</p><p>  N4=input('pleas

15、e input 零序短路編號(零序導納矩陣中短路點對應的編號) :');</p><p>  disp('fault表示輸入短路類型 1表示單相a短路 2表示兩相bc短路 3表示兩相bc短路接地')</p><p>  fault=input('please input fault的值 :');</p><p>  %求正

16、序和負序因子表</p><p>  YY1=zeros(N1,N1);</p><p>  for m=1:N1</p><p>  for n=m :N1 </p><p><b>  if n==m</b></p><p>  YY1(m,m)=Y1(m,m);</p><

17、p>  for k=1:m-1 </p><p>  YY1(m,m)=YY1(m,m)-YY1(k,m)*YY1(k,m)*YY1(k,k);</p><p><b>  end</b></p><p><b>  else</b></p><p>  YY1(m,n)=Y1(m,n);&l

18、t;/p><p>  for k=1:m-1</p><p>  YY1(m,n)=YY1(m,n)-YY1(k,m)*YY1(k,n)*YY1(k,k); </p><p><b>  end</b></p><p>  YY1(m,n)=YY1(m,n)/YY1(m,m);</p><p><

19、;b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p>  for m=1:N1</p><p>  YY1(m,m)=1/YY1(m,m);</p><p><b>  e

20、nd</b></p><p><b>  YY2=YY1;</b></p><p>  disp('正序因子表');</p><p><b>  YY1 </b></p><p>  disp('負序因子表');</p><p>&

21、lt;b>  YY2</b></p><p><b>  %求零序因子表</b></p><p>  YY0=zeros(N2,N2);</p><p>  for m=1:N2</p><p>  for n=m :N2 </p><p><b>  if n==m&

22、lt;/b></p><p>  YY0(m,m)=Y0(m,m);</p><p>  for k=1:m-1 </p><p>  YY0(m,m)=YY0(m,m)-YY0(k,m)*YY0(k,m)*YY0(k,k);</p><p><b>  end</b></p><p>&l

23、t;b>  else</b></p><p>  YY0(m,n)=Y0(m,n);</p><p>  for k=1:m-1</p><p>  YY0(m,n)=YY0(m,n)-YY0(k,m)*YY0(k,n)*YY0(k,k); </p><p><b>  end</b></p>

24、;<p>  YY0(m,n)=YY0(m,n)/YY0(m,m);</p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p>  for m=1:N2</p>

25、<p>  YY0(m,m)=1/YY0(m,m);</p><p><b>  end</b></p><p>  disp('零序因子表');</p><p><b>  YY0</b></p><p>  %正負序阻抗矩陣第N3列元素</p><

26、p>  f1=zeros(1,N1);</p><p>  h1=zeros(1,N1);</p><p>  z1=zeros(N1,N1);</p><p>  for k=1:N1</p><p><b>  if k==N3</b></p><p>  f1(1,k)=1;</

27、p><p>  h1(1,k)=f1(1,k)*YY1(k,k);</p><p><b>  end</b></p><p><b>  if k>N3</b></p><p>  for m=N3:k-1 %將k寫成N1了</p>

28、<p>  f1(1,k)=f1(1,k)-YY1(m,k)*f1(1,m);</p><p><b>  end</b></p><p>  h1(1,k)=f1(1,k)*YY1(k,k);</p><p><b>  end</b></p><p><b>  end<

29、;/b></p><p>  disp('f1的值');</p><p><b>  f1</b></p><p>  disp('h1的值');</p><p><b>  h1</b></p><p>  for m=N1:-1:1&

30、lt;/p><p>  z1(m,N3)=h1(1,m);</p><p>  for k=m+1:N1</p><p>  z1(m,N3)=z1(m,N3)-YY1(m,k)*z1(k,N3);</p><p><b>  end</b></p><p><b>  end</b&

31、gt;</p><p><b>  z2=z1;</b></p><p>  disp('正序阻抗矩陣第N3列元素');</p><p>  disp('z1=');</p><p>  disp(z1(:,N3));</p><p>  disp('負序阻

32、抗矩陣第N3列元素');</p><p>  disp('z2=');</p><p>  disp(z2(:,N3));</p><p>  %求零序阻抗矩陣第N4列元素</p><p>  f0=zeros(1,N2);</p><p>  h0=zeros(1,N2);</p>

33、<p>  z0=zeros(N1,N1);</p><p>  for k=1:N2</p><p>  if k==N4 % 這里的4表示零序導納矩陣中的短路節(jié)點號(實際上是6但是矩陣中為4)</p><p>  f0(1,k)=1;</p><p>  h0(1,k)=f0(1,k)*YY0(k,k);<

34、/p><p><b>  end</b></p><p><b>  if k>N4</b></p><p>  for m=N4:k-1</p><p>  f0(1,k)=f0(1,k)-YY0(m,k)*f0(1,m); %將k寫成N1了</p><p><b&

35、gt;  end</b></p><p>  h0(1,k)=f0(1,k)*YY0(k,k);</p><p><b>  end</b></p><p><b>  end</b></p><p>  disp('f0的值');</p><p>

36、;<b>  f0</b></p><p>  disp('h0的值');</p><p><b>  h0</b></p><p>  for m=N2:-1:1</p><p>  z0(m,N4)=h0(1,m);</p><p>  for k=m+1

37、:N2</p><p>  z0(m,N4)=z0(m,N4)-YY0(m,k)*z0(k,N4);</p><p><b>  end</b></p><p><b>  end</b></p><p>  disp('零序阻抗矩陣第N4列元素');</p><

38、p>  disp('z0=');</p><p>  disp(z0(:,N4));</p><p><b>  %求短路電流</b></p><p>  ZFF1=z1(N3,N3);</p><p><b>  ZFF2=ZFF1</b></p><p&

39、gt;  ZFF0=z0(N4,N4);</p><p><b>  IF1=0;</b></p><p><b>  IF2=0;</b></p><p><b>  IF0=0;</b></p><p>  if fault==1;</p><p>

40、  IF1=1/(ZFF0+ZFF1+ZFF2);</p><p><b>  IF2=IF1;</b></p><p><b>  IF0=IF1;</b></p><p><b>  else</b></p><p><b>  fault==2;</b&g

41、t;</p><p>  IF1=1/(ZFF1+ZFF2);</p><p><b>  IF2=-IF1;</b></p><p>  if fault==3</p><p>  IF1=1/(ZFF1+ZFF2*ZFF0/(ZFF0+ZFF2));</p><p>  IF2=-IF1*ZF

42、F0/(ZFF2+ZFF0);</p><p>  IF0=-IF1*ZFF2/(ZFF2+ZFF0);</p><p><b>  end</b></p><p><b>  end</b></p><p>  disp('ZFF1 ZFF2 ZFF0的值分別為:');</p

43、><p><b>  ZFF1</b></p><p><b>  ZFF2</b></p><p><b>  ZFF0</b></p><p>  disp('正序短路電流');</p><p>  IF1 %正序短路電流</

44、p><p>  disp('負序短路電流');</p><p>  IF2 %負序短路電流</p><p>  disp('零序短路電流');</p><p>  IF0 %零序短路電流</p><p><b>  %求節(jié)點電壓</b></p>

45、<p>  V1=zeros(N1,1);</p><p>  V2=zeros(N1,1);</p><p>  V0=zeros(N1,1);</p><p>  %由于z0矩陣中的編號和實際中的編號不一樣所以要進行轉換</p><p>  %用于求節(jié)點電壓的零序阻抗矩陣第N3列元素</p><p> 

46、 for k=1:N1</p><p>  V1(k,1)=1-IF1*z1(k,N3);</p><p>  V2(k,1)=-IF2*z2(k,N3);</p><p>  V0(k,1)=-IF0*z0(k,N4);</p><p><b>  end</b></p><p>  disp(

47、'節(jié)點正序電壓');</p><p><b>  V1</b></p><p>  disp('節(jié)點負序電壓');</p><p><b>  V2</b></p><p>  disp('節(jié)點零序電壓');</p><p>&

48、lt;b>  V0</b></p><p>  %計算短路點的線路電壓</p><p>  a=-1/2+sqrt(3)/2*j;</p><p>  if fault==1</p><p>  Vfb=j*[(a*a-a)*ZFF2+(a*a-1)*ZFF0]*IF1;</p><p>  Vfc=

49、j*[(a-a*a)*ZFF2+(a-1)*ZFF0]*IF1;</p><p><b>  Vfa=0;</b></p><p><b>  else</b></p><p>  if fault==2</p><p>  Vfa=2*j*ZFF2*IF1;</p><p>

50、;  Vfb=-1/2*Vfa;</p><p>  Vfc=-1/2*Vfa;</p><p><b>  else</b></p><p>  if fault==3</p><p>  Vfa=j*3*ZFF2*ZFF0/(ZFF2+ZFF0)*IF1;</p><p><b> 

51、 Vfb=0;</b></p><p><b>  Vfc=0;</b></p><p><b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p>

52、;  disp('短路點的線路電壓');</p><p><b>  Vfa</b></p><p><b>  Vfb</b></p><p><b>  Vfc</b></p><p><b>  Y1 =</b></p>

53、<p>  Columns 1 through 3</p><p>  0 -2.361941700000000e+001i 0 0 </p><p>  0

54、 0 -3.523809500000000e+001i 0 </p><p>  0 0

55、 0 -2.780952400000000e+001i</p><p>  0 +4.444444000000000e+000i 0 +2.857142900000000e+001i 0 </p><p

56、>  0 +9.523809999999999e+000i 0 0 +2.380952400000000e+001i</p><p>  0 +4.651163000000000e+000i 0

57、 0 </p><p>  Columns 4 through 6</p><p>  0 +4.444444000000000e+000i 0 +9.523809999999999e+000i

58、 0 +4.651163000000000e+000i</p><p>  0 +2.857142900000000e+001i 0 0 </p><p>  

59、0 0 +2.380524000000000e+001i 0 </p><p>  0 -4.101587300000000e+001i 0

60、 0 +8.000000000000000e+000i</p><p>  0 0 -3.904761900000000e+001i 0 +5.714286000000000e+000i<

61、/p><p>  0 +8.000000000000000e+000i 0 +5.714286000000000e+000i 0 -1.836544900000000e+001i</p><p><b>  Y0 =</b></p><p>  Co

62、lumns 1 through 3</p><p>  0 -3.137360000000000e+001i 0 0 </p><p>  0

63、 0 -3.383460000000000e+001i 0 </p><p>  0 0

64、 0 -9.208300000000000e+000i</p><p>  0 +2.500000000000000e+000i 0 +2.857140000000000e+001i 0 </p><p>  0 +5.000

65、000000000000e+000i 0 0 +5.208300000000000e+000i</p><p>  0 0

66、 0 </p><p>  Columns 4 through 6</p><p>  0 +2.500000000000000e+000i 0 +5.000000000000000e+000i

67、 0 +2.564103000000000e+000i</p><p>  0 +2.857140000000000e+001i 0 0 </p><p>  0

68、 0 +5.208300000000000e+000i 0 </p><p>  0 -3.507140000000000e+001i 0

69、 0 +4.000000000000000e+000i</p><p>  0 0 -5.208300000000000e+000i 0 +2.857100000000000e+000i</p><p&

70、gt;  0 +4.000000000000000e+000i 0 +2.857100000000000e+000i 0 -3.278400000000000e+000i</p><p><b>  正序導納矩陣</b></p><p><b>  Y1 =&

71、lt;/b></p><p>  Columns 1 through 3</p><p>  0 -2.361941700000000e+001i 0 0 </p><

72、p>  0 0 -3.523809500000000e+001i 0 </p><p>  0 0

73、 0 -2.780952400000000e+001i</p><p>  0 +4.444444000000000e+000i 0 +2.857142900000000e+001i 0

74、 </p><p>  0 +9.523809999999999e+000i 0 0 +2.380952400000000e+001i</p><p>  0 +4.651163000000000e+000i

75、 0 0 </p><p>  Columns 4 through 6</p><p>  0 +4.444444000000000e+000i 0 +9.523809999

76、999999e+000i 0 +4.651163000000000e+000i</p><p>  0 +2.857142900000000e+001i 0 0 <

77、;/p><p>  0 0 +2.380524000000000e+001i 0 </p><p>  0 -4.101587300000000e+001i

78、 0 0 +8.000000000000000e+000i</p><p>  0 0 -3.904761900000000e+001i 0 +5.71428

79、6000000000e+000i</p><p>  0 +8.000000000000000e+000i 0 +5.714286000000000e+000i 0 -1.836544900000000e+001i</p><p><b>  負序導納矩陣</b><

80、;/p><p><b>  Y0 =</b></p><p>  Columns 1 through 3</p><p>  0 -3.137360000000000e+001i 0 0

81、 </p><p>  0 0 -3.383460000000000e+001i 0 </p><p>  0

82、 0 0 -9.208300000000000e+000i</p><p>  0 +2.500000000000000e+000i 0 +2.857140000000000e+001i

83、 0 </p><p>  0 +5.000000000000000e+000i 0 0 +5.208300000000000e+000i</p><p>  0

84、 0 0 </p><p>  Columns 4 through 6</p><p>  0 +2.500000000000000e+000i

85、 0 +5.000000000000000e+000i 0 +2.564103000000000e+000i</p><p>  0 +2.857140000000000e+001i 0

86、 0 </p><p>  0 0 +5.208300000000000e+000i 0 </p><p>  0 -3.50714000

87、0000000e+001i 0 0 +4.000000000000000e+000i</p><p>  0 0 -5.208300000000000e+000i

88、 0 +2.857100000000000e+000i</p><p>  0 +4.000000000000000e+000i 0 +2.857100000000000e+000i 0 -3.278400000000000e+000i</p><p

89、>  please input 網絡節(jié)點數 :6</p><p>  please input 零序網絡節(jié)點數 :6</p><p>  please input 正負序短路編號(正負序導納矩陣中短路點對應的編號) :6</p><p>  please input 零序短路編號(零序導納矩陣中短路點對應的編號) :6</p><p&g

90、t;  fault表示輸入短路類型 1表示單相a短路 2表示兩相bc短路 3表示兩相bc短路接地</p><p>  please input fault的值 :3</p><p><b>  正序因子表</b></p><p><b>  YY1 =</b></p><p>  Columns

91、1 through 3</p><p>  0 +4.233804754791366e-002i 0 0 </p><p>  0

92、 0 +2.837837857012418e-002i 0 </p><p>  0 0

93、 0 +3.595890386329518e-002i</p><p>  0 0 0 </p><p>  0

94、 0 0 </p><p>  0 0

95、 0 </p><p>  Columns 4 through 6</p><p>  -1.881690813960396e-001 -4.032195206172955e-001 -1

96、.969211602470967e-001 </p><p>  -8.108108284514245e-001 0 0

97、 </p><p>  0 -8.560103366026690e-001 0 </p><p>  0 +5.877670844705451e-002i -1.0533295037

98、10766e-001 -5.216553407562652e-001 </p><p>  0 0 +6.830071284011730e-002i -5.822347028529298e-001

99、 </p><p>  0 0 0 +1.272842145325835e-001i</p><p><b>  負序因子表</b>

100、</p><p><b>  YY2 =</b></p><p>  Columns 1 through 3</p><p>  0 +4.233804754791366e-002i 0 0

101、 </p><p>  0 0 +2.837837857012418e-002i 0 </p><p>  0

102、 0 0 +3.595890386329518e-002i</p><p>  0 0

103、 0 </p><p>  0 0 0 </p><p>  0

104、 0 0 </p><p>  Columns 4 through 6</p><p>  -1.881690813960396e-001

105、 -4.032195206172955e-001 -1.969211602470967e-001 </p><p>  -8.108108284514245e-001

106、 0 0 </p><p>  0 -8.560103366026690e-001 0

107、 </p><p>  0 +5.877670844705451e-002i -1.053329503710766e-001 -5.216553407562652e-001 </p><p>  0

108、 0 +6.830071284011730e-002i -5.822347028529298e-001 </p><p>  0 0

109、 0 +1.272842145325835e-001i</p><p><b>  零序因子表</b></p><p><b>  YY0 =</b></p><p>  Columns 1 through 3</p><p>  0 +3.187393222327052e-002i

110、 0 0 </p><p>  0 0 +2.955554373333806e-002i

111、 0 </p><p>  0 0 0 +1.085976781816405e-001i</p><p>  0

112、 0 0 </p><p>  0 0

113、 0 </p><p>  0 0 0 </p><

114、;p>  Columns 4 through 6</p><p>  -7.968483055817630e-002 -1.593696611163526e-001 -8.172804523548460e-002 </p><p&

115、gt;  -8.444432622226951e-001 0 0 </p><p>  0 -5.6560928

116、72734382e-001 0 </p><p>  0 +9.306427593808672e-002i -3.707905529572902e-002 -3.912720071365359e-001

117、 </p><p>  0 0 +6.892674464886308e-001i -2.358419857589230e+000 </p><p>  0

118、 0 0 -1.504701183754390e-001i</p><p>  f1的值f1 = 0 0 0 0 0 1</p><p>  h1的值h1 = Columns 1 t

溫馨提示

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

評論

0/150

提交評論