完整版前推回代法计算流程.doc

前推回代法计算流程 要看懂前推回代法计算程序， 报告叙述计算原理及计算流程。 绘制计算流程框图。 确定前推回代 支路次序（广度优先，或深度优先） ，编写前推回代计算输入文件。进行潮流计算。 下列为节点配电网结构图及系统支路参数和系统负荷参数表。 1 7 6 2 8 9 10 3 11 12 4 5  13 1-2 节点配电网结构图 1 系统支路参数 支路 R(Ω) X （ Ω.） B/2(T K ) （ S） 1~2 3.367 3.685 0.0 2~3 2.356 2.541 0.0 3~4 1.145 1.28 0.0 4~5 4.524 5.04 0.0 2~6 0.856 1.14 0.0 6~7 2.745 2.965 0.0 2~8 3.743 4.251 0.0 8~9 2.237 2.756 0.0 3~10 4.14 4.696 0.0 3~11 1.328 1.763 0.0 11~12 2.436 2.866 0.0 4~13 3.521 3.966 0.0 表 2 系统负荷参数 节点编号 节点类型 节 点 初 始 电 压 Pi（ MV A ） Qi（ MV A ） kV ） 1 1 根节点 10.4 0 0 2 中间节点 10.0 0.0342 0.0301 3 中间节点 10.0 0.0693 0.0642 4 中间节点 10.0 0.0845 0.0763 5 叶节点 10.0 0.0295 0.0261 6 中间节点 10.0 0.0474 0.0409 7 叶节点 10.0 0.1176 0.0957 8 中间节点 10.0 0.0946 0.0857 9 叶节点 10.0 0.0916 0.0859 10 叶节点 10.0 0.0271 0.0229 11 中间节点 10.0 0.0696 0.0643 12 叶节点 10.0 0.0676 0.0579 13 叶节点 10.0 0.0298 0.0242 主程序清单： [PQ,FT,RX]=case114(); %调用数据文件 NN=size(PQ,1); %节点数 NB=size(FT,1); %支路数数 V=PQ(:,1); %V 初始电压相量 maxd=1 k=1 while maxd0.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 2 PL=PL+RX(i,1)*x; end angle(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)); %计算支路 end maxd=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 end maxd k=k+1 PQ1 %潮流分布 即支路首端潮流 MVA V=V2 % 节点电压模计算结果 kV angle %节点电压角度计算结果单位度 PL %网损单位 MW end clear 输入文件清单： function [PQ,FT,RX]=case114() PQ=[ %节点电压 有功 无功 10.4 0 0 3 10.0 0.0342 0.0301 10.0 0.0693 0.0642 10.0 0.0845 0.0763 10.0 0.0295 0.0261