電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)_第1頁
已閱讀1頁,還剩17頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

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

文檔簡介

1、<p><b>  目 錄</b></p><p><b>  1 概述 1</b></p><p>  1.1 課程設(shè)計(jì)目的 1</p><p>  1.2 課程設(shè)計(jì)任務(wù)和要求 1</p><p>  1.3 牛頓拉夫遜的基本計(jì)算原理 1</p><

2、;p>  1.4 課程題目 2</p><p>  2 手工計(jì)算 3</p><p>  3 系統(tǒng)程序 14</p><p>  3.1 MATLAB簡介 14</p><p>  3.2 系統(tǒng)牛頓-拉夫遜法潮流計(jì)算流程圖 14</p><p>  3.3 系統(tǒng)程序 15</p>

3、;<p><b>  總結(jié) 28</b></p><p><b>  參考文獻(xiàn) 29</b></p><p><b>  1 概述</b></p><p>  1.1 課程設(shè)計(jì)目的</p><p> ?。?)掌握電力牛頓-拉夫遜法潮流計(jì)算的基本原理;</

4、p><p> ?。?)掌握并能熟練運(yùn)用一門計(jì)算機(jī)語言(MATLAB語言、FORTRAN、C語言、C++語言);</p><p> ?。?)采用計(jì)算機(jī)語言對短路計(jì)算進(jìn)行計(jì)算機(jī)編程計(jì)算。</p><p>  通過課程設(shè)計(jì),使學(xué)生鞏固電力系統(tǒng)三相短路計(jì)算的基本原理與方法,掌握短路電流的數(shù)值求解方法(節(jié)點(diǎn)導(dǎo)納矩陣,修正方程),開發(fā)系統(tǒng)牛頓拉夫遜法的計(jì)算程序。讓學(xué)生掌握用計(jì)算機(jī)仿

5、真分析電力系統(tǒng)的方法。同時,通過軟件開發(fā),也有助于計(jì)算機(jī)操作能力和軟件開發(fā)能力的提高。</p><p>  1.2 課程設(shè)計(jì)任務(wù)和要求</p><p> ?。?)熟悉PSCAD軟件;</p><p> ?。?)編寫潮流計(jì)算流程圖;</p><p> ?。?)建立系統(tǒng)接線圖的仿真過程; </p><p> ?。?)得出仿

6、真結(jié)果。</p><p>  1.3 牛頓拉夫遜的基本計(jì)算原理</p><p>  設(shè)r是f(x)=0的根,選取x0作為r初始近似值,過點(diǎn)(x0,f(x0))做曲線y=f(x)的切線L,L的方程為y=f(x0)f'(x0)(x-x0),求出L與x軸交點(diǎn)的橫坐標(biāo)x1=x0-f(x0)/f'(x0),稱x1為r的一次近似值。過點(diǎn)(x1,f(x1))做曲線y=f(x)的切線,并求

7、該切線與x軸的橫坐標(biāo)x2=x1-f(x1)/f'(x1),稱x2為r的二次近似值。重復(fù)以上過程,得r的近似值序列,其中x(n+1)=x(n)-f(x(n))/f'(x(n)),稱為r的n+1次近似值,上式稱為牛頓迭代公式。</p><p>  解非線性方程f(x)=0的牛頓法是把非線性方程線性化的一種近似方法。把f(x)在x0點(diǎn)附近展開成泰勒級數(shù)f(x)=f(x0)+(x-x0)f'(x0

8、)+(x-x0)^2*f''(x0)/2!+…取其線性部分,作為非線性方程f(x)=0的近似方程,即泰勒展開的前兩項(xiàng),則有f(x0)+f'(x0)(x-x0)=f(x)=0,設(shè)f'(x0)≠0則其解為x1=x0-f(x0)/f'(x0)這樣,得到牛頓法的一個迭代序列:x(n+1)=x(n)-f(x(n))/f'(x(n))。</p><p><b>  1

