Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 905791, 9 pages
Research Article

Transformation Matrix for Time Discretization Based on Tustin’s Method

College of Precision Instrument & Optoelectronics Engineering, Tianjin University, Tianjin 300072, China

Received 16 April 2014; Revised 27 June 2014; Accepted 10 July 2014; Published 18 August 2014

Academic Editor: Ivanka Stamova

Copyright © 2014 Yiming Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


This paper studies rules in transformation of transfer function through time discretization. A method of using transformation matrix to realize bilinear transform (also known as Tustin’s method) is presented. This method can be described as the conversion between the coefficients of transfer functions, which are expressed as transform by certain matrix. For a polynomial of degree n, the corresponding transformation matrix of order n exists and is unique. Furthermore, the transformation matrix can be decomposed into an upper triangular matrix multiplied with another lower triangular matrix. And both have obvious regularity. The proposed method can achieve rapid bilinear transform used in automatic design of digital filter. The result of numerical simulation verifies the correctness of the theoretical results. Moreover, it also can be extended to other similar problems. Example in the last throws light on this point.

1. Introduction

In this paper we proposed a new method of using transformation matrix to realize bilinear transform. We describe the principles and characteristics of the method in Section 2. Also we have examinations and simulations on the method in Section 3.

In Section 2, Section 2.1 starts with exploration of how to express derivative of a function with discrete points under the method of bilinear transform. By studying the rules of that conversion, Section 2.2 reveals regularity in the arithmetic of bilinear transform operated on transfer function, so that it can be replaced by systematic mathematical derivations. Finally, the transformation matrix to realize Tustin transform is obtained by recursion. Also in the process, some inner rules of the Tustin transform can be revealed.

In Section 3, Section 3.1 provides a simple application to verify the mathematical correctness of the method and uses the obtained transformation matrix as an example. Section 3.2 explains how it works when the parameters of functions are uncertain. Section 3.3 shows an example about its expanded application as it can be used to construct transformation matrixes in other situations. The purpose is to enrich the method and to make it easier to understand. Of course they are all simplified.

1.1. About Tustin’s Method

Tustin’s method is also called the bilinear transform. As a convenient tool for time discretization, it is widely used in designing digital filters to obtain expected discrete time system with desired response characteristics. The transform, as a method to convert domain function to domain function, is a first-order approximation of the -to- mapping. Under the Laplace transform, it can be regarded as a discrete-time approximation [1]: Equivalently, where is the numerical integration step size of the trapezoidal rule used in the bilinear transform derivation.

Compared with transform based on the pulse response method, the Tustin transform eliminates the aliasing effect generated by overlapping of frequency spectrum, relying on the linear transformation from domain to domain as a one-to-one mapping. Consequently, it is suitable for the design of high-pass and band-stop filters. Furthermore, because of the simple algebraic relationship between and , the transfer function of digital filter can be obtained conveniently by directing the -to- mapping into the corresponding transfer function of analog filter [13].

1.2. Problems to Be Solved

Because the bilinear transform is a kind of fractional linear transformation, the amount of calculation increases sharply with the increase of the order of system. The problem of complex calculation appears when operating bilinear transform for high-order system. Particularly, in many situations the parameters of system are changeable; thus then traditional method fails to work. For these reasons, programs to realize automatic operation are required. An efficient way is to express the transform with matrix.

Many scholars have proposed methods to obtain the matrix for bilinear transform, of which the literatures [46] were to derive coefficient matrix of different order under the discussion of relationships between the matrix elements, while [7, 8] used method of induction to derive coefficient matrix of order n from coefficient matrix of order zero. Other relative reference as [914] also studied quick algorithms for different forms of linear -to- transform.

The purpose of this paper is to introduce an easier way to obtain the transformation matrix. The method applies when the order of transfer function is large, using computer as a direct operation, or when the system function is uncertain, using the method to get the dynamic differential equation.

2. Derivations

Bilinear transform is a bridge between domain and domain. Laplace transform can change a time domain function with t into a frequency domain function with , and it has a special quality that it can deform differential equations into algebraic equations to simplify mathematical operations. transform is similar to Laplace transform, but it is used to deal with discrete series of . It can change a time domain function with into domain. The result of every mapping from domain to domain is replacement of continuous function by discrete points. Accordingly, when processing a continuous function, the result of mapping is expression of relations between numerical values, integration, and differentiation by discrete points. Because transfer functions can be converted to differential equations, expressions of differential forms by discrete points should be studied.

2.1. Definition for Derivatives with Discrete Points

In order to discuss distinctions between deformations to differential forms before and after -to- mapping by bilinear transform, such relations below will be necessary to deduct.

Suppose are values at sampling points of discrete digital signal ( is an integer less than or equal to ). If the sampling step size is , then , .

