Non-Iterative Simulation Methods for Virtual Analog Modelling

被引:3
作者
Ducceschi, Michele [1 ]
Bilbao, Stefan [2 ]
机构
[1] Univ Bologna, Dept Ind Engn, I-40136 Bologna, Italy
[2] Univ Edinburgh, Acoust & Audio Grp, Edinburgh EH8 9DF, Midlothian, Scotland
基金
欧洲研究理事会;
关键词
Mathematical models; Numerical models; Stability analysis; Numerical stability; Thermal stability; Iterative methods; Standards; Virtual analog modeling; finite difference schemes; ordinary differential equations; SCHEMES; DISCRETIZATION;
D O I
10.1109/TASLP.2022.3209934
中图分类号
O42 [声学];
学科分类号
070206 ; 082403 ;
摘要
The simulation of nonlinear components is central to virtual analog simulation. In audio effects, circuits often include devices such as diodes and transistors, mostly operating in a strongly nonlinear regime. Mathematical models are in the form of systems of nonlinear ordinary differential equations (ODEs), and traditional integrators, such as the trapezoid and midpoint methods, can be employed as solvers. These methods are fully implicit and require the solution of a nonlinear algebraic system at each time step, introducing further complications regarding the existence and uniqueness of the solution, as well as the choice of halting conditions for the iterative root finder. On the other hand, fast explicit methods such as Forward Euler, are prone to unstable behaviour at standard audio sample rates. For these reasons, in this work, a family of linearly-implicit schemes is presented. These schemes take the form of a perturbation expansion, making the construction of higher-order schemes possible. Compared with classic implicit designs, the proposed methods have the advantage of efficiency, since the update is computed in a single iteration. Furthermore, the existence and uniqueness of the update are proven by simple inspection of the update matrix. Compared to classic explicit designs, the proposed schemes display stable behaviour at standard audio sample rates. In the case of a single scalar ODE, sufficient conditions for numerical stability can be derived, imposing constraints on the choice of the sampling rate. Several theoretical results are provided, as well as numerical examples for typical stiff equations used in virtual analog modelling.
引用
收藏
页码:3189 / 3198
页数:10
相关论文
共 38 条
[1]  
[Anonymous], 2010, P DIG AUD EFF GRAZ A
[2]  
[Anonymous], 2006, P INT C DIG AUD EFF
[3]  
Bernardini A, 2020, ADV WAVE DIGITAL MOD, P3
[4]   A Wave Digital Newton-Raphson Method for Virtual Analog Modeling of Audio Circuits with Multiple One-Port Nonlinearities [J].
Bernardini, Alberto ;
Bozzo, Enrico ;
Fontana, Federico ;
Sarti, Augusto .
IEEE-ACM TRANSACTIONS ON AUDIO SPEECH AND LANGUAGE PROCESSING, 2021, 29 :2162-2173
[5]   Linear Multistep Discretization Methods With Variable Step-Size in Nonlinear Wave Digital Structures for Virtual Analog Modeling [J].
Bernardini, Alberto ;
Maffezzoni, Paolo ;
Sarti, Augusto .
IEEE-ACM TRANSACTIONS ON AUDIO SPEECH AND LANGUAGE PROCESSING, 2019, 27 (11) :1763-1776
[6]   Higher-order accurate two-step finite difference schemes for the many-dimensional wave equation [J].
Bilbao, Stefan ;
Hamilton, Brian .
JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 367 :134-165
[7]  
Butcher J. C., 2003, NUMERICAL METHODS OR, DOI [DOI 10.1002/0470868279, 10.1002/0470868279]
[8]  
Danish M, 2021, INT CONF DIGIT AUDIO, P33, DOI 10.23919/DAFx51585.2021.9768301
[9]  
Dempwolf K., 2011, PROC INT C DIGIT AUD
[10]  
Ducceschi M, 2021, INT CONF DIGIT AUDIO, P25, DOI 10.23919/DAFx51585.2021.9768254