9、.4 課程題目</b></p><p>  題目五:如圖2所示,一個5節(jié)點(diǎn)系統(tǒng),已知節(jié)點(diǎn)5為平衡節(jié)點(diǎn),節(jié)點(diǎn)1為節(jié)點(diǎn),其余為節(jié)點(diǎn)。以為基準(zhǔn)的標(biāo)幺值支路數(shù)據(jù)由表1給出。,,,,,,給定電壓的初始值如表2所示,收斂系數(shù)。試求:</p><p>  圖2 5節(jié)點(diǎn)簡單電力系統(tǒng)</p><p>  表1 網(wǎng)絡(luò)各元件參數(shù)的標(biāo)幺值</p><p&

10、gt;  表2 各節(jié)點(diǎn)電壓(初值)標(biāo)幺值參數(shù)</p><p> ?。?)利用牛頓-拉夫遜法計(jì)算圖2網(wǎng)絡(luò)的潮流分布。</p><p><b>  2 系統(tǒng)程序</b></p><p>  2.1PSCAD軟件簡介</p><p>  2.2 系統(tǒng)牛頓-拉夫遜法潮流計(jì)算流程圖</p><p>  參

11、照《電力系統(tǒng)分析》一書的計(jì)算過程及原理,可以畫出以下程序流程圖。</p><p>  圖2.1 程序流程圖</p><p>  3.3 系統(tǒng)程序及結(jié)果</p><p>  根據(jù)以上的程序流程圖可編寫出系統(tǒng)的計(jì)算程序,具體源代碼如下。</p><p><b>  clear</b></p><p>

12、;  NN=input('輸入總結(jié)點(diǎn)數(shù):');%輸入原始信息</p><p>  PV=input('輸入PV節(jié)點(diǎn)數(shù):');</p><p>  biao=input('PQ、PV、平衡排列按1,PV、PQ、平衡排列按2:');</p><p>  Z=input('輸入阻抗矩陣:');&

13、lt;/p><p>  Y0=input('輸入對地導(dǎo)納:');</p><p>  U=input('輸入電壓初值:');</p><p>  S=input('輸入節(jié)點(diǎn)功率初值(流出為負(fù),流入為正):');</p><p>  R=input('輸入要求精度:');</p&

14、gt;<p>  GE=input('輸入變壓器個數(shù):');</p><p>  if GE~=0%當(dāng)有變壓器時</p><p>  for i=1:GE</p><p>  C=input('歸高按1,歸低按2:');%求有變壓器支路的參數(shù)</p><p><b

15、>  if C==1</b></p><p>  k1=input('輸入變比k:');%輸入信息</p><p>  n1=input('輸入高節(jié)點(diǎn)號n:');</p><p>  m1=input('輸入低節(jié)點(diǎn)號m:'); </p><p>  B1=Z(n1,m1)

16、;</p><p>  Z(m1,n1)=Z(n1,m1)/k1;</p><p>  Z(n1,m1)=Z(n1,m1)/k1;</p><p>  Y0(n1,m1)=(1-k1)/B1;</p><p>  Y0(m1,n1)=(k1*(k1-1))/B1;</p><p><b>  end</

17、b></p><p><b>  if C==2 </b></p><p>  k2=input('輸入變比k:');%輸入信息</p><p>  n2=input('輸入低節(jié)點(diǎn)號n:');</p><p>  m2=input('輸入高節(jié)點(diǎn)號m:');<

18、;/p><p>  B2=Z(n2,m2);</p><p>  Z(m2,n2)=Z(n2,m2)*k2;</p><p>  Z(n2,m2)=Z(n2,m2)*k2;</p><p>  Y0(n2,m2)=(k2-1)/(B2*k2);</p><p>  Y0(m2,n2)=(1-k2)/(B2*k2*k2);&

19、lt;/p><p><b>  end</b></p><p><b>  end</b></p><p><b>  else</b></p><p><b>  ;</b></p><p><b>  end</b&

20、gt;</p><p>  for i=1:NN%求y(i,j)</p><p>  for j=1:NN</p><p>  if Z(i,j)==0</p><p><b>  y(i,j)=0;</b></p><p><b>  else </b>

21、</p><p>  y(i,j)=1/Z(i,j);</p><p><b>  end </b></p><p><b>  end</b></p><p><b>  end</b></p><p>  for i=1:NN%求

