電力系統(tǒng)計(jì)算機(jī)輔助潮流計(jì)算--實(shí)驗(yàn)報(bào)告_第1頁(yè)
電力系統(tǒng)計(jì)算機(jī)輔助潮流計(jì)算--實(shí)驗(yàn)報(bào)告_第2頁(yè)
電力系統(tǒng)計(jì)算機(jī)輔助潮流計(jì)算--實(shí)驗(yàn)報(bào)告_第3頁(yè)
電力系統(tǒng)計(jì)算機(jī)輔助潮流計(jì)算--實(shí)驗(yàn)報(bào)告_第4頁(yè)
電力系統(tǒng)計(jì)算機(jī)輔助潮流計(jì)算--實(shí)驗(yàn)報(bào)告_第5頁(yè)
已閱讀5頁(yè)，還剩5頁(yè)未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

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

文檔簡(jiǎn)介

1、學(xué) 號(hào)： 7 電力系統(tǒng)計(jì)算機(jī)輔助潮流計(jì)算 實(shí)驗(yàn)報(bào)告（學(xué)生姓名：學(xué) 院：電力學(xué)院系 別：電力系 專 業(yè)：電氣工程及其自動(dòng)化 班 級(jí)：授課教師：二 一 一 年 十 二 月（2011-2012學(xué)年第一學(xué)期）1、實(shí)驗(yàn)?zāi)康模毫私庥?jì)算機(jī)潮流分析的基本原理、主要步驟；掌握節(jié)點(diǎn)導(dǎo)納矩陣形成和修改的方法，掌握數(shù)據(jù)處理的基本方法；熟悉Matlab運(yùn)行環(huán)境，了解Matlab基本編程語(yǔ)句和語(yǔ)法；運(yùn)用潮流分析程序?qū)o定網(wǎng)絡(luò)的運(yùn)行方式做潮流分析，并初步分析計(jì)算結(jié)果2、實(shí)驗(yàn)要求：通過(guò)預(yù)習(xí)，對(duì)計(jì)算機(jī)潮流分析基本理論有深入了解；為程序準(zhǔn)備必要的、準(zhǔn)確的原始數(shù)據(jù)；熟悉Matlab運(yùn)行環(huán)境，輸入潮流程序，上機(jī)獨(dú)立完成程序的調(diào)試，給

2、出潮流分析的結(jié)果并按要求繪制潮流分布圖3、實(shí)驗(yàn)內(nèi)容：輸入網(wǎng)絡(luò)參數(shù)，包括節(jié)點(diǎn)號(hào)、節(jié)點(diǎn)導(dǎo)納矩陣、節(jié)點(diǎn)功率等；輸入潮流程序、調(diào)試并輸出結(jié)果，繪制潮流分布圖 4、實(shí)驗(yàn)步驟：1、熟悉原始資料：根據(jù)計(jì)算要求，整理數(shù)據(jù)，包括：計(jì)算網(wǎng)絡(luò)中線路、變壓器的參數(shù)、形成節(jié)點(diǎn)導(dǎo)納矩陣；表示各節(jié)點(diǎn)的注入功率。（以上數(shù)據(jù)均采用有名值計(jì)算）2、讀通潮流程序：完成程序的解釋和說(shuō)明，必要時(shí)附加對(duì)應(yīng)的公式和程序語(yǔ)言的說(shuō)明3、上機(jī)調(diào)試：熟悉Matlab的運(yùn)行環(huán)境，準(zhǔn)確輸入原始數(shù)據(jù)、節(jié)點(diǎn)編號(hào)、節(jié)點(diǎn)注入功率等信息4、整理計(jì)算結(jié)果：根據(jù)計(jì)算結(jié)果作電網(wǎng)潮流分布圖原始網(wǎng)絡(luò)：5、實(shí)驗(yàn)數(shù)據(jù)及處理：一、實(shí)驗(yàn)程序：clearG(1,1)=3.75;B

3、(1,1)=-11.25;G(1,2)=-2.5;B(1,2)=7.5;G(1,3)=0;B(1,3)=0;G(1,4)=-1.25;B(1,4)=3.75;G(1,5)=0;B(1,5)=0;G(2,1)=-2.5;B(2,1)=7.5;G(2,2)=10.834;B(2,2)=-32.5;G(2,3)=-1.667;B(2,3)=5;G(2,4)=-1.667;B(2,4)=5;G(2,5)=-5;B(2,5)=15;G(3,1)=0;B(3,1)=0;G(3,2)=-1.667;B(3,2)=5;G(3,3)=12.917;B(3,3)=-38.75;G(3,4)=-10;B(3,4)=

