function [Y_B]=YmatrixConstruct(para_line,para_trans2,para_trans3,para_capac,SB,NodeNum)
% 自動生成節(jié)點(diǎn)導(dǎo)納矩陣程序
% SB為功率基準(zhǔn)值，NodeNum是電力系統(tǒng)網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)(含三繞組變壓器增加的節(jié)點(diǎn)數(shù))

Y_B=zeros(NodeNum,NodeNum);

% 對線路
for i=1:size(para_line,1)
    [Y,Z]=LineModel(para_line(i,:),SB);
    act_node=para_line(i,2:3);
    Y_B(act_node(1),act_node(1))=Y_B(act_node(1),act_node(1))+1/Z+Y;
    Y_B(act_node(2),act_node(2))=Y_B(act_node(2),act_node(2))+1/Z+Y;
    Y_B(act_node(1),act_node(2))=Y_B(act_node(1),act_node(2))-1/Z;
end

% 對兩繞組變壓器
for i=1:size(para_trans2,1)
    [Y,Z,Gt,Bt]=Trans2Model(para_trans2(i,:),SB);
    act_node=para_trans2(i,2:3);
    Y_B(act_node(1),act_node(1))=Y_B(act_node(1),act_node(1))+1/Z+Y(1);
    Y_B(act_node(2),act_node(2))=Y_B(act_node(2),act_node(2))+1/Z+Y(2);
    Y_B(act_node(1),act_node(2))=Y_B(act_node(1),act_node(2))-1/Z;
end

% 對三繞組變壓器
for i=1:size(para_trans3,1)
    [Y_out,Z_out]=Trans3Model(para_trans3(i,:),SB);
    Y_B(para_trans3(5),para_trans3(5))=Y_B(para_trans3(5),para_trans3(5))+sum(1./Z_out)+Y_out(4);
    act_node=para_trans3(i,2:4);
    M_sequence=1;
    for j=act_node
        Y_B(j,j)=Y_B(j,j)+1/Z_out(M_sequence)+Y_out(M_sequence);
        Y_B(j,para_trans3(5))=Y_B(j,para_trans3(5))-1/Z_out(M_sequence);
        M_sequence=M_sequence+1;
    end
end

% construct Y_B  
temp_triu=triu(Y_B,1)
Y_B=Y_B+(triu(Y_B,1)).';
Y_B=Y_B+(tril(Y_B,-1)).';
Y_B=Y_B-temp_triu;