clear%清除內(nèi)存變量
format?long
jc=0;
%--------------------------節(jié)點(diǎn)原始數(shù)據(jù)模塊--------------------------------%
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];%輸入節(jié)點(diǎn)數(shù)據(jù)
JieDianZongShu=size（JieDianShuJu1）;%節(jié)點(diǎn)總數(shù)
PVJieDianShu=5;%PV節(jié)點(diǎn)數(shù)
PQJieDianShu=JieDianZongShu-PVJieDianShu-1;%PQ節(jié)點(diǎn)數(shù)
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];%標(biāo)示某節(jié)點(diǎn)是否為PV節(jié)點(diǎn)的矩陣
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];%標(biāo)示某節(jié)點(diǎn)是否為PQ節(jié)點(diǎn)的矩陣
PVR=[1?2?5?8?11?13];%松弛節(jié)點(diǎn)與PV節(jié)點(diǎn)編號
PV=[2?5?8?11?13];%PV節(jié)點(diǎn)編號
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節(jié)點(diǎn)編號
SC=[1?4?7?10?12];%NLM中需刪除的行與列
U=JieDianShuJu（:2）;%電壓有效值
Theta=JieDianShuJu（:3）;%電壓相位
PG=JieDianShuJu（:4）;%發(fā)電機(jī)有功
QG=JieDianShuJu（:5）;%發(fā)電機(jī)無功
PL=JieDianShuJu（:6）;%負(fù)載有功
QL=JieDianShuJu（:7）;%負(fù)載無功
DetP=zeros（JieDianZongShu1）;%節(jié)點(diǎn)有功功率誤差向量
DetQ=zeros（JieDianZongShu1）;%節(jié)點(diǎn)無功功率誤差向量
sump=0;%臨時變量
sumq=0;%臨時變量
DetTheta=0;%臨時變量
%==========================節(jié)點(diǎn)原始數(shù)據(jù)模塊================================%
%
%--------------------------定義雅可比陣和修正向量---------------------------%
XiuZhengXiangLiang=zeros（JieDianZongShu+PQJieDianShu-11）;%定義修正向量
ThetaXiuZheng=XiuZhengXiangLiang（1:JieDianZongShu-1）;%定義幅角修正向量
UXiuZheng=XiuZhengXiangLiang（JieD