14节点潮流 牛顿法

clear%清除内存变量
format long
jc=0;
%--------------------------节点原始数据模块--------------------------------%
% 1-节点编号 2-电压初值 3-电压相位初值  4-发电机有功 5-发电机无功 6-负载有功 7-负载无功
JieDianShuJu=[1 1.0600 0.0000 0.0000 0.0000 0.0000 0.0000;
              2 1.0450 0.0000 0.4000 0.4000 0.2170 0.1270;
              3 1.0100 0.0000 0.0000 0.2000 0.9420 0.1900;
              4 1.0000 0.0000 0.0000 0.0000 0.4780 -0.0390;
              5 1.0000 0.0000 0.0000 0.0000 0.0760 0.0160;
              6 1.0700 0.0000 0.0000 0.3000 0.1120 0.0750;
              7 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000;
              8 1.0900 0.0000 0.0000 0.3000 0.0000 0.0000;
              9 1.0000 0.0000 0.0000 0.0000 0.2950 0.1660;
              10 1.0000 0.0000 0.0000 0.0000 0.0900 0.0580;
              11 1.0000 0.0000 0.0000 0.0000 0.0350 0.0180;
              12 1.0000 0.0000 0.0000 0.0000 0.0610 0.0160;
              13 1.0000 0.0000 0.0000 0.0000 0.1350 0.0580;
              14 1.0000 0.0000 0.0000 0.0000 0.1490 0.0500];%输入节点数据
JieDianZongShu=size（JieDianShuJu1）;%节点总数
PVJieDianShu=4;%PV节点数
PQJieDianShu=JieDianZongShu-PVJieDianShu-1;%PQ节点数
PVI=[0 1 1 0 0 1 0 1 0 0 0 0 0 0];%标示某节点是否为PV节点的矩阵
PQI=[0 0 0 1 1 0 1 0 1 1 1 1 1 1];%标示某节点是否为PQ节点的矩阵
PVR=[1 2 3 6 8];%松弛节点与PV节点编号
PV=[2 3 6 8];%PV节点编号
PQ=[4 5 7 9 10 11 12 13 14];%PQ节点编号
SC=[1 2 5 7];%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（JieDianZongShu:end）;%定义电压修正向量
J=zeros（JieDianZongShu*2-2）;%定义雅可比矩阵
H=zeros（JieDianZongShu-1）;%初始化H阵
N=zeros（JieDianZongShu-1）;%初始化N阵
L=zeros（JieDianZongShu-1）;%初始化L阵
M=zeros（JieDianZongShu-1）;%初始化M阵
H1=H;N1=N;M1=M;L1=L;
%==========================定义雅可比阵和修正向量===========================%
%
%--------------------------支路原始数据模块--------------------------------%
% 1-序号 2-支路首节点 3-支路末节点 4-电阻 5-电抗 6-1/2对地电纳 7-变压器变比
ZhiLuShuJu=[-1 4 7 0.00000 0.20912 0.00000 0.978;
            -2 4 9 0.00000 0.55678 0.00000 0.969;
            -3 5 6 0.00000 0.25202 0.00000 0.932;
             0 0 9 0.00000 0.00000 0.01900 0.000
             1 1 2 0.01938 0.05917 0.02640 0.000;
             2 2 3 0.04699 0.19797 0.02190 0.000;
             3 2 4 0.05811 0.17632 0.01870 0.000;
             4 1 5 0.05403 0.22304 0.02460 0.000;
             5 2 5 0.05695 0.17388 0.01700 0.000;
             6 3 4 0.06701 0.17103 0.01730 0.000;
             7 4 5 0.01335 0.04211 0.00640 0.000;
             8 7 8 0.00000 0.17615 0.00000 0.000;