基于IEEE33节点前推回代 matlab

%*********************%潮流计算*****************************************
function flow（）
global X Y Z  Ploss Num V DG
format long
Branch_Data=dlmread（‘branchdata.txt‘）;   %读取线路参数矩阵
Node_Data=dlmread（‘nodedata.txt‘）;       %读取节点参数矩阵
UB = 12.66;                             % 电压基准 kV
SB = 10;                                % 功率基准 MVA
ZB = UB^2/SB;                           % 阻抗基准 ohm
Branch_Data（:[34]） = Branch_Data（:[34]） / ZB;     % 阻抗标幺值
Node_Data（:[23]） = Node_Data（:[23]） / SB / 1000;  % 功率标幺值
for j=1:DG
Node_Data（X（j）-12） = Node_Data（X（j）-12） - Y（j）;        %节点编号从1开始
Node_Data（X（j）-13） = Node_Data（X（j）-13） - Z（j）;
end
A = zeros（Num）;
for n = 1:Num-1
    A（Branch_Data（n1）Branch_Data（n2）） = 1;
end                                     % 形成邻接矩阵A
AT = A‘;                                        % 形成A的转置矩阵
S = [0;Node_Data（:2） + i*Node_Data（:3）];           % 形成节点负荷矩阵，根节点赋值0
ZL = [0;Branch_Data（:3） + i*Branch_Data（:4）];       % 形成线路阻抗矩阵
 
V = ones（Num1）;                         % 各个节点电压赋初值
Delta = 1;                              % 收敛判据赋初值
TempV = V;