22、Y(i,j)</p><p>  for j=1:NN</p><p><b>  if i~=j</b></p><p>  Y(i,j)=-y(i,j);</p><p><b>  else </b></p><p>  Y(i,j)=sum(y(i,:))+sum(Y

23、0(i,:));</p><p><b>  end </b></p><p><b>  end</b></p><p><b>  end</b></p><p>  for i=1:NN%求yi0</p><p>  y00(i

24、)=sum(Y0(i,:));</p><p><b>  end</b></p><p>  E=real(U);</p><p>  F=imag(U);</p><p>  for i=1:PV%保存Vis</p><p>  if biao==2</p>

25、<p>  bao(i)=F(i)^2+E(i)^2;</p><p><b>  else</b></p><p>  bao(NN-i)=F(NN-i)^2+E(NN-i)^2;</p><p><b>  end</b></p><p><b>  end</b>

26、;</p><p>  G=real(Y);</p><p>  B=imag(Y);</p><p>  P=real(S);</p><p>  Q=imag(S); </p><p>  Error_Power=ones(2*(NN-1),1); </p><p>  count=0;

27、%迭代次數(shù) </p><p>  while max(abs(Error_Power))>R %迭代條件</p><p>  for i=1:NN-1 %求</p><p>  Error_P(i)=0;</p><p>  Error_Q(i)=0;</p><p><

28、b>  end </b></p><p>  for i=1:NN-1 </p><p>  for j=1:NN </p><p>  Error_P(i)=Error_P(i)-(E(i)*(E(j)*G(i,j)-F(j)*B(i,j))+F(i)*(F(j)*G(i,j)+E(j)*B(i,j))); </p&g

29、t;<p>  if biao==2%當(dāng)PV、PQ、平衡排列時</p><p><b>  if i>PV</b></p><p>  Error_Q(i)=Error_Q(i)-(F(i)*(E(j)*G(i,j)-F(j)*B(i,j))-E(i)*(F(j)*G(i,j)+E(j)*B(i,j))); </p>&

30、lt;p><b>  else </b></p><p>  Error_Q(i)=-(E(i)^2+F(i)^2);</p><p><b>  end</b></p><p>  else%當(dāng)PQ、PV、平衡排列時</p><p>  if i<NN-PV</p

31、><p>  Error_Q(i)=Error_Q(i)-(F(i)*(E(j)*G(i,j)-F(j)*B(i,j))-E(i)*(F(j)*G(i,j)+E(j)*B(i,j))); </p><p><b>  else </b></p><p>  Error_Q(i)=-(E(i)^2+F(i)^2);</p><p&

32、gt;<b>  end</b></p><p><b>  end</b></p><p><b>  end</b></p><p>  if biao==2%當(dāng)PV、PQ、平衡排列時</p><p><b>  if i>PV</b>

33、;</p><p>  Error_Q(i)=Error_Q(i)+Q(i); </p><p><b>  else</b></p><p>  Error_Q(i)=Error_Q(i)+bao(i); </p><p><b>  end</b></p><p>  e

34、lse %當(dāng)PQ、PV、平衡排列時</p><p>  if i<NN-PV</p><p>  Error_Q(i)=Error_Q(i)+Q(i); </p><p><b>  else</b></p><p>  Error_Q(i)=Error_Q(i)+bao(i); </p>

35、;<p><b>  end</b></p><p><b>  end</b></p><p>  Error_P(i)=Error_P(i)+P(i);</p><p><b>  end </b></p><p>  for i=1:(NN-1)<

36、/p><p>  Error_Power(2*(i-1)+1)=Error_P(i);%按序排列形成ΔW</p><p>  Error_Power(2*i)=Error_Q(i);</p><p><b>  end</b></p><p>  for i=1:NN%求雅可比矩陣</p>

37、;<p>  for j=1:NN</p><p>  T_P1Err(i,j)=0;%累計(jì)初值清0</p><p>  T_Q1Err(i,j)=0;</p><p>  T_P2Err(i,j)=0;</p><p>  T_Q2Err(i,j)=0;</p><p><b> 

