It is common to think and write in time domain. creating the mathematical description of the. Continuous systems- using Laplace or s-

Size: px
Start display at page:

Download "It is common to think and write in time domain. creating the mathematical description of the. Continuous systems- using Laplace or s-"

Transcription

1

2 It is common to think and write in time domain quantities, but this is not the best thing to do in creating the mathematical description of the system we are dealing with. Continuous systems- using Laplace or s- transform to create the mathematical description of the system will simplified the process. When dealing with discrete systems, the - transform will be used to simplify the process. 2

3 The transform is a mathematical tool commonly used for the analysis and synthesis of discrete-time control systems. The transform in discrete-time systems play a similar role as the Laplace transform in continuous-time systems 3

4 In a linear discrete-time control system, a linear difference equation characteries the dynamics of the system. To determine the system's response to a given input, such a difference equation must be solved. transformation transforms linear time- invariant difference equations into algebraic equations in. 4

5 To present definitions of the transform, basic theorems associated with the transform, methods for finding the inverse transform. Discussed how to solved the difference equations by using the transform method. 5

6 Discrete-time signals arise if the system involves a sampling operation of continuous-time signals. The sampled signal is x(0,x(t, x(2t,..., T = sampling period. Such a sequence of values arising from the sampling operation is usually written as x(kt. If the system involves an iterative process carried out by a digital computer, the signal involved is a number sequence x(0, x(, x(2... The sequence of numbers is usually written as x(k, k indicates the order in which the number occurs in the sequence, for example, x(0,x(,x(2... Although x(k is a number sequence, it can be considered as a sampled signal of x(t when the sampling period T is sec. 6

7 continuous-time signal x, sampled signal x(, and the number sequence x(. 7

8 The transform of a time function x(t, where t is nonnegative, or of a sequence of values x(kt where k takes ero or positive integers and T is the sampling period, is defined by the following equation: k X( = Z [x(t] = Z [x(kt] = (2. 2 k 0 x( kt For a sequence of numbers x(k, the transform is defined by k X( = Z [x(k] = x( k (2. 3 k 0 The transform defined by Equation (2. 2 or (2. 3 is referred to as the one-sided transform The symbol Z denotes "the transform of." 8

9 In the one-sided x(t = 0 for t < 0 or x(k = 0 for k < 0. transform, we assume 9

10 The transform of x(t, where - <t<, or of x(k, where k takes integer values (k = 0, ±, ±2,, is defined by k X( = Z [x(t] = Z [x(kt] = x ( kt (2. 4 k Or X( = Z [x(k] = x ( k (2. 5 k The transform defined by Equation (2. 4 or (2.5 is referred to as the two-sided transform k 0

11 In the two-sided transform, the time function x(t is assumed to be nonero for t < 0 and the sequence x(k is considered to have nonero values for k < 0

12 Both the one-sided and two-sided transforms are series in powers of -. (The latter involves both positive and negative powers of -. In our discussion, only the one-sided transform is considered in detail. 2

13 Expansion of the right-hand side of Equation (2. 2 gives: X( = x(0 + x(t - + x(2t x(kt -k + (2. 6 3

14 4

15 Let us find the transform of the unit-step function ( t 0 t x( t 0, t 0 In sampling a unit-step function we assume that this function is continuous from the right; that is, (0 =. Then, referring to Equation (2. 2, we have X( = Z [(t] = = k k 0 k 0 k It is noted that (k as defined by ( k, k 0, k 0 0,,2,... is commonly called a unit-step sequence. 5

16 Consider the unit-ramp function x( t t, 0 0, t Notice that x(kt = kt, k = 0,,2,... t 0 Figure sampled unitramp signal. The magnitudes of the sampled values are proporti- onal to the sampling period T 6

17 transform of the unit-ramp function: X( = Z [t] = x( kt k kt k T k k k 0 k 0 k 0 T ( T ( T 2 ( 2 7

18 k Obtain the transform of x(k as defined by x( k k a, k 0,,2,... 0, k 0 a = constant 8

19 k Referring to the definition of the transform given by Equation (2. 3, we obtain X( = Z [a k ] = k 0 x( k k = + a - + a a a a k 0 a k k 9

20 Find the transform of at e, 0 x( t 0, t t 0 Since x(kt = e -akt, k = 0,, 2, we have X( = Z [e -at ] = x ( kt k 0 k e akt = + e -at - + e -2aT -2 + e -3aT -3 + e at k 0 e at k 20

21 Consider the sinusoidal function Noting that Sin t, x ( t 0, e j t = cos t + j sin t e -j t = cos t - j sin t we have 0 t t 0 Sin t j t j t 2 j ( e e Cos t j t j t ( e 2 e 2

22 Since the transform of the exponential function is Z [e -at ] = at Z [e -at ] = X( = Z [Sin t] = Z e at ( t j t j e e X( = Z [Sin t] = Z ( 2 e e j 2 e e j T j T j 2 ( ( 2 e e e e j T j T j T j T j 2 sin T sin 2 T 22 2 cos 2 T cos 2 2 T

23 EXERCISE Obtain the transform of the cosine function x t Cos t t t 23

24 Solution X( = Z [cos t] = Z ( t j t j e e X( = Z [cos t] = Z ( e e T j T j 2 e e T j T j ( 2 e e T j T j 2 ( ( 2 2 e e e e T j T j 2 cos 2 cos T T cos 2 cos 2 2 T T 24 cos 2 2 T

25 EXERCISE 2 (example 2-2 Obtain the transform of X (s s s( s Hint Whenever a function in s is given, one approach for finding the corresponding transform is to convert X(s into x(t and then find the transform of x(t. Another approach is to expand X(s into partial fractions and use a transform table to find the transforms of the expanded terms. 25

26 Solution The inverse Laplace transform of X(s is x(t = e -t, 0 t Hence, X( = Z [ e -t ] = T e ( ( e T T ( e ( ( e ( T e T 26

27 Just as in working with the Laplace transformation, a table of transforms of commonly encountered functions is very useful for solving problems in the field of discrete-time systems. Table 2- and table 2.(ctnd is such a table. 27

28 IMPORTANT PROPERTIES AND THEOREMS OF THE TRANSFORM Will discussed important properties and useful theorems of the transform. Assumption time function x(t is -transformable and x(t = 0 for t < 0. 28

29 Will discussed properties and theorems:. Multiplication by a Constant 2. Linearity of the Transform 3. Multiplication by a k 4. Shifting Theorem 5. Complex Translation Theorem 6. Initial Value Theorem 7. Final Value Theorem 29

30 Multiplication by a Constant. If X( is the transform of x(t, then Z [ax(t] = az [x(t] = ax( where a is a constant. To prove this, note that by definition Z [ax(t] = k 0 ax ( kt k a k 0 x ( kt k ax ( 30

31 Linearity of the Transform. The transform possesses an important property: linearity. This means that, if f(k and g(k are -transformable and and are scalars, then x(k formed by a linear combination x(k = f(k + g(k has the transform X( = F( + G( where F( and G( are the transforms of f(k and g(k, respectively. Prove this property!!! 3

32 Multiplication by a k If X( is the transform of x(k, then the transform of a k x(k can be given by X(a - : Z [a k x(k] = X(a - This can be proved as follows: Z [a k x(k] = k a k x( k 0 k k 0 x( k( a k = X(a - 32

33 Shifting Theorem Also call real translation theorem. If x(t = 0 for t < 0 and x(t has the transform X(, then Z [x(t-nt] = -n X( (2.8 and n k Z [x(t+nt] = X ( x ( kt (2.9 n = ero or a positive integer Prove eqn (2.8 and eqn (2.9 n k 0 33

34 Example 2-3 Find the transforms of unit-step functions that are delayed by sampling period and 4 sampling periods, respectively, as shown in figure (a and (b below 34

35 Solution Using the shifting theorem given by Equation (2. 8, we have Z [(t-t] = - Z [(t] = Also, Z [(t-4t] = -4 Z [(t] = 4 4 (Note that - represents a delay of sampling period T, regardless of the value of T. 35

36 Example 2-4 Obtain the transform of f ( a a k 0,, k,2,3,... Solution: Referring to Equation (2. 8, we have Z [x(k-] = - X( The transform of a k is Z [a k ] = a and so Z [f(a] = Z [a k- ] = where k =,2,3,... k 0 a a 36

37 Example 2-5 Consider the function y(k, which is a sum of functions x(h, where h = 0,, 2,..., k, such that y( k 0 x( h, where y(k = 0 for k < 0. Obtain the transform of y(k. h k k 0,,2,... 37

38 Solution First note that y(k = x(0 + x( + y(k- = x(0 + x( +. + x(k- + x(k. + x(k- Hence, y(k y(k- = x(k, k = 0,, 2, Therefore, Z [y(k y(k-] = Z[x(k] or Y( - Y( = X( Which yield, Y ( X ( Where, X( = Z [x(k] 38

39 Complex Translation Theorem If x(t has the transform X(, then the transform of e -at x(t can be given by X(e at. To prove Z [e -at x(t] x( kt k 0 e akt k k 0 x( kt ( e at k X ( e at Thus, we see that replacing in X( by e at gives the transform of e -at x(t. 39

40 Example 2-6 By using the complex translation theorem, obtain the transforms of: e -at sin t and e -at cos t, respectively, 40

41 Solution Noting that Z [sin t ] = sin T 2 cos T Using the complex translation theorem e sin cos T T e at at Z [ e sin t ] at 2aT 2 2e 2 4

42 Solution Z [cos t ] cos T 2 cos T 2 Using the complex translation theorem Z [e -at cos t ] at e cos T at 2 at 2e cos T e 2 42

43 Exercise By using the complex translation theorem, Obtain the transform of te -at. Hint: T Z [t] = X ( 2 ( 43

44 Initial Value Theorem. If x(t has the transform X( and if lim X ( exists, then the initial value x(0 of x(t or x(k is given by x ( 0 lim X ( 44

45 Initial Value Theorem (ctnd. The initial value theorem is convenient for checking transform calculations for possible errors. Since x(0 is usually known, a check of the initial value by lim X ( can easily spot errors in X(, if any exist. 45

46 Example 2-8 Determine the initial value x(0 if the transform of x(t is given by T ( e X ( T ( ( e By using the initial value theorem, we find T ( e x( 0 lim 0 T ( ( e Referring to Example 2-2, notice that this X( was the transform of t x ( t e and thus x(0 = 0, which agrees with the result obtained earlier. 46

47 Final Value Theorem. The final value of x(k, that is, the value of x(k as k approaches infinity, can be given by lim x( k k lim[( X ( ] (2. 27 Proved? 47

48 Example 2-9 Determine the final value x( of X (, a at e by using the final value theorem. By applying the final value theorem to the given X(, we obtain x( lim ( X ( lim ( 0 at e lim e at 48

49 In this section we have presented important properties and theorems of the transform that will prove to be useful in solving many transform problems. For the purpose of convenient reference, these important properties and theorems are summaried in Table

50 The transformation serves the same role for discrete-time control systems that the Laplace transformation serves for continuous-time control systems. The notation for the inverse transform is Z -. The inverse transform of X( yields the corresponding time sequence x(k. 50

51 An obvious method for finding the inverse transform is to refer to a transform table. However, unless we refer to an extensive transform table, we may not be able to find the inverse transform of a complicated function of. 5

52 Other methods: Direct division method Computational method Partial-fraction-expansion method Inversion integral method 52

53 In obtaining the inverse transform, we assume, as usual, that the time sequence x(kt or x(k is ero for k 0. Before we present the four methods, however, a few comments on poles and eros of the pulse transfer function are in order. 53

54 Poles and Zeros in the Plane. In engineering applications of the transform method, X( may have the form or m m b0 b... bm X ( ( m n n n a... a b0 ( ( 2...( m X ( ( m n ( p ( p...( p ( 2 pn n (2.28 where the p i s (i =,2,...,n are the poles of X( and the j s (j =,2,...,m the eros of X(. 54

55 Note that in control engineering and signal processing X( is frequently expressed as a ratio of polynomials in -, as follows: X b b ( n m ( n m 0 ( a a b a n m n n (2.29 where - is interpreted as the unit delay operator. 55

56 In finding the poles and eros of (, it is convenient to express ( as a ratio of polynomials in. For example, X 2 2 Clearly, ( has poles at = - and = -2 and eros at = 0 and =

57 If X( is written as a ratio of polynomials in -, however, the preceding X( can be written as X ( ( ( 2 Although poles at = -l and = -2 and a ero at = -0.5 are clearly seen from the expression, a ero at = 0 is not explicitly shown, and so the beginner may fail to see the existence of a ero at = 0. Therefore, in dealing with the poles and eros of X(, it is preferable to express X( as a ratio of polynomials in, rather than polynomials in - 57

58 Direct Division Method. Obtain the inverse transform by expanding ( into an infinite power series in -. Useful when it is difficult to obtain the closed-form expression for the inverse transform or it is desired to find only the first several terms of x(k. 58

59 The direct division method stems from the fact that if X( is expanded into a power series in -, that is, if or X ( k 0 x( kt k x(0 x( T x(2t 2... x( kt k... X ( k 0 x( k k x(0 x( x( x( k k... then x(kt or x(k is the coefficient of the -k term. Hence, the values of x(kt or x(k for k = 0,,2,... can be determined by inspection. 59

60 If X( is given in the form of a rational function, the expansion into an infinite power series in increasing powers of - can be accomplished by simply dividing the numerator by the denominator, where both the numerator and denominator of X( are written in increasing powers of -. If the resulting series is convergent, the coefficients of the -k term in the series are the values x(kt of the time sequence or the values of x(k of the number sequence. 60

61 Although the present method gives the values of x(0, x(t, x(2t,... or the values of x(0, x(, x(2,... in a sequential manner, it is usually difficult to obtain an expression for the general term from a set of values of x(kt or x(k. 6

62 Example 2-0 Find x(k for k = 0,,2,3,4 when X( is given by X ( 0 ( ( First, rewrite X( as a ratio of polynomials in -, as follows: X (

63 Dividing the numerator by the denominator, we have have

64 Thus, X ( By comparing this infinite series expansion of X( with X( = k 0 x( k x(0 = 0 x( = 0 x(2 = 7 x(3 = 8.4 x(4 = 8.68 k, we obtain

65 As seen from this example, the direct division method may be carried out by hand calculations if only the first several terms of the sequence are desired. In general, the method does not yield a closed-form expression for x(k 65

66 Example 2- Find x(k when X( is given by X ( Solution 66

67 Example 2-2 Obtain the inverse transform of X( = The transform X( is already in the form of a power series in -. Since X( has a finite number of terms, it corresponds to a signal of finite length. By inspection, we find x(0 = x( = 2 x(2 = 3 x(3 = 4 All other x(k values are ero. 67

68 Two computational approaches to obtain the inverse transform. MATLAB approach Difference equation approach Consider a system G( defined by G( (

69 In finding the inverse transform, we utilie the Kronecker delta function 0 (kt where, 0(kT =, for k = 0 = 0 for k 0 Assume that x(k, the input to the system G(, is the Kronecker delta input, or x(k =, for k = 0 = 0, for k 0 The transform of the Kronecker delta input is X( = 69

70 Using the Kronecker delta input, Equation (2. 30 can be rewritten as 2 Y ( G ( 2 X ( (2.3 70

71 MATLAB Approach. MATLAB can be used for finding the inverse transform. Referring to Equation (2. 3, the input X( is the transform of the Kronecker delta input. In MATLAB the Kronecker delta input is given by x = [ eros(,n] where N corresponds to the end of the discretetime duration of the process considered. 7

72 Since the transform of the Kronecker delta input X( is equal to unity, the response of the system to this input is Y ( G( Hence the inverse transform of G( is given by y(0,y(l,y(2,... Let us obtain y(k up to k =

73 To obtain the inverse transform of G( with MATLAB, we proceed as follows: Enter the numerator and denominator as follows: num = [ ] den = [ ] Enter the Kronecker delta input. x = [ eros(,40] Then enter the command y = filter(num,den,x to obtain the response y(k from k = 0 to k =

74 MATLAB Program 2- % Finding inverse transform % ***** Finding the inverse transform of C( is the same as % finding the response of the system Y(/X( = G( to the % Kronecker delta input ***** % ***** Enter the numerator and denominator of C( ***** num = [ ]; den = [ ]; % ***** Enter the Kronecker delta input x and filter command % y = filter(num,den,x ***** x=[ eros(,40]; y = filter(num,den,x 74

75 Computational Method MATLAB (ctnd If this program is executed, the screen will show the output y(k from k = 0 to 40 as follows: y = Columns through Columns 8 through Columns 5 through Columns 22 through Columns 29 through Columns 36 through

76 Computational Method-MATLAB (ctnd (Note that MATLAB computations begin from column and end at column 4, rather than from column 0 to column 40. These values give the inverse transform of G(. That is, y(0 = 0 y( = y(2 = y(3 = y(40 =

77 To plot the values of the inverse transform of G(, follow the procedure given in the following. Plotting Response to the Kronecker Delta Input. Consider the system given by Equation (2. 3. A possible MATLAB program to obtain the response of this system to the Kronecker delta input is shown in MATLAB Program 2-2. The corresponding plot is shown in Figure

78 MATLAB Program 2-2 % Response to Kronecker delta input num = [ ]; den = [ ]; x = [ eros(,40]; k = 0:40; y = filter(num,den,x; plot(k,y, o v=[ ]; axis(v; grid title ( Response to Kronecker Delta Input xlabel( k ylabel( y(k 78

79 Figure 2-2 Response of the system defined by Equation (2. 3 to the Kronecker delta input. 79

80 Difference Equation Approach Noting that Equation (2. 3 can be written as ( Y( = ( X( Convert this equation into the difference equation as follows: y(k y(k y(k = x(k x(k ---- (2. 32 where x(0 = and x(k = 0 for k 0, and y(k = 0 for k < 0. [x(k is the Kronecker delta input.] 80

81 The initial data y(0 and y( can be determined as follows: By substituting k = -2 into Equation (2.32, we find y( y(-l y(-2 = x( x(-2 from which we get y(0 = 0 Next, by substituting k = - into Equation (2. 32, we obtain y( y( y(- = x( x(- from which we get: y( =

82 To find the inverse transform of Y( - solve the following difference equation for y(k: y(k y(k y(k = x(k x(k ---- (2. 33 with the initial data y(0 = 0, y( = , x(0 =, and x(k = 0 for k 0. Equation (2. 33 can be solved easily by hand, or by use of BASIC, FORTRAN, or others. 82

83 The partial-fraction expansion method presented here, which is parallel to the partial-fraction- expansion method used in Laplace transformation, is widely used in routine problems involving transforms. The method requires that all terms in the partial fraction expansion be easily recogniable in the table of transform pairs. 83

84 To find the inverse transform, if X( has one or more eros at the origin ( = 0, then X(/ or X( is expanded into a sum of simple first or second-order terms by partial fraction expansion, and a transform table is used to find the corresponding time function of each expanded term. It is noted that the only reason that we expand X(/ into partial fractions is that each expanded term has a form that may easily be found from commonly available transform tables 84

85 Example 2-3 Before we discuss the partial-fraction-expansion method, we shall review the shifting theorem. Consider the following X(: X ( a By writing X( as Y(, we obtain X ( Y ( a 85

86 Referring to Table 2-, the inverse transform of Y( can be obtained as follows: Z - [ Y ( ] y( k k a Hence, the inverse transform of X( = - Y( is given by Z - [ X ( ] x( k y( k Since y(k is assumed to be ero for all k < 0, we have k y( k a, k,2,3,... x( k 0, k 0 86

87 Consider X( as given by m m b 0 b... b m b m X (, m n n n a... a a To expand X( into partial fractions, we first factor the denominator polynomial of X( and find the poles of X(: m m b b... bm X ( ( p ( p...( p 0 m 2 n n n b 87

88 Expand X(/ into partial fractions so that each term is easily recogniable in a table of transforms. If the shifting theorem is utilied in taking inverse transforms, however, X( instead of X(/, may be expanded into partial fractions. The inverse transform of X( is obtained as the sum of the inverse transforms of the partial fractions. 88

89 A commonly used procedure for the case where all the poles are of simple order and there is at least one ero at the origin (that is, b m = 0 is to divide both sides of X( by and then expand X(/ into partial fractions. Once X(/ is expanded, it will be of the form X ( a p a 2 p 2... a n p n 89

90 The coefficient a i can be determined by multiplying both sides of this last equation by - p i and setting = p i. This will result in ero for all the terms on the right-hand side except the a i term, in which the multiplicative factor p i has been cancelled by the denominator. Hence, we have a ( pi i X ( Note that such determination of a i is valid only for simple poles. p i 90

91 If X(/ involves a multiple pole, for example, a double pole at = p and no other poles, then X(/ will have the form X ( c c2 2 ( p ( p The coefficients c and c 2 are determined from 2 X ( c ( p It is noted that if X(/ involves a triple pole at = p, then the partial fractions must include a term ( + p /( p 3. p 9

92 Example 2-4 Given the transform at ( e X ( ( e ( at where a is a constant and T is the sampling period, determine the inverse transform x(kt by use of the partial-fraction-expansion method. 92

93 Partial-Fraction-Expansion Method (ctnd ( at B A e X ( ( ( at at at e B A e e X ( ( ( ( ( at at e e X A ( ( ( ( ( at at at e e X e B e ( ( ( ( ( at at e at e e e e e e B at at ( ( at at at at e e e e at e X ( 93 e

94 The partial fraction expansion of X( is found to be X ( at e From Table 2- we find - Z Z - e at e akt Hence, the inverse transform of X( is akt x ( kt e, k 0,,2,... 94

95 Difference equations can be solved easily by use of a digital computer, provided the numerical values of all coefficients and parameters are given. However, closed-form expressions for x(k cannot be obtained from the computer solution, except for very special cases. The usefulness of the transform method is that it enables us to obtain the closed-form expression for x(k. 95

96 Consider the linear time-invariant discrete-time system characteried by the following linear difference equation: x ( k a x ( k... a x ( k n b u( k b u( k... 0 n n b u( k n (2.34 where u(k and x(k are the system's input and output, respectively, at the kth iteration. 96

97 Let us define Z [x(k] = X( Then x( k, x( k 2, x( k 3,... and x ( k, x ( k 2, x ( k 3,... can be expressed in terms of X( and the initial conditions. Their transforms were summaried in Table

98 Table 2.3 transform of x(k+m and x(k-m 98

99 Example 2-8 Solve the following difference equation by use of the transform method: x(k x(k + + 2x(k = 0, x(0 = 0, x( = Solution First note that the transforms of x(k + 2, x(k +, and x(k are given, respectively, by Z [x(k + 2] = 2 X( - 2 x(0 - x( Z [x(k + ] = X( - x(0 Z [x(k] = X( 99

100 Taking the transforms of both sides of the given difference equation, we obtain 2 X( - 2 x(0 - x( + 3X( - 3x(0 + 2X( = 0 Substituting the initial data and simplifying gives X ( ( (

101 Noting that k Z - (, Z - 2 ( k k we have x( k ( ( 2, k 0,,2,... 2 k 0

102 Example 2-9 Obtain the solution of the following difference equation in terms of x(0 and x(: x(k (a + bx(k + + abx(k = 0 where a and b are constants and k = 0,,2,... The transform of this difference equation can be given by [ 2 X( - 2 x(0 - x(l] + (a + b[x( - x(0] + abx( = 0 02

103 Solving this last equation for X( gives [ 2 ( a b ] x (0 x ( X ( 2 ( a b ab Notice that constants a and b are the negatives of the two roots of the characteristic equation. We shall now consider separately two cases: (a a b and (b a = b. 03

104 (a For the case where a b, expanding X(/ in to partial fractions, we obtain X ( bx (0 x ( ax (0 x (, a b a a a b b from which we get bx(0 x( ax(0 x( X ( b a a a b b The inverse transform of X( gives bx (0 x ( k ax (0 x ( k x ( k ( a ( b, b a a b where k = 0,,2,.. b a b 04

105 (b For the case where a = b, the transform X( becomes 2 ( 2a x(0 x( X ( a a x(0 [ ax(0 x(] 2 a ( a Entry 8 x(0 [ ax(0 a ( a x(] 2 The inverse transform of X( gives x(k = x(0(-a k + [ax(0 + x(l]k(-a k-, a = b where k = 0,,2,... Entry 9 05

106 RECONSTRUCTING ORIGINAL SIGNALS FROM SAMPLED SIGNALS If the sampling frequency is sufficiently high compared with the highest-frequency component involved in the continuous-time signal, the amplitude characteristics of the continuous-time signal may be preserved in the envelope of the sampled signal. To reconstruct the original signal from a sampled signal, there is a certain minimum frequency that the sampling operation must satisfy. Such a minimum frequency is specified in the sampling theorem. We shall assume that a continuous-time signal x(t has a frequency spectrum as shown in Figure 2-3. This signal x(t does not contain any frequency components above radians per second. 06

107 RECONSTRUCTING ORIGINAL SIGNALS FROM SAMPLED SIGNALS (ctnd Figure 2-3 A frequency spectrum. 07

108 RECONSTRUCTING ORIGINAL SIGNALS FROM SAMPLED SIGNALS (ctnd Sampling Theorem. If s, defined as 2 /T, where T is the sampling period, is greater than 2 ; or s>2 is the highest-frequency component present in the continuous-time signal x(t, then the signal x(t can be reconstructed completely from the sampled signal x*(t. The theorem implies that if s>2 then, from the knowledge of the sampled signal, it is theoretically possible to reconstruct exactly the original continuous-time signal. 08

109 Sampling Theorem (ctnd. To show the validity of this sampling theorem, we need to find the frequency spectrum of the sampled signal x*(t. The Laplace transform of x*(t is given by X *( s T k X ( s j sk (2.35 X ( X ( s j s k T k s ( / T ln (

110 Sampling Theorem (ctnd. To obtain the frequency spectrum, we substitute j for s in Equation ( Thus, X *( j T k X ( j j s k T X ( j( s T X ( j T X ( j( s (2.37 0

111 Sampling Theorem (ctnd. Equation (2. 37 gives the frequency spectrum of the sampled signal x*(t. The frequency spectrum of the impulse-sampled signal is reproduced an infinite number of times and is attenuated by the factor /T. Thus, the process of impulse modulation of the continuous-time signal produces a series of sidebands. Since X*(s is periodic with period 2 / s or X*(s = X*(s ± j sk, s k = 0,,2,.. if a function X(s has a pole at s = s, then X*(s has poles at s = s ± j sk (k= 0,,2,...

112 Sampling Theorem (ctnd. Figure 2-4(a and (b show plots of the frequency spectra X*(j versus for two values of the sampling frequency s. Figure 2-4(a corresponds to s > 2 while Figure 2-4(b corresponds to s < 2. Each plot of X*(j versus consists of X(j /T repeated every s = 2 /T rad/sec. In the frequency spectrum of X*(j the component X(j /T is called the primary component, and the other components, X(j( ± s k /T, are called complementary components. 2

113 Figure 2-4 Plots of the frequency spectra X*(j versus for two values of sampling frequency s: (a s > 2 ; (b s< 2 3

114 Sampling Theorem (ctnd. If s>2, no two components of X*(j will overlap, and the sampled frequency spectrum will be repeated every s rad/sec. If s < 2, the original shape of X(j no longer appears in the plot of X*(j versus because of the superposition of the spectra. Therefore, we see that the continuous-time signal x(t can be reconstructed from the impulse-sampled signal x*(t by filtering if and only if s > 2 4

115 Sampling Theorem (ctnd. It is noted that although the requirement on the minimum sampling frequency is specified by the sampling theorem as s>2 where is the highest-frequency component present in the signal, practical considerations on the stability of the closed-loop system and other design considerations may make it necessary to sample at a frequency much higher than this theoretical minimum. (Frequently, s is chosen to be 0 to 20. 5

116 Ideal Low-Pass Filter. The amplitude frequency spectrum of the ideal lowpass filter G I (j is shown in Figure 2-5. The magnitude of the ideal filter is unity over the frequency range and is ero outside this frequency range. 6

117 Ideal Low-Pass Filter (ctnd. Figure 2-5 Amplitude frequency spectrum of the ideal low-pass filter. 7

118 Ideal Low-Pass Filter (ctnd. The sampling process introduces an infinite number of complementary components (sideband components in addition to the primary component. The ideal filter will attenuate all complementary components to ero and will pass only the primary component, provided the sampling frequency s is greater than twice the highest-frequency component of the continuous-time signal. Ideal filter reconstructs the continuous-time signal represented by the samples. 8

119 Ideal Low-Pass Filter (ctnd. Figure 2-6 shows the frequency spectra of the signals before and after ideal filtering. Figure 2-6 Frequency spectra of the signals before and after ideal filtering. 9

120 The frequency spectrum at the output of the ideal filter is /T times the frequency spectrum of the original continuous-time signal x(t. Since the ideal filter has constant-magnitude characteristics for the frequency region, 2 s 2 s, there is no distortion at any frequency within this frequency range. there is no phase shift in the frequency spectrum of the ideal filter. (The phase shift of the ideal filter is ero. 20

121 Ideal Low-Pass Filter (ctnd. It is noted that if the sampling frequency is less than twice the highest-frequency component of the original continuous-time signal, then because of the frequency spectrum overlap of the primary component and complementary components, even the ideal filter cannot reconstruct the original continuoustime signal. (In practice, the frequency spectrum of the continuous-time signal in a control system may extend beyond 2 s even though the amplitudes at the higher frequencies are small. 2

122 Ideal Low-Pass Filter (ctnd. Let find the impulse response function of the ideal filter. Since the frequency spectrum of the ideal filter is given by 2 s 2 s GI ( j 0 elsewhere the inverse Fourier transform of the frequency spectrum gives g I ( t sin( T s t t / 2 / 2 s (

123 Ideal Low-Pass Filter (ctnd. Equation (2. 38 gives the unit-impulse response of the ideal filter. Figure 2-7 shows a plot of g I (t versus t Notice that the response extends from t = - to t =. This implies that there is a response for t < 0 to a unit impulse applied at t = 0. (That is, the time response begins before an input is applied. This cannot be true in the physical world. Hence, such an ideal filter is physically unrealiable. 23

124 Figure 2-7 Impulse response g I (t of ideal filter. 24

125 Ideal Low-Pass Filter (ctnd. In practice ideal low-pass filter is not physically realiable. Because of the ideal filter is unrealiable and signals in practical control systems generally have higher-frequency components and are not ideally band limited, it is not possible, in practice, to exactly reconstruct a continuous-time signal from the sampled signal, no matter what sampling frequency is chosen. In other words, it is not possible to reconstruct exactly a continuous-time signal in a practical control system once it is sampled 25

126 Frequency-Response Characteristics of the Zero- Order Hold. The transfer function of a ero-order hold is Ts ( e Gh0 s (2.39 s To compare the ero-order hold with the ideal filter, obtain the frequency-response characteristics of the transfer function of the ero-order hold. 26

127 Frequency-Response Characteristics of the Zero- Order Hold (ctnd. frequency-response x-tics of the ZOH.(Figure 2-8(a Undesirable gain peaks at frequencies of 3 s /2, 5 s /2 and so on. The magnitude is more than 3 db down (0.637 = db at 2 frequency. s Because the magnitude decreases gradually as the frequency increases, the complementary components gradually attenuate to ero. Since the magnitude characteristics of the ero-order hold are not constant, if a system is connected to the sampler and eroorder hold, distortion of the frequency spectra occurs in the system. 27

128 Figure 2.8 (a Frequency-response curves for the eroorder hold; (b equivalent Bode diagram when T = sec. 28

129 Frequency-Response Characteristics of the Zero-Order Hold (ctnd. The phase-shift characteristics of the ZOH Note that sin ( T/2 alternates positive and negative values as increases from 0 to s, s to 2 s, 2 s to 3 s, and so on. Thus, the phase curve [Figure 2-8(a, bottom] is discontinuous at = k s = 2 k/t where k =,2,3,... Such a discontinuity or a switch from a positive value to a negative value, or vice versa, may be considered to be a phase shift of ±80. In Figure 2-8(a, phase shift is assumed to be

130 Frequency-Response Characteristics of the Zero- Order Hold (ctnd. Summary: The frequency spectrum of the output of the ZOH includes complementary components, since the magnitude characteristics show that the magnitude of G h0 (j is not ero for, except at points 2 s where = ± s, = ± 2 s, = ± 3 s,. In the phase curve there are phase discontinuities of ±80 at frequency points that are multiples of s. Except for these phase discontinuities, the phase characteristic is linear in. 30

131 Figure 2-9 shows the comparison of the ideal filter and the ero-order hold Figure 2-9 Comparison of the ideal filter and the ero-order hold 3

132 Folding. The phenomenon of the overlap in the frequency spectra is known as folding. Figure 2-20 shows the regions where folding error occurs. 2 The frequency s is called the folding frequency or Nyquist frequency N. That is, N 2 s T 32

133 Figure 2-20 Diagram showing the regions where folding error occurs 33

134 Folding(ctnd. In practice, signals in control systems have high- frequency components, and some folding effect will almost always exist. For example, in an electromechanical system some signal may be contaminated by noises. The frequency spectrum of the signal, therefore, may include low- frequency components as well as high-frequency noise components (that is, noises at 60 or 400 H. Since sampling at frequencies higher than 400 H is not practical, the high frequency will be folded in and will appear as a low frequency. Remember that all signals with frequencies higher than 2 s appear as signals of frequencies between 0 and. 2 s 34

135 Aliasing. In the frequency spectra of an impulse-sampled signal x* (t, where s < 2 as shown in Figure 2-2, consider an arbitrary frequency point 2 that falls in the region of the overlap of the frequency spectra. The frequency spectrum at = 2 comprises two components, X*(j 2 and X*(j( s- 2. The latter component comes from the frequency spectrum centered at = s. 35

136 Aliasing (ctnd. Figure 2-2 Frequency spectra of an impulsesampled signal x*(t. 36

137 Aliasing (ctnd. Thus, the frequency spectrum of the sampled signal at = 2 includes components not only at frequency 2 but also at frequency s - 2 (in general, at s ± 2, where n is an integer. When the composite spectrum is filtered by a lowpass filter, such as a ero-order hold, some higher harmonics will still be present in the output. The frequency component at = n s ± 2 (where n is an integer will appear in the output as if it were a frequency component at = 2 It is not possible to distinguish the frequency spectrum at = 2 from that at = n s ± 2. 37

138 Aliasing (ctnd. As shown in Figure 2-2, the phenomenon that the frequency component s - 2 (in general, n s ± 2, where n is an integer shows up at frequency 2 when the signal x(t is sampled is called aliasing. This frequency s - 2 (in general, n s ± 2 is called an alias of 2. It is important to remember that the sampled signals are the same if the two frequencies differ by an integral multiple of the sampling frequency s. s If a signal is sampled at a slow frequency such that the sampling theorem is not satisfied, then high frequencies are "folded in" and appear as low frequencies. 38

139 Aliasing (ctnd. To avoid aliasing, choose the sampling frequency high enough s > 2, where is the highest-frequency component present in the signal or use a prefilter ahead of the sampler to reshape the frequency spectrum of the signal (so that the frequency spectrum for 2 s is negligible before the signal is sampled. 39

140 MAPPING BETWEEN THE s PLANE AND THE PLANE The absolute stability and relative stability of the linear time-invariant continuous-time closed-loop control system are determined by the locations of the closed-loop poles in the s plane. For example, complex closed-loop poles in the left half of the s plane near the j axis will exhibit oscillatory behavior, and closed-loop poles on the negative real axis will exhibit exponential decay. 40

141 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd Complex variables and s relationship = e Ts pole and ero locations in the plane are related to the pole and ero locations in the s plane. The stability of the linear time-invariant discrete-time closed-loop system can be determined in terms of the locations of the poles of the closed-loop pulse transfer function. 4

142 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd Dynamic behavior of the discrete-time control system depends on the sampling period T. Locations of poles and eros in the plane depend on the sampling period T. Therefore, a change in the sampling period T modifies the pole and ero locations in the plane and causes the response behavior to change. 42

143 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd Mapping of the Left Half of the s Plane into the Plane. In the design of a continuous-time control system, the locations of the poles and eros in the s plane are very important in predicting the dynamic behavior of the system. In designing discrete-time control systems, the locations of the poles and eros in the plane are very important. Will investigate how the locations of the poles and eros in the s plane compare with the locations of the poles and eros in the plane. 43

144 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd Mapping of the Left Half of the s Plane into the Plane (ctnd complex variable s has real part and imaginary part, we have s = + j and e T ( j T jt T j( T 2 k e e e e 44

145 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd Mapping of the Left Half of the s Plane into the Plane (ctnd Since is negative in the left half of the s plane, the left half of the s plane corresponds to = et < The j axis in the s plane corresponds to =. That is, the imaginary axis in the s plane (the line = 0 corresponds to the unit circle in the plane, The left half of the s plane corresponds to the interior of the unit circle 45

146 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd Conclusion. The j axis in the s plane corresponds to = 2. The left half of the s plane corresponds to the interior of the unit circle 3. The right half of the s plane corresponds to the exterior of the unit circle 46

147 MAPPING BETWEEN THE s PLANE AND THE PLANE (ctnd j Unit circle s plane Figure : Mapping s plane plane plane 47

148 48

149 OPEN-LOOP DISCRETE-TIME SYSTEMS The Relationship Between E( and E*(s -The -transform of the number sequence e(k was defined as: Z [e(k] = E( = e(0 + e( - + e(2-2 + (4- -The starred transform for the time function e(t was defined as: E*(s = e(0 + e(te -Ts + e(2te -2Ts +... (4-2 -If e(k of (4- is equal to e(kt of (4-2, and if e st = in (4-2, then (4-2 becomes the -transform. Hence in this case E ( E *( s e st...(4.3 49

150 The Pulse Transfer Function To develop an expression for the -transform of the output of open-loop sampled-data systems. Consider the open-loop system of Figure 4-la, Gp(s - plant transfer fn. G(s - product of the plant transfer fn and the ero-order hold transfer fn, that is Ts e G( s Gp ( s s Hence this system can be represented like in Figure 4-lb. ZOH Figure 4- Open-loop sampled-data system 50

151 The Pulse Transfer Function (cont In Figure 4-, C(s = G(sE*(s.(4.5 Assume that c(t is continuous at all sampling instants. From defination, C * ( s C( s jn s [ G( s E * ( s]*...(4.6 T n Thus, from (4-5 and (4-6, C * ( s G( s jn s E * ( s jn s T n...(4.7 From property of E*(s, E*(s+jn s = E*(s 5

152 The Pulse Transfer Function (cont Thus (4-7 becomes C *( s E *( s G( s jn s E *( s G *( s T n...(4.8 and then from (4-3, C( = E(G(.(4-9 G( is called the pulse transfer function and is the transfer function between the sampled input and the output at the sampling instants. 52

153 The derivation above is completely general. Thus given any function that can be expressed as A (s = B(sF*(s... (4-0 where F*(s must be expressible as F*(s = f 0 + f e -Ts + f 2 e -2Ts +. then, from the preceding development, A* (s = B*(sF*(s. (4- Hence, from (4-3, A( = B(F(.(4-2 where B(s is a function of s and F*(s is a function of e Ts ; that is, in F*(s, s appears only in the form e Ts. Then, in (4-2, B( = Z [B(s], F( = F*(s Ts e = (4-3 The following examples illustrate this procedure. 53

154 Example 4.2 Find the -transform of e A ( s s( s Solution Ts From (4-0, we consider : s( s ( e Ts and B( s s(s( s F *( s e Ts F ( F *( s e Ts 54

155 The Pulse Transfer Function (cont Then, from the -transform tables, The Pulse Transfer Function (cont Then, from the -transform tables, ( ( T e Z B ( ( ( ( ( T e e s s Z B and T T T T T T e e e e F B A ( ( ( ( ( ( e e ( ( 55

156 Example 4.3 Given the system shown in Figure 4-2, with input e(t a unit step function, determine the output function C(. Figure

157 Solution C( s Ts e E *( s s(s( s G( s E *( s In Example 4.2 it was shown that G( Z e s( s Ts In addition, from the table E ( Z ( t Thus C( ENTRANCE E 8 G( E( ( ( e T e e T T ( e T e T 57

158 and the inverse -transform of this function yields c(nt = e -nt Entrance 8 This response is plotted in Figure 4-2b. 58

159 The Pulse Transfer Function (cont Investigate open-loop systems of other configurations in Fig 4.3. Figure 4.3 Open-loop sampled-data systems 59

160 The Pulse Transfer Function (cont Figure 4-3a - two plants, both G (s and G 2 (s contain the transfer functions of the data holds. Then, C( s G2( s A*( s C( G2( A( A ( s G ( s E *( s A( G ( E( from (4.5 and (4.6 (4.5 (4.6 C ( G ( G 2 ( E ( (4.7 and the total transfer function is the product of the pulse transfer functions. 60

161 The Pulse Transfer Function (cont Figure 4-3b - G 2 (s not contain a data-hold transfer function. Then C( s G ( s G2( s E *( s C( GG 2( E( where GG 2( Z G ( s G2( s (4.8 bar above a product term indicates that the product must be performed in the s-domain before the - transform is taken. GG 2( G ( G2( ( 4.9 6

162 The Pulse Transfer Function (cont For the system of Figure 4-3c, Thus C( s G2( s A* ( s G2( s GE *( s C( G2( GE ( (

163 63

164 Stability Analysis of a Closed-Loop System. Will discuss the stability of linear time-invariant single-input-single-output discrete-time control systems. Consider the following closed-loop pulse-transfer function system: C( R( G( GH ( (

165 Stability Analysis of a Closed-Loop System (ctnd The stability of the system defined by Equation (2. 42, as well as of other types of discrete-time control systems, may be determined from the locations of the closed- loop poles in the plane, or the roots of the characteristic equation as follows: P( = + GH( = 0 65

166 Stability Analysis of a Closed-Loop System (ctnd. For the system to be stable, the closed-loop poles or the roots of the characteristic equation must lie within the unit circle in the plane. Any closed- loop pole outside the unit circle makes the system unstable. 2. If a simple pole lies at =, then the system becomes critically stable. Also, the system becomes critically stable if a single pair of conjugate complex poles lies on the unit circle in the plane. Any multiple closed-loop pole on the unit circle makes the system unstable. 3. Closed-loop eros do not affect the absolute stability and therefore maybe located anywhere in the plane. 66

167 Stability Analysis of a Closed-Loop System (ctnd Thus, a linear time-invariant single-input-single- output discrete-time closed-loop control system becomes unstable if any of the closed-loop poles lies outside the unit circle and/or any multiple closed-loop pole lies on the unit circle in the plane. 67

168 Stability Analysis of a Closed-Loop System (ctnd Example 4-2 Consider the closed-loop control system shown in Figure Determine the stability of the system when K =. The open-loop transfer function G(s of the system is G( s s e s s( s 68

169 Figure 2-34 Closed-loop control system of Example

170 Stability Analysis of a Closed-Loop System (ctnd Solution The transform of G(s is G( = Z e s s s( s ( ( (2.44 Since the closed-loop pulse transfer function for the system is C( G( R ( G ( the characteristic equation is + G( = 0 70

171 Stability Analysis of a Closed-Loop System (ctnd which becomes ( ( = 0 Or = 0 The roots of the characteristic equation are found to be = j0.68, 2 = 0.5 j0.68 Since = 2 < the system is stable. 7

172 Methods for Testing Absolute Stability. Three stability tests can be applied directly to the characteristic equation P( = 0 without solving for the roots. Schur-Cohn stability test Jury stability test. Routh stability criterion Will only discussed Jury stability test 72

173 The Jury Stability Test. In applying the Jury stability test to characteristic equation P( = 0, we construct a table whose elements are based on the coefficients of P(. Assume that the characteristic equation P( is a polynomial in as follows: n n P( a a... a n 0 n a (2.45 where a 0 > 0. Then the Jury table becomes as given in Table

174 Table 2. 4 General Form Of The Jury Stability Table 74

175 First row consist of the coefficients in P( arranged First row consist of the coefficients in P( arranged in the ascending order of powers of. Second row consist of the coefficients of P( arranged in the descending order of powers of. The elements for rows 3 through 2n - 3 are given by the following determinants: 0,,2,...,, 0 n k a a a a b k k n n k 2 0,,2,...,, 0 2 n k b b b b c k k n n k p p b b k k 75 0,,2, k p p p p q k k k

176 Stability Criterion by the Jury Test. A system with the characteristic equation P( = 0 given by Equation (2. 44, rewritten as n n P( a0 a... an an where a 0 > 0, is stable if the following conditions are all satisfied:. a n < a 0 2. P( 0 ( 0 for n even P ( 3. 0 for n odd 76

177 b n b 4. n 0 c n 2 c 0 q 2 q 0 77

178 Example 4-3 Construct the Jury stability table for the following characteristic equation: P ( a0 a a2 a3 a4 where a 0 > 0. Write the stability conditions. 78

179 Solution Referring to the general case of the Jury stability table given by Table 2. 4, a Jury stability table for the fourth-order system may be constructed as shown in Table This table is slightly modified from the standard form and is convenient for the computations of the b's and c's. The determinant given in middle of each row gives the value of b or c written on the right-hand side of the same row. 79

180 Solution The stability conditions are as follows:. a 4 < a 0 2. P( a0 a a2 a3 a P( a0 a a2 a3 a4 4. b 3 b 0 0, n 4 even c 2 c 0 It is noted that the value of c (or, in the case of the nth-order system, the value of q is not used in the stability test, and therefore the computation of c (or q may be omitted. 80

181 Table 2. 5 Jury Stability Table For The Fourth-Order System 8

182 Example 4-4 Examine the stability of the following characteristic equation: P( Solution Notice that for this characteristic equation a 0 a a a 2 3 a

183 Solution (ctnd Examine st condition for stability, a 4 < a 0, < satisfied. 2nd condition for stability P ( 0 : ( P(l = = 0.09 > 0 satisfied. 3rd condition for stability (P(- > 0, n= even P(- = =.89 > 0, satisfied. 83

184 Solution (ctnd 4 th condition, construct the Jury stability table. Referring to Example 4-3, we compute the values of b 3, b 2, b, and b 0 and c 2 and c 0. The result is shown in Table 2.6. From table 2.6, we get b 3 = > = b 0 c 2 = > 0.35 = c 0 Thus both parts of the fourth condition given in Example 4-3 are satisfied. Since all conditions for stability are satisfied, the given characteristic equation is stable, or all roots lie inside the unit circle in the plane. 84

185 Table 2. 6 JURY STABILITY TABLE FOR THE SYSTEM OF EXAMPLE

186 Solution (ctnd As a matter of fact, the given characteristic equation P( can be factored as follows: P( = ( - 0.8( + 0.5( - 0.5( As a matter of course, the result obtained above agrees with the fact that all roots are within the unit circle in the plane. 86

187 87

188 This document was created with Win2PDF available at The unregistered version of Win2PDF is for evaluation or non-commercial use only. This page will not be added after purchasing Win2PDF.

Let H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )

Let H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 ) Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:

More information

C(s) R(s) 1 C(s) C(s) C(s) = s - T. Ts + 1 = 1 s - 1. s + (1 T) Taking the inverse Laplace transform of Equation (5 2), we obtain

C(s) R(s) 1 C(s) C(s) C(s) = s - T. Ts + 1 = 1 s - 1. s + (1 T) Taking the inverse Laplace transform of Equation (5 2), we obtain analyses of the step response, ramp response, and impulse response of the second-order systems are presented. Section 5 4 discusses the transient-response analysis of higherorder systems. Section 5 5 gives

More information

Use: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)}

Use: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)} 1 VI. Z Transform Ch 24 Use: Analysis of systems, simple convolution, shorthand for e jw, stability. A. Definition: X(z) = x(n) z z - transforms Motivation easier to write Or Note if X(z) = Z {x(n)} z

More information

A system that is both linear and time-invariant is called linear time-invariant (LTI).

A system that is both linear and time-invariant is called linear time-invariant (LTI). The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped

More information

Frequency Response and Continuous-time Fourier Series

Frequency Response and Continuous-time Fourier Series Frequency Response and Continuous-time Fourier Series Recall course objectives Main Course Objective: Fundamentals of systems/signals interaction (we d like to understand how systems transform or affect

More information

Lecture 27 Frequency Response 2

Lecture 27 Frequency Response 2 Lecture 27 Frequency Response 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/12 1 Application of Ideal Filters Suppose we can generate a square wave with a fundamental period

More information

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:

More information

Chapter 7. Digital Control Systems

Chapter 7. Digital Control Systems Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,

More information

Automatic Control Systems theory overview (discrete time systems)

Automatic Control Systems theory overview (discrete time systems) Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations

More information

Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain)

Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain) 1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear systems (frequency domain) 2 Motivations Consider an LTI system Thanks to the Lagrange s formula we can compute the motion of

More information

3.2 Complex Sinusoids and Frequency Response of LTI Systems

3.2 Complex Sinusoids and Frequency Response of LTI Systems 3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex

More information

Laplace Transform Analysis of Signals and Systems

Laplace Transform Analysis of Signals and Systems Laplace Transform Analysis of Signals and Systems Transfer Functions Transfer functions of CT systems can be found from analysis of Differential Equations Block Diagrams Circuit Diagrams 5/10/04 M. J.

More information

z Transform System Analysis

z Transform System Analysis z Transform System Analysis Block Diagrams and Transfer Functions Just as with continuous-time systems, discrete-time systems are conveniently described by block diagrams and transfer functions can be

More information

Methods for analysis and control of. Lecture 6: Introduction to digital control

Methods for analysis and control of. Lecture 6: Introduction to digital control Methods for analysis and of Lecture 6: to digital O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr www.lag.ensieg.inpg.fr/sename 6th May 2009 Outline Some interesting books:

More information

How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches

How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches Difference Equations (an LTI system) x[n]: input, y[n]: output That is, building a system that maes use of the current and previous

More information

Recursive, Infinite Impulse Response (IIR) Digital Filters:

Recursive, Infinite Impulse Response (IIR) Digital Filters: Recursive, Infinite Impulse Response (IIR) Digital Filters: Filters defined by Laplace Domain transfer functions (analog devices) can be easily converted to Z domain transfer functions (digital, sampled

More information

Digital Signal Processing

Digital Signal Processing COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even

More information

STABILITY ANALYSIS TECHNIQUES

STABILITY ANALYSIS TECHNIQUES ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response

More information

INTRODUCTION TO DIGITAL CONTROL

INTRODUCTION TO DIGITAL CONTROL ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant

More information

Chapter 2 z-transform. X(z) = Z[x(t)] = Z[x(kT)] = Z[x(k)] x(kt)z k = x(kt)z k = x(k)z k. X(z)z k 1 dz 2πj c = 1

Chapter 2 z-transform. X(z) = Z[x(t)] = Z[x(kT)] = Z[x(k)] x(kt)z k = x(kt)z k = x(k)z k. X(z)z k 1 dz 2πj c = 1 One-sided -transform Two-sided -transform Chapter 2 -TRANSFORM X() Z[x(t)] Z[x(kT)] Z[x(k)] x(kt) k x(k) k X() k0 k x(kt) k k0 k x(k) k Note that X() x(0) + x(t ) + x(2t ) 2 + + x(kt) k + Inverse -transform

More information

ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization

ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)

More information

Chapter 6: The Laplace Transform. Chih-Wei Liu

Chapter 6: The Laplace Transform. Chih-Wei Liu Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace

More information

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27 1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system

More information

Identification Methods for Structural Systems

Identification Methods for Structural Systems Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from

More information

ECE 301 Fall 2010 Division 2 Homework 10 Solutions. { 1, if 2n t < 2n + 1, for any integer n, x(t) = 0, if 2n 1 t < 2n, for any integer n.

ECE 301 Fall 2010 Division 2 Homework 10 Solutions. { 1, if 2n t < 2n + 1, for any integer n, x(t) = 0, if 2n 1 t < 2n, for any integer n. ECE 3 Fall Division Homework Solutions Problem. Reconstruction of a continuous-time signal from its samples. Consider the following periodic signal, depicted below: {, if n t < n +, for any integer n,

More information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014 Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can

More information

Department of Electronics and Instrumentation Engineering M. E- CONTROL AND INSTRUMENTATION ENGINEERING CL7101 CONTROL SYSTEM DESIGN Unit I- BASICS AND ROOT-LOCUS DESIGN PART-A (2 marks) 1. What are the

More information

Z-Transform. x (n) Sampler

Z-Transform. x (n) Sampler Chapter Two A- Discrete Time Signals: The discrete time signal x(n) is obtained by taking samples of the analog signal xa (t) every Ts seconds as shown in Figure below. Analog signal Discrete time signal

More information

UNIT-II Z-TRANSFORM. This expression is also called a one sided z-transform. This non causal sequence produces positive powers of z in X (z).

UNIT-II Z-TRANSFORM. This expression is also called a one sided z-transform. This non causal sequence produces positive powers of z in X (z). Page no: 1 UNIT-II Z-TRANSFORM The Z-Transform The direct -transform, properties of the -transform, rational -transforms, inversion of the transform, analysis of linear time-invariant systems in the -

More information

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},

More information

Solution of ODEs using Laplace Transforms. Process Dynamics and Control

Solution of ODEs using Laplace Transforms. Process Dynamics and Control Solution of ODEs using Laplace Transforms Process Dynamics and Control 1 Linear ODEs For linear ODEs, we can solve without integrating by using Laplace transforms Integrate out time and transform to Laplace

More information

ESS Finite Impulse Response Filters and the Z-transform

ESS Finite Impulse Response Filters and the Z-transform 9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of

More information

Chapter 5 THE APPLICATION OF THE Z TRANSFORM. 5.3 Stability

Chapter 5 THE APPLICATION OF THE Z TRANSFORM. 5.3 Stability Chapter 5 THE APPLICATION OF THE Z TRANSFORM 5.3 Stability Copyright c 2005- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org February 13, 2008 Frame # 1 Slide # 1 A. Antoniou Digital Signal

More information

EE Experiment 11 The Laplace Transform and Control System Characteristics

EE Experiment 11 The Laplace Transform and Control System Characteristics EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal

More information

EE451/551: Digital Control. Chapter 3: Modeling of Digital Control Systems

EE451/551: Digital Control. Chapter 3: Modeling of Digital Control Systems EE451/551: Digital Control Chapter 3: Modeling of Digital Control Systems Common Digital Control Configurations AsnotedinCh1 commondigitalcontrolconfigurations As noted in Ch 1, common digital control

More information

Laplace Transforms Chapter 3

Laplace Transforms Chapter 3 Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important

More information

Section 7.4: Inverse Laplace Transform

Section 7.4: Inverse Laplace Transform Section 74: Inverse Laplace Transform A natural question to ask about any function is whether it has an inverse function We now ask this question about the Laplace transform: given a function F (s), will

More information

V. IIR Digital Filters

V. IIR Digital Filters Digital Signal Processing 5 March 5, V. IIR Digital Filters (Deleted in 7 Syllabus). (dded in 7 Syllabus). 7 Syllabus: nalog filter approximations Butterworth and Chebyshev, Design of IIR digital filters

More information

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation A. Hably 1 1 Gipsa-lab, Grenoble-INP ahmad.hably@grenoble-inp.fr Commande d un robot mobile (Gipsa-lab (DA)) ASI 1 / 25 Outline 1 2 (Gipsa-lab (DA)) ASI 2 / 25 of controllers Signals must be sampled and

More information

E2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)

E2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2) E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,

More information

Some of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e

Some of the different forms of a signal, obtained by transformations, are shown in the figure. jwt e z. jwt z e Transform methods Some of the different forms of a signal, obtained by transformations, are shown in the figure. X(s) X(t) L - L F - F jw s s jw X(jw) X*(t) F - F X*(jw) jwt e z jwt z e X(nT) Z - Z X(z)

More information

The Laplace Transform

The Laplace Transform The Laplace Transform Syllabus ECE 316, Spring 2015 Final Grades Homework (6 problems per week): 25% Exams (midterm and final): 50% (25:25) Random Quiz: 25% Textbook M. Roberts, Signals and Systems, 2nd

More information

Root Locus Methods. The root locus procedure

Root Locus Methods. The root locus procedure Root Locus Methods Design of a position control system using the root locus method Design of a phase lag compensator using the root locus method The root locus procedure To determine the value of the gain

More information

Representing a Signal

Representing a Signal The Fourier Series Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity and timeinvariance of the system and represents the

More information

The Discrete-Time Fourier

The Discrete-Time Fourier Chapter 3 The Discrete-Time Fourier Transform 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 3-1-1 Continuous-Time Fourier Transform Definition The CTFT of

More information

37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.

37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40. 28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t

More information

ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin

ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency

More information

Chapter Intended Learning Outcomes: (i) Understanding the relationship between transform and the Fourier transform for discrete-time signals

Chapter Intended Learning Outcomes: (i) Understanding the relationship between transform and the Fourier transform for discrete-time signals z Transform Chapter Intended Learning Outcomes: (i) Understanding the relationship between transform and the Fourier transform for discrete-time signals (ii) Understanding the characteristics and properties

More information

(i) Represent discrete-time signals using transform. (ii) Understand the relationship between transform and discrete-time Fourier transform

(i) Represent discrete-time signals using transform. (ii) Understand the relationship between transform and discrete-time Fourier transform z Transform Chapter Intended Learning Outcomes: (i) Represent discrete-time signals using transform (ii) Understand the relationship between transform and discrete-time Fourier transform (iii) Understand

More information

Digital Control Systems

Digital Control Systems Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist

More information

The Laplace Transform

The Laplace Transform C H A P T E R 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive forcing terms. For such problems the methods described

More information

MAS107 Control Theory Exam Solutions 2008

MAS107 Control Theory Exam Solutions 2008 MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve

More information

Control of Electromechanical Systems

Control of Electromechanical Systems Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance

More information

Discrete-Time David Johns and Ken Martin University of Toronto

Discrete-Time David Johns and Ken Martin University of Toronto Discrete-Time David Johns and Ken Martin University of Toronto (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) University of Toronto 1 of 40 Overview of Some Signal Spectra x c () t st () x s () t xn

More information

Quadrature-Mirror Filter Bank

Quadrature-Mirror Filter Bank Quadrature-Mirror Filter Bank In many applications, a discrete-time signal x[n] is split into a number of subband signals { v k [ n]} by means of an analysis filter bank The subband signals are then processed

More information

The Continuous-time Fourier

The Continuous-time Fourier The Continuous-time Fourier Transform Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline Representation of Aperiodic signals:

More information

Sampling and Discrete Time. Discrete-Time Signal Description. Sinusoids. Sampling and Discrete Time. Sinusoids An Aperiodic Sinusoid.

Sampling and Discrete Time. Discrete-Time Signal Description. Sinusoids. Sampling and Discrete Time. Sinusoids An Aperiodic Sinusoid. Sampling and Discrete Time Discrete-Time Signal Description Sampling is the acquisition of the values of a continuous-time signal at discrete points in time. x t discrete-time signal. ( ) is a continuous-time

More information

EE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL

EE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL EE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL EE 213 Fall 2009 LABORATORY #1 INTRODUCTION TO MATLAB INTRODUCTION The purpose of this laboratory is to introduce you to Matlab and to illustrate some of its circuit

More information

X. Chen More on Sampling

X. Chen More on Sampling X. Chen More on Sampling 9 More on Sampling 9.1 Notations denotes the sampling time in second. Ω s = 2π/ and Ω s /2 are, respectively, the sampling frequency and Nyquist frequency in rad/sec. Ω and ω denote,

More information

EC6405 - CONTROL SYSTEM ENGINEERING Questions and Answers Unit - I Control System Modeling Two marks 1. What is control system? A system consists of a number of components connected together to perform

More information

8.5 Taylor Polynomials and Taylor Series

8.5 Taylor Polynomials and Taylor Series 8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:

More information

Filter Banks II. Prof. Dr.-Ing. G. Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany

Filter Banks II. Prof. Dr.-Ing. G. Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany Filter Banks II Prof. Dr.-Ing. G. Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany Page Modulated Filter Banks Extending the DCT The DCT IV transform can be seen as modulated

More information

Design of IIR filters

Design of IIR filters Design of IIR filters Standard methods of design of digital infinite impulse response (IIR) filters usually consist of three steps, namely: 1 design of a continuous-time (CT) prototype low-pass filter;

More information

Multimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2

Multimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December

More information

Analysis of Discrete-Time Systems

Analysis of Discrete-Time Systems TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time

More information

E : Lecture 1 Introduction

E : Lecture 1 Introduction E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation

More information

EC Signals and Systems

EC Signals and Systems UNIT I CLASSIFICATION OF SIGNALS AND SYSTEMS Continuous time signals (CT signals), discrete time signals (DT signals) Step, Ramp, Pulse, Impulse, Exponential 1. Define Unit Impulse Signal [M/J 1], [M/J

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information

UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH *

UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH * 4.4 UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH 19 Discussion Problems 59. Two roots of a cubic auxiliary equation with real coeffi cients are m 1 1 and m i. What is the corresponding homogeneous

More information

GATE EE Topic wise Questions SIGNALS & SYSTEMS

GATE EE Topic wise Questions SIGNALS & SYSTEMS www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)

More information

Dynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.

Dynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology. Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control

More information

7.17. Determine the z-transform and ROC for the following time signals: Sketch the ROC, poles, and zeros in the z-plane. X(z) = x[n]z n.

7.17. Determine the z-transform and ROC for the following time signals: Sketch the ROC, poles, and zeros in the z-plane. X(z) = x[n]z n. Solutions to Additional Problems 7.7. Determine the -transform and ROC for the following time signals: Sketch the ROC, poles, and eros in the -plane. (a) x[n] δ[n k], k > 0 X() x[n] n n k, 0 Im k multiple

More information

MEM 355 Performance Enhancement of Dynamical Systems

MEM 355 Performance Enhancement of Dynamical Systems MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Intro Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /5/27 Outline Closed Loop Transfer

More information

The z-transform Part 2

The z-transform Part 2 http://faculty.kfupm.edu.sa/ee/muqaibel/ The z-transform Part 2 Dr. Ali Hussein Muqaibel The material to be covered in this lecture is as follows: Properties of the z-transform Linearity Initial and final

More information

Lecture 28 Continuous-Time Fourier Transform 2

Lecture 28 Continuous-Time Fourier Transform 2 Lecture 28 Continuous-Time Fourier Transform 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/14 1 Limit of the Fourier Series Rewrite (11.9) and (11.10) as As, the fundamental

More information

# FIR. [ ] = b k. # [ ]x[ n " k] [ ] = h k. x[ n] = Ae j" e j# ˆ n Complex exponential input. [ ]Ae j" e j ˆ. ˆ )Ae j# e j ˆ. y n. y n.

# FIR. [ ] = b k. # [ ]x[ n  k] [ ] = h k. x[ n] = Ae j e j# ˆ n Complex exponential input. [ ]Ae j e j ˆ. ˆ )Ae j# e j ˆ. y n. y n. [ ] = h k M [ ] = b k x[ n " k] FIR k= M [ ]x[ n " k] convolution k= x[ n] = Ae j" e j ˆ n Complex exponential input [ ] = h k M % k= [ ]Ae j" e j ˆ % M = ' h[ k]e " j ˆ & k= k = H (" ˆ )Ae j e j ˆ ( )

More information

Lecture 3 - Design of Digital Filters

Lecture 3 - Design of Digital Filters Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can

More information

UNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?

UNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter? UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics

More information

Module 1: Signals & System

Module 1: Signals & System Module 1: Signals & System Lecture 6: Basic Signals in Detail Basic Signals in detail We now introduce formally some of the basic signals namely 1) The Unit Impulse function. 2) The Unit Step function

More information

Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN

Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN Pusat Pengajian Kejuruteraan Komputer Dan Perhubungan Universiti Malaysia Perlis Discrete-Time

More information

Stability of Feedback Control Systems: Absolute and Relative

Stability of Feedback Control Systems: Absolute and Relative Stability of Feedback Control Systems: Absolute and Relative Dr. Kevin Craig Greenheck Chair in Engineering Design & Professor of Mechanical Engineering Marquette University Stability: Absolute and Relative

More information

1 Introduction & Objective

1 Introduction & Objective Signal Processing First Lab 13: Numerical Evaluation of Fourier Series Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in

More information

Chapter 3 Data Acquisition and Manipulation

Chapter 3 Data Acquisition and Manipulation 1 Chapter 3 Data Acquisition and Manipulation In this chapter we introduce z transf orm, or the discrete Laplace Transform, to solve linear recursions. Section 3.1 z-transform Given a data stream x {x

More information

Table of Laplacetransform

Table of Laplacetransform Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun

More information

DIGITAL CONTROL OF POWER CONVERTERS. 3 Digital controller design

DIGITAL CONTROL OF POWER CONVERTERS. 3 Digital controller design DIGITAL CONTROL OF POWER CONVERTERS 3 Digital controller design Frequency response of discrete systems H(z) Properties: z e j T s 1 DC Gain z=1 H(1)=DC 2 Periodic nature j Ts z e jt e s cos( jt ) j sin(

More information

Need for transformation?

Need for transformation? Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations

More information

Theory and Problems of Signals and Systems

Theory and Problems of Signals and Systems SCHAUM'S OUTLINES OF Theory and Problems of Signals and Systems HWEI P. HSU is Professor of Electrical Engineering at Fairleigh Dickinson University. He received his B.S. from National Taiwan University

More information

Experiment 13 Poles and zeros in the z plane: IIR systems

Experiment 13 Poles and zeros in the z plane: IIR systems Experiment 13 Poles and zeros in the z plane: IIR systems Achievements in this experiment You will be able to interpret the poles and zeros of the transfer function of discrete-time filters to visualize

More information

CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation

CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is

More information

Lecture 11. Frequency Response in Discrete Time Control Systems

Lecture 11. Frequency Response in Discrete Time Control Systems EE42 - Discrete Time Systems Spring 28 Lecturer: Asst. Prof. M. Mert Ankarali Lecture.. Frequency Response in Discrete Time Control Systems Let s assume u[k], y[k], and G(z) represents the input, output,

More information

APPLICATIONS FOR ROBOTICS

APPLICATIONS FOR ROBOTICS Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table

More information

Chapter 13 Digital Control

Chapter 13 Digital Control Chapter 13 Digital Control Chapter 12 was concerned with building models for systems acting under digital control. We next turn to the question of control itself. Topics to be covered include: why one

More information

SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB

SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB INTRODUCTION SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB Laplace transform pairs are very useful tools for solving ordinary differential equations. Most

More information

Notes for ECE-320. Winter by R. Throne

Notes for ECE-320. Winter by R. Throne Notes for ECE-3 Winter 4-5 by R. Throne Contents Table of Laplace Transforms 5 Laplace Transform Review 6. Poles and Zeros.................................... 6. Proper and Strictly Proper Transfer Functions...................

More information

CMPT 889: Lecture 3 Fundamentals of Digital Audio, Discrete-Time Signals

CMPT 889: Lecture 3 Fundamentals of Digital Audio, Discrete-Time Signals CMPT 889: Lecture 3 Fundamentals of Digital Audio, Discrete-Time Signals Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 6, 2005 1 Sound Sound waves are longitudinal

More information

CONTROL OF DIGITAL SYSTEMS

CONTROL OF DIGITAL SYSTEMS AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it

More information

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,

More information

Numeric Matlab for Laplace Transforms

Numeric Matlab for Laplace Transforms EE 213 Spring 2008 LABORATORY # 4 Numeric Matlab for Laplace Transforms Objective: The objective of this laboratory is to introduce some numeric Matlab commands that are useful with Laplace transforms

More information

Filter Analysis and Design

Filter Analysis and Design Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved

More information

Analysis of Discrete-Time Systems

Analysis of Discrete-Time Systems TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time

More information