前推回代法计算流程.docx

1、1前推回代法计算流程要看懂前推回代法计算程序，报告叙述计算原理及计算流程。绘制计算流程框图。确定前推回代支路次序（广度优先，或深度优先） ，编写前推回代计算输入文件。进行潮流计算。下列为节点配电网结构图及系统支路参数和系统负荷参数表。 1234567891012图 1-2 节点配电网结构图表 1 系统支路参数支路 R() X（.） B/2(TK) （ S）12 3.367 3.685 0.023 2.356 2.541 0.034 1.145 1.28 0.045 4.524 5.04 0.026 0.856 1.14 0.067 2.745 2.965 0.028 3.743 4.251 0

2、.089 2.237 2.756 0.0310 4.14 4.696 0.0311 1.328 1.763 0.01112 2.436 2.866 0.0413 3.521 3.966 0.0表 2 系统负荷参数节点编号 节点类型 节点初始电压（kV）Pi（MVA ） Qi（MVA）21 根节点 10.4 0 02 中间节点 10.0 0.0342 0.03013 中间节点 10.0 0.0693 0.06424 中间节点 10.0 0.0845 0.07635 叶节点 10.0 0.0295 0.02616 中间节点 10.0 0.0474 0.04097 叶节点 10.0 0.1176 0

3、.09578 中间节点 10.0 0.0946 0.08579 叶节点 10.0 0.0916 0.085910 叶节点 10.0 0.0271 0.022911 中间节点 10.0 0.0696 0.064312 叶节点 10.0 0.0676 0.057913 叶节点 10.0 0.0298 0.0242主程序清单：PQ,FT,RX=case114(); %调用数据文件NN=size(PQ,1); %节点数NB=size(FT,1); %支路数数V=PQ(:,1); %V 初始电压相量maxd=1k=1while maxd0.0001PQ2=PQ; %每一次迭代各节点的注入有功和无功相同P

4、L=0.0;for i=1:NBkf=FT(i,1); %前推始节点号kt=FT(i,2); %前推终节点号x=(PQ2(kf,2)2+PQ2(kf,3)2)/V(kf)/V(kf); %计算沿线电流平方 APQ1(i,1)=PQ2(kf,2)+RX(i,1)*x; %计算支路首端有功 /MW RX(i,1)RPQ1(i,2)=PQ2(kf,3)+RX(i,2)*x; %计算沿支路的无功损耗 /Mvar RX(i,2)XPQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用 PQ1 去修正支路末端节点的有功 P 单位 MWPQ2(kt,3)= PQ2(kt,3)+PQ1(i,2);

