Laplace Transfer Function

Overview

The Laplace transfer function device implements a linear device defined in the frequency domain by a Laplace transform. For example the Laplace transform ???MATH???\frac{1}{s+1}???MATH??? defines a first order low pass filter while ???MATH???exp(-s)???MATH??? defines a 1 second delay.

The SIMetrix Laplace transfer function device features two different methods of implementation, namely "lumped network" and "convolution". The lumped network implementation is fast and efficient and may be used for transfer functions that can be represented in polynomial form. The convolution method is more general purpose and can model a wide range of transfer functions. High performance is achieved by using a fast FFT-based convolution method.

User Interface

Selecting the menu Place > Analog Behavioural > Laplace Transfer Function brings up the following dialog The operation of the various controls is described below.

The Advanced Options tab displays the following:

Laplace Expression

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 in the Define Laplace Transfer Function dialog box in the Definition section - see the Delay entry. 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.

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

Parameters

Parameters (defined with .PARAM) may be used in Laplace expressions but currently only if convolution implementation is selected.

For the lumped network implementation you can use parameters in the underlying simulator device if you define the expression in terms of the numerator and denominator coefficients. Refer to the S-domain Transfer Function Block in the Simulator Reference Manual for further details.

Device type

Implementation

Lumped Network

The diagram below shows the configuration used for a third order network. This method can implement a quotient of polynomials in the s variable and any expression that can be reduced to a quotient of polynomials. For such expressions it is fast and efficient.

This method cannot implement expressions containing arbitrary functions such as square-root or exponentials. Such Laplace expressions must be implemented using the convolution method.

Convolution

This is a general purpose method and can implement a wide range of Laplace transforms.

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 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 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\SIMetrix810\LaplaceCache

where <login-name> is the name used to log in to your system. 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.

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.