CN113326671A

CN113326671A - 一种对电路中线性电阻计算电路方程节点电流的方法

Info

Publication number: CN113326671A
Application number: CN202110661913.8A
Authority: CN
Inventors: 陶雄; 周振亚; 吴大可; 程明厚; 刘强
Original assignee: Shenzhen Huada Jiutian Technology Co ltd
Current assignee: Shenzhen Huada Jiutian Technology Co ltd
Priority date: 2021-06-15
Filing date: 2021-06-15
Publication date: 2021-08-31

Abstract

一种对电路中线性电阻计算电路方程节点电流的方法，包括以下步骤：1）根据线性电阻的端口位置建立第一电导矩阵；2）对第一电导矩阵进行去掉空行空列处理，获得第二电导矩阵；3）用电导值填充第二电导矩阵；4）利用第二电导矩阵计算电阻电流向量；5）计算总的电流向量。本发明利用电路中器件的连接关系不会变化同时线性电阻的电导值不变的特点，用线性电阻的电导值形成一个固定的电导矩阵，通过电导矩阵和电压向量的乘法得到线性电阻的电流向量，避免遍历每个线性电阻，使计算能够得到加速。

Description

一种对电路中线性电阻计算电路方程节点电流的方法

技术领域

本发明涉及模拟电路数值仿真技术领域，特别涉及一种对电路中线性电阻计算电路方程节点电流的方法。

背景技术

电路中的电流遵循基尔霍夫电流定律(KCL方程)，任意节点连接的所有支路的总电流和为零。在对电路进行SPICE精度的仿真过程中，需要计算出连到同一个节点上所有支路的电流和，并用作电路方程的右端项，通过解方程使电路状态向满足基尔霍夫电流定律的方向迭代。因此计算节点电流是整个仿真过程的重要部分。

一般计算节点电流的做法是，遍历所有器件，由器件根据电路状态计算出端口电流，添加到对应的节点上。对于大部分的电路，尤其是后仿真电路，线性的电阻电容占器件的大多数。对线性电阻，可以简单的用电导乘以电压差，得到端口电流。如果按照普通器件的处理方式，遍历每一个线性电阻，对电阻器件的访问会造成大量的内存随机访问。

因而，亟需一种能够高效计算线性电阻节点电流的方法。

发明内容

为了解决现有技术存在的不足，本发明的目的在于提供对电路中线性电阻计算电路方程节点电流的方法，以提高节点电荷的计算效率。

为实现上述目的，本发明提供的对电路中线性电阻计算电路方程节点电流的方法，包括以下步骤：

1)根据线性电阻的端口位置建立第一电导矩阵；

2)对第一电导矩阵进行去掉空行空列处理，获得第二电导矩阵；

3)用电导值填充第二电导矩阵；

4)利用第二电导矩阵计算电阻电流向量；

5)计算总的电流向量。

进一步地，所述步骤2)进一步包括以下步骤：

21)记录第一电导矩阵各元素的行号和列号；

22)去掉空行和空列，得到第二电导矩阵；

23)记录第一电导矩阵和第二电导矩阵之间各元素的行号和列号的映射关系。

进一步地，所述第一电导矩阵和所述第二电导矩阵为对称矩阵。

进一步地，所述步骤4)进一步包括以下步骤：

41)计算第二电导矩阵中各个节点对应的电压值，构成电压向量；

42)将第二电导矩阵与电压向量相乘，获得第二电导矩阵各个节点对应的第二电流向量。

进一步地，所述步骤5)进一步包括，将第二电流向量按照第一电导矩阵和第二电导矩阵之间的映射关系加至各个节点，再与其他器件的等效电流叠加，得到总的节点电流向量。

更进一步地，对电路结构连接关系不改变的电路中线性电阻计算电路方程节点电流时，步骤1)-步骤3)只需进行一次，步骤4)-步骤5)可以根据电压值改变重复多次。

