redj0.m

可以计算电力系统概率潮流

function J0red=redJ0(J0, nswing)

% remove nswing row and column
glo_var

J0red=zeros(nbus-1, nbus-1);

J0red(1:nswing-1, 1:nswing-1)=J0(1:nswing-1, 1:nswing-1);
J0red(nswing:nbus-1, 1:nswing-1)=J0(nswing+1:nbus, 1:nswing-1);
J0red(1:nswing-1, nswing:nbus-1)=J0(1:nswing-1, nswing+1:nbus);
J0red(nswing:nbus-1, nswing:nbus-1)=J0(nswing+1:nbus, nswing+1:nbus);


return