matlab求节点导纳矩阵

节点导纳矩阵可以通过电路的基本元件参数计算得到，下面是一个示例代码：

假设电路有n个节点，m个支路，邻接矩阵为G，其中G(i,j)表示第i个节点和第j个节点之间的支路的电导。则节点导纳矩阵为Y，其中Y(i,i)为第i个节点的总导纳，Y(i,j)为第i个节点和第j个节点之间的互导纳，有以下代码实现：

n = 3;  % 节点数
m = 4;  % 支路数

% 构造邻接矩阵G，假设电路为以下情况：
%    1---R1---2
%     \      |
%      C     R2
%     /      |
%    3---L---4
G = zeros(n, n);
G(1, 2) = 1/R1;
G(1, 3) = 1/(1i*w*C);
G(2, 1) = 1/R1;
G(2, 4) = 1/R2;
G(3, 1) = 1/(1i*w*C);
G(3, 4) = -1/(1i*w*L);
G(4, 2) = 1/R2;
G(4, 3) = -1/(1i*w*L);

% 构造节点导纳矩阵Y
Y = zeros(n, n);
for i = 1:n
    for j = 1:n
        if i == j
            % 计算节点i的总导纳
            Y(i, i) = sum(G(i, :));
        else
            % 计算节点i和节点j之间的互导纳
            Y(i, j) = -G(i, j);
        end
    end
end

其中，R1、R2、C和L为电路中的元件参数，w为角频率。