为实现上述目的，本发明还提供一种对电路中线性电阻计算电路方程节点电流的装置，包括存储器和处理器，所述存储器上储存有在所述处理器上运行的程序，所述处理器运行所述程序时执行上述的对电路中线性电阻计算电路方程节点电流的方法的步骤。

为实现上述目的，本发明还提供一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行时执行上述的对电路中线性电阻计算电路方程节点电流的方法的步骤。

有益效果：本发明利用电路中器件的连接关系不会变化，同时线性电阻的电导值不变的特点，用线性电阻的电导值形成一个固定的矩阵。然后用电导矩阵和电压向量相乘算出线性电阻的端口电流，代替了遍历每个电阻单独计算电流的过程。剩余的器件，可以遍历器件来计算节点电流。将所有的电流合并，就可以得到方程的右端项。在做矩阵向量乘法时，依据线性电导矩阵的对称性，可以只对上三角或下三角矩阵做乘法并只存储对应的元素，提高计算效率。本发明通过矩阵和向量的乘法得到线性电阻的电流向量，避免遍历每个线性电阻，使计算能够得到加速。

本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。

附图说明

附图用来提供对本发明的进一步理解，并且构成说明书的一部分，并与本发明的实施例一起，用于解释本发明，并不构成对本发明的限制。在附图中：

图1为本发明的对电路中线性电阻计算电路方程节点电流的方法的流程图；

图2为根据本发明实施方式的示例电路图；

图3为一个线性电阻往电导矩阵中填入电导的实施方法示意图；

图4为根据图2中的电路得到的电导矩阵，其中保留了空行空列；

图5为整体矩阵中行号到电导矩阵的行号对应关系示意图，其中-1代表没有对应关系；

图6为电导矩阵行号到整体矩阵行号的对应关系示意图；

图7为由图2中电路得到的电导矩阵(去除了空行空列)和电压向量做乘法的过程和结果示意图；

图8为得到的电流向量加至节点总电流向量的示意图。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

图1为根据本发明的对电路中线性电阻计算电路方程节点电流的方法流程图，下面将参考图1，对本发明的对电路中线性电阻计算电路方程节点电流的方法进行详细描述。

在步骤11，根据线性电阻的端口位置建立整体电导矩阵。

在该步骤中，遍历所有线性电阻，根据线性电阻的端口位置建立电导矩阵。

图2为根据本发明实施方式的示例电路图，在图中已经标出了节点位置以及线性电阻，其中共包含2个线性电阻和4个电路节点以及一个接地节点。对图2所示的电路中的矩阵建立电导矩阵首先需要对各个电阻分别建立电导矩阵。图3为一个线性电阻往电导矩阵中填入电导的实施方法示意图，其中假设线性电阻的两个端口编号为i和j，电导值为G。在对图2中的G1电阻进行分析后可知，其邻近电路节点为节点1和节点2，根据基尔霍夫电流定律以及电阻电流计算公式，可以根据图3所示的矩阵确定G1电阻的电导矩阵，其中行号和列号分别对应节点号。

每个电阻根据图3的描述在电导矩阵中添加非零元的位置(由所有非零元定下矩阵的稀疏结构，再申请计算机内存，以方便填入电导)。图4为得到的整个电导矩阵，在里面添加了电导的数据，以表明非零元的来源。

建立的整体电导矩阵结构一般是稀疏矩阵，并且结构不会变化，采用压缩方式存储。

在步骤22，对整体电导矩阵进行处理，获得去除空行空列的电导矩阵。

在该步骤中，对整体电导矩阵中的元素位置进行标记，遍历整体电导矩阵，去掉空行空列，对行号和列号重新编号，记录和原始矩阵行列号映射关系。图5记录了整体矩阵行号到电导矩阵行号的映射，图6记录了电导矩阵行号到整体矩阵行号的映射。列号的映射和行号一致。

在步骤33，用电导值填充电导矩阵。

在该步骤中，遍历所有线性电阻，将电导值加入电导矩阵。由于电导矩阵已经去除了空行空列，所以要使用图5中的映射关系以填入电导矩阵正确的位置。