If represents first-order derivative with discrete values , of two adjacent points, such recursion can be obtained as represents the second-order derivative with three, represents the third-order derivative with four, and so on. Here we define where and are -dimensional row and column vectors and is equal the number of samples. Let ( are nonnegative integers, and ):

Thus is the expression of the th order derivative for continuous time signal with values. In addition, each sampling point can be expressed by the linear combination of as below: That is,

To facilitate the expression, let

is an upper triangular matrix of and satisfies the following properties:

Extensionally, if is defined as an -dimensional column vector by moving forward the elements of for steps, we get another property (10) as

For example, when , , where, , , , , and , , and .

2.2. Principles and Steps to Obtain Transformation Matrix

With the definitions above, the transformation matrix for bilinear transform will be presented and the relations between coefficients before and after the transform will be deducted.

Since the digital signal is not continuous, is used as the th order derivative of noncontinuous signal . For a general transfer function,

is equal to the differential equation as follows: where, and are polynomials of , and . The general method of using bilinear transformation is substituting the -to- mapping (2) into the transfer function. The domain transfer function can be shown as follows:

According to the definition of transform, (14) is equal to the difference equation as follows:

Here define , , , and . Combined with the relation between and obtained by the transfer function , equation of and can be obtained. The following are deduced on the relationship between and , , and .

Consider the simple case when , , and ; then

According to the Laplace transform, (16) is equal to (17):

Take input sequence from signal and output sequence from signal, and the step size is . Meanwhile, the average values are and weighed by .

Replace the first-order derivative of by ; according to (17), equals , as which can be expressed as

Thus it can be concluded that when taking input sequence and output sequence into related operations, as a kind of approximation, takes the position of ( and ), and takes the position of the first-order derivative of y. The mapping can be shown as below

Take input sequence and output sequence sequentially. According to (17),

Subtract (22) from (21) and then multiply by ; we get as

Because is the expression of second-order derivative, is the average first-order derivative weighed by . Similar to (17), (24) can be expressed by the differential equation as follows:

After the Laplace transform, (26) can be obtained which can be deduced by (16):

Add (22) to (21) and we get as

Because is the expression of first-order derivative , is the average value of weighed by . Equation (28) can be expressed by the differential equation as (16).

Thus it can be concluded that, when taking input sequence and output sequence into related operations, takes the position of , takes the position of the first-order derivative of , and takes the position of the second-order derivative of . The mapping can be shown as follows:

Now consider when , , and . The corresponding differential equation is

With conclusions above, taking input sequence and output sequence , (30) transforms into the following:

The coefficients of and in (31) are the same as coefficients in the difference equation after bilinear transform.

By this method, take input sequence and output sequence with points, and the mapping can be expressed as

Such rules can be used to deal with a general transfer function (12), for (13) can be expressed as

As a conclusion, for a discrete digital signal whose interval is , in order to get its th order derivative, at least values are needed to take into operations. If the mth order derivative of is expressed by values, such mapping can be shown as follows:

As changes, the mapping shown by (34) satisfies the distribution rule of binomial coefficients [4, 5], shown in Figure 1.

Figure 1: Using values to express th order derivative of .

The tree structure in Figure 1 can be expressed by Figure 2, in which is the element in th row and th column. Meanwhile, (34) can also be expressed as .

Figure 2: Relations between elements of Figure 1 [4].

For example, according to that theory, while taking from signal , according to (34), represents the th order derivative of ; represents the th order derivative of (). And according to (10), such equation exists as follows:

For example, when , ,

Back to the deduction on the relationship between and , and and , first multiply both sides by ; referring to the conclusions above, the right side of (12) can be transformed by such mapping:

According to (34), the equation can be shown as follows: where is lower triangular matrix of ; is defined as are -dimensional column vectors, and where , is partitioned matrix of and is zero matrix of .

According to (9), the right side of (15) transforms as follows:

Obviously is the coefficient vector of . Since (13) equals (15) and (38) equals (41), consequently,

In the same way,

According to (8), , the coefficient of (15) can be expressed as follows: where and express the coefficients of the numerator and denominator polynomials of . The transformation matrix for bilinear transformation is obtained by multiplying and , , and exist and are unique when is set to a fixed value.

3. Examples

3.1. An Example of Digital Filter Design

For example, the design of a 4-order Butterworth low-pass filter with cut-off frequency .

Set  s; then the cut-off frequency of the corresponding analog filter can be obtained [4, 5]:

According to 4-order Butterworth normalized system function [2],

Take into ; thus So and ; according to (41), where, according to (7), (39), and , For , each column can be calculated by such rules: where Take and , and into (47); and can be obtained as Thus,

could be expressed by difference equation:

The amplitude-frequency response of Butterworth digital low-pass filter obtained by Tustin transform is shown in Figure 3.

Figure 3: Simulation of amplitude frequency response.
3.2. An Example of Computing -Domain Transfer Function

