Laplace Transfer Function - Convolution Implementation

Netlist Entry

Uxxxx p n cp cn modelname [M=m IC=ic VOLTAGE_MODE=vm]

Model Syntax

.model modelname LAPLACE ( parameters )

Model Parameters

Laplace transfer function

The Laplace expression defines the behaviour of the device in the frequency domain. For example '1/(s+1)', defines a simple single pole low-pass filter. The expression may contain arithmetic operators and a number of functions as described in the following sections. Be aware that not any expression is physically realisable; for example 'exp(s)' defines a negative delay, that is a device whose output responds to an input in the future.

Operators

+ - * / ^

Constants

Any decimal number following normal rules. SPICE style engineering suffixes are accepted.

s Variable

This can be raised to a power using ^, for example s^2. If the power is an integer between 0 and 9 the ^ may be omitted. For example: s2 is the same as s^2.

Functions

Filter response functions

Filter response functions may be used in both lumped network implementation and convolution implementation.

These are described in the following table:

Lookup tables

The frequency response of a system may be defined in tabular form using lookup tables. The lookup table consists of a sequence of values arranged in triplets. Each triplet is in the form ???MATH???frequency, value1, value2???MATH??? where ???MATH???value1???MATH??? and ???MATH???value2???MATH??? define the magnitude and phase in various ways as described below.

There are five variants of the lookup table as described below:

Example

Table_M(0.01,1,30,  1,1,30,  10,0.1,30,  100,0.01,-30)

The above defines a response with a gain of 1 and a phase shift of 30 degrees from 0.01Hz to 1Hz. From 1Hz to 10Hz the gain falls from 1 to 0.1 with the phase shift remaining at 30 degrees. From 10Hz to 100Hz the gain falls from 0.1 to 0.01 and the phase changes from 30 to -30 degrees.

In real life it is not possible to implement a characteristic with the above behaviour as it is non-causal. However, it can be implemented if a suitable delay is added to the characteristic. The delay may be added using the DELAY parameter. Alternatively it may be added to the Laplace expression using ???MATH???exp(-delay.s)???MATH???. E.g:

Table_M(0.01,1,30,  1,1,30,  10,0.1,30,  100,0.01,-30) * exp(-s*2)

The above adds a 2 second delay.

For ease of reading, each triplet in the table may be placed on a separate line by prefixing each line with a '+' character. E.g.

Table_M(
+ 0.01,		1,		30,
+ 1,		1,		30,
+ 10,		0.1,		30,
+ 100,		0.01,		-30)

Interpolation Between defined frequency points, the Laplace transfer function finds the magnitude and phase using interpolation. The interpolation is performed logarithmically on the magnitude data. Values outside the table frequency range are defined by their respective end points.

If the table contains zero frequency terms or zero magnitude terms, the values are replaced with the values of parameters TABLE_ZERO_FREQ and TABLE_ZERO_MAG parameters respectively. This is because the interpolation is performed logarithmically and this requires all values to be > 0.

For the above example the magnitude and phase are shown below:

Parameters

Parameters (defined with .PARAM) may be used in Laplace expressions.

Implementation

The Laplace transfer function is implemented using a convolution method.

In the frequency domain, a linear system may be represented by: \[ Vout(s) = Vin(s).f(s) \] where ???MATH???f(s)???MATH??? is the transfer function.

In the time domain the same system may be represented by: \[ Vout(t) = Vin(t) \ast f(t) \] Where ???MATH???f(t)???MATH??? is the impulse response of ???MATH???f(s)???MATH??? and ???MATH???\ast???MATH??? is the convolution operator.

This method is challenging to implement as simple convolution applied to simulation data in its raw form is an ???MATH???\mathcal{O}(n^2)???MATH??? algorithm. This means that the number of computations required at each step is proportional to the square of the number of steps. In practice this becomes unacceptably slow when there are more than about 10000 time steps.

SIMetrix overcomes this performance limitation by using an FFT based fast convolution method with a computation speed of ???MATH???\mathcal{O}(n\log{}^2n)???MATH???. Although dramatically faster, this algorithm requires data to be presented to it at fixed-interval steps which has the effect of placing an upper frequency limit. The actual upper frequency limit can be arbitrarily increased by increasing the number of steps (the value of ???MATH???n???MATH???) at the expense of speed and the memory consumed. In practice the default settings work well in most applications and give good performance and accuracy.

Impulse Response

Problems with Impulse Response Extraction

For the convolution method an impulse response of the Laplace transfer function must be extracted. As mentioned above, three methods are attempted. The inverse FFT method and Stehfest methods are approximate in nature and for this reason, the result from each method is tested. If the test fails the impulse response is rejected and another method is attempted. If all methods fail the simulation will abort.

