版權(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)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)
- 電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)
- 電力系統(tǒng)穩(wěn)態(tài)潮流計(jì)算課程設(shè)計(jì)
- 電力系統(tǒng)穩(wěn)態(tài)課程設(shè)計(jì)--潮流計(jì)算
- 電力系統(tǒng)課程設(shè)計(jì)--電力系統(tǒng)規(guī)劃設(shè)計(jì)
- 電力系統(tǒng)課程設(shè)計(jì)--電力系統(tǒng)穩(wěn)定分析和計(jì)算
- 電力系統(tǒng)課程設(shè)計(jì)--電力系統(tǒng)穩(wěn)定分析和計(jì)算
- 電力系統(tǒng)暫態(tài)分析課程設(shè)計(jì)
- 電力系統(tǒng)暫態(tài)分析課程設(shè)計(jì)
- 電力網(wǎng)損耗分析論文-《電力系統(tǒng)穩(wěn)態(tài)分析》 _課程設(shè)計(jì)
- 電力系統(tǒng)課程設(shè)計(jì)--電力系統(tǒng)穩(wěn)定分析和計(jì)算設(shè)計(jì)
- 電力系統(tǒng)課程設(shè)計(jì)
- 電力系統(tǒng)分析課程設(shè)計(jì)電力系統(tǒng)短路計(jì)算
- 電力系統(tǒng)課程設(shè)計(jì)
- 電力系統(tǒng)分析課程設(shè)計(jì)---電力系統(tǒng)短路計(jì)算
- 電力系統(tǒng)分析課程設(shè)計(jì)-- 電力系統(tǒng)潮流計(jì)算
- 電力系統(tǒng)分析課程設(shè)計(jì)--電力系統(tǒng)潮流計(jì)算
- 電力系統(tǒng)分析課程設(shè)計(jì)電力系統(tǒng)短路計(jì)算
- 答案 電力系統(tǒng)穩(wěn)態(tài)分析
- 課程設(shè)計(jì)---電力系統(tǒng)分析
評論
0/150
提交評論