CN113255275A - 一种基于不光滑波形的时间离散格式切换方法 - Google Patents
一种基于不光滑波形的时间离散格式切换方法 Download PDFInfo
- Publication number
- CN113255275A CN113255275A CN202110557125.4A CN202110557125A CN113255275A CN 113255275 A CN113255275 A CN 113255275A CN 202110557125 A CN202110557125 A CN 202110557125A CN 113255275 A CN113255275 A CN 113255275A
- Authority
- CN
- China
- Prior art keywords
- time
- unsmooth
- point
- format
- discrete
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于不光滑波形的时间离散格式切换方法,包括以下步骤:1)确定仿真过程中的时域不光滑点;2)将时间离散格式切换为隐式Euler格式计算下一步的时间离散。本发明的基于不光滑波形的时间离散格式切换方法,能够判断在仿真过程中遇到波形不光滑的情形并切换时间离散格式,减小时间离散的误差。
Description
技术领域
本发明涉及EDA模拟电路时域仿真技术领域,尤其涉及时域仿真中的时间离散格式切换。
背景技术
在EDA模拟电路时域仿真中,需要对通过电路网络建模得到的微分代数方程(DAE)进行时间离散。离散的常见格式有隐式Euler格式、梯形格式、Gear格式等,其中Euler格式为1阶精度,梯形格式和Gear格式为二阶精度。
上述离散格式的精度仅对解是光滑的情形有效,然而在实际电路仿真过程中,需要针对非光滑的情形进行处理。例如,对于梯形脉冲的电压源,如果在不光滑的拐点处仍采用高精度格式进行离散,就会引入额外的误差。事实上,对于波形不光滑的情形,DAE解的定义本身存在问题。目前的处理方案只是假定在不光滑点附近做了以时间步长为尺度的磨光处理,很难找到精确的结果。
发明内容
为了解决现有技术存在的不足,本发明的目的在于提供一种基于不光滑波形的时间离散格式切换方法,以判断在仿真过程中遇到波形不光滑的情形并切换时间离散格式,减小时间离散的误差。
为实现上述目的,本发明提供的一种基于不光滑波形的时间离散格式切换方法,包括以下步骤:
1)确定仿真过程中的时域不光滑点;
2)将时间离散格式切换为隐式Euler格式计算下一步的时间离散。
进一步地,所述步骤1)还包括,根据元件本身波形的不光滑特性确定不光滑点。
进一步地,所述步骤1)还包括,利用一阶向前差商计算出当前仿真的时间点上所有网格点的电荷的数值导数:如果某个网格点上的电荷的数值导数的绝对值大于该网格点前6个时间步的数值导数绝对值的平均值的20倍,则将该时间点标记为不光滑点。
进一步地,所述步骤1)还包括,
计算出当前点与前一个时间点电荷数值导数的差值的绝对值D,计算出当前点前5个时间步的电荷数值导数的绝对值的算数平均值ED;
如果网格点满足D/ED>0.7,则将该时间点标记为不光滑点。
进一步地,所述步骤2)还包括,切换为隐式Euler格式的时间步数等于原有时间离散格式的精度阶。
更进一步地,如果当前不光滑点的原有时间离散格式是隐式Euler格式,则无需进行切换。
为实现上述目的,本发明还提供一种电子设备,包括存储器和处理器,所述存储器上储存有在所述处理器上运行的程序,所述处理器运行所述程序时执行上述基于不光滑波形的时间离散格式切换方法的步骤。
为实现上述目的,本发明还提供一种计算机可读存储介质,其上存储有计算机指令,所述计算机指令运行时执行上述基于不光滑波形的时间离散格式切换方法的步骤。
有益效果:本发明和不切换时间离散格式相比,其时间格式切换策略能够明显减少仿真结果中不光滑点附近出现破坏守恒性等非物理结果的情形。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,并与本发明的实施例一起,用于解释本发明,并不构成对本发明的限制。在附图中:
图1为根据本发明的基于不光滑波形的时间离散格式切换方法流程图;
图2为仿真过程中的梯形脉冲波形示意图;
图3为根据本发明实施方式的另一仿真波形部分示意图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
本发明实施例中,在进行仿真时采用的是自适应时间步长策略,时间步长会随时间而改变。但本发明本身不涉及自适应时间步长策略,因此假定:当仿真进行到当前时间点时,所有之前的时间步长均已经确定。
实施例1
图1为根据本发明的基于不光滑波形的时间离散格式切换方法流程图,下面将参考图1,对本发明的基于不光滑波形的时间离散格式切换方法进行详细描述。
在步骤101,识别仿真过程中的时域不光滑点。
本发明实施例中,在对不光滑点进行处理之前,首先需要识别出仿真过程中的时域不光滑点。本发明通过判断仿真时间点是否满足以下三个条件动态识别仿真过程中的不光滑点:
a)根据元件本身波形的不光滑特性直接确定不光滑点,例如如果网络中包含梯形脉冲的电压源,则梯形的角点则被主动识别为不光滑点;
b)利用一阶向前差商计算出当前仿真的时间点上所有网格点的电荷Q的数值导数。如果某个网格点上的数值导数的绝对值大于该网格点前6个时间步的数值导数绝对值的平均值的20倍,则将该时间点标记为不光滑点。
c)计算出当前点与前一个时间点的电荷数值导数的差值的绝对值D和当前点前5个时间步的电荷数值导数的绝对值的算数平均值ED,如果网格点满足D/ED>0.7,则将该时间点标记为不光滑点。
本发明实施例中,只要仿真时间点满足上述a)、b)和c)之一,就被标记为不光滑点;如果仿真时间点不满足上述任何一个条件,则该时间点被标记为光滑点。
步骤102,切换时间离散格式。
当通过步骤101判断某仿真时间点为不光滑点时,计算下一步的时间离散时,时间离散格式会被临时切换为隐式Euler格式。临时切换为隐式Euler格式的时间步数等于原有格式的精度阶。例如,如果原有格式为2阶精度的Gear格式,则切换为隐式Euler格式执行两步时间推进,往后的时间推进再回到原有的2阶Gear格式。如果原有格式就是隐式Euler格式,则无需进行切换。
实施例2
图2为仿真过程中的梯形脉冲波形示意图,如果电路网络中存在如图2所示的梯形脉冲电压源,则根据条件a),根据波形的不光滑特性,可在图2中所示的A、B、C、D点均会被标记为不光滑点。
实施例3
如果在仿真过程中的波形出现如图3所示的部分波形,则图中所示的A、B、C点均会被标记为不光滑点。其中A点的判定由条件b)所确定,B、C两点的判定由条件c)所确定。
另外,本发明的针对不光滑波形切换时间离散格式的方法仅限于求解电路时域仿真的DAE中的时间离散方法。不适用于求解其他DAE问题的时间离散方法或求解电路时域仿真的DAE中与时间无关的其他方法。
实施例4
本发明还提供针一种电子设备,包括存储器和处理器,存储器上储存有在处理器上运行的程序,处理器运行所述程序时执行上述一种基于不光滑波形的时间离散格式切换方法的步骤。
实施例5
本发明还提供了一种计算机可读存储介质,其上存储有计算机指令,所述计算机指令运行时执行上述的一种基于不光滑波形的时间离散格式切换方法的步骤,所述一种基于不光滑波形的时间离散格式切换方法参见前述部分的介绍,不再赘述。
本领域普通技术人员可以理解:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于不光滑波形的时间离散格式切换方法,包括以下步骤:
1)确定仿真过程中的时域不光滑点;
2)将时间离散格式切换为隐式Euler格式计算下一步的时间离散。
2.根据权利要求1所述的基于不光滑波形的时间离散格式切换方法,其特征在于,所述步骤1)还包括,根据元件本身波形的不光滑特性确定不光滑点。
3.根据权利要求1所述的基于不光滑波形的时间离散格式切换方法,其特征在于,所述步骤1)还包括,利用一阶向前差商计算出当前仿真的时间点上所有网格点的电荷的数值导数:如果某个网格点上的电荷的数值导数的绝对值大于该网格点前6个时间步的数值导数绝对值的平均值的20倍,则将该时间点标记为不光滑点。
4.根据权利要求1所述的基于不光滑波形的时间离散格式切换方法,其特征在于,所述步骤1)还包括,
计算出当前点与前一个时间点电荷数值导数的差值的绝对值D,计算出当前点前5个时间步的电荷数值导数的绝对值的算数平均值ED;
如果网格点满足D/ED>0.7,则将该时间点标记为不光滑点。
5.根据权利要求1所述的基于不光滑波形的时间离散格式切换方法,其特征在于,所述步骤2)还包括,切换为隐式Euler格式的时间步数等于原有时间离散格式的精度阶。
6.根据权利要求1所述的基于不光滑波形的时间离散格式切换方法,其特征在于,如果当前不光滑点的时间离散格式是隐式Euler格式,则无需进行切换。
7.一种针对不光滑波形切换时间离散格式的装置,其特征在于,包括存储器和处理器,所述存储器上储存有在所述处理器上运行的程序,所述处理器运行所述程序时执行权利要求1-6任一项所述的基于不光滑波形的时间离散格式切换方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机指令,其特征在于,所述计算机指令运行时执行权利要求1-6任一项所述的基于不光滑波形的时间离散格式切换方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110557125.4A CN113255275B (zh) | 2021-05-21 | 2021-05-21 | 一种基于不光滑波形的时间离散格式切换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110557125.4A CN113255275B (zh) | 2021-05-21 | 2021-05-21 | 一种基于不光滑波形的时间离散格式切换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113255275A true CN113255275A (zh) | 2021-08-13 |
CN113255275B CN113255275B (zh) | 2022-05-24 |
Family
ID=77183507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110557125.4A Active CN113255275B (zh) | 2021-05-21 | 2021-05-21 | 一种基于不光滑波形的时间离散格式切换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113255275B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050262454A1 (en) * | 1999-09-27 | 2005-11-24 | International Business Machines Corporation | Method for VLSI system debug and timing analysis |
US20110071812A1 (en) * | 2009-09-24 | 2011-03-24 | Texas Instruments Incorporated | Small-Signal Stability Analysis at Transient Time Points of Integrated Circuit Simulation |
CN105740528A (zh) * | 2016-01-27 | 2016-07-06 | 东南大学 | 一种自适应同步方式的电力通信联合仿真平台及其同步方法 |
CN107193036A (zh) * | 2017-06-26 | 2017-09-22 | 成都理工大学 | 一种改进型核信号梯形脉冲成形方法和装置 |
US20190208348A1 (en) * | 2016-09-01 | 2019-07-04 | Universiteit Antwerpen | Method of determining a personalized head-related transfer function and interaural time difference function, and computer program product for performing same |
-
2021
- 2021-05-21 CN CN202110557125.4A patent/CN113255275B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050262454A1 (en) * | 1999-09-27 | 2005-11-24 | International Business Machines Corporation | Method for VLSI system debug and timing analysis |
US20110071812A1 (en) * | 2009-09-24 | 2011-03-24 | Texas Instruments Incorporated | Small-Signal Stability Analysis at Transient Time Points of Integrated Circuit Simulation |
CN105740528A (zh) * | 2016-01-27 | 2016-07-06 | 东南大学 | 一种自适应同步方式的电力通信联合仿真平台及其同步方法 |
US20190208348A1 (en) * | 2016-09-01 | 2019-07-04 | Universiteit Antwerpen | Method of determining a personalized head-related transfer function and interaural time difference function, and computer program product for performing same |
CN107193036A (zh) * | 2017-06-26 | 2017-09-22 | 成都理工大学 | 一种改进型核信号梯形脉冲成形方法和装置 |
Non-Patent Citations (4)
Title |
---|
CLAUDE BARDOS 等: "The Euler equations in planar nonsmooth convex domains", 《JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS》 * |
刘勇刚: "基于时域分析的电缆线路故障测距", 《电气开关》 * |
孙明轩 等: "采用干扰差分补偿的无切换吸引离散时间控制方法", 《控制与决策》 * |
张伟伟 等: "非定常流动模拟的时间离散方法", 《力学进展》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113255275B (zh) | 2022-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9443041B2 (en) | Simulation system and method for testing a simulation of a device against one or more violation rules | |
CN107025335B (zh) | 基于状态变量离散化的仿真计算方法和仿真系统 | |
TW201602819A (zh) | 調整用於靜態時序分析的時序降額的方法 | |
CN102436533A (zh) | 标准单元库模型的时序验证方法 | |
März et al. | Recent results in solving index-2 differential-algebraic equations in circuit simulation | |
CN103699013A (zh) | 一种有效提高电机硬件在环实时仿真速度和精度的方法 | |
CN112115676B (zh) | 静态电压降修复方法、装置、设备及存储介质 | |
US20140096100A1 (en) | Method of sharing and re-using timing models in a chip across multiple voltage domains | |
CN113255275B (zh) | 一种基于不光滑波形的时间离散格式切换方法 | |
CN109190210B (zh) | 基于Saber平台建模仿真的电路性能可靠性分析方法 | |
CN110765584A (zh) | 一种含多开关元件的电磁暂态仿真算法、系统、介质及设备 | |
EP3579126A1 (en) | Co-simulation method and device | |
KR101449212B1 (ko) | Can 버스 시뮬레이터의 제어기 모델링 장치 및 방법 | |
CN113033128B (zh) | 一种选取电路仿真中牛顿迭代的初值的方法 | |
JP2013025467A (ja) | 対数関数の近似演算回路 | |
CN110968938B (zh) | 一种用于电磁暂态仿真的理想开关过程分析方法及系统 | |
JP2009122732A (ja) | 動作タイミング検証装置及びプログラム | |
US8151233B1 (en) | Circuit design with incremental simultaneous switching noise analysis | |
CN112989739B (zh) | 一种Trap-Gear时间离散格式的时间步长设定方法 | |
TW201721419A (zh) | 系統模擬方法 | |
CN113255268B (zh) | 一种电路仿真中瞬态分析不收敛的检测及修复方法 | |
Loeckx et al. | Efficient identification of major contributions to EMI-induced rectification effects in analog automotive circuits | |
TWI608372B (zh) | 電力狀態覆蓋度量標準及其估計方法 | |
CN112989738B (zh) | 一种电路仿真中牛顿迭代的收敛判定的改进方法 | |
US7516429B2 (en) | System for analyzing an electronic circuit described by characterization data |
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 |