前推回代法計(jì)算流程蒼松書屋

1、 前推回代法計(jì)算流程要看懂前推回代法計(jì)算程序，報(bào)告敘述計(jì)算原理及計(jì)算流程。繪制計(jì)算流程框圖。確定前推回代支路次序（廣度優(yōu)先，或深度優(yōu)先），編寫前推回代計(jì)算輸入文件。進(jìn)行潮流計(jì)算。下列為節(jié)點(diǎn)配電網(wǎng)結(jié)構(gòu)圖及系統(tǒng)支路參數(shù)和系統(tǒng)負(fù)荷參數(shù)表。圖1-2 節(jié)點(diǎn)配電網(wǎng)結(jié)構(gòu)圖表1 系統(tǒng)支路參數(shù)支路r()x（.）b/2(tk) （s）123.3673.6850.0232.3562.5410.0341.1451.280.0454.5245.040.0260.8561.140.0672.7452.9650.0283.7434.2510.0892.2372.7560.03104.144.6960.03111.3281.

2、7630.011122.4362.8660.04133.5213.9660.0表2 系統(tǒng)負(fù)荷參數(shù)節(jié)點(diǎn)編號(hào)節(jié)點(diǎn)類型節(jié)點(diǎn)初始電壓（kv）pi（mva）qi（mva）1根節(jié)點(diǎn)10.4002中間節(jié)點(diǎn)10.00.03420.03013中間節(jié)點(diǎn)10.00.06930.06424中間節(jié)點(diǎn)10.00.08450.07635葉節(jié)點(diǎn)10.00.02950.02616中間節(jié)點(diǎn)10.00.04740.04097葉節(jié)點(diǎn)10.00.11760.09578中間節(jié)點(diǎn)10.00.09460.08579葉節(jié)點(diǎn)10.00.09160.085910葉節(jié)點(diǎn)10.00.02710.022911中間節(jié)點(diǎn)10.00.06960.064312

3、葉節(jié)點(diǎn)10.00.06760.057913葉節(jié)點(diǎn)10.00.02980.0242主程序清單：pq,ft,rx=case114(); %調(diào)用數(shù)據(jù)文件nn=size(pq,1); %節(jié)點(diǎn)數(shù)nb=size(ft,1); %支路數(shù)數(shù)v=pq(:,1); %v初始電壓相量maxd=1k=1while maxd>0.0001 pq2=pq; %每一次迭代各節(jié)點(diǎn)的注入有功和無(wú)功相同 pl=0.0;for i=1:nb kf=ft(i,1); %前推始節(jié)點(diǎn)號(hào) kt=ft(i,2); %前推終節(jié)點(diǎn)號(hào) x=(pq2(kf,2)2+pq2(kf,3)2)/v(kf)/v(kf); %計(jì)算沿線電流平方a pq1

4、(i,1)=pq2(kf,2)+rx(i,1)*x; %計(jì)算支路首端有功/mw rx(i,1)r pq1(i,2)=pq2(kf,3)+rx(i,2)*x; %計(jì)算沿支路的無(wú)功損耗/mvar rx(i,2)x pq2(kt,2)= pq2(kt,2)+pq1(i,1); %用pq1去修正支路末端節(jié)點(diǎn)的有功p 單位mw pq2(kt,3)= pq2(kt,3)+pq1(i,2); %用pq1去修正支路末端節(jié)點(diǎn)的有功q 單位mvar pl=pl+rx(i,1)*x;endangle(1)=0.0;for i=nb:-1:1 kf=ft(i,2); %回代始節(jié)點(diǎn)號(hào) kt=ft(i,1); %回代終節(jié)

