1、精選優質文檔-傾情為你奉上學 號： 7 電力系統計算機輔助潮流計算 實驗報告（學生姓名：學 院：電力學院系 別：電力系 專 業：電氣工程及其自動化 班 級：授課教師：二 一 一 年 十 二 月（2011-2012學年第一學期）1、實驗目的：了解計算機潮流分析的基本原理、主要步驟；掌握節點導納矩陣形成和修改的方法，掌握數據處理的基本方法；熟悉Matlab運行環境，了解Matlab基本編程語句和語法；運用潮流分析程序對給定網絡的運行方式做潮流分析，并初步分析計算結果2、實驗要求：通過預習，對計算機潮流分析基本理論有深入了解；為程序準備必要的、準確的原始數據；熟悉Matlab運行環境，輸入潮流程序，

2、上機獨立完成程序的調試，給出潮流分析的結果并按要求繪制潮流分布圖3、實驗內容：輸入網絡參數，包括節點號、節點導納矩陣、節點功率等；輸入潮流程序、調試并輸出結果，繪制潮流分布圖 4、實驗步驟：1、熟悉原始資料：根據計算要求，整理數據，包括：計算網絡中線路、變壓器的參數、形成節點導納矩陣；表示各節點的注入功率。（以上數據均采用有名值計算）2、讀通潮流程序：完成程序的解釋和說明，必要時附加對應的公式和程序語言的說明3、上機調試：熟悉Matlab的運行環境，準確輸入原始數據、節點編號、節點注入功率等信息4、整理計算結果：根據計算結果作電網潮流分布圖原始網絡：5、實驗數據及處理：一、實驗程序：clear

3、G(1,1)=3.75;B(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

4、)=-10;B(3,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 %形成節點導納矩陣delt(1)=0;delt(2)=0

5、;delt(3)=0; delt(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; %設迭代初值k=1;precision=1 %設迭代次數和精度N1=4; %PQ節點數while precision>0. %判斷是否滿足精度要求delt(5)=0;u(5)=1.06; %給定平衡節點編號for m=1:N1 for n=1:N1+1pt(n)=u(m)*u(n)*(G

6、(m,n)*cos(delt(m)-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); %計算PQ節點的注入功率 dp(m)=ps(m)-pi(m); dq(m)=qs(m)-qi(m); %計算PQ節點的功率不平衡量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)

7、2*G(m,m); J(m,m)=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

8、); N(m,n)=-J(m,n);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 %計算jocbi各項,并放入統一矩陣JJ中，對JJ下標統一編號JJ for m=1:N1 PP(2*m-1)=dp(m); PP(2*m)=dq(m); End %按統一矩陣形成功率不平衡uu=inv(JJ)*PP'precision=max(abs(uu); %判斷是否收斂for n=1:N1delt(n)=delt(n)+u

9、u(2*n-1);u(n)=u(n)+uu(2*n)*u(n); %將結果分解為電壓幅值和角度end %求解修正方程，得電壓幅值變化量（標幺值）和角度變化量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); %求平衡節點的注入電流endS5=U(5)*sum(conj(I) %求平衡節點的注入功率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);% endenden

10、d %求節點i,j節點之間的功率，方向為由i指向j, S %顯示支路功率二、實驗結果：1、節點導納矩陣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

11、 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i 0 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 6.2500 -18.7500i節點導納矩陣特點：1.節點導納矩陣的對角元就等于各該節點所連接導納的總和2.節點導納矩陣是稀疏矩陣3.節點導納矩陣一般是對稱矩陣4.節點導納矩陣的非對角元Yij等于連接節點i、j支路導納的負值5. 節點導納矩陣是方陣2、迭代過程數據：l 電壓變化量du = 0.0094 - 0.0010i 0.0448 - 0.0021i 0.0156 - 0.0013i 0.0142 - 0.001

12、3idu =-0.0077 - 0.0000i -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.40

13、00 - 0.0500idS = 0.0208 - 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

14、.7500 -1.2500 -3.7500 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.

15、7500 -1.2500 -5.0000 -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.151

16、1 -8.0396 -11.5380 35.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

17、.3861 3.1619 -7.6217 -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

18、.5953 39.2852 9.9876 -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.

19、2899 -1.5439 -5.2907 -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、收斂后數據：l 支路功率S = 0 -0