13节点配电网流计算—上机.docVIP

配网前推后代潮流计算, 对给定的13节点配电网，熟悉MATLAB语言前腿后代潮流计算程序。 13节点配电网结构如图1所示，支路参数见表1；节点类型及迭代初始电压及各节点所接负荷见表2。 要看懂前推回代法计算程序，报告叙述计算原理及计算流程。绘制计算流程框图。确定前推回代支路次序（广度优先，或深度优先），编写前推回代计算输入文件。进行潮流计算。 整理专利计算结果。 内容包括节点电压；线路潮流分布，网损 总结配电网潮流计算的要点 图1 13节点配电网结构图 表1 系统支路参数 支路 R（Ω.）Ω.） ））） 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 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]=case113() PQ=[ %节点电压 有功 无功 10.5 0 0 10. 0.6 0.45 10 0.4 0.3 10 0.4 0.28 10 0.6 0.4 10 0.4 0.3 10 0.5 0.35 10 0.5 0.4 ]; FT=[ %首端 末端 4 3 6 5 7 5 3 2 5 2 8 2 2 1 ]; RX=[ % R X 0.6 0.35 1.0 0.55 0.65 0.35 0.62 0.5 0.72 0.75 0.90 0.5 0.54 0.65 ]; 计算结果清单： k =7 PQ1 = 0.4014 0.2808 0.4025 0.3014 0.5025 0.3513 0.8074 0.5857 1.5294 1.0782 0.5036 0.4020 3.5356 2.6303 V =10.5000 10.1557 10.0776 10.0439 9.9677 9.9107 9.9226 10.0913 angle = 0 0.4716 0.4944 0.4785