38、 end </b></p><p><b>  end</b></p><p>  if biao==2%當(dāng)PV、PQ、平衡排列時</p><p>  for t=1:PV%對于PV節(jié)點(diǎn)</p><p>  for j=1:NN-1 </p><p>

39、  if j>PV%當(dāng)i~=j時</p><p>  T_P1(t,j)=-(G(t,j)*E(t)+B(t,j)*F(t));</p><p>  T_Q1(t,j)=B(t,j)*E(t)-G(t,j)*F(t);</p><p>  T_P2(t,j)=0;</p><p>  T_Q2(t,j)=0;<

40、/p><p>  else %當(dāng)i=j時</p><p>  for k=1:NN</p><p>  T_P1Err(t,j)=T_P1Err(t,j)-(G(t,k)*E(k)-B(t,k)*F(k));</p><p>  T_Q1Err(t,j)=T_Q1Err(t,j)-(G(t,k)*F(k)+B(t,k)*E(k));

41、</p><p><b>  end</b></p><p>  T_P1(t,j)=T_P1Err(t,j)-G(t,j)*E(t)-B(t,j)*F(t);</p><p>  T_Q1(t,j)=T_Q1Err(t,j)+B(t,j)*E(t)-G(t,j)*F(t);</p><p>  T_P2(t,j)=

42、-2*E(t);</p><p>  T_Q2(t,j)=-2*F(t);</p><p><b>  end</b></p><p>  J(2*t-1,2*j-1)=T_P1(t,j);%放入雅可比矩陣相應(yīng)處</p><p>  J(2*t-1,2*j)=T_Q1(t,j);</p><p&

43、gt;  J(2*t,2*j-1)=T_P2(t,j);</p><p>  J(2*t,2*j)=T_Q2(t,j);</p><p><b>  end</b></p><p><b>  end</b></p><p>  for i=PV+1:NN-1%對于PQ節(jié)點(diǎn)</p

44、><p>  for j=1:NN-1</p><p>  if i~=j%當(dāng)i~=j時</p><p>  T_P1(i,j)=-(G(i,j)*E(i)+B(i,j)*F(i));</p><p>  T_Q1(i,j)=B(i,j)*E(i)-G(i,j)*F(i);</p><p>  T_P2

45、(i,j)=B(i,j)*E(i)-G(i,j)*F(i);</p><p>  T_Q2(i,j)=G(i,j)*E(i)+B(i,j)*F(i);</p><p>  else%當(dāng)i=j時</p><p>  for k=1:NN</p><p>  T_P1Err(i,j)=T_P1Err(i,j)-(G(i,k)*E

46、(k)-B(i,k)*F(k));</p><p>  T_Q1Err(i,j)=T_Q1Err(i,j)-(G(i,k)*F(k)+B(i,k)*E(k));</p><p>  T_P2Err(i,j)=T_P2Err(i,j)+(G(i,k)*F(k)+B(i,k)*E(k));</p><p>  T_Q2Err(i,j)=T_Q2Err(i,j)-(G

47、(i,k)*E(k)-B(i,k)*F(k));</p><p><b>  end </b></p><p>  T_P1(i,j)=T_P1Err(i,j)-G(i,j)*E(i)-B(i,j)*F(i);</p><p>  T_Q1(i,j)=T_Q1Err(i,j)+B(i,j)*E(i)-G(i,j)*F(i);</p&

48、gt;<p>  T_P2(i,j)=T_P2Err(i,j)+B(i,j)*E(i)-G(i,j)*F(i);</p><p>  T_Q2(i,j)=T_Q2Err(i,j)+G(i,j)*E(i)+B(i,j)*F(i);</p><p><b>  end</b></p><p>  J(2*i-1,2*j-1)=T_P1

49、(i,j);%放入雅可比矩陣相應(yīng)處</p><p>  J(2*i-1,2*j)=T_Q1(i,j);</p><p>  J(2*i,2*j-1)=T_P2(i,j);</p><p>  J(2*i,2*j)=T_Q2(i,j);</p><p><b>  end</b></p><

50、p><b>  end</b></p><p>  else%當(dāng)PQ、PV、平衡排列時</p><p>  for t=NN-PV:NN-1%對于PV節(jié)點(diǎn)</p><p>  for j=1:NN-1</p><p>  if j<NN-PV%當(dāng)i~=j時<

