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;
???????????