在步骤44，利用电导矩阵计算电流向量。

在该步骤中，将原始电路的电压向量映射成电导矩阵对应的电压向量，并用电导矩阵乘以电压向量，得到电流向量。图7示出了计算过程。

在做矩阵向量乘法时，依据线性电导矩阵的对称性，可以只对上三角或下三角矩阵做乘法并只存储对应的元素，提高计算效率。

在步骤55，计算总的电流向量。

在该步骤中，使用了图6中的映射关系，将利用电导矩阵计算得到的电流向量加至整体节点电流，再叠加其他器件的电流，得到总的电流向量。图8示出了合并的做法过程，其中对于其他剩余的器件，可以通过遍历器件来计算节点电流。

另外，对电路结构连接关系不改变的电路中线性电阻计算电路方程节点电流时，步骤11-步骤33只需进行一次，步骤44-步骤55可以根据电压值改变重复多次。

本发明依据线性电导矩阵结构和值都对称的特点，改进矩阵和向量的乘法的实现，减少乘法次数和存储空间。

本发明还提供对电路中线性电阻计算电路方程节点电流的装置，包括存储器和处理器，存储器上储存有在处理器上运行的程序，处理器运行所述程序时执行上述对电路中线性电阻计算电路方程节点电流的方法的步骤。

本发明还提供了一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行时执行上述的对电路中线性电阻计算电路方程节点电流的方法的步骤，所述对电路中线性电阻计算电路方程节点电流的方法参见前述部分的介绍，不再赘述。

本领域普通技术人员可以理解：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims

1.一种对电路中线性电阻计算电路方程节点电流的方法，其特征在于，包括以下步骤：

1）根据线性电阻的端口位置建立第一电导矩阵；

2）对第一电导矩阵进行去掉空行空列处理，获得第二电导矩阵；

3）用电导值填充第二电导矩阵；

4）利用第二电导矩阵计算电阻电流向量；

5）计算总的电流向量。

2.根据权利要求1所述的对电路中线性电阻计算电路方程节点电流的方法，其特征在于，所述步骤2）进一步包括以下步骤：

21）记录第一电导矩阵各元素的行号和列号；

22）去掉空行和空列，得到第二电导矩阵；

23）记录第一电导矩阵和第二电导矩阵之间各元素的行号和列号的映射关系。

3.根据权利要求2所述的对电路中线性电阻计算电路方程节点电流方法，其特征在于，所述第一电导矩阵和所述第二电导矩阵为对称矩阵。

4.根据权利要求3所述的对电路中线性电阻计算电路方程节点电流的方法，其特征在于，所述步骤4）进一步包括以下步骤：

41）计算第二电导矩阵中各个节点对应的电压值，构成电压向量；

42）将第二电导矩阵与电压向量相乘，获得第二电导矩阵各个节点对应的第二电流向量。

5.根据权利要求4所述的对电路中线性电阻计算电路方程节点电流的方法，其特征在于，所述步骤5）进一步包括，将第二电流向量按照第一电导矩阵和第二电导矩阵之间的映射关系加至各个节点，再与其他器件的等效电流叠加，得到总的节点电流向量。

6.根据权利要求1所述的对电路中线性电阻计算电路方程节点电流的方法，其特征在于，对电路结构连接关系不改变的电路中线性电阻计算电路方程节点电流时，步骤1）-步骤3）只需进行一次，步骤4）-步骤5）可以根据电压值改变重复多次。

7.一种对电路中线性电阻计算电路方程节点电流的装置，其特征在于，包括存储器和处理器，所述存储器上储存有在所述处理器上运行的程序，所述处理器运行所述程序时执行权利要求1-6任一项所述的对电路中线性电阻计算电路方程节点电流的方法的步骤。

8.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述计算机指令运行时执行权利要求1-6任一项所述的对电路中线性电阻计算电路方程节点电流的方法的步骤。