Transfer functions that decay slowly are likely to fail the inverse FFT method as this method requires that the impulse response decay to zero or nearly to zero. If the Stehfest method also fails, the problem can sometimes be resolved by increasing the inverse FFT size. This is controlled by the IMPULSE_IFFT_N parameter. This increases the duration over which the inverse FFT is evaluated so providing a longer time to decay.

In many situations the impulse response extraction fails because the transfer function is non-causal. This means that it has a response in negative time. Non-causal transfer functions can be made causal by adding a delay.

Impulse Response Testing

If the inverse FFT method or Stehfest methods are used to extract the impulse response, the result is tested before allowing it to be used for the simulation.

The Stehfest method is tested by performing a spot frequency check. That is the impulse response is used to measure the gain of a sequence of sine waves in the time domain to check that the magnitudes match to that predicted in the frequency domain. The parameters IMPULSE_TEST_FMAX and IMPULSE_TEST_FMIN determine the range of frequencies used.

The inverse FFT method may be tested by one of two methods determined by the IMPULSE_IFFT_TEST_METHOD parameter. If set to 0, the spot frequency method as described above will be used. If set to 1, the default, the impulse response is tested by measuring its decay to zero. The inverse FFT method is accurate if the response decays to zero and this is used to test the result.

The tolerance required for the test is set by the IMPULSE_TEST_TOL parameter.

Impulse Response Cache

The impulse response is cached to speed up subsequent simulations with the same response. The cache is located here:

C:\Users\login-name\AppData\Roaming\SIMetrix Technologies\SIMetrixvvv\LaplaceCache

where login-name is the name used to log in to your system and vvv is the version, e.g. 830 for version 8.3. This can be deleted at any time. The cache is size limited to 250Mbytes by default. This can be changed using the option variable LaplaceCacheSizeMBytes. Type 'Set LaplaceCacheSizeMBytes=nnn' at the command line to set a new value.

The cache may be disabled by setting the IMPULSE_ENABLE_CACHE parameter to 0.

Impulse Response Analytical Extraction

The impulse response analytical extraction attempts to match the transfer function to one of the following patterns. The parameters ???MATH???a???MATH???, ???MATH???b???MATH???, ???MATH???c???MATH??? etc represent constant values.

The analytic matching algorithm will attempt to match partial expressions as well as the whole expression if they are combined by addition, subtraction or multiplication. For example:

???MATH???(1/(s+1))*exp(-s)???MATH???

will be matched to the product of ???MATH???1/(s+1)???MATH??? and ???MATH???exp(-s)???MATH???. Respectively these have impulse responses of ???MATH???e^{-t}???MATH??? and ???MATH???\delta(1+t)???MATH??? and those two impulse responses will be convolved to yield the final result. Note that this only succeeds if all partial expressions can be matched analytically.

To exploit this method, ensure that the expression is entered in a manner that keeps each partial expression distinct. For example, the following would be matched correctly:

???MATH???1/(s+1)*1/(s+2)???MATH???

as this would be seen as the two partial expressions ???MATH???1/(s+1)???MATH??? and ???MATH???1/(s+2)???MATH??? each of which can be resolved analytically. The final result would be obtained by convolving the two impulse responses.

However the following mathematically identical expression would not be matched:

???MATH???1/((s+1)*(s+2))???MATH???

The analysis algorithm does not attempt to decompose denominator products so this will not be recognised analytically. The above would be extracted using inverse FFT which would usually be successful but might not be if the run time is considerably shorter than the time constant.

Run Time Error Control

During the simulation run the time step is controlled using a truncation error method. The tolerances used for this error control are set using the TRUNC_RELTOL and TRUNC_ABSTOL parameters. The truncation error test may be disabled altogether by setting the TRUNC_TEST parameter to zero.

PSpice LAPLACE and FREQ compatibility

SIMetrix is compatible with PSpice Laplace transfer functions at the netlist level entered in one of these forms

Exxx n1 n2 LAPLACE { input_expression } { laplace_expression }
Gxxx n1 n2 LAPLACE { input_expression } { laplace_expression }

Exxx n1 n2 FREQ {input_expression} [keyword1] [keyword2] table_values DELAY=delay
Gxxx n1 n2 FREQ {input_expression} [keyword1] [keyword2] table_values DELAY=delay

SIMetrix will implement the PSpice expression in the best way possible. If the Laplace expression can be expressed as a quotient of polynomials, the lumped model (see Laplace Transfer Function - Lumped Implementation) will be used instead of the convolution model as long as it does not contain any parameters.