内点法解最优潮流

clear; %clc
%数据加载
c=input（‘请输入要计算的节点系统（6）:‘）
  load Node6.txt;%节点数据
  load Branch6.txt;%支路数据
  load Generator6.txt;%发电机数据
  Node=Node6;
  Branch=Branch6;
  Generator=Generator6;

  
%节点数据处理
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（:1）;%发电机序号
Nbus=Generator（:2）;%所在母线号
Pumax=Generator（:3）;%发电机有功出力上界
Qumax=Generato