redy.m

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

function Yred=redY(Y, nswing)

% remove nswing row and column
glo_var

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

Yred(1:nswing-1, 1:nswing-1)=Y(1:nswing-1, 1:nswing-1);
Yred(nswing:nbus-1, 1:nswing-1)=Y(nswing+1:nbus, 1:nswing-1);
Yred(1:nswing-1, nswing:nbus-1)=Y(1:nswing-1, nswing+1:nbus);
Yred(nswing:nbus-1, nswing:nbus-1)=Y(nswing+1:nbus, nswing+1:nbus);


return