51、/p><p>  T_P1(t,j)=-(G(t,j)*E(t)+B(t,j)*F(t));</p><p>  T_Q1(t,j)=B(t,j)*E(t)-G(t,j)*F(t);</p><p>  T_P2(t,j)=0;</p><p>  T_Q2(t,j)=0;</p><p>  else

52、%當(dāng)i=j時</p><p>  for k=1:NN</p><p>  T_P1Err(t,j)=T_P1Err(t,j)-(G(t,k)*E(k)-B(t,k)*F(k));</p><p>  T_Q1Err(t,j)=T_Q1Err(t,j)-(G(t,k)*F(k)+B(t,k)*E(k));</p><p><b>

53、  end</b></p><p>  T_P1(t,j)=T_P1Err(t,j)-G(t,j)*E(t)-B(t,j)*F(t);</p><p>  T_Q1(t,j)=T_Q1Err(t,j)+B(t,j)*E(t)-G(t,j)*F(t);</p><p>  T_P2(t,j)=-2*E(t);</p><p> 

54、 T_Q2(t,j)=-2*F(t);</p><p><b>  end</b></p><p>  J(2*t-1,2*j-1)=T_P1(t,j);%放入雅可比矩陣相應(yīng)處</p><p>  J(2*t-1,2*j)=T_Q1(t,j);</p><p>  J(2*t,2*j-1)=T_P2(t,j)

55、;</p><p>  J(2*t,2*j)=T_Q2(t,j);</p><p><b>  end</b></p><p><b>  end</b></p><p>  for i=1:NN-PV-1%對于PQ節(jié)點(diǎn)</p><p>  for j=1:NN-

56、1 </p><p>  if i~=j%當(dāng)i~=j時</p><p>  T_P1(i,j)=-(G(i,j)*E(i)+B(i,j)*F(i));</p><p>  T_Q1(i,j)=B(i,j)*E(i)-G(i,j)*F(i);</p><p>  T_P2(i,j)=B(i,j)*E(i

57、)-G(i,j)*F(i);</p><p>  T_Q2(i,j)=G(i,j)*E(i)+B(i,j)*F(i);</p><p>  Else%當(dāng)i=j時</p><p>  for k=1:NN</p><p>  T_P1Err(i,j)=T_P1Err(i,j)-(G(i,k)*E(k)-B(i,k)*F(k));

58、</p><p>  T_Q1Err(i,j)=T_Q1Err(i,j)-(G(i,k)*F(k)+B(i,k)*E(k));</p><p>  T_P2Err(i,j)=T_P2Err(i,j)+(G(i,k)*F(k)+B(i,k)*E(k));</p><p>  T_Q2Err(i,j)=T_Q2Err(i,j)-(G(i,k)*E(k)-B(i,k)

59、*F(k));</p><p><b>  end </b></p><p>  T_P1(i,j)=T_P1Err(i,j)-G(i,j)*E(i)-B(i,j)*F(i);</p><p>  T_Q1(i,j)=T_Q1Err(i,j)+B(i,j)*E(i)-G(i,j)*F(i);</p><p>  T

60、_P2(i,j)=T_P2Err(i,j)+B(i,j)*E(i)-G(i,j)*F(i);</p><p>  T_Q2(i,j)=T_Q2Err(i,j)+G(i,j)*E(i)+B(i,j)*F(i);</p><p><b>  end</b></p><p>  J(2*i-1,2*j-1)=T_P1(i,j);%放入雅可比矩

61、陣相應(yīng)處</p><p>  J(2*i-1,2*j)=T_Q1(i,j);</p><p>  J(2*i,2*j-1)=T_P2(i,j);</p><p>  J(2*i,2*j)=T_Q2(i,j);</p><p><b>  end</b></p><p><b>  e

62、nd</b></p><p><b>  end</b></p><p>  J;%雅可比矩陣</p><p>  Error_U=-(inv(J)*Error_Power); %求</p><p>  Error_E=zeros(1,NN-1); </p>&l

