matlab内点法

%数据加载 n=input('请输入要计算的节点系统(5):') load Node5.txt;%节点数据 load Branch5.txt;%支路数据 load Generator5.txt;%发电机数据 Node=Node5; Branch=Branch5; Generator=Generator5;

clear;
%clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%初始化%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%数据加载
n=input（‘请输入要计算的节点系统（5）:‘）
    load Node5.txt;%节点数据
    load Branch5.txt;%支路数据
    load Generator5.txt;%发电机数据
    Node=Node5;
    Branch=Branch5;
    Generator=Generator5;

%节点数据处理
N=Node（:1）;%节点号
Type=Node（:2）;%节点类型
Uamp=Node（:3）;%节点电压幅值
Dlta=Node（:4）;%节点电压相角
Pd=Node（:5）;%节点负荷有功
Qd=Node（:6）;%节点负荷无功
Pg=Node（:7）;%节点出力有功
Qg=Node（:8）;%节点出力无功
Umax=Node（:9）;%节点电压幅值上限
Umin=Node（:10）;%节点电压幅值下限
Bc=Node（:11）;%节点补偿电容电纳值
%支路数据处理
Nbr=Branch（:1）;%支路号
Nl=Branch（:2）;%支路首节点
Nr=Branch（:3）;%支路末节点
R=Branch（:4）;%支路电阻
X=Branch（:5）;%支路电抗
Z=R+1i*X;%支路阻抗=支路电阻+支路电抗
Bn=Branch（:6）;%支路对地电纳
K=Branch（:7）;%支路变压器变比，0表示无变压器
Ptmax=Branch（:8）;%线路传输功率上限

%发电机数据处理
Ng=Generator（: