[电路的数值仿真--节点列表法]

最新推荐文章于 2024-09-29 23:16:27 发布

Mai_Zuning

最新推荐文章于 2024-09-29 23:16:27 发布

阅读量164

点赞数

文章标签： matlab

本文链接：https://blog.csdn.net/mai_zuning/article/details/132200506

版权

节点列表法

WHY

回路电流法不允许存在无伴电流源支路,且规定的复合支路不允许存在受控电流源;结点电压法和割集电压法不允许存在无伴电压源支路,且规定的复合支路不允许存在受控电压源。这就使上述几种电路分析法有一定的局限性。列表法对支路类型无过多限制,适应性强,但方程数较多。

列表法的矩阵形式

KCL、KVL、支路方程
$A I = 0$ ; $U-A^{T}U_{n}=0$ ; $FU+HI=U_{s}+I_{s}$
$\left[ \begin{matrix} 0& 0& A\\ -A^T& I_b& 0\\ 0& F& H\\ \end{matrix} \right] \left[ \begin{array}{c} U_n\\ U\\ I\\ \end{array} \right] =\left[ \begin{array}{c} 0\\ 0\\ U_S+I_S\\ \end{array} \right]$

电路图的关联矩阵

设一条支路连接于某两个结点,则称该支路与这两个结点相关联。支路与结点的关联性质可以用所谓关联矩阵描述。设有向图的结点数为n ,
支路数为b,且所有结点与支路均加以编号。于是,该有向图的关联矩阵为一个(n×b)的矩阵,用A_a表示。它的行对应结点,列对应支路,它的
任一元素 $a_{jk}$ 定义如下:
$a_{jk}=1$ ,表示支路k与结点j关联并且它的方向背离结点(支路方向流出节点）
$a_{jk}=-1$ ,表示支路k与结点j关联并且它指向结点;(支路方向流入节点）
$a_{jk}=0$ ,表示支路é与结点j无关联。
$A_{a}$ 将j=0的行删除，构成独立关联矩阵A.

节点列表法，表格各列的说明。
当列2或列3出现数值0,说明是参考节点。max{第二列最大值，第三列最大值}=独立总节点数；总节点数=独立总节点数+1；
支路方向：表格每一行的列2指向列3.
元件类型：1->电阻R;2->电感L;3->电容C;4->电压源Us;5->电流源Is;6->VCVS;7->VCCS;8->CCVS;9->CCCS

支路号	始结点号	终结点号	元件类型	元件值	受控源控制的支路号

有向图可以用表格来表示。那么程序求解关联矩阵就水到渠成。

function A = guanlianju(table)
[m,n] = size (table);
node_num = max (max (table (:,2:3)))+1;
edge_num = m;
A = zeros (node_num, edge_num);
for i = 1:edge_num
    start = table (i,2);
    endl = table (i,3);
    if start == 0
        start=node_num;
    end
    if endl==0
        endl=node_num;
    end
    A (start,i) = 1;
    A (endl,i) = -1;
end
A(end,:)=[];
end