4、30;G(3,5)=-1.25;B(3,5)=3.75;G(4,1)=-1.25;B(4,1)=3.75;G(4,2)=-1.667;B(4,2)=5;G(4,3)=-10;B(4,3)=30;G(4,4)=12.917;B(4,4)=-38.750;G(4,5)=0;B(4,5)=0;G(5,1)=0;B(5,1)=0;G(5,2)=-5;B(5,2)=15;G(5,3)=-1.25;B(5,3)=3.75;G(5,4)=0;B(5,4)=0;G(5,5)=6.25;B(5,5)=-18.75;Y=G+j*B %形成節(jié)點(diǎn)導(dǎo)納矩陣delt(1)=0;delt(2)=0;delt(3)=0; d

5、elt(4)=0; u(1)=1.0;u(2)=1.0;u(3)=1.0;u(4)=1.0; ps(1)=-0.6;qs(1)=-0.10;ps(2)=0.2;qs(2)=0.2; ps(3)=-0.45;qs(3)=-0.15; ps(4)=-0.4;qs(4)=-0.05; %設(shè)迭代初值k=1;precision=1 %設(shè)迭代次數(shù)和精度N1=4; %PQ節(jié)點(diǎn)數(shù)while precision0. %判斷是否滿足精度要求delt(5)=0;u(5)=1.06; %給定平衡節(jié)點(diǎn)編號(hào)for m=1:N1 for n=1:N1+1pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)

6、-delt(n)+B(m,n)*sin(delt(m)-delt(n);qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n); end pi(m)=sum(pt);qi(m)=sum(qt); %計(jì)算PQ節(jié)點(diǎn)的注入功率 dp(m)=ps(m)-pi(m); dq(m)=qs(m)-qi(m); %計(jì)算PQ節(jié)點(diǎn)的功率不平衡量endfor m=1:N1 for n=1:N1 if m=n H(m,m)=-qi(m)-u(m)2*B(m,m); N(m,m)=pi(m)+u(m)2*G(m,m); J(m,m)=

