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.

Non-Iterative Simulation Methods for Virtual Analog Modelling

Ducceschi, M
Primo
Conceptualization
;
2022

Abstract

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.
File in questo prodotto:
File Dimensione Formato  
Non-Iterative_Simulation_Methods_for_Virtual_Analog_Modelling.pdf

accesso aperto

Tipo: Versione (PDF) editoriale
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione (CCBY)
Dimensione 1.76 MB
Formato Adobe PDF
1.76 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/901519
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact