基于Matpower编写的静态电压稳定评估仿真程序-IEEE39系统仿真，连续潮流法

基于Matpower编写了静态电压稳定评估仿真程序，用于负荷极限求解，静态电压稳定裕度求解，仿真算例为IEEE39系统
程序如下,安装Matpower后，可直接在matlab中运行：
define_constants;
mpopt = mpoption(‘out.all’, 0, ‘verbose’, 2);
mpopt = mpoption(mpopt, ‘cpf.stop_at’, ‘NOSE’, ‘cpf.step’, 0.2);
mpopt = mpoption(mpopt, ‘cpf.plot.level’, 2);
mpcb = loadcase(‘case39’); % load base case
allsum=[];
for i=1:15000
mpcb1=mpcb;
plall=0;
qlall=0;
pgall=0;
qgall=0;
Ptotal=0;
%flag=round(rand(1,1)*2);
%
k1p=(rand(1,1)*2-1)*0.3;
k1q=(rand(1,1)*2-1)*0.3;

k3p=(rand(1,1)*2-1)*0.3;
k3q=(rand(1,1)*2-1)*0.3;

k4p=(rand(1,1)*2-1)*0.3;
k4q=(rand(1,1)*2-1)*0.3;

k7p=(rand(1,1)*2-1)*0.3;
k7q=(rand(1,1)*2-1)*0.3;

k8p=(rand(1,1)*2-1)*0.3;
k8q=(rand(1,1)*2-1)*0.3;

k9p=(rand(1,1)*2-1)*0.3;
k9q=(rand(1,1)*2-1)*0.3;

k12p=(rand(1,1)*2-1)*0.3;
k12q=(rand(1,1)*2-1)*0.3;

k15p=(rand(1,1)*2-1)*0.3;
k15q=(rand(1,1)*2-1)*0.3;

k16p=(rand(1,1)*2-1)*0.3;
k16q=(rand(1,1)*2-1)*0.3;

k18p=(rand(1,1)*2-1)*0.3;
k18q=(rand(1,1)*2-1)*0.3;

k20p=(rand(1,1)*2-1)*0.3;
k20q=(rand(1,1)*2-1)*0.3;

k21p=(rand(1,1)*2-1)*0.3;
k21q=(rand(1,1)*2-1)*0.3;

k23p=(rand(1,1)*2-1)*0.3;
k23q=(rand(1,1)*2-1)*0.3;

k24p=(rand(1,1)*2-1)*0.3;
k24q=(rand(1,1)*2-1)*0.3;

k25p=(rand(1,1)*2-1)*0.3;
k25q=(rand(1,1)*2-1)*0.3;

k26p=(rand(1,1)*2-1)*0.3;
k26q=(rand(1,1)*2-1)*0.3;

k27p=(rand(1,1)*2-1)*0.3;
k27q=(rand(1,1)*2-1)*0.3;

k28p=(rand(1,1)*2-1)*0.3;
k28q=(rand(1,1)*2-1)*0.3;

k29p=(rand(1,1)*2-1)*0.3;
k29q=(rand(1,1)*2-1)*0.3;

k31p=(rand(1,1)*2-1)*0.3;
k31q=(rand(1,1)*2-1)*0.3;

k39p=(rand(1,1)*2-1)*0.3;
k39q=(rand(1,1)*2-1)*0.3;

g30p=(rand(1,1)*2-1)*0.3;