任务：

（1）能读懂并理解前推回代法的Matlab程序；

（2）根据所给的辐射式配电网参数，要求能进行程序的修改并执行，得出结果，和分析。

8节点7支路辐射式配电网潮流计算算例的Matlab源代码【部分语句需同学们自行补齐】

%%%%%%%%%%

%主程序清单：

[PQ,FT,RX]=data();      %调用数据文件，data()为一子函数，功能是原始数据

NN=size(PQ,1);      %节点数

NB=size(FT,1);      %支路数数

V=PQ(:,1);          %V初始电压相量

maxd=1

k=1

while maxd>0.0001

    PQ2=PQ;     %每一次迭代各节点的注入有功和无功相同

    PL=0.0;

for i=1:NB

    kf=FT(i,1);     %前推始节点号

    kt=FT(i,2);     %前推终节点号

    x=***;   %计算沿线电流平方A

    PQ1(i,1)=PQ2(kf,2)+****;        %计算支路首端有功/MW  RX(i,1)~R

    PQ1(i,2)=PQ2(kf,3)+*****;        %计算沿支路的无功损耗/Mvar  RX(i,2)~X

    PQ2(kt,2)= PQ2(kt,2)+***;       %用PQ1去修正支路末端节点的有功P 单位MW

    PQ2(kt,3)= PQ2(kt,3)+***;       %用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=******;    %计算支路电压损耗的纵分量dv1

    dv2=******;    %计算支路电压损耗的横分量dv2

    V2(kt)=****;                 %计算支路末端电压/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

%输入文件清单,data()子函数如下：

function [PQ,FT,RX]=data()

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

];