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

配网前推后代潮流计算, 对给定的13节点配电网，熟悉MATLAB语言前腿后代潮流计算程序。 13节点配电网结构如图1所示，支路参数见表1；节点类型及迭代初始电压及各节点所接负荷见表2。 要看懂前推回代法计算程序，报告叙述计算原理及计算流程。绘制计算流程框图。确定前推回代支路次序（广度优先，或深度优先），编写前推回代计算输入文件。进行潮流计算。 整理专利计算结果。 内容包括节点电压；线路潮流分布，网损 总结配电网潮流计算的要点 图1 13节点配电网结构图表1 系统支路参数 支路 R（Ω.）Ω.） ））） maxd=1 k=1 while maxd0.0001PQ2=PQ;%每一次迭代各节点的注入有功和无功相同PL=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);%用PQ1去修正支路末端节点的有功Q 单位MvarPL=PL+RX(i,1)*x; end angle(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+dv2^2);%计算支路末端电压/kVangle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1));%计算支路 end maxd=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));endend 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=[ % RX 0.60.35 1.00.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.28080.4025 0.30140.5025 0.35130.8074 0.58571.5294 1.07820.5036 0.40203.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 0.6815 0.6351 0.6511 0.4101 PL =0.1356MW