7、pi(m)-u(m)2*G(m,m); L(m,m)=qi(m)-u(m)2*B(m,m); JJ(2*m-1,2*m-1)=H(m,m); JJ(2*m-1,2*m)=N(m,m); JJ(2*m,2*m-1)=J(m,m); JJ(2*m,2*m)=L(m,m); else H(m,n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n)-B(m,n)*cos(delt(m)-delt(n);J(m,n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n)+B(m,n)*sin(delt(m)-delt(n); N(m,n)=-J(m,n)

8、;L(m,n)=H(m,n); JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n); JJ(2*m,2*n-1)=J(m,n); JJ(2*m,2*n)=L(m,n); End endend %計(jì)算jocbi各項(xiàng),并放入統(tǒng)一矩陣JJ中，對(duì)JJ下標(biāo)統(tǒng)一編號(hào)JJ for m=1:N1 PP(2*m-1)=dp(m); PP(2*m)=dq(m); End %按統(tǒng)一矩陣形成功率不平衡uu=inv(JJ)*PP;precision=max(abs(uu); %判斷是否收斂for n=1:N1delt(n)=delt(n)+uu(2*n-1);u(n)=u(n)+uu

9、(2*n)*u(n); %將結(jié)果分解為電壓幅值和角度end %求解修正方程，得電壓幅值變化量（標(biāo)幺值）和角度變化量k=k+1; end for n=1:N1+1U(n)=u(n)*(cos(delt(n)+j*sin(delt(n); endfor m=1:N1+1 I(m)=Y(5,m)*U(m); %求平衡節(jié)點(diǎn)的注入電流endS5=U(5)*sum(conj(I) %求平衡節(jié)點(diǎn)的注入功率for m=1:N1+1 for n=1:N1+1 S(m,n)=U(m)*(conj(U(m)-conj(U(n)*conj(-Y(m,n);% endendend %求節(jié)點(diǎn)i,j節(jié)點(diǎn)之間的功率，方向?yàn)橛?/p>

10、i指向j, S %顯示支路功率二、實(shí)驗(yàn)結(jié)果：1、節(jié)點(diǎn)導(dǎo)納矩陣Y =3.7500 -11.2500i -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 0 -2.5000 + 7.5000i 10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -5.0000 +15.0000i 0 -1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i -1.2500 + 3.7500i-1.2500 + 3.7500i -1.6670 + 5.0000i -10.0000 +

11、30.0000i 12.9170 -38.7500i 0 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 6.2500 -18.7500i節(jié)點(diǎn)導(dǎo)納矩陣特點(diǎn)：1.節(jié)點(diǎn)導(dǎo)納矩陣的對(duì)角元就等于各該節(jié)點(diǎn)所連接導(dǎo)納的總和2.節(jié)點(diǎn)導(dǎo)納矩陣是稀疏矩陣3.節(jié)點(diǎn)導(dǎo)納矩陣一般是對(duì)稱矩陣4.節(jié)點(diǎn)導(dǎo)納矩陣的非對(duì)角元Yij等于連接節(jié)點(diǎn)i、j支路導(dǎo)納的負(fù)值5. 節(jié)點(diǎn)導(dǎo)納矩陣是方陣2、迭代過(guò)程數(shù)據(jù)：l 電壓變化量du = 0.0094 - 0.0010i 0.0448 - 0.0021i 0.0156 - 0.0013i 0.0142 - 0.0013idu =-0.0077 - 0.000

12、0i -0.0064 - 0.0000i -0.0065 - 0.0000i -0.0068 - 0.0000idu = -0.6060 - 0.0000i -0.4686 - 0.0000i -0.5040 - 0.0000i -0.5213 - 0.0000idu = -0.3934 - 0.0000i -0.2741 + 0.0000i -0.3107 - 0.0000i -0.3238 - 0.0000il 功率不平衡量dS = -0.6000 - 0.1000i 0.5000 + 1.1000i -0.3750 + 0.0750i -0.4000 - 0.0500idS = 0.02

13、08 - 0.0206i -0.0468 - 0.0876i 0.0047 - 0.0226i 0.0155 - 0.0096idS= 1.0e-003 *0.2233 - 0.2193i -0.4990 - 0.5087i -0.0052 - 0.1658i 0.1556 - 0.0997idS =1.0e-007 *0.2011 - 0.1958i -0.4104 - 0.2076i -0.0402 - 0.0900i 0.1273 - 0.0836il 雅可比矩陣JJ = 11.2500 3.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500 -3.7500

14、 11.2500 2.5000 -7.5000 0 0 1.2500 -3.7500 -7.5000 -2.5000 33.4000 10.5340 -5.0000 -1.6670 -5.0000 -1.6670 2.5000 -7.5000 -11.1340 31.6000 1.6670 -5.0000 1.6670 -5.0000 0 0 -5.0000 -1.6670 38.9750 12.8420 -30.0000 -10.0000 0 0 1.6670 -5.0000 -12.9920 38.5250 10.0000 -30.0000 -3.7500 -1.2500 -5.0000

15、-1.6670 -30.0000 -10.0000 38.7500 12.91701.2500 -3.7500 1.6670 -5.0000 10.0000 -30.0000 -12.9170 38.7500JJ = 11.5406 3.1996 -7.7223 -3.1029 0 0 -3.8183 -1.3384 -4.4412 11.3818 3.1029 -7.7223 0 0 1.3384 -3.8183 -8.0396 -2.1511 35.0648 12.0317 -5.3599 -1.5576 -5.3622 -1.5238 2.1511 -8.0396 -11.5380 35

16、.6400 1.5576 -5.3599 1.5238 -5.3622 0 0 -5.2222 -1.9705 40.0793 12.8630 -30.9519 -10.1136 0 0 1.9705 -5.2222 -13.7724 39.8246 10.1136 -30.9519 -3.8576 -1.2203 -5.2038 -1.9989 -30.8297 -10.4802 39.8912 12.8685 1.2203 -3.8576 1.9989 -5.2038 10.4802 -30.8297 -13.6994 39.8104JJ = 11.3861 3.1619 -7.6217

17、-3.0453 0 0 -3.7644 -1.3171 -4.3623 11.1866 3.0453 -7.6217 0 0 1.3171 -3.7644 -7.9246 -2.1368 34.7163 11.8401 -5.2904 -1.5440 -5.2913 -1.5116 2.1368 -7.9246 -11.4391 35.1173 1.5440 -5.2904 1.5116 -5.2913 0 0 -5.1585 -1.9397 39.5849 12.6953 -30.5425 -9.9876 0 0 1.9397 -5.1585 -13.5953 39.2852 9.9876

18、-30.5425 -3.8018 -1.2050 -5.1397 -1.9661 -30.4266 -10.3354 39.3681 12.7062 1.2050 -3.8018 1.9661 -5.1397 10.3354 -30.4266 -13.5065 39.2683JJ = 11.3850 3.1617 -7.6210 -3.0448 0 0 -3.7640 -1.3169 -4.3617 11.1850 3.0448 -7.6210 0 0 1.3169 -3.7640 -7.9237 -2.1368 34.7136 11.8386 -5.2899 -1.5439 -5.2907

19、-1.5115 2.1368 -7.9237 -11.4386 35.1136 1.5439 -5.2899 1.5115 -5.2907 0 0 -5.1580 -1.9394 39.5811 12.6940 -30.5394 -9.9867 0 0 1.9394 -5.1580 -13.5940 39.2811 9.9867 -30.5394 -3.8013 -1.2049 -5.1392 -1.9659 -30.4235 -10.3343 39.3641 12.7050 1.2049 -3.8013 1.9659 -5.1392 10.3343 -30.4235 -13.5050 39.26413、收斂后數(shù)據(jù)：l 支路功率S =

溫馨提示

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

評(píng)論

0/150

提交評(píng)論