ieee30节点

clear%清除内存变量
format long
jc=0;
%--------------------------节点原始数据模块--------------------------------%
JieDianShuJu=[1 1.0500 0.0000 0.0000 0.0000 0.0000 0.0000;
              2 1.0338 0.0000 0.5756 0.0000 0.2170 0.1270;
              3 1.0000 0.0000 0.0000 0.0000 0.0240 0.0120;
              4 1.0000 0.0000 0.0000 0.0000 0.0760 0.0160;
              5 1.0058 0.0000 0.2456 0.0000 0.9420 0.1900;
              6 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              7 1.0000 0.0000 0.0000 0.0000 0.2280 0.1090;
              8 1.0230 0.0000 0.3500 0.0000 0.3000 0.3000;
              9 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              10 1.0000 0.0000 0.0000 0.0000 0.0580 0.0200;
              11 1.0913 0.0000 0.1793 0.0000 0.0000 0.0000;
              12 1.0000 0.0000 0.0000 0.0000 0.1120 0.0750;
              13 1.0883 0.0000 0.1691 0.0000 0.0000 0.0000;
              14 1.0000 0.0000 0.0000 0.0000 0.0620 0.0160;
              15 1.0000 0.0000 0.0000 0.0000 0.0820 0.0250;
              16 1.0000 0.0000 0.0000 0.0000 0.0350 0.0180;
              17 1.0000 0.0000 0.0000 0.0000 0.0900 0.0580;
              18 1.0000 0.0000 0.0000 0.0000 0.0320 0.0090;
              19 1.0000 0.0000 0.0000 0.0000 0.0950 0.0340;
              20 1.0000 0.0000 0.0000 0.0000 0.0220 0.0070;
              21 1.0000 0.0000 0.0000 0.0000 0.1750 0.1120;
              22 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              23 1.0000 0.0000 0.0000 0.0000 0.0320 0.0160;
              24 1.0000 0.0000 0.0000 0.0000 0.0870 0.0670;
              25 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              26 1.0000 0.0000 0.0000 0.0000 0.0350 0.0230;
              27 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              28 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              29 1.0000 0.0000 0.0000 0.0000 0.0240 0.0090;
              30 1.0000 0.0000 0.0000 0.0000 0.1060 0.0190];%输入节点数据
JieDianZongShu=size（JieDianShuJu1）;%节点总数
PVJieDianShu=5;%PV节点数
PQJieDianShu=JieDianZongShu-PVJieDianShu-1;%PQ节点数
PVI=[0 1 0 0 1 0 0 1 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];%标示某节点是否为PV节点的矩阵
PQI=[0 0 1 1 0 1 1 0 1 1 0 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1];%标示某节点是否为PQ节点的矩阵
PVR=[1 2 5 8 11 13];%松弛节点与PV节点编号
PV=[2 5 8 11 13];%PV节点编号
PQ=[3 4 6 7 9 10 12 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30];%PQ节点编号
SC=[1 4 7 10 12];%NLM中需删除的行与列
U=JieDianShuJu（:2）;%电压有效值
Theta=JieDianShuJu（:3）;%电压相位
PG=JieDianShuJu（:4）;%发电机有功
QG=JieDianShuJu（:5）;%发电机无功
PL=JieDianShuJu（:6）;%负载有功
QL=JieDianShuJu（:7）;%负载无功
DetP=zeros（JieDianZongShu1）;%节点有功功率误差向量
DetQ=zeros（JieDianZongShu1）;%节点无功功率误差向量
sump=0;%临时变量
sumq=0;%临时变量
DetTheta=0;%临时变量
%==========================节点原始数据模块================================%
%
%--------------------------定义雅可比阵和修正向量---------------------------%
XiuZhengXiangLiang=zeros（JieDianZongShu+PQJieDianShu-11）;%定义修正向量
ThetaXiuZheng=XiuZhengXiangLiang（1:JieDianZongShu-1）;%定义幅角修正向量
UXiuZheng=XiuZhengXiangLiang（JieD