直流潮流程序，有用的

function [load_shed EXIT] = load_shedding1（canshu）

n=10;%抽样次数
n0=30;%系统节点数
n1=6;%发电机台数
n2=41;%线路数
generatorgzl=0.04;
% linegzl=[0.5 1.54 1.45 0.32 1.65 1.47 0.35 0.97 0.68 0.35 1.73 0.92 1.17 2.13 1.09 1.66 1.66 1.61 1.82 1.08 0.57 1.74 0.70 0.62 1.25 0.2 1.68 1.49 2.25 2.74 3.17 1.74 3.46 5.02 3.78
%         1.67 0.5 1.73 4.63 2.13 3.3];
R=rand（n1+n2n）;%生成随机数
lie=[0.004;0.004;0.004;0.004;0.004;0.004;0.004;0.002;0.002;0.002;0.002;0.002;0.002;0.002;0.002;0.002;0.002;0.002;0.002;0.02;0.002;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;0.02;
        0.002;0.002;0.002;0.002;0.002;0.002];
[lengthwide]=size（R）;
for i=1:wide
    R（find（R（:i）<=lie（:1））i）=0;
    R（find（R（:i）>lie（:1））i）=1;
   
    
end
T=（unique（R‘‘rows‘））;%取出相同列
%%%判断是否有孤岛产生%%%%%

M=T‘

[lengthwide]=size（M）;
for i=1:wide
    mpc=case30;
    for j=1:n1
        if M（ji）==0
            mpc.gen（j8）=0;
        end
    end
        for j=n1+1:n1+n2
            if M（ji）==0
               mpc.branch（j-n111）=0;
            end
        end
        canshu{i}=mpc%保存系统发电机、线路参数
        jieguo{i}=rundcpf（mpc）%对每次系统状态的潮流结果进行保存
end
%   canshu{1}=mpc;%保存系统发电机、线路参数
%         jieguo{1}=rundcpf（mpc）;%对每次系统状态的潮流结果进行保存

%   canshu{1}=mpc;%保存系统发电机、线路参数
%         jieguo{1}=rundcpf（mpc）;%对每次系统状态的潮流结果进行保存

%该函数根据输入的系统数据mpc，计算得到各节点切负荷情况
%mpc：系统数据，其中包括线路、发电机状态
%load_shed：两列的矩阵，第一列为节点编号，第二列为对应切负荷量，单位MW
%EXIT：优化算法收敛标志，如果为1则表明找到最优解


%形成节点、线路、发电机参数，定义基准功率
mpc=canshu{i};
bus = mpc.bus;
gen = mpc.gen;
branch = mpc.branch;
baseMVA = mpc.baseMVA;
%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%修改线路，去掉停运线路
index = find（branch（:11） == 0）;
%%%%%%%%%%%%%%%%%%%%%%%%%
%形成B0矩阵
B0 = makeBdc（mpc.bas