Non-Iterative Simulation Methods for Virtual Analog Modelling

doi:10.1109/TASLP.2022.3209934

2022-01-01

Michele Ducceschi *

Stefan Bilbao

DOI：10.1109/TASLP.2022.3209934

原文链接

原文求助

摘要

En 中文

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.

Keyword:

Mathematical models

Numerical models

Stability analysis

Numerical stability

Thermal stability

Iterative methods

Standards

Virtual analog modeling

finite difference schemes

ordinary differential equations