5、 %用 PQ1 去修正支路末端节点的有功 Q 单位 Mvar3PL=PL+RX(i,1)*x;endangle(1)=0.0;for i=NB:-1:1kf=FT(i,2); %回代始节点号 kt=FT(i,1); %回代终节点号dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf); %计算支路电压损耗的纵分量 dv1dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf); %计算支路电压损耗的横分量 dv2V2(kt)=sqrt(V(kf)-dv1)2+dv22); %计算支路末端电压/kVangle(kt)=angle(

6、kf)+atand(dv2/(V(kf)-dv1); %计算支路endmaxd=abs(V2(2)-V(2);V2(1)=V(1);for i=3:1:NNif abs(V2(i)-V(i)maxd;maxd=abs(V2(i)-V(i); end endmaxdk=k+1PQ1 %潮流分布 即支路首端潮流 MVAV=V2 %节点电压模计算结果 kVangle %节点电压角度计算结果单位度PL %网损单位 MWendclear输入文件清单：function PQ,FT,RX=case114()PQ=%节点电压 有功 无功10.4 0 0410.0 0.0342 0.030110.0 0.069

7、3 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 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.8

8、661.328 1.7632.745 2.9650.856 1.142.237 2.7563.743 4.2512.356 2.5413.367 3.685;计算过程maxd =15k =1maxd =0.1780k =2PQ1 =0.0296 0.02620.0299 0.02430.1443 0.12720.0272 0.02300.0678 0.05810.1378 0.12300.1182 0.09640.1660 0.13780.0920 0.08630.1890 0.17480.3847 0.34390.8099 0.7260V =Columns 1 through 810.400

9、0 9.8807 9.8220 9.9672 9.9734 9.9701 9.9390 9.8550Columns 9 through 139.9556 9.9780 9.9600 9.9668 9.9799angle =6Columns 1 through 80 0.3011 0.3986 0.4211 0.4387 0.3421 0.3916 0.3878Columns 9 through 130.4225 0.4173 0.4444 0.4747 0.4400PL =0.0471maxd =0.1787k =3PQ1 =0.0296 0.02620.0299 0.02430.1443 0

10、.12720.0272 0.02300.0678 0.05810.1378 0.12300.1182 0.09640.1660 0.13780.0920 0.08630.1890 0.17490.3849 0.34420.8112 0.7274V =Columns 1 through 810.4000 9.8798 9.7004 9.7886 9.9405 9.8504 9.9089 9.73387Columns 9 through 139.8100 9.7996 9.7813 9.9267 9.9470angle =Columns 1 through 80 0.3011 0.4011 0.4

11、244 0.4421 0.3431 0.3929 0.3899Columns 9 through 130.4257 0.4204 0.4486 0.4791 0.4434PL =0.0484maxd =0.1793k =4PQ1 =0.0296 0.02620.0299 0.02430.1444 0.12720.0272 0.02300.0678 0.05810.1379 0.12310.1182 0.09640.1661 0.13780.0920 0.08640.1891 0.17490.3851 0.34440.8115 0.72778V =Columns 1 through 810.40

12、00 9.8796 9.6994 9.6666 9.7614 9.8495 9.7884 9.7329Columns 9 through 139.6883 9.6777 9.6591 9.7474 9.7680angle =Columns 1 through 80 0.3011 0.4011 0.4250 0.4433 0.3431 0.3942 0.3899Columns 9 through 130.4266 0.4209 0.4498 0.4814 0.4447PL =0.0487maxd =0.1226k =5PQ1 =0.0296 0.02620.0299 0.02430.1444 0

13、.12720.0272 0.023090.0678 0.05810.1379 0.12310.1183 0.09640.1661 0.13790.0920 0.08640.1891 0.17500.3852 0.34440.8115 0.7278V =Columns 1 through 810.4000 9.8795 9.6991 9.6656 9.6391 9.8493 9.7875 9.7326Columns 9 through 139.6873 9.6767 9.6581 9.6248 9.6457angle =Columns 1 through 80 0.3011 0.4011 0.4

14、250 0.4438 0.3431 0.3942 0.3899Columns 9 through 130.4266 0.4209 0.4498 0.4822 0.4452PL =0.0487maxd =0.0010k =610PQ1 =0.0296 0.02620.0299 0.02430.1444 0.12730.0272 0.02300.0678 0.05810.1379 0.12310.1183 0.09640.1661 0.13790.0920 0.08640.1891 0.17500.3852 0.34450.8116 0.7278V =Columns 1 through 810.4

15、000 9.8795 9.6991 9.6653 9.6381 9.8492 9.7873 9.7326Columns 9 through 139.6870 9.6764 9.6579 9.6238 9.6447angle =Columns 1 through 80 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899Columns 9 through 130.4266 0.4209 0.4498 0.4823 0.4452PL =0.048811maxd =2.6021e-004k =7PQ1 =0.0296 0.02620.0299 0.0243

16、0.1444 0.12730.0272 0.02300.0678 0.05810.1379 0.12310.1183 0.09640.1661 0.13790.0920 0.08640.1891 0.17500.3852 0.34450.8116 0.7278V =Columns 1 through 810.4000 9.8795 9.6991 9.6652 9.6378 9.8492 9.7872 9.7326Columns 9 through 139.6870 9.6764 9.6578 9.6235 9.6445angle =Columns 1 through 80 0.3011 0.4

17、011 0.4250 0.4438 0.3431 0.3942 0.3899Columns 9 through 13120.4266 0.4209 0.4498 0.4823 0.4452PL =0.0488maxd =6.1046e-005k =8PQ1 =0.0296 0.02620.0299 0.02430.1444 0.12730.0272 0.02300.0678 0.05810.1379 0.12310.1183 0.09640.1661 0.13790.0920 0.08640.1891 0.17500.3852 0.34450.8116 0.7278V =Columns 1 t

18、hrough 810.4000 9.8795 9.6991 9.6652 9.6377 9.8492 9.7872 9.7326Columns 9 through 139.6870 9.6764 9.6578 9.6235 9.6444angle =13Columns 1 through 80 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899Columns 9 through 130.4266 0.4209 0.4498 0.4823 0.4452PL =0.0488计算结果清单：maxd = 1k = 1maxd = 0.1780k = 2ma

19、xd = 0.1787k = 3maxd = 0.1793k = 4maxd = 0.1226k = 5maxd = 0.0010k = 6maxd = 2.6021e-04k = 7maxd = 6.1046e-05k = 8PQ1 =0.0296 0.02620.0299 0.02430.1444 0.12730.0272 0.02300.0678 0.05810.1379 0.12310.1183 0.09640.1661 0.13790.0920 0.08640.1891 0.17500.3852 0.34450.8116 0.7278V =1410.4000 9.8795 9.6991 9.6652 9.6377 9.8492 9.7872 9.7326 9.6870 9.6764 9.6578 9.6235 9.6444angle =0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 0.4266 0.4209 0.4498 0.4823 0.4452PL = 0.0488参考文献1 何仰赞 温增银 电力系统分析 华中科技大学出版社.2 李维波. MATLAB 在电气工程中应用 中国电力出版社.2007