5、點(diǎn)號(hào) dv1=(pq1(i,1)*rx(i,1)+pq1(i,2)*rx(i,2)/v(kf); %計(jì)算支路電壓損耗的縱分量dv1 dv2=(pq1(i,1)*rx(i,2)-pq1(i,2)*rx(i,1)/v(kf); %計(jì)算支路電壓損耗的橫分量dv2 v2(kt)=sqrt(v(kf)-dv1)2+dv22); %計(jì)算支路末端電壓/kv angle(kt)=angle(kf)+atand(dv2/(v(kf)-dv1); %計(jì)算支路endmaxd=abs(v2(2)-v(2);v2(1)=v(1);for i=3:1:nn if abs(v2(i)-v(i)>maxd; maxd=

6、abs(v2(i)-v(i); end endmaxdk=k+1pq1 %潮流分布 即支路首端潮流mvav=v2 %節(jié)點(diǎn)電壓模計(jì)算結(jié)果kvangle %節(jié)點(diǎn)電壓角度計(jì)算結(jié)果單位度pl %網(wǎng)損單位mwendclear輸入文件清單：function pq,ft,rx=case114()pq=%節(jié)點(diǎn)電壓 有功 無(wú)功10.4 0 010.0 0.0342 0.030110.0 0.0693 0.064210.0 0.0845 0.076310.0 0.0295 0.026110.0 0.0474 0.040910.0 0.1176 0.095710.0 0.0946 0.085710.0 0.0916

7、 0.085910.0 0.0271 0.022910.0 0.0696 0.064310.0 0.0676 0.057910.0 0.0298 0.0242;ft=%首端 末端 5 413 44 310 312 1111 37 66 29 88 23 22 1;rx=% r x 4.524 5.043.521 3.9661.145 1.284.14 4.6962.436 2.8661.328 1.7632.745 2.9650.856 1.142.237 2.7563.743 4.2512.356 2.5413.367 3.685;計(jì)算過(guò)程maxd = 1k = 1maxd = 0.1780

8、k = 2pq1 = 0.0296 0.0262 0.0299 0.0243 0.1443 0.1272 0.0272 0.0230 0.0678 0.0581 0.1378 0.1230 0.1182 0.0964 0.1660 0.1378 0.0920 0.0863 0.1890 0.1748 0.3847 0.3439 0.8099 0.7260v = columns 1 through 8 10.4000 9.8807 9.8220 9.9672 9.9734 9.9701 9.9390 9.8550 columns 9 through 13 9.9556 9.9780 9.9600

9、 9.9668 9.9799angle = columns 1 through 8 0 0.3011 0.3986 0.4211 0.4387 0.3421 0.3916 0.3878 columns 9 through 13 0.4225 0.4173 0.4444 0.4747 0.4400pl = 0.0471maxd = 0.1787k = 3pq1 = 0.0296 0.0262 0.0299 0.0243 0.1443 0.1272 0.0272 0.0230 0.0678 0.0581 0.1378 0.1230 0.1182 0.0964 0.1660 0.1378 0.092

10、0 0.0863 0.1890 0.1749 0.3849 0.3442 0.8112 0.7274v = columns 1 through 8 10.4000 9.8798 9.7004 9.7886 9.9405 9.8504 9.9089 9.7338 columns 9 through 13 9.8100 9.7996 9.7813 9.9267 9.9470angle = columns 1 through 8 0 0.3011 0.4011 0.4244 0.4421 0.3431 0.3929 0.3899 columns 9 through 13 0.4257 0.4204

11、0.4486 0.4791 0.4434pl = 0.0484maxd = 0.1793k = 4pq1 = 0.0296 0.0262 0.0299 0.0243 0.1444 0.1272 0.0272 0.0230 0.0678 0.0581 0.1379 0.1231 0.1182 0.0964 0.1661 0.1378 0.0920 0.0864 0.1891 0.1749 0.3851 0.3444 0.8115 0.7277v = columns 1 through 8 10.4000 9.8796 9.6994 9.6666 9.7614 9.8495 9.7884 9.73

12、29 columns 9 through 13 9.6883 9.6777 9.6591 9.7474 9.7680angle = columns 1 through 8 0 0.3011 0.4011 0.4250 0.4433 0.3431 0.3942 0.3899 columns 9 through 13 0.4266 0.4209 0.4498 0.4814 0.4447pl = 0.0487maxd = 0.1226k = 5pq1 = 0.0296 0.0262 0.0299 0.0243 0.1444 0.1272 0.0272 0.0230 0.0678 0.0581 0.1

13、379 0.1231 0.1183 0.0964 0.1661 0.1379 0.0920 0.0864 0.1891 0.1750 0.3852 0.3444 0.8115 0.7278v = columns 1 through 8 10.4000 9.8795 9.6991 9.6656 9.6391 9.8493 9.7875 9.7326 columns 9 through 13 9.6873 9.6767 9.6581 9.6248 9.6457angle = columns 1 through 8 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.394

14、2 0.3899 columns 9 through 13 0.4266 0.4209 0.4498 0.4822 0.4452pl = 0.0487maxd = 0.0010k = 6pq1 = 0.0296 0.0262 0.0299 0.0243 0.1444 0.1273 0.0272 0.0230 0.0678 0.0581 0.1379 0.1231 0.1183 0.0964 0.1661 0.1379 0.0920 0.0864 0.1891 0.1750 0.3852 0.3445 0.8116 0.7278v = columns 1 through 8 10.4000 9.

15、8795 9.6991 9.6653 9.6381 9.8492 9.7873 9.7326 columns 9 through 13 9.6870 9.6764 9.6579 9.6238 9.6447angle = columns 1 through 8 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 columns 9 through 13 0.4266 0.4209 0.4498 0.4823 0.4452pl = 0.0488maxd = 2.6021e-004k = 7pq1 = 0.0296 0.0262 0.0299 0.0

16、243 0.1444 0.1273 0.0272 0.0230 0.0678 0.0581 0.1379 0.1231 0.1183 0.0964 0.1661 0.1379 0.0920 0.0864 0.1891 0.1750 0.3852 0.3445 0.8116 0.7278v = columns 1 through 8 10.4000 9.8795 9.6991 9.6652 9.6378 9.8492 9.7872 9.7326 columns 9 through 13 9.6870 9.6764 9.6578 9.6235 9.6445angle = columns 1 thr

17、ough 8 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 columns 9 through 13 0.4266 0.4209 0.4498 0.4823 0.4452pl = 0.0488maxd = 6.1046e-005k = 8pq1 = 0.0296 0.0262 0.0299 0.0243 0.1444 0.1273 0.0272 0.0230 0.0678 0.0581 0.1379 0.1231 0.1183 0.0964 0.1661 0.1379 0.0920 0.0864 0.1891 0.1750 0.3852

18、0.3445 0.8116 0.7278v = columns 1 through 8 10.4000 9.8795 9.6991 9.6652 9.6377 9.8492 9.7872 9.7326 columns 9 through 13 9.6870 9.6764 9.6578 9.6235 9.6444angle = columns 1 through 8 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 columns 9 through 13 0.4266 0.4209 0.4498 0.4823 0.4452pl = 0.0488計(jì)算結(jié)果清單：maxd = 1k = 1maxd = 0.1780k = 2maxd = 0.1787k = 3maxd = 0.1793k = 4maxd = 0.1226k = 5maxd = 0.0010k = 6maxd = 2.6021e-04k = 7maxd = 6.1046e-05k = 8pq1 = 0.0296 0.0262 0.0299 0.0243 0.1444 0.1273 0.0272 0.0230 0.0678 0.0581 0.1379 0.1231 0.1183 0.0964 0.1661 0.1379 0.0920 0.0864 0.1891 0.