CN112417794B - 一种散射参数计算方法 - Google Patents
一种散射参数计算方法 Download PDFInfo
- Publication number
- CN112417794B CN112417794B CN202011414046.XA CN202011414046A CN112417794B CN 112417794 B CN112417794 B CN 112417794B CN 202011414046 A CN202011414046 A CN 202011414046A CN 112417794 B CN112417794 B CN 112417794B
- Authority
- CN
- China
- Prior art keywords
- parameter
- circuit
- calculating
- port
- scattering
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 230000005284 excitation Effects 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 12
- 230000004044 response Effects 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 235000013599 spices Nutrition 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/31—Design entry, e.g. editors specifically adapted for circuit design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种散射参数计算方法，包括以下步骤：根据设置的频率，基于AC交流分析，计算出电路网络的导纳矩阵；根据给定的交流激励信号，逐个计算出每一个端口的响应电压；根据所述导纳矩阵和所述端口的响应电压，计算电路散射参数。本发明的散射参数计算方法，能够处理不同参考阻抗的散射参数计算，对于大规模的电路，特别是对于多端口电路，大大提高了计算效率。
Description
技术领域
本发明涉及集成电路计算机辅助设计(Integrated Circuit/Computer AidedDesign)技术领域，特别是涉及电路仿真中快速准确求解散射参数的方法。
背景技术
散射参数用来模拟电子电路在不同频率下的行为，可以将电路描述成一个黑盒子。散射参数分析(也可以简称为S参数分析)是spice仿真器的一种分析类型，目的是使电路在直流工作点附近线性化，并且把原始电路等效成一个包含N个端口的等效电路，并输出在各个离散的频率点处的S参数值。
现有技术常用Y参数通过一系列的矩阵变换得到S参数，该方法可以很直观的根据端口的电流电压关系得到电路各个端口间的导纳系数矩阵即Y参数，然后通过一些列矩阵运算得到散射参数矩阵即S参数。还有一种直接计算S参数的方法，但是该方法只给出了相同参考阻抗的S参数计算方法。
发明内容
为了解决现有技术存在的不足，本发明的目的在于提供一种散射参数计算方法，能够处理不同参考阻抗的散射参数计算，对于大规模的电路，特别是对于多端口电路，大大提高了计算效率。
为实现上述目的，本发明提供的一种散射参数计算方法，包括以下步骤：
根据设置的频率，基于AC交流分析，计算出电路网络的导纳矩阵；
根据给定的交流激励信号，逐个计算出每一个端口的响应电压；
根据所述导纳矩阵和所述端口的响应电压，计算电路散射参数。
进一步地，所述根据设置的频率，基于AC交流分析，计算出电路网络的导纳矩阵的步骤，进一步包括，根据端口对应的特征阻抗建立AC交流分析的电路方程。
进一步地，所述根据给定的交流激励信号，逐个计算出每一个端口的响应电压的步骤，进一步包括，
求解对应的AC交流分析的线性方程组，得到解向量；
将所述解向量中对应于AC交流激励信号源位置的分量端口位置减去1；
所述解向量中其他位置的分量乘以参考阻抗比例因子。
更进一步地，所述交流激励信号，其强度为2伏。
为实现上述目的，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序运行时执行如上文所述的散射参数计算方法步骤。
为实现上述目的，本发明还提供一种散射参数计算设备，包括存储器和处理器，所述存储器上储存有在所述处理器上运行的计算机程序，所述处理器运行所述计算机程序时执行如上文所述的散射参数计算方法步骤。
本发明的一种散射参数计算方法，具有以下有益效果：
1)可以处理不同参考阻抗的散射参数计算。
2)一方面避免了从Y参数到S参数的复杂转换过程，另一方面也避免了在转换过程中的数值误差积累。
3)对于大规模的电路，特别是对于多端口电路，大大提高了计算效率。
本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，并与本发明的实施例一起，用于解释本发明，并不构成对本发明的限制。在附图中：
图1为根据本发明的散射参数计算方法流程图；
图2为根据本发明的实施例一两端口网络示意图；
图3为根据本发明的实施例一端口电路结构示意图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
图1为根据本发明的散射参数计算方法流程图，下面将参考图1，对本发明的散射参数计算方法进行详细描述。
首先，在步骤101，对原始多端口网络进行AC分析。
优选地，设置一个频率，计算AC分析(交流分析)的导纳矩阵。
本实施例中，基于AC分析的结果，以及不同port处的参考阻抗来直接计算散射参数。假定电路有n个port，每个port的参考阻抗是Ri(其中i＝1,...,n)。
在步骤102，分别在每个端口逐次添加2伏的AC激励信号源，并根据AC分析结果计算其余端口相应的AC响应。
优选地，逐个在每个端口添加2伏的激励信号，通过仿真器求解AC分析的结果，得到其余各端口相应的AC激励源的响应。
在步骤103，建立AC分析的线性方程组，得到解向量，即每个端口的电压。
优选地，建立满足Y参数的电路方程(公式1)，其中Y指Y参数，I是单位矩阵，Z0是对角矩阵，对角元素分别是端口对应的特征阻抗。
建立S参数与S1之间的等价关系(公式2)。
根据公式1中Y参数的方程求解AC分析得到的v和i之后，计算S参数的公式(公式3)，将公式1和公式2代入公式3计算原始电路的S参数。
在步骤104，将解向量对应于AC交流激励信号源位置的分量端口位置减去1。
在步骤105，将解分量里其他位置的分量乘以参考阻抗比例因子
下面结合一具体实施例对本发明的布线方法做进一步的说明。
图2为根据本发明的实施例一两端口网络示意图，如图2所示，一个简单的两端口电路网络，端口特征阻抗是不同的，指定需要模拟的端口，分别对两端口轮流指定一个2伏的ac激励源，通过仿真器求解ac分析的结果，对其中一个端口指定AC激励信号时，得到两个端口对该端口的响应；根据S参数计算公式，求得系统的S参数，仿真出该待分析电路网络的S参数特性。
图3为根据本发明的实施例一端口电路结构示意图，如图3所示，带有特征阻抗的0伏电压源，每个端口的特征阻抗可以不同。
本发明提供一种散射参数计算方法，由于电路矩阵的条件数较差，通过复杂的矩阵计算会产生较大的数值误差，因此我们采用直接计算S参数的方法，将已有的相同参考阻抗的S参数计算方法推广到不同参考阻抗的情形，从而满足设计工程师不同的设计需求。
本发明的一个实施例中，还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序运行时执行如上文所述的散射参数计算方法的步骤。
本发明的一个实施例中，还提供一种散射参数计算设备，包括存储器和处理器，所述存储器上储存有在所述处理器上运行的计算机程序，所述处理器运行所述计算机程序时执行如上文所述的散射参数计算方法的步骤。
本领域普通技术人员可以理解：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种散射参数计算方法，其特征在于，包括以下步骤：
根据设置的频率，基于AC交流分析，计算出电路网络的导纳矩阵；
根据给定的交流激励信号，逐个计算出每一个端口的响应电压；
根据所述导纳矩阵和所述端口的响应电压，计算电路散射参数；
建立满足Y参数的电路方程其中，Y指Y参数，I是单位矩阵，Z0是对角矩阵，对角元素分别是端口对应的特征阻抗；
建立S参数与S1之间的等价关系式
根据Y参数的电路方程求解AC分析得到的v和i之后，计算S参数的公式将所述Y参数的电路方程和S参数与S1之间的等价关系式代入S参数的公式，计算原始电路的S参数。
2.根据权利要求1所述的散射参数计算方法，其特征在于，所述根据设置的频率，基于AC交流分析，计算出电路网络的导纳矩阵的步骤，进一步包括，根据端口对应的特征阻抗建立AC交流分析的电路方程。
3.根据权利要求1所述的散射参数计算方法，其特征在于，所述根据给定的交流激励信号，逐个计算出每一个端口的响应电压的步骤，进一步包括，
求解对应的AC交流分析的线性方程组，得到解向量；
将所述解向量中对应于AC交流激励信号源端口位置的分量位置减去1；
所述解向量中其他位置的分量乘以参考阻抗比例因子。
4.根据权利要求1所述的散射参数计算方法，其特征在于，所述交流激励信号，其强度为2伏。
5.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序运行时执行权利要求1至4任一项所述的散射参数计算方法步骤。
6.一种散射参数计算设备，其特征在于，包括存储器和处理器，所述存储器上储存有在所述处理器上运行的计算机程序，所述处理器运行所述计算机程序时执行权利要求1至4任一项所述的散射参数计算方法步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011414046.XA CN112417794B (zh) | 2020-12-04 | 2020-12-04 | 一种散射参数计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011414046.XA CN112417794B (zh) | 2020-12-04 | 2020-12-04 | 一种散射参数计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112417794A CN112417794A (zh) | 2021-02-26 |
| CN112417794B true CN112417794B (zh) | 2024-03-08 |
Family
ID=74774934
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011414046.XA Active CN112417794B (zh) | 2020-12-04 | 2020-12-04 | 一种散射参数计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112417794B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004333320A (ja) * | 2003-05-08 | 2004-11-25 | Matsushita Electric Ind Co Ltd | 高周波電圧印加システム、高周波電圧印加方法、及び高周波電圧印加プログラム |
| CN101629980A (zh) * | 2009-09-10 | 2010-01-20 | 南京师范大学 | 基于散射参数的emi滤波器性能测试方法 |
| CN105184032A (zh) * | 2015-11-04 | 2015-12-23 | 河南科技大学 | 一种改进的hbt小信号等效电路模型建立方法 |
| CN106199212A (zh) * | 2016-07-07 | 2016-12-07 | 中国人民解放军国防科学技术大学 | 一种多导体传输线网络的时频分析方法 |
| CN107247685A (zh) * | 2017-05-26 | 2017-10-13 | 京信通信系统（中国）有限公司 | Mems器件端口特性参数提取方法和装置 |
| CN110750929A (zh) * | 2019-10-15 | 2020-02-04 | 电子科技大学 | 一种基于特征模理论的天线辐射和散射特性的分析方法 |
| CN110781609A (zh) * | 2019-11-11 | 2020-02-11 | 中国电子科技集团公司第二十九研究所 | 一种滤波器s参数去相位加载方法、系统、介质和设备 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030065498A1 (en) * | 2001-07-06 | 2003-04-03 | Bois Karl J. | Software tool for generation of scattering parameter models of N-port lumped element circuits for use in spice simulators |
| WO2015090478A1 (en) * | 2013-12-20 | 2015-06-25 | Advantest Corporation | Multi-port measurement technique for determining s-parameters |
-
2020
- 2020-12-04 CN CN202011414046.XA patent/CN112417794B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004333320A (ja) * | 2003-05-08 | 2004-11-25 | Matsushita Electric Ind Co Ltd | 高周波電圧印加システム、高周波電圧印加方法、及び高周波電圧印加プログラム |
| CN101629980A (zh) * | 2009-09-10 | 2010-01-20 | 南京师范大学 | 基于散射参数的emi滤波器性能测试方法 |
| CN105184032A (zh) * | 2015-11-04 | 2015-12-23 | 河南科技大学 | 一种改进的hbt小信号等效电路模型建立方法 |
| CN106199212A (zh) * | 2016-07-07 | 2016-12-07 | 中国人民解放军国防科学技术大学 | 一种多导体传输线网络的时频分析方法 |
| CN107247685A (zh) * | 2017-05-26 | 2017-10-13 | 京信通信系统（中国）有限公司 | Mems器件端口特性参数提取方法和装置 |
| CN110750929A (zh) * | 2019-10-15 | 2020-02-04 | 电子科技大学 | 一种基于特征模理论的天线辐射和散射特性的分析方法 |
| CN110781609A (zh) * | 2019-11-11 | 2020-02-11 | 中国电子科技集团公司第二十九研究所 | 一种滤波器s参数去相位加载方法、系统、介质和设备 |
Non-Patent Citations (4)
| Title |
|---|
| PSPICE优化技术及其在微波电路计算机辅助设计中的应用;桂垣等;《河北建筑工程学院学报》;20030630;第21卷(第2期);全文 * |
| S散射参数的普遍化理论;梁昌洪;《西北电讯工程学院学报》;19841231(第4期);全文 * |
| Voltage controlled scattering from porous silicon Mie-particles in liquid crystals;Lakshmi Madhuri Pappu等;《Journal of Molecular Liquids》;20190501;全文 * |
| 复杂谐振系统散射参数的计算;李润旗，高葆新;电子学报;19940625(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112417794A (zh) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5946482A (en) | Method and apparatus for using parameters to simulate an electronic circuit | |
| CN106777827A (zh) | 一种机电‑电磁混合仿真方法及系统 | |
| CN111859835A (zh) | 一种电路互连网络模型的降阶方法、降阶装置及降阶设备 | |
| Dolin et al. | Modified methods of circuit simulation of radio engineering devices in the time domain | |
| US7363208B2 (en) | Power consumption estimation | |
| CN106372348A (zh) | 一种线性系统中基于误差控制的矢量拟合模型降阶方法 | |
| Yao et al. | Contingency analysis based on partitioned and parallel holomorphic embedding | |
| CN112417794B (zh) | 一种散射参数计算方法 | |
| CN104881508B (zh) | 一种基于查表法的半导体器件的建模方法和系统 | |
| CN110703038B (zh) | 一种适用于风机接入配电网的谐波阻抗估算方法 | |
| CN117277350A (zh) | 阻抗计算、电网稳定性分析方法及存储介质和终端设备 | |
| CN109033660B (zh) | 不对称故障分析方法及装置 | |
| CN110135090A (zh) | 一种基于响应面法的电路系统容差建模与分析方法 | |
| Rauh et al. | Verification and reachability analysis of fractional-order differential equations using interval analysis | |
| CN115864502A (zh) | 交直流电力系统小信号稳定性分析与致稳设计方法及装置 | |
| CN109670254B (zh) | 机电暂态与电磁暂态混合仿真的接口位置选择方法及系统 | |
| CN112818537A (zh) | 一种光伏并网系统稳定性分析方法及装置 | |
| CN110350507A (zh) | 直流微网等效电容确定方法、装置及计算机可读存储介质 | |
| CN114692903A (zh) | 用于设备故障检测的方法及终端设备 | |
| CN109388911A (zh) | 一种半导体器件的等效模型建立方法、装置及终端设备 | |
| CN110502818B (zh) | 一种多功率接口数模混合仿真系统准确性评估方法 | |
| CN112632883B (zh) | 器件模型的仿真结果的测试方法、装置、设备和介质 | |
| Wang et al. | Gramian‐based model order reduction of parameterized time‐delay systems | |
| CN112257323B (zh) | 基于电磁暂态并行计算的自动分网方法、装置及终端设备 | |
| Elfadel | Convergence of transverse waveform relaxation for the electrical analysis of very wide transmission line buses |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |