写出一个matlab潮流计算的程序

下面是一个简单的MATLAB潮流计算程序，用于计算一个4节点、3支支路的电力系统中的节点电压和线路电流。程序主要步骤包括输入系统数据、构建总线节点入射矩阵和线路阻抗矩阵、计算负荷功率、计算节点电压和线路电流等。

% 输入系统数据
n = 4; % 节点数
l = 3; % 支路数
Pload = [0; 0; 0; 1]; % 负荷功率
Zline = [0.1+0.2i, 0.05+0.1i, 0; 0.05+0.1i, 0, 0.2+0.4i; 0, 0.2+0.4i, 0.1+0.2i]; % 线路阻抗

% 构建总线节点入射矩阵
B = [1 0 0; -1 1 0; 0 -1 1; 0 0 -1];

% 计算节点电压和线路电流
V = ones(n, 1); % 初始猜测节点电压
tol = 1e-6; % 收敛精度
iter = 0; % 迭代次数
while true
    Iline = inv(Zline) * B * (V + Pload); % 计算线路电流
    Vnew = Zline * Iline + Pload; % 计算节点电压
    if max(abs(Vnew - V)) < tol % 判断是否收敛
        break;
    end
    V = Vnew;
    iter = iter + 1;
    if iter > 100 % 判断是否达到最大迭代次数
        warning('潮流计算未收敛！');
        break;
    end
end

% 输出结果
disp(['节点电压：']);
disp(V);
disp(['线路电流：']);
disp(Iline);

其中，tol为收敛精度，iter为迭代次数，当节点电压的改变量小于tol时，认为计算收敛；若迭代次数超过100次仍未收敛，则输出警告信息。