Reference [15] mentioned design of switch capacitor circuit by Tustin transform. The -domain transfer function of first-order and second-order are as follows: where , , , , and are real.

In order to get -domain transfer function, operate by Tustin transformation as

then operate by Tustin transformation as

Thus, and are the -domain transfer functions.

3.3. An Example of Extensional Application

Reference [9] presented a matrix method for a biquadratic transformation. By the algorithm recommended in this paper, we get the matrix with little deformation.

For an th order discrete-time polynomial where is an -dimensional column vector, .

Take the transformation form (61) into : then transformed into a continuous-time polynomial as where is an 2-dimensional column vector, .

The relationship between coefficients in (60) and (62) can be expressed as where is the transformation matrix. can be obtained by this method.

Equal to (61), the transformation form can be deformed as

Take (64) into , then transforms into as follows:

According to (44), the relation between and can be expressed as (66); here consider the step as :

Last, take the equation between and into , as follows:

Similar to the deduction of (44), the relation between and can be expressed as where is a matrix:

Finally, can be expressed as

while : thus can be obtained by (70), which is the same as shown in [8]:

4. Conclusions

The method of using corresponding matrix for Tustin transform representation reflects the inherent properties; on the other hand, it allows computer programs to deal with the transfer functions of order automatically. In this paper, the transformation matrix is obtained by multiplying and , which dramatically reduces the complexity of the program without any tabulation. Furthermore, this method can also be modified properly to suit similar problems about coefficients after linear transform.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


  1. Z. Y. Wu, Digital Signal Processing, China Higher Education Press, 2004.
  2. Z. Y. Song, “Design of digital filter using bilinear Z transform,” China Computer & Network, vol. 4, no. 1, pp. 33–66, 1978. View at Google Scholar
  3. V. Havu and J. Malinen, “The Cayley transform as a time discretization scheme,” Numerical Functional Analysis and Optimization, vol. 28, no. 7-8, pp. 825–851, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  4. Q. Y. Shen, X. Y. Miao, and A. L. Song, “A transformation matrix algorithm for computing bilinear transformation,” Information and Control, vol. 11, no. 4, pp. 25–30, 1982. View at Google Scholar
  5. B. Pšenička, F. García-Ugalde, and A. Herrera-Camacho, “The bilinear Z transform by Pascal matrix and its application in the design of digital filters,” IEEE Signal Processing Letters, vol. 9, no. 11, pp. 368–370, 2002. View at Publisher · View at Google Scholar · View at Scopus
  6. S. Chivapreecha, S. Sriyapong, S. Junnapiya, and K. Dejhan, “Bilinear s-z with frequency transformation using pascal matrix operation,” in Proceedings of the International Symposium on Communications and Information Technologies (ISCIT '05), pp. 764–767, October 2005. View at Publisher · View at Google Scholar · View at Scopus
  7. G. L. Liu and P. K. Li, “A simple recursive method for computing bilinear,” Journal of North China Institute of Electric Power, vol. 14, no. 1, pp. 75–82, 1988. View at Google Scholar
  8. M. Jiang and D. Z. Cao, “A simple algorithm for bilinear transforming,” Journal of Southeast University, vol. 21, no. 3, pp. 124–127, 1991. View at Google Scholar
  9. D. Kim and Y. Choo, “Computing biquadratic transformation matrix,” Journal of Institute of Control, Robotics and Systems, vol. 15, no. 2, pp. 128–131, 2009. View at Publisher · View at Google Scholar · View at Scopus
  10. C. Zhu and Y. Zou, “Improved recursive algorithm for fractional-order system solution based on PSE and Tustin transform,” Systems Engineering and Electronics, vol. 31, no. 11, pp. 2736–2741, 2009. View at Google Scholar · View at Scopus
  11. G. Zhang, X. Chen, and T. Chen, “Digital redesign via the generalised bilinear transformation,” International Journal of Control, vol. 82, no. 4, pp. 741–754, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus
  12. S. Yan, N. Shiratori, and L. Xu, “Simple state-space formulations of 2-D frequency transformation and double bilinear transformation,” Multidimensional Systems and Signal Processing, vol. 21, no. 1, pp. 3–23, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  13. S. Yan, N. Shiratori, H. Shieh, and L. Xu, “A general state-space representation of n-variable bilinear transformation,” Signal Processing, vol. 91, no. 2, pp. 185–190, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus
  14. P. Karimaghaee and N. Noroozi, “Frequency weighted discrete-time controller order reduction using bilinear transformation,” Journal of Electrical Engineering, vol. 62, no. 1, pp. 44–48, 2011. View at Google Scholar · View at Scopus
  15. X. J. Yang, “Design of switched-capacitor one-order and two-order sections using bilinear transformation,” Acta Electronica Sinica, vol. 21, no. 4, pp. 31–40, 1982. View at Google Scholar