Optimizing Differentiated Discretization for Audio Circuits Beyond Driving Point Transfer Functions

One goal of Virtual Analog modeling of audio circuits is to produce digital models whose behavior matches analog prototypes as closely as possible. Discretization methods provide a systematic approach to generate those models but they introduce frequency response error, such as the well-known frequency warping for the trapezoidal method. Recent work showed how using different discretization methods for each reactive elements could reduce such error. It further provided a procedure that optimizes that error according to a chosen metric through joint selection of the discretization parameters. In this work, we extend that approach to the general case of transfer functions with one input and an arbitrary number of out- puts expressed as linear combinations of the network variables. To demonstrate the validity of our approach, we apply the optimization procedure for the response of a Hammond organ vibrato/chorus ladder filter, a 19-output, 36th order filter, where each output frequency response presents many features spread across its passband.