63、t;p>  Error_F=zeros(1,NN-1); </p><p>  for i=1:NN-1 </p><p>  Error_E(i)=Error_U(i*2-1); %求ΔE</p><p><b>  end </b></p><p>  for i=1:NN-1</p>

64、<p>  Error_F(i)=Error_U(i*2);%求ΔF</p><p><b>  end</b></p><p>  for i=1:NN-1 %求E</p><p>  E(i)=E(i)+Error_E(i);</p><p><b>  end

65、</b></p><p>  for i=1:NN-1%求F</p><p>  F(i)=F(i)+Error_F(i); </p><p><b>  end</b></p><p>  count=count+1; %count加1</p><p&g

66、t;<b>  end </b></p><p>  for j=1:NN %求從平衡節(jié)點(diǎn)的注入功率 </p><p>  P(NN)=E(NN)*(E(j)*G(NN,j)-F(j)*B(NN,j))+F(NN)*(F(j)*G(NN,j)+E(j)*B(NN,j))+P(NN); %平衡節(jié)點(diǎn)功率實(shí)部</p&

67、gt;<p>  Q(NN)=F(NN)*(E(j)*G(NN,j)-F(j)*B(NN,j))-E(NN)*(F(j)*G(NN,j)+E(j)*B(NN,j))+Q(NN); %平衡節(jié)點(diǎn)功率虛部</p><p><b>  end </b></p><p>  S(NN)=P(NN)+Q(NN)*sqrt(-1); %平衡節(jié)

68、點(diǎn)功率用復(fù)數(shù)表示 </p><p>  for i=1:NN %節(jié)點(diǎn)電壓用復(fù)數(shù)表示 </p><p>  Complex_U(i)=E(i)+F(i)*sqrt(-1); </p><p><b>  end </b></p><p>  for i=1:NN %線路功率</p

69、><p>  for j=1:NN </p><p>  if y(i,j)~=0 </p><p>  Line_S(i,j)=(Complex_U(i))^2*conj(y00(i))+Complex_U(i)*(conj(Complex_U(i))-conj(Complex_U(j)))*conj(y(i,j)); </p><

70、p><b>  else</b></p><p>  Line_S(i,j)=0;</p><p><b>  end </b></p><p><b>  end</b></p><p><b>  end </b></p>&l

71、t;p>  Y,Y0,count,Complex_U,Error_U,S,J,Line_S,Error_Power%輸出結(jié)果</p><p><b>  到此全部程序結(jié)束。</b></p><p><b>  以下是計(jì)算結(jié)果。</b></p><p><b>  Y=</b></p>

72、;<p>  10.7136-32.4317i -1.6667+5.0000i -1.7794+5.6940i -2.3810+7.1429i -5.0000 +15.0000i</p><p>  -1.6667+5.0000i 8.7990-32.2244i -5.8824 +23.5294i 0 -1.2500 + 3.7500i</p><p> 

73、 -1.7794+5.6940i -5.8824+23.5294i 8.9117-32.9184i -1.2500+3.7500i 0 </p><p>  -2.3810+7.1429i 0 -1.2500+3.7500i 3.7500-11.2250i 0 </p><p>  -5.0000+15.0000i -1.

74、2500+3.7500i 0 0 6.2500-18.7000i</p><p><b>  Y0=</b></p><p>  0 0+0.0200i 0+0.0200i -0.1134+0.3401i 0+0.0250i</p><p>  0+0

75、.0200i 0 0+0.0100i 0 0+0.0250i</p><p>  0+0.0200i 0+0.0100i 0 0+0.0250i 0 </p><p>  0.1190-0.3571i 0 0+0.

76、0250i 0 0 </p><p>  0+0.0250i 0+0.0250i 0 0 0 </p><p><b>  count=4</b></p><p>  Complex_U=</p

77、><p>  0.9992-0.0398i 0.9756-0.0806i 0.9701-0.0873i 0.9288-0.1007i 1.0500 </p><p><b>  Error_U=</b></p><p><b>  1.0e-007*</b></p><p>

78、<b>  -0.0055</b></p><p><b>  -0.0372</b></p><p><b>  -0.0515</b></p><p><b>  -0.0146</b></p><p><b>  -0.0640</b

