配电网前推回代潮流算法

clc;
clear;
%配电网编号从零开始，电源点为0，电源点为平衡点
z=input（‘请输入配电网参数z=‘）;    %第一列为支路首节点号，第二列为支路末节点号，第三列为支路参数，第四列为末节点负荷第五列为末节点类型（1为PQ节点，2为PV节点）
n=input（‘请输入配电网节点数n=‘）;
e=input（‘请输入收敛精度e=‘）;
Sb=1;
Vb=10.5;
Zb=Vb^2/Sb;
Z（:3）=Z（:3）/Zb;
Z（:4）=Z（:4）/Sb/1000;           %修改配电网参数，使之统一至标幺值
V=ones（n-11）;                   %为各节点赋电压初值（电源点除外）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%开始计算
for i=1:20                       %i为迭代次数，最大迭代次数为20，超过20则收敛性较差
    I=（z（:4）.‘）*（diag（1./V））;   %求各节点注入电流
    I=I‘;
    I1=I;
   for k=n-1:-1:1                %前推支路电流及功率分布
      A=（find（z（:1）==k））;       %判断矩阵Z第1列是否存在值为k的元素，并把行值赋给矩阵A，亦即找出网络的末节点
      if size（A）==0              %k为末节点
         I（k1）;
      else 
         I（k1）=sum（I（A1））+I（k1）;%利用KCL求支路电流