节点导纳矩阵的Matlab算法

function Y = Admittance_matrix（nbAB）
%形成导纳矩阵（节点数，支路数，支路参数矩阵A，电容对地导纳矩阵B）
%函数说明
%{
    1.支路信息不可重复输入，即已经输入i-j支路的信息，就不可输入j-i支路的信息
    2.含变压器的支路第一个节点编号默认为变压器一次侧（即变压器等值电路中的阻抗归算侧）
    3.变比为1表示线路；变比不为1表示变压器
   *4.考虑到课本例题4-1，节点可能接有电容，输入矩阵B作为节点电容对地导纳参数。如果没有可以置为：
       B=[1   0.00;
          2   0.00;               
          3   0.00;
          ... ... ];
%}

    Y=zeros（n）;%生成导纳矩阵
    y1=zeros（n）;%串联支路导纳矩阵
    y2=zeros（n）;%节点处并联支路导纳矩阵

    %考虑到节点可能接有电容（课本例题4-1），计算节点电容对地导纳参数
    for p=1:n
        if B（p2）~=0;
            t=B（p1）;
            y2（tt）=B（p2）*i;
        end
    end
    %计算串联支路导纳
    for p=1:n
        if A（p6）~=1    %变比不为1时
            y1（A（p1）A（p2））=（1/（A（p4）*i））/A（p6）;
            y1（A（p2）A（p1））=y1（A（p