79、></p><p><b>  -0.0145</b></p><p><b>  -0.1815</b></p><p><b>  0.0112</b></p><p><b>  S =</b></p><p>  0.

80、2000 -0.4500-0.1500i -0.4000-0.0500i -0.6000 - 0.1000i 1.3090+0.7229i</p><p><b>  J=</b></p><p>  -12.2160 -32.4620 1.8645 4.9297 2.0047 5.6186 2.6635 7.0424</p>

81、<p>  -1.9984 0.0796 0 0 0 0 0 0</p><p>  2.0292 4.7436 -10.7376 -30.9194 7.6364 22.4811 0 0</p><p>  4.7436 -2.0292 -30.5383 11.6286 22.4811 -7.6364 0 0</p&

82、gt;<p>  2.2231 5.3682 7.7602 22.3118 -11.1140 -31.2432 1.5399 3.5287</p><p>  5.3682 -2.2231 22.3118 -7.7602-31.0674 11.9228 3.5287 -1.5399</p><p>  2.9305 6.3946 0 0 1.5385 3.

83、3571 -3.9861 -10.2238</p><p>  6.3946 -2.9305 0 0 3.3571 -1.5385 -9.8726 5.2400</p><p><b>  Line_S=</b></p><p>  0 0.1000-0.3546i 0.1798-0.3263i 0.

84、4708-0.0610i -1.0183-0.9233i</p><p>  -0.2502-0.0811i 0 0.1830+0.0217i 0 -0.4088 - 0.1939i</p><p>  -0.3289-0.1022i-0.2005-0.1233i 0 0.1015+0.0737i 0

85、 </p><p>  -0.4319-0.0110i 0 0.0552+0.1429i 0 0 </p><p>  0.8939+0.5358i 0.4152+0.1319i 0 0 0 </p><p&

86、gt;  Error_Power=</p><p><b>  1.0e-006*</b></p><p><b>  -0.0492</b></p><p><b>  -0.0008</b></p><p><b>  -0.0003</b></

87、p><p><b>  -0.0125</b></p><p><b>  0.0013</b></p><p><b>  -0.0174</b></p><p><b>  -0.0208</b></p><p><b>

88、  -0.1732</b></p><p><b>  總結(jié)</b></p><p><b> ?。?)系統(tǒng)編程總結(jié)</b></p><p>  我組的任務(wù)是用牛頓—拉夫遜法去對五個節(jié)點(diǎn)的系統(tǒng)進(jìn)行潮流計(jì)算。具體來講,我是被分配在了程序設(shè)計(jì)這一塊,主要是進(jìn)行程序設(shè)計(jì)??梢哉f,首先接到這種題目,要用程序去進(jìn)行計(jì)算我

89、實(shí)在是不知從何下手,更何況是用自己不是很熟悉的MATLAB語言進(jìn)行。當(dāng)然,也可以用我們比較熟悉的C語言,但用C語言很多函數(shù)程序都得自己一一編寫,工作量要大得多,限于時間有限,加上MATLAB語言與C語言有很多相似之處,所以我最后還是選擇用MATLAB語言來編寫。鑒于之前對此語言沒有任何了解,所以一開始我將工作的重心放在了對MATLAB語言的認(rèn)識和了解。這一點(diǎn)我是通過去圖書館借了一些有關(guān)MATLAB的書及到網(wǎng)上下一些資料來完成這些初步工作

90、的。接下來就是畫出一個很概要的程序流程圖,然后根據(jù)一些資料來對自己認(rèn)為可行的一些程序進(jìn)行修改或者根據(jù)自己的思路來寫。當(dāng)我開始對這種語言有了更進(jìn)一步的了解時,我就對程序進(jìn)行自己的修改和編寫了。對于計(jì)算的過程的編寫,我是完全參照書本上的計(jì)算方法、步驟及公式來進(jìn)行的。這時,一個最令我頭痛的事來了,那就是好像所有資料及書本上的計(jì)算對于節(jié)點(diǎn)的分布都是先PQ,再PV最后是平衡節(jié)點(diǎn),而題目是所給的卻是先P</p><p>&l

溫馨提示

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

最新文檔

評論

0/150

提交評論