matlab 前推回代,前推回代法计算流程.docx

前推回代法计算流程

前推回代法计算流程要看懂前推回代法计算程序，报告叙述计算原理及计算流程。绘制计算流程框图。确定前推回代支路次序(广度优先，或深度优先)，编写前推回代计算输入文件。进行潮流计算。下列为节点配电网结构图及系统支路参数和系统负荷参数表。图1-2 节点配电网结构图表1 系统支路参数支路R(Ω)X(Ω.)B/2(TK) (S)1~23.3673.6850.02~32.3562.5410.03~41.1451.280.04~54.5245.040.02~60.8561.140.06~72.7452.9650.02~83.7434.2510.08~92.2372.7560.03~104.144.6960.03~111.3281.7630.011~122.4362.8660.04~133.5213.9660.0表2 系统负荷参数节点编号节点类型节点初始电压(kV)Pi(MVA)Qi(MVA)1根节点10.4002中间节点10.00.03420.03013中间节点10.00.06930.06424中间节点10.00.08450.07635叶节点10.00.02950.02616中间节点10.00.04740.04097叶节点10.00.11760.09578中间节点10.00.09460.08579叶节点10.00.09160.085910叶节点10.00.02710.022911中间节点10.00.06960.064312叶节点10.00.06760.057913叶节点10.00.02980.0242主程序清单：[PQ,FT,RX]=case114(); %调用数据文件NN=size(PQ,1); %节点数NB=size(FT,1); %支路数数V=PQ(:,1); %V初始电压相量maxd=1k=1while maxd>0.0001 PQ2=PQ; %每一次迭代各节点的注入有功和无功相同 PL=0.0;for i=1:NB kf=FT(i,1); %前推始节点号 kt=FT(i,2); %前推终节点号 x=(PQ2(kf,2)^2+PQ2(kf,3)^2)/V(kf)/V(kf); %计算沿线电流平方A PQ1(i,1)=PQ2(kf,2)+RX(i,1)*x; %计算支路首端有功/MW RX(i,1)~R PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x; %计算沿支路的无功损耗/Mvar RX(i,2)~X PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用PQ1去修正支路末端节点的有功P 单位MW PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用PQ1去修正支路末端节点的有功Q 单位Mvar PL=PL+RX(i,1)*x;endangle(1)=0.0;for i=NB:-1:1 kf=FT(i,2); %回代始节点号 kt=FT(i,1); %回代终节点号 dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2))/V(kf); %计算支路电压损耗的纵分量dv1 dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1))/V(kf); %计算支路电压损耗的横分量dv2 V2(kt)=sqrt((V(kf)-dv1)^2+dv2^2); %计算支路末端电压/kV angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1)); %计算支路endmaxd=abs(V2(2)-V(2));V2(1)=V(1);for i=3:1:NN if 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 010.0 0.0342 0.030110.0 0.0693 0.064210.0 0.0845 0.076310.0