电力系统中的Kron简化，含MATLAB代码（全网唯一）

1、关于Kron降阶的叙述

电力系统中的Kron简化（Kron Reduction）_shizk_cn的博客-CSDN博客_kron简化电力系统中的Kron简化简单介绍一下电力系统导纳矩阵简化方法Kron简化（Kron Reduction），如有不对的地方还请大家指正。本文参考Power Systems Analysis, Second Edition by Arthur R. Bergen and Vijay Vittal. 机械工业出版社有影印版。...https://blog.csdn.net/weixin_42217751/article/details/1069872212、直接上代码，大家自行复制即可

%% 组织：东北电力大学
%% 作者：Yang Zheng
%% 日期：2022-11-6
function Y_new=kron(Y,Reduction_index)
%Y为需要降阶的矩阵
%Reduction_index为降阶节点
    [raw,column]=size(Y);
    Y_new=zeros(raw,column);
    Reduction_mat=[];
    for n = 1:max(size(Reduction_index))
        for i = 1:raw
            if(ismember(i,Reduction_mat))
                continue;
            end
            for j = 1:column
                if(ismember(j,Reduction_mat))
                    continue;
                end
                Y_new(i,j)=Y(i,j)-Y(i,Reduction_index(n))*Y(Reduction_index(n),j)/Y(Reduction_index(n),Reduction_index(n));
            end
        end
        Reduction_mat=[Reduction_mat Reduction_index(n)];
        Y=Y_new;
    end
    Y_new=Y;
end
 
%test
% Y=[
%     5+6i, -1-2j,0,   -4-4i;
%     -1-2i,6+9i ,-3-4i,-2-3i;
%     0,    -3-4i,7+8i, -4-4i;
%     -4-4i,-2-3i,-4-4i,10+11i;];
% Reduction_index=[2 4];