Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition

Copyright 2012 Orchard Publications. All rights reserved. Printed in the United States of America. No part of this
publication may be reproduced or distributed in any form or by any means, or stored in a data base or retrieval system,
without the prior written permission of the publisher.

Direct all inquiries to Orchard Publications, [email protected]

Product and corporate names are trademarks or registered trademarks of the Microsoft Corporation and The
MathWorks Inc. They are used only for identification and explanation, without intent to infringe.

Library of Congress Cataloging-in-Publication Data

Catalog record is available from the Library of Congress

Library of Congress Control Number: 2011939678

ISBN13: 9781934404232
ISBN10: 1934404233

Copyright TXu 1778061

 Preface

This text contains a comprehensive discussion on continuous and discrete time signals and systems
with many MATLAB and several Simulink examples. It is written for junior and senior
electrical and computer engineering students, and for selfstudy by working professionals. The
prerequisites are a basic course in differential and integral calculus, and basic electric circuit theory.
This book can be used in a twoquarter, or one semester course. This author has taught the subject
material for many years and was able to cover all material in 16 weeks, with 2 lecture hours per
week.
To get the most out of this text, it is highly recommended that Appendix A is thoroughly reviewed.
This appendix serves as an introduction to MATLAB, and is intended for those who are not
familiar with it. The Student Edition of MATLAB is an inexpensive, and yet a very powerful
software package; it can be found in many college bookstores, or can be obtained directly from
The MathWorks Inc., 3 Apple Hill Drive, Natick, MA 017602098
Phone: 508 6477000, Fax: 508 6477001
http://www.mathworks.com
email: [email protected]
The elementary signals are reviewed in Chapter 1, and several examples are given. The purpose of
this chapter is to enable the reader to express any waveform in terms of the unit step function, and
subsequently the derivation of the Laplace transform of it. Chapters 2 through 4 are devoted to
Laplace transformation and circuit analysis using this transform. Chapter 5 is an introduction to
statespace and contains many illustrative examples. Chapter 6 discusses the impulse response.
Chapters 7 and 8 are devoted to Fourier series and transform respectively. Chapter 9 introduces
discretetime signals and the Z transform. Considerable time was spent on Chapter 10 to present
the Discrete Fourier transform and FFT with the simplest possible explanations. Chapter 11
contains a thorough discussion to analog and digital filters analysis and design procedures. As
mentioned above, Appendix A is an introduction to MATLAB. Appendix B is an introduction to
Simulink, Appendix C contains a review of complex numbers, Appendix D is an introduction to
matrix theory, Appendix E contains a comprehensive discussion on window functions, Appendix F
describes the cross correlation and autocorrelation functions, and Appendix G presents an example
of a nonlinear system and derives its describing function.
New to the Fifth Edition
The most notable change is the addition of Appendixes F and G. All chapters and appendixes are
rewritten and the MATLAB scripts and Simulink models are based on Release R2011b (MATLAB
Version 7.13, Simulink Version 7.8.)
The author wishes to express his gratitude to the staff of The MathWorks, the developers of
MATLAB and Simulink, especially to The MathWorks Book Program Team, for the
encouragement and unlimited support they have provided me with during the production of this
and all other texts by this publisher.
Our heartfelt thanks also to Ms. Sally Wright, P.E., of Renewable Energy Research Laboratory
University of Massachusetts, Amherst, for bringing some errors and suggestions on a previous
edition to our attention.

Orchard Publications
www.orchardpublications.com
[email protected]

2
Table of Contents
1 Elementary Signals 11
1.1 Signals Described in Math Form ............................................................................... 11
1.2 The Unit Step Function ............................................................................................ 12
1.3 The Unit Ramp Function.......................................................................................... 19
1.4 The Delta Function.................................................................................................. 111
1.4.1 The Sampling Property of the Delta Function .............................................. 111
1.4.2 The Sifting Property of the Delta Function .................................................. 112
1.5 Higher Order Delta Functions ................................................................................ 113
1.6 Summary .................................................................................................................. 123
1.7 Exercises ................................................................................................................... 124
1.8 Solutions to EndofChapter Exercises.................................................................. 125
MATLAB Computing
Pages 119 through 122
Simulink Modeling
Page 117

2 The Laplace Transformation 21

2.1 Definition of the Laplace Transformation ................................................................21
2.2 Properties and Theorems of the Laplace Transform.................................................22
2.2.1 Linearity Property ..........................................................................................23
2.2.2 Time Shifting Property ...................................................................................23
2.2.3 Frequency Shifting Property ..........................................................................24
2.2.4 Scaling Property .............................................................................................24
2.2.5 Differentiation in Time Domain Property ....................................................24
2.2.6 Differentiation in Complex Frequency Domain Property............................26
2.2.7 Integration in Time Domain Property ..........................................................26
2.2.8 Integration in Complex Frequency Domain Property ..................................27
2.2.9 Time Periodicity Property ..............................................................................28
2.2.10 Initial Value Theorem ...................................................................................29
2.2.11 Final Value Theorem ...................................................................................210
2.2.12 Convolution in Time Domain Property......................................................211
2.2.13 Convolution in Complex Frequency Domain Property .............................211
2.3 The Laplace Transform of Common Functions of Time........................................212
2.3.1 The Laplace Transform of the Unit Step Function u 0 ( t ) ...........................212
2.3.2 The Laplace Transform of the Ramp Function u 1 ( t ) .................................212
2.3.3 The Laplace Transform of t n u 0 ( t ) ...............................................................214
2.3.4 The Laplace Transform of the Delta Function ( t ) ....................................217
2.3.5 The Laplace Transform of the Delayed Delta Function ( t a ) .................217

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition i
Copyright Orchard Publications
 2.3.6 The Laplace Transform of e at u0 ( t ) ............................................................ 218
at
2.3.7 The Laplace Transform of t n e u 0 ( t ) ..........................................................218
2.3.8 The Laplace Transform of sin t u 0 t ............................................................219
2.3.9 The Laplace Transform of cos t u 0 t .......................................................... 219
2.3.10 The Laplace Transform of e at sin t u0 ( t ) .................................................. 220
2.3.11 The Laplace Transform of e at cos t u 0 ( t ) .................................................. 220
2.4 The Laplace Transform of Common Waveforms ................................................... 221
2.4.1 The Laplace Transform of a Pulse................................................................. 221
2.4.2 The Laplace Transform of a Linear Segment................................................ 222
2.4.3 The Laplace Transform of a Triangular Waveform ...................................... 223
2.4.4 The Laplace Transform of a Rectangular Periodic Waveform ..................... 224
2.4.5 The Laplace Transform of a HalfRectified Sine Waveform ....................... 225
2.5 Using MATLAB for Finding the Laplace Transforms of Time Functions ............. 226
2.6 Summary................................................................................................................... 227
2.7 Exercises ................................................................................................................... 229
The Laplace Transform of a Sawtooth Periodic Waveform ................................. 230
The Laplace Transform of a FullRectified Sine Waveform................................ 230
2.8 Solutions to EndofChapter Exercises .................................................................. 231

3 The Inverse Laplace Transform 31

3.1 The Inverse Laplace Transform Integral.....................................................................31
3.2 Partial Fraction Expansion..........................................................................................31
3.2.1 Distinct Poles....................................................................................................32
3.2.2 Complex Poles ..................................................................................................35
3.2.3 Multiple (Repeated) Poles ................................................................................37
3.3 Case where F(s) is Improper Rational Function.......................................................312
3.4 Alternate Method of Partial Fraction Expansion .....................................................314
3.5 Summary ...................................................................................................................317
3.6 Exercises ....................................................................................................................319
3.7 Solutions to EndofChapter Exercises...................................................................320
MATLAB Computing
Pages 33, 34, 35, 36, 37, 310, 311, 312, 313, 320

4 Circuit Analysis with Laplace Transforms 41

4.1 Circuit Transformation from Time to Complex Frequency......................................41
4.1.1 Resistive Network Transformation ..................................................................41
4.1.2 Inductive Network Transformation .................................................................41
4.1.3 Capacitive Network Transformation ...............................................................41
4.2 Complex Impedance Z(s)............................................................................................48
4.3 Complex Admittance Y(s) ........................................................................................411
4.4 Transfer Functions....................................................................................................413

ii Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 4.5 Using the Simulink Transfer Fcn Block .................................................................. 417
4.6 Summary................................................................................................................... 420
4.7 Exercises.................................................................................................................... 421
4.8 Solutions to EndofChapter Exercises .................................................................. 424
MATLAB Computing
Pages 46, 48, 412, 416, 417, 418, 426, 427, 428, 429, 434
Simulink Modeling
Page 418

5 State Variables and State Equations 51

5.1 Expressing Differential Equations in State Equation Form ....................................51
5.2 Solution of Single State Equations...........................................................................56
5.3 The State Transition Matrix .....................................................................................58
5.4 Computation of the State Transition Matrix.........................................................510
5.4.1 Distinct Eigenvalues ....................................................................................511
5.4.2 Multiple (Repeated) Eigenvalues .................................................................515
5.5 Eigenvectors ............................................................................................................517
5.6 Circuit Analysis with State Variables .....................................................................521
5.7 Relationship between State Equations and Laplace Transform ............................528
5.8 Summary .................................................................................................................536
5.9 Exercises ..................................................................................................................539
5.10 Solutions to EndofChapter Exercises.................................................................541
MATLAB Computing
Pages 514, 517, 525, 534, 545, 546, 548
Simulink Modeling
Pages 526, 535, 543

6 The Impulse Response and Convolution 61

6.1 The Impulse Response in Time Domain ..................................................................61
6.2 Even and Odd Functions of Time.............................................................................64
6.3 Convolution...............................................................................................................66
6.4 Graphical Evaluation of the Convolution Integral ...................................................68
6.5 Circuit Analysis with the Convolution Integral ......................................................617
6.6 Summary ..................................................................................................................620
6.7 Exercises ...................................................................................................................621
6.8 Solutions to EndofChapter Exercises..................................................................623
MATLAB Applications
Pages 611, 614, 615, 628

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition iii
Copyright Orchard Publications
7 Fourier Series 71
7.1 Wave Analysis.............................................................................................................71
7.2 Evaluation of the Coefficients....................................................................................72
7.3 Symmetry in Trigonometric Fourier Series................................................................76
7.3.1 Symmetry in Square Waveform........................................................................78
7.3.2 Symmetry in Square Waveform with Ordinate Axis Shifted...........................78
7.3.3 Symmetry in Sawtooth Waveform....................................................................79
7.3.4 Symmetry in Triangular Waveform..................................................................79
7.3.5 Symmetry in Fundamental, Second, and Third Harmonics............................79
7.4 Trigonometric Form of Fourier Series for Common Waveforms...........................710
7.4.1 Trigonometric Fourier Series for Square Waveform .....................................710
7.4.2 Trigonometric Fourier Series for Sawtooth Waveform .................................714
7.4.3 Trigonometric Fourier Series for Triangular Waveform ...............................716
7.4.4 Trigonometric Fourier Series for HalfWave Rectifier Waveform ...............717
7.4.5 Trigonometric Fourier Series for FullWave Rectifier Waveform ................720
7.5 Gibbs Phenomenon .................................................................................................723
7.6 Alternate Forms of the Trigonometric Fourier Series .............................................724
7.7 Circuit Analysis with Trigonometric Fourier Series .............................................. 727
7.8 The Exponential Form of the Fourier Series ......................................................... 730
7.9 Symmetry in Exponential Fourier Series................................................................ 732
7.9.1 Even Functions ............................................................................................. 732
7.9.2 Odd Functions.............................................................................................. 733
7.9.3 Half-Wave Symmetry .................................................................................... 733
7.9.4 No Symmetry ................................................................................................ 733
7.9.5 Relation of C n to C n .................................................................................. 733
7.10 Line Spectra ............................................................................................................ 735
7.11 Computation of RMS Values from Fourier Series ................................................ 739
7.12 Computation of Average Power from Fourier Series............................................. 741
7.13 Evaluation of Fourier Coefficients Using Excel.................................................. 743
7.14 Evaluation of Fourier Coefficients Using MATLAB .......................................... 744
7.15 Summary ................................................................................................................. 748
7.16 Exercises .................................................................................................................. 751
7.17 Solutions to EndofChapter Exercises................................................................. 753
MATLAB Computing
Pages 737, 745
Simulink Modeling
Page 730

8 The Fourier Transform 81

8.1 Definition and Special Forms .................................................................................. 81
8.2 Special Forms of the Fourier Transform.................................................................. 82

iv Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 8.2.1 Real Time Functions.....................................................................................83
8.2.2 Imaginary Time Functions............................................................................86
8.3 Properties and Theorems of the Fourier Transform ................................................89
8.3.1 Linearity.........................................................................................................89
8.3.2 Symmetry .......................................................................................................89
8.3.3 Time Scaling ................................................................................................810
8.3.4 Time Shifting...............................................................................................810
8.3.5 Frequency Shifting ......................................................................................811
8.3.6 Time Differentiation ...................................................................................812
8.3.7 Frequency Differentiation...........................................................................812
8.3.8 Time Integration .........................................................................................813
8.3.9 Conjugate Time and Frequency Functions ................................................813
8.3.10 Time Convolution.......................................................................................814
8.3.11 Frequency Convolution ..............................................................................814
8.3.12 Area Under f ( t ) ...........................................................................................815
8.3.13 Area Under F ( ) ........................................................................................815
8.3.14 Parsevals Theorem......................................................................................815
8.4 Fourier Transform Pairs of Common Functions ...................................................816
8.4.1 The Delta Function Pair .............................................................................816
8.4.2 The Constant Function Pair.......................................................................818
8.4.3 The Cosine Function Pair...........................................................................818
8.4.4 The Sine Function Pair................................................................................819
8.4.5 The Signum Function Pair ..........................................................................819
8.4.6 The Unit Step Function Pair .......................................................................821
j 0 t
8.4.7 The e u0 ( t ) Function Pair.......................................................................823
8.4.8 The ( cos 0 t ) ( u 0 t ) Function Pair..................................................................823
8.4.9 The ( sin 0 t ) ( u 0 t ) Function Pair ..................................................................823
8.5 Derivation of the Fourier Transform from the Laplace Transform ......................824
8.6 Fourier Transforms of Common Waveforms........................................................826
8.6.1 The Transform of f ( t ) = A [ u 0 ( t + T ) u0 ( t T ) ] .........................................826
8.6.2 The Transform of f ( t ) = A [ u 0 ( t ) u 0 ( t 2T ) ] .............................................827
8.6.3 The Transform of f ( t ) = A [ u 0 ( t + T ) + u 0 ( t ) u 0 ( t T ) u 0 ( t 2T ) ] ............828
8.6.4 The Transform of f ( t ) = A cos 0 t [ u0 ( t + T ) u 0 ( t T ) ] ...............................829
8.6.5 The Transform of a Periodic Time Function with Period T ......................829

8.6.6 The Transform of the Periodic Time Function f ( t ) = A
n =
( t nT ) .....830

8.7 Using MATLAB for Finding the Fourier Transform of Time Functions .............831
8.8 The System Function and Applications to Circuit Analysis..................................834

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition v
Copyright Orchard Publications
 8.9 Summary ................................................................................................................ 840
8.10 Exercises ................................................................................................................. 845
8.11 Solutions to EndofChapter Exercises................................................................ 847
MATLAB Computing
Pages 831, 832, 848, 852, 853, 854, 855, 857, 858

9 DiscreteTime Systems and the Z Transform 91

9.1 Definition and Special Forms of the Z Transform.................................................91
9.2 Properties and Theorems of the Z Transform ........................................................93
9.2.1 Linearity ........................................................................................................93
9.2.2 Shift of f [ n ]u 0 [ n ] in the DiscreteTime Domain .......................................93
9.2.3 Right Shift in the DiscreteTime Domain...................................................94
9.2.4 Left Shift in the DiscreteTime Domain .....................................................95
n
9.2.5 Multiplication by a in the DiscreteTime Domain...................................96
naT
9.2.6 Multiplication by e in the DiscreteTime Domain..............................96
9.2.7 Multiplication by n and n2 in the DiscreteTime Domain .......................96
9.2.8 Summation in the DiscreteTime Domain .................................................97
9.2.9 Convolution in the DiscreteTime Domain................................................98
9.2.10 Convolution in the DiscreteFrequency Domain .......................................99
9.2.11 Initial Value Theorem ..................................................................................99
9.2.12 Final Value Theorem................................................................................... 99
9.3 The Z Transform of Common DiscreteTime Functions....................................911
9.3.1 The Transform of the Geometric Sequence ..............................................912
9.3.2 The Transform of the DiscreteTime Unit Step Function .......................914
9.3.3 The Transform of the DiscreteTime Exponential Sequence ...................915
9.3.4 The Transform of the DiscreteTime Cosine and Sine Functions ...........915
9.3.5 The Transform of the DiscreteTime Unit Ramp Function.....................916
9.4 Computation of the Z Transform with Contour Integration .............................917
9.5 Transformation Between s and zDomains ........................................................920
9.6 The Inverse Z Transform......................................................................................923
9.6.1 Partial Fraction Expansion .......................................................................923
9.6.2 The Inversion Integral ..............................................................................929
9.6.3 Long Division of Polynomials..................................................................933
9.7 The Transfer Function of DiscreteTime Systems...............................................935
9.8 State Equations for DiscreteTime Systems .........................................................940
9.9 Summary ...............................................................................................................944
9.10 Exercises ................................................................................................................949
9.11 Solutions to EndofChapter Exercises...............................................................951
MATLAB Computing
Pages 933, 934, 937, 938, 943, 955, 957

vi Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Simulink Modeling
Page 940
Excel Plots
Pages 932, 940

10 The DFT and the FFT Algorithm 101

10.1 The Discrete Fourier Transform (DFT) ...............................................................101
10.2 Even and Odd Properties of the DFT..................................................................109
10.3 Common Properties and Theorems of the DFT ............................................... 1010
10.3.1 Linearity .................................................................................................1011
10.3.2 Time Shift ..............................................................................................1011
10.3.3 Frequency Shift......................................................................................1012
10.3.4 Time Convolution .................................................................................1012
10.3.5 Frequency Convolution.........................................................................1013
10.4 The Sampling Theorem .....................................................................................1013
10.5 Number of Operations Required to Compute the DFT...................................1016
10.6 The Fast Fourier Transform (FFT).....................................................................1017
10.7 Summary.............................................................................................................1028
10.8 Exercises..............................................................................................................1031
10.9 Solutions to EndofChapter Exercises ............................................................1033
MATLAB Computing
Pages 105, 107, 1034
Simulink Modeling
Pages 105, 106
Excel Analysis ToolPak
Pages 107, 108

11 Analog and Digital Filters

11.1 Filter Types and Classifications............................................................................111
11.2 Basic Analog Filters ..............................................................................................112
11.2.1 RC LowPass Filter..................................................................................112
11.2.2 RC HighPass Filter ................................................................................114
11.2.3 RLC BandPass Filter.............................................................................. 116
11.2.4 RLC BandElimination Filter ................................................................. 118
11.3 LowPass Analog Filter Prototypes .................................................................... 1110
11.3.1 Butterworth Analog LowPass Filter Design......................................... 1113
11.3.2 Chebyshev Type I Analog LowPass Filter Design ............................... 1124
11.3.3 Chebyshev Type II Analog LowPass Filter Design .............................. 1136
11.3.4 Elliptic Analog LowPass Filter Design................................................. 1137
11.4 HighPass, BandPass, and BandElimination Filter Design .......................... 1139

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition vii
Copyright Orchard Publications
 11.5 Digital Filters ......................................................................................................1149
11.6 Digital Filter Design with Simulink ...................................................................1167
11.6.1 The Direct Form I Realization of a Digital Filter ..................................1167
11.6.2 The Direct Form II Realization of a Digital Filter.................................1168
11.6.3 The Series Form Realization of a Digital Filter .....................................1170
11.6.4 The Parallel Form Realization of a Digital Filter...................................1172
11.6.5 The Digital Filter Design Block .............................................................1174
11.7 Summary .............................................................................................................1183
11.8 Exercises..............................................................................................................1187
11.9 Solutions to EndofChapter Exercises ............................................................1192
MATLAB Computing
Pages 113 through 119, 1113, 1115, 1118, 1122, 1123, 1125, 1129,
1133, 1134, 1135, 1136, 1137, 1139, 1141, 1143, 1144, 1146,
1148, 1152, 1154, 1155, 1156, 1157, 1158, 1159, 1161 through 1166,
1176, and 1192 through 11100
Simulink Modeling
Pages 1169, 1171, 1174, 1177, 1178, 1179, 1180, 1181

A Introduction to MATLAB A1
A.1 MATLAB and Simulink............................................................................... A1
A.2 Command Window............................................................................................ A1
A.3 Roots of Polynomials.......................................................................................... A3
A.4 Polynomial Construction from Known Roots................................................... A4
A.5 Evaluation of a Polynomial at Specified Values................................................. A5
A.6 Rational Polynomials.......................................................................................... A7
A.7 Using MATLAB to Make Plots .......................................................................... A9
A.8 Subplots ............................................................................................................ A16
A.9 Multiplication, Division, and Exponentiation................................................. A17
A.10 Script and Function Files ................................................................................. A24
A.11 Display Formats ................................................................................................ A29

MATLAB Computing
Pages A3 through A9, A11, A13, A15, A16,
A20, A21, A23, A26, A27

B Introduction to Simulink B 1
B.1 Simulink and its Relation to MATLAB ............................................................... B1
B.2 Simulink Demos ................................................................................................. B20
MATLAB Computing
Page B4

viii Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Simulink Modeling
Pages B7, B12, B14, B18

C A Review of Complex Numbers C1

C.1 Definition of a Complex Number........................................................................ C1
C.2 Addition and Subtraction of Complex Numbers ................................................ C2
C.3 Multiplication of Complex Numbers................................................................... C2
C.4 Division of Complex Numbers ............................................................................ C4
C.5 Exponential and Polar Forms of Complex Numbers .......................................... C4
MATLAB Computing
Pages C6, C7, C8
Simulink Modeling
Page C7

D Matrices and Determinants D1

D.1 Matrix Definition............................................................................................... D1
D.2 Matrix Operations ............................................................................................. D2
D.3 Special Forms of Matrices.................................................................................. D5
D.4 Determinants ..................................................................................................... D9
D.5 Minors and Cofactors......................................................................................D11
D.6 Cramers Rule ..................................................................................................D15
D.7 Gaussian Elimination Method ........................................................................D17
D.8 The Adjoint of a Matrix...................................................................................D19
D.9 Singular and NonSingular Matrices ..............................................................D19
D.10 The Inverse of a Matrix ...................................................................................D20
D.11 Solution of Simultaneous Equations with Matrices .......................................D22

MATLAB Computing
Pages D2, D3, D5, D6, D7, D8, D10,
D13 D14, D17, D21, D23, D24, D27
Simulink Modeling
Page D3, D28
Excel Spreadsheet
Page D25

E Window Functions E1
E.1 Window Function Defined.....................................................................................E1
E.2 Common Window Functions.................................................................................E1
E.2.1 Rectangular Window Function....................................................................E3
E.2.2 Triangular Window Function......................................................................E5

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition ix
Copyright Orchard Publications
 E.2.3 Hanning Window Function........................................................................ E7
E.2.4 Hamming Window Function...................................................................... E9
E.2.5 Blackman Window Function .................................................................... E11
E.2.6 Kaiser Family of Window Functions......................................................... E13
E.3 Other Window Functions.................................................................................... E14
E.4 Fourier Series Method for Approximating an FIR Amplitude Response ........... E15
MATLAB Computing
Pages E3 through E34

F Correlation Functions F -1
F.1 Cross Correlation .....................................................................................................F-1
F.2 Autocorrelation ........................................................................................................F-6
MATLAB Computing
Pages F8, F9
Simulink Modeling
Pages F3, F5, F7, F9,

G Nonlinear Systems G 1
E.1 Describing Functions .............................................................................................G1
Simulink Modeling
Page G4

References R1

Index IN1

x Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 1 Elementary Signals

Other forms of the unit step function are shown in Figure 1.8.

t t t
0 0 0
(a) (b) (c)
A A A
A u0 ( t ) A u0 ( t T ) A u0 ( t + T )

Au 0 ( t ) Au 0 ( t + T ) Au 0 ( t T )
A A A

t t 0
t
0 0 (e) (f)
(d)

t t 0 t
0 (g) 0 (h) (i)
A A A
A u0 ( t ) A u0 ( t + T ) A u0 ( t T )

Figure 1.8. Other forms of the unit step function

Unit step functions can be used to represent other timevarying functions such as the rectangular
pulse shown in Figure 1.9.
u0 ( t )
1
1
t t t
0 1 0 0
(c)
(a) (b)
u0 ( t 1 )
Figure 1.9. A rectangular pulse expressed as the sum of two unit step functions
Thus, the pulse of Figure 1.9(a) is the sum of the unit step functions of Figures 1.9(b) and 1.9(c) and
it is represented as u 0 ( t ) u 0 ( t 1 ) .

The unit step function offers a convenient method of describing the sudden application of a voltage
or current source. For example, a constant voltage source of 24 V applied at t = 0 , can be denoted
as 24u 0 ( t ) V . Likewise, a sinusoidal voltage source v ( t ) = V m cos t V that is applied to a circuit at
t = t 0 , can be described as v ( t ) = ( V m cos t )u 0 ( t t 0 ) V . Also, if the excitation in a circuit is a
rectangular, or triangular, or sawtooth, or any other recurring pulse, it can be represented as a sum
(difference) of unit step functions.

14 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Higher Order Delta Functions

8. To drag a point along the xaxis, we select that point, and we hold down the Shift key while drag-
ging that point.
9. When we select a line segment on the time axis (xaxis) we observe that at the lower end of the
waveform display window the Left Point and Right Point fields become visible. We can then
reshape the given waveform by specifying the Time (T) and Amplitude (Y) points.

Figure 1.26. Waveforms for the Simulink model of Figure 1.24

The two positive spikes that occur at t = 2 , and t = 7 , are clearly shown in Figure 1.26.
MATLAB* has built-in functions for the unit step, and the delta functions. These are denoted by
the names of the mathematicians who used them in their work. The unit step function u 0 ( t ) is
referred to as Heaviside(t), and the delta function ( t ) is referred to as Dirac(t). Their use is illus-
trated with the examples below.
syms k a t;% Define symbolic variables
u=k*sym('Heaviside(t-a)')% Create unit step function at t = a
u=
k*Heaviside(t-a)
d=diff(u)% Compute the derivative of the unit step function
d=
k*Dirac(t-a)
int(d)% Integrate the delta function
ans =
Heaviside(t-a)*k
The MATLAB heaviside function can be used to plot the unit step, unit impulse, and unit ramp
functions as illustrated in Figures 1.27 through 1.29, the unit impulse in Figure 1.30, and the unit
ramp in Figure 1.31.
The plot in Figure 1.27 was generated with the MATLAB script
ezplot('heaviside(t-1)',[0 10]); grid

* An introduction to MATLAB is presented in Appendix A.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 119
Copyright Orchard Publications
Chapter 2 The Laplace Transformation

2.2.3 Frequency Shifting Property

The frequency shifting property states that if we multiply a time domain function f ( t ) by an exponen-
at
tial function e where a is an arbitrary positive constant, this multiplication will produce a shift
of the s variable in the complex frequency domain by a units. Thus,
at
e f(t) F(s + a ) (2.14)
Proof:

( s + a )t
0 0 f ( t ) e
at at st
L {e f( t) } = e f(t )e dt = dt = F ( s + a )

Note 2:
A change of scale is represented by multiplication of the time variable t by a positive scaling factor
a . Thus, the function f ( t ) after scaling the time axis, becomes f ( at ) .

2.2.4 Scaling Property

Let a be an arbitrary positive constant; then, the scaling property states that

1 s
f ( at ) --- F -- (2.15)
a a
Proof:

0 f ( at ) e
st
L { f ( at ) } = dt

and letting t = a , we obtain

1 1 s
d -- = --- d ( ) = --- F --
s ( a ) ( s a )
L { f ( at ) } = 0 f( )e
a a 0 f ( ) e a a
Note 3:
Generally, the initial value of f ( t ) is taken at t = 0 to include any discontinuity that may be pres-
ent at t = 0 . If it is known that no such discontinuity exists at t = 0 , we simply interpret f ( 0 ) as
f( 0) .

2.2.5 Differentiation in Time Domain Property

The differentiation in time domain property states that differentiation in the time domain corre-
sponds to multiplication by s in the complex frequency domain, minus the initial value of f ( t ) at

t = 0 . Thus,
d
f ' ( t ) = ----- f ( t ) sF ( s ) f ( 0 ) (2.16)
dt

24 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Laplace Transforms of Common Functions of Time

TABLE 2.1 Summary of Laplace Transform Properties and Theorems

Property/Theorem Time Domain Complex Frequency Domain

1 Linearity c1 f1 ( t ) + c2 f2 ( t ) c1 F1 ( s ) + c2 F2 ( s )
+ + cn fn ( t ) + + cn Fn ( s )

2 Time Shifting f ( t a )u 0 ( t a ) as
e F(s)
3 Frequency Shifting as F( s + a)
e f(t)
4 Time Scaling f ( at ) 1
--- F -s-
a a

5 Time Differentiation d- sF ( s ) f ( 0 )

---- f(t)
See also (2.18) through (2.20) dt
6 Frequency Differentiation tf ( t ) d
See also (2.22) ----- F ( s )
ds
7 Time Integration t
f ( ) d F( s) f ( 0 )
---------- + -------------
s s
8 Frequency Integration f(t)
--------
t s F ( s ) ds
9 Time Periodicity f ( t + nT ) T
st
0 f ( t ) e dt
------------------------------
sT
1e
10 Initial Value Theorem lim f ( t ) lim sF ( s ) = f ( 0 )

t0 s
11 Final Value Theorem lim f ( t ) lim sF ( s ) = f ( )
t s0
12 Time Convolution f 1 ( t )*f 2 ( t ) F 1 ( s )F 2 ( s )

13 Frequency Convolution f 1 ( t )f 2 ( t ) 1-
------- F ( s )*F 2 ( s )
2j 1

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 213
Copyright Orchard Publications
Chapter 3
The Inverse Laplace Transformation

T
his chapter is a continuation to the Laplace transformation topic of the previous chapter and
presents several methods of finding the Inverse Laplace Transformation. The partial fraction
expansion method is explained thoroughly and it is illustrated with several examples.

3.1 The Inverse Laplace Transform Integral

The Inverse Laplace Transform Integral was stated in the previous chapter; it is repeated here for
convenience.
+ j
1
j
1 st
L { F ( s ) } = f ( t ) = -------- F ( s ) e ds (3.1)
2j

This integral is difficult to evaluate because it requires contour integration using complex variables
theory. Fortunately, for most engineering problems we can refer to Tables of Properties, and Com-
mon Laplace transform pairs to lookup the Inverse Laplace transform.

3.2 Partial Fraction Expansion

Quite often the Laplace transform expressions are not in recognizable form, but in most cases
appear in a rational form of s , that is,
N(s)
F ( s ) = ----------- (3.2)
D(s)

where N ( s ) and D ( s ) are polynomials, and thus (3.2) can be expressed as

m m1 m2
N(s) bm s + bm 1 s + bm 2 s + + b1 s + b0
F ( s ) = ----------- = -------------------------------------------------------------------------------------------------------------------
- (3.3)
D(s) n
an s + an 1 s
n 1
+ an 2 s
n 2
+ + a1 s + a0

The coefficients a k and b k are real numbers for k = 1, 2, , n , and if the highest power m of
N ( s ) is less than the highest power n of D ( s ) , i.e., m < n , F ( s ) is said to be expressed as a proper
rational function. If m n , F ( s ) is an improper rational function.
In a proper rational function, the roots of N ( s ) in (3.3) are found by setting N ( s ) = 0 ; these are
called the zeros of F ( s ) . The roots of D ( s ) , found by setting D ( s ) = 0 , are called the poles of F ( s ) .
We assume that F ( s ) in (3.3) is a proper rational function. Then, it is customary and very conve-
n
nient to make the coefficient of s unity; thus, we rewrite F ( s ) as

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 31
Copyright Orchard Publications
 Partial Fraction Expansion

or k1
1 d
- [ ( s p1 )m F ( s ) ]
r 1k = lim ------------------ ------------- (3.42)
s p 1 ( k 1 )! ds k 1

Example 3.4
Use the partial fraction expansion method to simplify F 4 ( s ) of (3.43) below, and find the time
domain function f 4 ( t ) corresponding to F 4 ( s ) .

s+3
F 4 ( s ) = -----------------------------------2 (3.43)
(s + 2)(s + 1)
Solution:
We observe that there is a pole of multiplicity 2 at s = 1 , and thus in partial fraction expansion
form, F 4 ( s ) is written as
s+3 r1 r 21 r 22
F 4 ( s ) = -----------------------------------2 = ---------------- + ------------------2 + ---------------- (3.44)
(s + 2)(s + 1) (s + 2) (s + 1) (s + 1)
The residues are
s+3
r 1 = ------------------2 = 1
(s + 1) s = 2

s+3
r 21 = ----------- = 2
s+2 s = 1

d s+3 ( s + 2 ) ( s + 3 -)
r 22 = ----- ----------- = -------------------------------------- = 1
ds s + 2 (s + 2)
2
s = 1 s = 1

The value of the residue r 22 can also be found without differentiation as follows:

Substitution of the already known values of r 1 and r 21 into (3.44), and letting s = 0 *, we obtain

s+3 1 2 r 22
----------------------------------- = ---------------- + ------------------2 + ----------------
(s + 1) (s + 2)
2 (s + 2) s=0 (s + 1) (s + 1) s=0
s=0 s=0
or
--- = 1
3 --- + 2 + r 22
2 2
from which r 22 = 1 as before. Finally,
s+3 1 2 1 2t t t
F 4 ( s ) = -----------------------------------2 = ---------------- + ------------------2 + ---------------- e + 2te e = f 4 ( t ) (3.45)
(s + 2)(s + 1) (s + 2) (s + 1) (s + 1)

* This is permissible since (3.44) is an identity.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 39
Copyright Orchard Publications
Chapter 3 The Inverse Laplace Transformation

Check with MATLAB:

syms s t; Fs=(s+3)/((s+2)*(s+1)^2); ft=ilaplace(Fs)
ft = exp(-2*t)+2*t*exp(-t)-exp(-t)
We can use the following script to check the partial fraction expansion.
syms s
Ns = [1 3]; % Coefficients of the numerator N(s) of F(s)
expand((s + 1)^2); % Expands (s + 1)^2 to s^2 + 2*s + 1;
d1 = [1 2 1]; % Coefficients of (s + 1)^2 = s^2 + 2*s + 1 term in D(s)
d2 = [0 1 2]; % Coefficients of (s + 2) term in D(s)
Ds=conv(d1,d2); % Multiplies polynomials d1 and d2 to express the
% denominator D(s) of F(s) as a polynomial
[r,p,k]=residue(Ns,Ds)
r =
1.0000
-1.0000
2.0000
p =
-2.0000
-1.0000
-1.0000
k =
[]

Example 3.5
Use the partial fraction expansion method to simplify F 5 ( s ) of (3.46) below, and find the time
domain function f 5 ( t ) corresponding to the given F 5 ( s ) .
2
s + 3s + 1
F 5 ( s ) = -------------------------------------
3
-
2
(3.46)
(s + 1) (s + 2)
Solution:
We observe that there is a pole of multiplicity 3 at s = 1 , and a pole of multiplicity 2 at s = 2 .
Then, in partial fraction expansion form, F 5 ( s ) is written as

r 11 r 12 r 13 r 21 r 22
F 5 ( s ) = ------------------ + ------------------ + ---------------
- + ------------------ + ---------------
- (3.47)
(s + 1)
3
(s + 1)
2 (s + 1) (s + 2) 2 (s + 2)

The residues are

310 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 4 Circuit Analysis with Laplace Transforms

t=0:0.01:10;...
Vout=1.36.*exp(6.57.*t)+0.64.*exp(0.715.*t).*cos(0.316.*t)1.84.*exp(0.715.*t).*sin(0.316.*t);...
plot(t,Vout); grid

1.5

0.5

-0.5
0 2 4 6 8 10

Figure 4.10. Plot of v out ( t ) for the circuit of Example 4.3

4.2 Complex Impedance Z(s)

Consider the s domain RLC series circuit of Figure 4.11, where the initial conditions are
assumed to be zero.

R sL
+
+ V out ( s )
I(s ) 1-
VS ( s ) -----
sC
Figure 4.11. Series RLC circuit in sdomain

1
For this circuit, the sum R + sL + ------ represents the total opposition to current flow. Then,
sC

VS ( s )
I ( s ) = ------------------------------------ (4.14)
R + sL + 1 sC

and defining the ratio V s ( s ) I ( s ) as Z ( s ) , we obtain

48 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Using the Simulink Transfer Fcn Block

w=1:10:10000; Gs=1./(2.5.*10.^(6).*w.^25.*j.*10.^(3).*w+5);...
semilogx(w,abs(Gs)); xlabel('Radian Frequency w'); ylabel('|Vout/Vin|');...
title('Magnitude Vout/Vin vs. Radian Frequency'); grid
The plot is shown in Figure 4.22. We observe that the given op amp circuit is a second order low
pass filter whose cutoff frequency ( 3 dB ) occurs at about 700 r s .
Magnitude Vout/Vin vs. Radian Frequency
0.2

0.15
|Vout/Vin|

0.1

0.05

0
0 1 2 3 4
10 10 10 10 10
Radian Frequency w

Figure 4.22. G ( j ) versus for the circuit of Example 4.7

4.5 Using the Simulink Transfer Fcn Block

The Simulink Transfer Fcn block shown above implements a transfer function where the input
V IN ( s ) and the output V OUT ( s ) can be expressed in transfer function form as

V OUT ( s )
G ( s ) = --------------------- (4.33)
V IN ( s )

Example 4.8
Let us reconsider the active lowpass filter op amp circuit of Figure 4.21, Page 4-15 where we
found that the transfer function is
V out ( s ) 1
G ( s ) = ------------------- = -------------------------------------------------------------------------------------------------------------------------------- (4.34)
V in ( s ) R 1 [ ( 1 R 1 + 1 R 2 + 1 R 3 + sC 1 ) ( sR 3 C 2 ) + 1 R 2 ]

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 417
Copyright Orchard Publications
Chapter 5 State Variables and State Equations

We observe that
x 1 = x 2
x 2 = x 3
x 3 = x 4 (5.17)
4
d y
--------- = x 4 = a 0 x 1 a 1 x 2 a 2 x 3 a 3 x 4 + u ( t )
4
dt
and in matrix form
x 1 0 1 0 0 x1 0
x 2 0 0 1 0 x2
= + 0 u(t) (5.18)
x 3 0 0 0 1 x3 0
x 4 a0 a1 a2 a3 x4 1

In compact form, (5.18) is written as

x = Ax + bu (5.19)
where
x 1 0 1 0 0 x1 0
x 2 0 0 1 0 x2
x = , A= , x= , b= 0, and u = u ( t )
x 3 0 0 0 1 x3 0
x 4 a0 a1 a2 a3 x4 1

We can also obtain the state equations directly from given circuits. We choose the state variables to
represent inductor currents and capacitor voltages. In other words, we assign state variables to
energy storing devices. The examples below illustrate the procedure.

Example 5.3
Write state equation(s) for the circuit of Figure 5.2, given that v C ( 0 ) = 0 , and u 0 ( t ) is the unit
step function.
R

+ +
v C ( t ) = v out ( t )
C

vS u0 ( t )

Figure 5.2. Circuit for Example 5.3

Solution:
This circuit contains only one energystoring device, the capacitor. Therefore, we need only one
state variable. We choose the state variable to denote the voltage across the capacitor as shown in

54 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Computation of the State Transition Matrix

2 n1 1 t
a0 + a1 1 + a2 1 + + an 1 1 = e
2 n1 2 t
a0 + a1 2 + a2 2 + + an 1 2 = e

2 n1 n t
a0 + a1 n + a2 n + + an 1 n = e

We use as many equations as the number of the eigenvalues, and we solve for the coefficients
ai .

4. We substitute the a i coefficients into the state transition matrix of (5.54), and we simplify.

Example 5.7
At
Compute the state transition matrix e given that

5 7 5
A = 0 4 1 (5.55)
2 8 3
Solution:
1. We first compute the eigenvalues from det [ A I ] = 0 . We obtain [ A I ] at once, by sub-
tracting from each of the main diagonal elements of A . Then,

5 7 5
det [ A I ] = det 0 4 1 = 0 (5.56)
2 8 3

and expansion of this determinant yields the polynomial

3 2
6 + 11 6 = 0 (5.57)
We will use MATLAB roots(p) function to obtain the roots of (5.57).
p=[1 6 11 6]; r=roots(p); fprintf(' \n'); fprintf('lambda1 = %5.2f \t', r(1));...
fprintf('lambda2 = %5.2f \t', r(2)); fprintf('lambda3 = %5.2f', r(3))
lambda1 = 3.00 lambda2 = 2.00 lambda3 = 1.00
and thus the eigenvalues are
1 = 1 2 = 2 3 = 3 (5.58)

2. Since A is a 3 3 matrix, we use the first 3 terms of (5.54), that is,

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 513
Copyright Orchard Publications
Chapter 5 State Variables and State Equations

We can obtain the plot of Figure 5.8 with the Simulink StateSpace block with the unit step func-
tion as the input using the Step block, and the capacitor voltage as the output displayed on the
Scope block as shown in the model of Figure 5.9 where for the StateSpace block Function Block
Parameters dialog box we have entered:
A: [-4 -4; 3/4 0]
B: [4 0]
C: [0 1]
D: [ 0 ]
Initial conditions: [0 1/2]

Figure 5.9. Simulink model for Example 5.10

The waveform for the capacitor voltage for the simulation time interval 0 t 10 seconds is shown
in Figure 5.9 where we observe that the initial condition v C ( 0 ) = 0.5 V is also displayed.

Example 5.11
A network is described by the state equation
x = Ax + bu (5.98)
where

A = 1 0 x0 = 1 b = 1 and u = ( t ) (5.99)
1 1 0 1
Compute the state vector
x1
x =
x2
Solution:
We compute the eigenvalues from

526 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 6 The Impulse Response and Convolution

A(t ) A
e
At
y( t) = e bu ( ) d = e bu ( ) d (6.21)

6.4 Graphical Evaluation of the Convolution Integral

The convolution integral is more conveniently evaluated by the graphical evaluation. The procedure
is best illustrated with the following examples.

Example 6.4
The signals h ( t ) and u ( t ) are as shown in Figure 6.6. Compute h ( t )*u ( t ) using the graphical evalu-
ation.
u ( t ) = u0 ( t ) u0 ( t 1 )
1 h(t) = t + 1 1

t t
0 1 0 1

Figure 6.6. Signals for Example 6.4

Solution:
The convolution integral states that

h ( t ) u ( t ) = u ( t )h ( ) d (6.22)

where is a dummy variable, that is, u ( ) and h ( ) , are considered to be the same as u ( t ) and
h ( t ) . We form u ( t ) by first constructing the image of u ( ) ; this is shown as u ( ) in Figure 6.7.

u ( )
1

1 0

Figure 6.7. Construction of u ( ) for Example 6.4
Next, we form u ( t ) by shifting u ( ) to the right by some value t as shown in Figure 6.8.

1 u ( t )

0 t

Figure 6.8. Formation of u ( t ) for Example 6.4

68 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 6 The Impulse Response and Convolution

Example 6.7
For the circuit of Figure 6.31, use the convolution integral to find the capacitor voltage when the
input is the unit step function u 0 ( t ) , and v C ( 0 ) = 0 .

1
C +
+ vC ( t )
1F
u0 ( t )

Figure 6.31. Circuit for Example 6.7

Solution:
Before we apply the convolution integral, we must know the impulse response h ( t ) of this circuit.
The circuit of Figure 6.31 was analyzed in Example 6.1, Page 62, where we found that
1 t RC
h ( t ) = -------- e u0 ( t ) (6.33)
RC
With the given values, (6.33) reduces to
t
h ( t ) = e u0 ( t ) (6.34)

Next, we use the graphical evaluation of the convolution integral as shown in Figures 6.32 through
6.34.
The formation of u 0 ( ) is shown in Figure 6.32.
u0 ( )
1

Figure 6.32. Formation of u 0 ( ) for Example 6.7

Figure 6.33 shows the formation of u 0 ( t ) .

u0 ( t ) 1

0 t

Figure 6.33. Formation of u 0 ( t ) for Example 6.7

618 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Evaluation of the Coefficients

sin x cos x

sin x cos x

2
Figure 7.2. Graphical proof of 0 ( sin mt ) ( cos nt ) dt = 0

Moreover, if m and n are different integers, then,

2

0 ( sin mt ) ( sin nt ) dt = 0 (7.6)

since
1
( sin x ) ( sin y ) = --- [ cos ( x y ) cos ( x y ) ]
2

The integral of (7.6) can also be confirmed graphically as shown in Figure 7.3, where m = 2 and
n = 3 . We observe that the net shaded area above and below the time axis is zero.

sin 2x sin 3x sin 2x sin 3x

2
Figure 7.3. Graphical proof of 0 ( sin mt ) ( sin nt ) dt = 0 for m = 2 and n = 3

Also, if m and n are different integers, then,

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 73
Copyright Orchard Publications
 Trigonometric Form of Fourier Series for Common Waveforms

ics will be present since this waveform has also halfwave symmetry. However, we will compute all
coefficients to verify this. Also, for brevity, we will assume that = 1

T
A

2
t
0

Figure 7.12. Square waveform as odd function

The a i coefficients are found from
2 2
1 1 A
0 0
2
a n = --- f ( t ) cos nt dt = --- A cos nt dt + ( A ) cos nt dt = ------ ( sin nt 0 sin nt )
n
(7.19)
A A
= ------ ( sin n 0 sin n2 + sin n ) = ------ ( 2 sin n sin n2 )
n n

and since n is an integer (positive or negative) or zero, the terms inside the parentheses on the sec-
ond line of (7.19) are zero and therefore, all a i coefficients are zero, as expected since the square
waveform has odd symmetry. Also, by inspection, the average ( DC ) value is zero, but if we attempt
to verify this using (7.19), we will obtain the indeterminate form 0 0 . To work around this prob-
lem, we will evaluate a 0 directly from (7.12), Page 76. Thus,
2
1 A
a 0 = ---
0 A dt + ( A ) dt = ---- ( 0 2 + ) = 0

(7.20)

The b i coefficients are found from (7.14), Page 76, that is,
2 2
1 1 A
0 0
2
b n = --- f ( t ) sin nt dt = --- A sin nt dt + ( A ) sin nt dt = ------ ( cos n t 0 + cos nt )
n
(7.21)
A A
= ------ ( cos n + 1 + cos 2n cos n ) = ------ ( 1 2 cos n + cos 2n )
n n

For n = even , (7.21) yields

A
b n = ------ ( 1 2 + 1 ) = 0
n
as expected, since the square waveform has halfwave symmetry.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 711
Copyright Orchard Publications
Chapter 7 Fourier Series

Figure 7.30. Waveform for relation (7.87)

The waveform of Figure 7.30 is a rudimentary presentation of the capacitor voltage for the circuit
of Figure 7.27. However, it will improve if we add a sufficient number of harmonics in (7.87).
We can obtain a more accurate waveform for the capacitor voltage in Figure 7.27 with the Simu-
link model of Figure 7.31.

Figure 7.31. Simulink model for the circuit of Figure 7.27

7.8 The Exponential Form of the Fourier Series

The Fourier series are often expressed in exponential form. The advantage of the exponential form
is that we only need to perform one integration rather than two, one for the a n , and another for
the b n coefficients in the trigonometric form of the series. Moreover, in most cases the integration
is simpler.
The exponential form is derived from the trigonometric form by substitution of
jt jt
e +e
cos t = ---------------------------- (7.88)
2
jt jt
e e
sin t = --------------------------- (7.89)
j2
into f ( t ) . Thus,

730 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 8 The Fourier Transform

TABLE 8.7 Time Domain and Frequency Domain Correspondence (Completed Table)

f(t) F()
Real Imaginary Complex Even Odd
Real
Real and Even
Real and Odd
Imaginary
Imaginary and Even
Imaginary and Odd

Table 8.7 is now complete and shows that if f ( t ) is real (even or odd), the real part of F ( ) is even,
and the imaginary part is odd. Then,
F Re ( ) = F Re ( ) f ( t ) = Real (8.31)
and
F Im ( ) = F Im ( ) f ( t ) = Real (8.32)
Since,
F ( ) = F Re ( ) + jF Im ( ) (8.33)
it follows that
F ( ) = F Re ( ) + jF Im ( ) = F Re ( ) jF Im ( )
or
F ( ) = F ( ) f ( t ) = Real (8.34)

Now, if F ( ) of some function of time f ( t ) is known, and F ( ) is such that F ( ) = F ( ) , can

we conclude that f ( t ) is real? The answer is yes; we can verify this with (8.14), Page 83, which is
repeated here for convenience.

1
f Im ( t ) = ------
2 [ FRe ( ) sin t + FIm ( ) cos t ] d (8.35)

We observe that the integrand of (8.35) is zero since it is an odd function with respect to because
both products inside the brackets are odd functions*.
Therefore, f Im ( t ) = 0 , that is, f ( t ) is real.
Accordingly, we can state that a necessary and sufficient condition for f ( t ) to be real, is that
F ( ) = F ( ) .

* In relations (8.31) and (8.32) above, we determined that F Re ( ) is even and F Im ( ) is odd.

88 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Fourier Transform Pairs of Common Functions

TABLE 8.8 Fourier Transform Properties and Theorems

Property f(t) F( )
Linearity a1 f1 ( t ) + a2 f2 ( t ) + a1 F1 ( ) + a2 F2 ( ) +
Symmetry F(t) 2f ( )
Time Scaling f ( at ) 1
----- F ----
a a
Time Shifting f ( t t0 ) jt 0
F ( )e
Frequency Shifting j 0 t F ( 0 )
e f(t)
Time Differentiation n n
d ( j ) F ( )
--------- f ( t )
n
dt
Frequency Differentiation n n
( jt ) f ( t ) d
----------- F ( )
n
d
Time Integration t F ( -)
----------- + F ( 0 ) ( )
f ( ) d j

Conjugate Functions f ( t ) F ( )
Time Convolution f 1 ( t ) f 2 ( t ) F1 ( ) F2 ( )
Frequency Convolution f1 ( t ) f2 ( t ) 1
------ F 1 ( ) F 2 ( )
2
Area under f ( t )
F( 0) = f ( t ) dt
Area under F ( ) 1

f ( 0 ) = ------
2
F ( ) d

Parsevals Theorem
1

2 2
f ( t ) dt = ------ F ( ) d
2

Likewise, the Fourier transform for the shifted delta function ( t t 0 ) is

jt 0
( t t0 ) e (8.61)

We will use the notation f ( t ) F ( ) to show the time domain to frequency domain correspon-
dence. Thus, (8.60) may also be denoted as in Figure 8.1.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 817
Copyright Orchard Publications
Chapter 8 The Fourier Transform

f(t) F() 1
(t)

0 t 0

Figure 8.1. The Fourier transform of the delta function

8.4.2 The Constant Function Pair

A 2A ( ) (8.62)
Proof:

1 1
F { 2A ( ) } = ------
2 2A ( )e
jt
d = A ( )e
jt
d = Ae
jt
=0
= A

and (8.62) follows.

The f ( t ) F ( ) correspondence is also shown in Figure 8.2.

f(t) F()
A
2A ( )

t
0 0
Figure 8.2. The Fourier transform of constant A
Also, by direct application of the Inverse Fourier transform, or the frequency shifting property and
(8.62), we derive the transform
j 0 t
e 2 ( 0 ) (8.63)

The transform pairs of (8.62) and (8.63) can also be derived from (8.60) and (8.61) by using the
symmetry property F ( t ) 2f ( )

8.4.3 The Cosine Function Pair

1 j0 t j 0 t
cos 0 t = --- ( e +e ) ( 0 ) + ( + 0 ) (8.64)
2

Proof:
This transform pair follows directly from (8.63). The f ( t ) F ( ) correspondence is also shown in
Figure 8.3.

818 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 The Z Transform of Common DiscreteTime Functions

where A ( z ) and B ( z ) are polynomials with real coefficients.

For convenience, we summarize the properties and theorems of the Z transform in Table 9.1.

TABLE 9.1 Properties and Theorems of the Z transform

Property / Theorem Time Domain Z transform
Linearity af 1 [ n ] + bf 2 [ n ] + aF 1 ( z ) + bF 2 ( z ) +
Shift of x [ n ]u 0 [ n ] f [ n m ]u 0 [ n m ] z
m
F(z )
Right Shift f [n m] m1

f [ n m ]z
m n
z F(z ) +
n=0

Left Shift f [n + m] 1

m n
z F(z) + f [ n + m ]z
n = m

Multiplication by a n n
z
a f [n] F --
a
naT naT aT
Multiplication by e e f [n] F(e z)
Multiplication by n nf [ n ] d
z ------ F ( z )
dz
Multiplication by n 2 2
n f [n] d 2d
2
z ------ F ( z ) + z -------2- F ( z )
dz dz
Summation in Time n
z -
----------
f [m]
m=0
z 1
F(z)

Time Convolution f 1 [ n ] f 2 [ n ] F1 ( z ) F2 ( z )

Frequency Convolution f1 [ n ] f2 [ n ] 1 z

xF ( v )F --v- v
1
-------- 1 2 dv
j 2
Initial Value Theorem f [ 0 ] = lim F ( z )
z
Final Value Theorem lim f [ n ] = lim ( z 1 ) F ( z )
n z1

9.3 The Z Transform of Common DiscreteTime Functions

In this section we will provide several examples to find the Z transform of some discretetime
functions. In this section, we will derive the Z transforms of the most common discretetime
functions in Subsections 9.3.1 through 9.3.5 below.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 911
Copyright Orchard Publications
 The Inverse Z Transform

syms n z
fn=2*(0.5)^n9*(0.75)^n+8; % This is the answer in (9.80)
Fz=ztrans(fn,n,z); simple(Fz) % Verify answer by first taking Z transform of f[n]
ans = 8*z^3/(2*z-1)/(4*z-3)/(z-1)
iztrans(Fz) % Now, verify that Inverse of F(z) gives back f[n]
ans = 2*(1/2)^n-9*(3/4)^n+8
We can use Microsoft Excel to obtain and plot the values of f [ n ] . The spreadsheet of Figure 9.6
shows the first 25 values of n but only part of the spreadsheet is shown.

n f[n] n n
Discrete Time Sequence f[n] = 2(0.5) 9(0.75) + 8
0.000 1.0000 8
1.000 2.2500
2.000 3.438 7
3.000 4.453 6
4.000 5.277
5.000 5.927 5
6.000 6.429 4
7.000 6.814
8.000 7.107 3
9.000 7.328
2
10.000 7.495
11.000 7.621 1
12.000 7.715
0
13.000 7.786
1

11

13

15

17

19

21

23

25
14.000 7.84
15.000 7.88
n n
Figure 9.6. The discretetime sequence f [ n ] = 2 ( 0.5 ) 9 ( 0.75 ) + 8 for Example 9.4

Example 9.5
Use the partial fraction expansion method to compute the Inverse Z transform of
12z
F ( z ) = -----------------------------------2- (9.81)
Solution: (z + 1)(z 1)
Division of both sides by z and partial expansion yields
F(z) 12 r1 r2 r3
----------- = ------------------------------------ = ---------------- + ------------------ + ---------------
-
z (z + 1)(z 1)
2 ( z + 1 ) (z 1)
2 ( z 1)
The residues are
12 12
r 1 = ------------------2 = ----------------------2- = 3
(z 1) z = 1
( 1 1)

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 925
Copyright Orchard Publications
 The Discrete Fourier Transform (DFT)

A B C D E
1 Input data x(n) are same as in Example 10.1
2 and are entered in cells A11 through A14
3
4 From the Tools drop down menu, we select
5 Data Analysis and from it, Fourier Analysis
6
7 The Input Range is A11 through A14 (A11:A14)
8 and the Output Range is B11 through B14 (B11:B14)
9
10 x(n) X(m)
11 1 6
12 2 -1-i
13 2 0
14 1 -1+i
15
16 To obtain the discrete time sequence, we select
17 Inverse from the Fourier Analysis menu
18
19 Input data are the same as in Example 10.2
20
21 The Input Range is A25 through A28 (A25:A28)
22 and the Output Range is B25 through B28 (B25:B28)
23
24 X(m) x(n)
25 6 1
26 -1-j 2
27 0 2
28 -1+j 1
Figure 10.3. Using Excel to find the DFT and Inverse DFT

Example 10.3
Use MATLAB to compute the magnitude of the frequency components of the following discrete
time function. Then, use Excel to display the time and frequency values.

x[0] x[1] x[2] x[3] x[4] x[5] x[6] x[7] x[8] x[9] x[10] x[11] x[12] x[13] x[14] x[15]
1.0 1.5 2.0 2.3 2.7 3.0 3.4 4.1 4.7 4.2 3.8 3.6 3.2 2.9 2.5 1.8

Solution:
We compute the magnitude of the frequency spectrum with the MATLAB script below.
xn=[1 1.5 2 2.3 2.7 3 3.4 4.1 4.7 4.2 3.8 3.6 3.2 2.9 2.5 1.8]; magXm=abs(fft(xn));...
fprintf(' \n'); fprintf('magXm1 = %4.2f \t', magXm(1)); fprintf('magXm2 = %4.2f \t', magXm(2));...

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 107
Copyright Orchard Publications
Chapter 10 The DFT and the FFT Algorithm

0 0 0 0
WN WN WN WN
X [0] 0 1 2 N1
x [0]
X [1] WN WN WN WN x [1]
X [2] = W0 W2 WN
4
WN
2(N 1) x [2] (10.49)
N N

X [N 1] 2 x [N 1]
0 N1 2(N 1) (N 1)
WN WN WN WN

This is a complex Vandermonde matrix and it is expressed in a more compact form as

X [ m ] = WN x [ n ] (10.50)

The algorithm that was developed by Cooley and Tukey, is based on matrix decomposition methods,
where the matrix W N in (10.50) is factored into L smaller matrices, that is,

WN = W1 W2 WL (10.51)

L
where L is chosen as L = log 2 N or N = 2 .

Each row of the matrices on the right side of (10.51) contains only two non-zero terms, unity and
k
W N . Then, the vector X [ m ] is obtained from

X [ m ] = W1 W2 WL x [ n ] (10.52)

The FFT computation begins with matrix W L in (10.52). This matrix operates on vector x [ n ]
producing a new vector, and each component of this new vector, is obtained by one multiplication
and one addition. This is because there are only two nonzero elements on a given row, and one of
these elements is unity. Since there are N components on x [ n ] , there will be N complex additions
and N complex multiplications. This new vector is then operated on by the [ W L 1 ] matrix, then
on [ W L 2 ] , and so on, until the entire computation is completed. It appears that the entire com-
putation would require NL = Nlog 2 N complex multiplications, and also Nlog 2 N additions for a
0 N2
total of 2Nlog 2 N arithmetic operations. However, since W N = 1 , W N = 1 , and other reduc-
tions, it is estimated that only about half of these, that is, Nlog 2 N total arithmetic operations are
2
required by the FFT versus the N computations required by the DFT.
Under those assumptions, we construct Table 10.3 to compare the percentage of computations
achieved by the use of FFT versus the DFT.

1018 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Basic Analog Filters

1
G ( j ) = --------------------------------------------- (11.5)
2 2 2
1 + 1 ( R C )
and the phase angle or argument, is
1
= arg { G ( j ) } = tan ( 1 ( RC ) ) (11.6)

We can obtain a quick sketch for the magnitude G ( j ) versus by evaluating (11.5) at = 0 ,
= 1 RC , and . Thus,

As 0 ,
G(j) 0
For = 1 RC ,
G ( j ) = 1 2 = 0.707
and as ,
G(j) 1

We will use the MATLAB script below to plot G ( j ) versus radian frequency . This is shown in
Figure 11.6 where, for convenience, we let RC = 1 .
w=0:0.02:100; RC=1; magGjw=1./sqrt(1+1./(w.*RC).^2); semilogx(w,magGjw);...
xlabel('Frequency in rad/sec log scale'); ylabel('Magnitude of Vout/Vin');...
title('Magnitude Characteristics of basic RC highpass filter'); grid

Halfpower point

Figure 11.6. Magnitude characteristics of the basic RC highpass filter with RC = 1

We can also obtain a quick sketch for the phase angle, i.e., = arg { G ( j ) } versus , by evaluat-
ing (11.6) at = 0 , = 1 RC , = 1 RC , , and . Thus,
As 0 , atan 0 0
For = 1 RC ,

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition 115
Copyright Orchard Publications
Chapter 11 Analog and Digital Filters

s
G ( s ) actual = G ---------------- (11.34)
actual
that is, we replace s with s actual
Quite often, we require that C , that is, in the stop band of the lowpass filter, the attenuation
to be larger than 20 dB decade , i.e., we require a sharper cutoff. As we have seen from the plots
of Figure 11.15, Page 1114, the Butterworth lowpass filter cutoff becomes sharper for larger val-
ues of k . Accordingly, we generate the plot for different values of k shown in Figure 11.17 using
the MATLAB script below.
w_w0=1:0.02:10; dBk1=20.*log10(sqrt(1./(w_w0.^2+1)));...
dBk2=20.*log10(sqrt(1./(w_w0.^4+1))); dBk3=20.*log10(sqrt(1./(w_w0.^6+1)));...
dBk4=20.*log10(sqrt(1./(w_w0.^8+1))); dBk5=20.*log10(sqrt(1./(w_w0.^10+1)));...
dBk6=20.*log10(sqrt(1./(w_w0.^12+1))); dBk7=20.*log10(sqrt(1./(w_w0.^14+1)));...
dBk8=20.*log10(sqrt(1./(w_w0.^16+1))); semilogx(w_w0,dBk1,w_w0,dBk2,w_w0,dBk3,...
w_w0,dBk4,w_w0,dBk5,w_w0,dBk6,w_w0,dBk7,w_w0,dBk8);...
xlabel('Normalized Frequency (rads/sec) log scale'); ylabel ('Magnitude Response (dB)');...
title('Magnitude Attenuation as a Function of Normalized Frequency');...
set(gca, 'XTick', [1 2 3 4 5 6 7 8 9 10]); grid

Figure 11.17. Attenuation for different values of k

Figure 11.17 indicates that for k = 1 the attenuation is 20 dB decade , for k = 2 the attenua-
tion is 40 dB decade , and so on.

1118 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 11 Analog and Digital Filters

Figure 11.22. Bode plots for example 11.6 using MATLABs bode function in Hz

11.2.2 Chebyshev Type I Analog LowPass Filter Design

The Chebyshev Type I filters are based on approximations derived from the Chebyshev polynomials
C k ( x ) which constitute a set of orthogonal functions.* The coefficients of these polynomials are
tabulated in math tables. See, for example, the Handbook of Mathematical Functions, Dover Publica-
tions. These polynomials are derived from the equations
1
C k ( x ) = cos ( kcos x ) ( x 1) (11.43)
and 1
C k ( x ) = cos h ( kcosh x ) x >1 (11.44)

From (11.43), with k = 0 , we obtain

1
C 0 ( x ) = cos ( 0cos x ) = 1 (11.45)
With k = 1 ,
C 1 ( x ) = cos ( 1cos x ) = x
1
(11.46)

* Two functions are said to be orthogonal if, when multiplied together and integrated over the domain of interest, the
integral becomes zero. The property of orthogonality is usually applied to a class of functions that differ by one or
more variables.
1
We recall that if x = cos y , then y = cos x , and cos y = x .

1124 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 11 Analog and Digital Filters

11.8 Solutions to EndofChapter Exercises

1. We use MATLAB for all computations.
% PART I Find Resistor values for second order Butterworth filter, a=sqrt(2), b = 1
%
a=sqrt(2); b =1; C1=10^(8); C2=C1; fc=1000; wc=2*pi*fc; K=2;
R2=(4*b)/(C1*sqrt(a^2+8*b*(K1))*wc);
R1=b/(C1^2*R2*wc^2); R3=(K*R2)/(K1); R4=K*R2; fprintf(' \n');...
fprintf('R1 = %5.0f Ohms \t',R1); fprintf('R2 = %5.0f Ohms \t',R2);...
fprintf('R3 = %5.0f Ohms \t',R3); fprintf('R4 = %5.0f Ohms \t',R4)
R1=12582 Ohms R2=20132 Ohms R3=40263 Ohms R4=40263 Ohms
We choose standard resistors as close as possible to those found above. These are shown in the
MATLAB script below.
% PART II Plot with standard resistors R1=12.7 K, R2=20.0 K, R3=40.2 K, R4= R3
%
f=10:10:20000; w=2*pi*f; R1=12700; R2=20000; R3=40200; R4=R3; K=1+R4/R3;...
wc=(4*b)/(C1*sqrt(a^2+8*b*(K1))*R2); s=w*j; Gw=(K.*s.^2)./(s.^2+a.*wc.*s./b+wc.^2./b);...
semilogx(f,abs(Gw)); xlabel('Frequency, Hz log scale), ylabel('|Vout/Vin| absolute values');...
title('2nd Order Butterworth HighPass Filter Response'); grid

2. We use MATLAB for all computations.

% PART I Find Resistor values for second order bandpass filter f0 = 1 KHz
%
Q=10; K=10; C1=10^(7); C2=C1; f0=1000; w0=2*pi*f0; R1=(2*Q)/(C1*w0*K);...
R2=(2*Q)/(C1*w0*(1+sqrt((K1)^2+8*Q^2))); R3=(1/(C1^2*w0^2))*(1/R1+1/R2); R4=2*R3;...
R5=R4; fprintf(' \n'); fprintf('R1 = %5.0f Ohms \t',R1); fprintf('R2 = %5.0f Ohms \t',R2);...
fprintf('R3 = %5.0f Ohms \t',R3); fprintf('R4 = %5.0f Ohms \t',R4);...

1192 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Chapter 11 Analog and Digital Filters

Next, we compute G 2 ( z ) without prewarping using the following MATLAB script:

N=2; % # of poles
Rp=3; % Pass band ripple in dB
Ts=0.25; % Sampling period
wc=4; % Analog cutoff frequency
wd=(wc*Ts)/pi;
[Nz,Dz]=cheby1(N,Rp,wd,'high');
%
fprintf('The numerator N(z) coefficients in descending powers of z are: \n\n');
fprintf('%8.4f \t',[Nz]); fprintf(' \n');
fprintf('The denominator D(z) coefficients in descending powers of z are: \n\n');
fprintf('%8.4f \t',[Dz]); fprintf(' \n');
%
fprintf('Press any key to see the plot \n');
pause;
%
w=0:2*pi/300:pi; Gz=freqz(Nz,Dz,w); plot(w,abs(Gz)); grid; xlabel('Frequency (rads/sec)');
ylabel('|H|'); title('HighPass Digital Filter without prewarping')
The numerator N(z) coefficients in descending powers of z are:
0.3689 -0.7377 0.3689
The denominator D(z) coefficients in descending powers of z
are:
1.0000 -0.6028 0.4814
and thus the transfer function and the plot without prewarping are as shown below.
2
0.3689z 0.7377z + 0.3689-
G 2 ( z ) = ---------------------------------------------------------------------
2
z 0.6028z + 0.4814

11100 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Appendix A
Introduction to MATLAB

T
his appendix serves as an introduction to the basic MATLAB commands and functions, pro-
cedures for naming and saving the user generated files, comment lines, access to MATLABs
Editor / Debugger, finding the roots of a polynomial, and making plots. Several examples are
provided with detailed explanations.

A.1 MATLAB and Simulink

MATLAB and Simulink are products of The MathWorks, Inc. These are two outstanding soft-
ware packages for scientific and engineering computations and are used in educational institutions
and in industries including automotive, aerospace, electronics, telecommunications, and environ-
mental applications. MATLAB enables us to solve many advanced numerical problems rapidly and
efficiently.

A.2 Command Window

To distinguish the screen displays from the user commands, important terms, and MATLAB func-
tions, we will use the following conventions:
Click: Click the left button of the mouse
Courier Font: Screen displays
Helvetica Font: User inputs at MATLABs command window prompt >> or EDU>>*
Helvetica Bold: MATLAB functions
Normal Font Bold Italic: Important terms and facts, notes and file names
When we first start MATLAB, we see various help topics and other information. Initially, we are
interested in the command screen which can be selected from the Window drop menu. When the
command screen, we see the prompt >> or EDU>>. This prompt is displayed also after execution of
a command; MATLAB now waits for a new command from the user. It is highly recommended that
we use the Editor/Debugger to write our program, save it, and return to the command screen to exe-
cute the program as explained below.
To use the Editor/Debugger:
1. From the File menu on the toolbar, we choose New and click on MFile. This takes us to the Edi-
tor Window where we can type our script (list of statements) for a new file, or open a previously
saved file. We must save our program with a file name which starts with a letter. Important!
MATLAB is case sensitive, that is, it distinguishes between upper and lowercase letters. Thus, t

* EDU>> is the MATLAB prompt in the Student Version

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition A1
Copyright Orchard Publications
Appendix A Introduction to MATLAB

R1
A

R2
C
V L

Figure A.1. Electric circuit for Example A.10

The ammeter readings were then recorded for each frequency. The magnitude of the impedance
|Z| was computed as Z = V A and the data were tabulated on Table A.1.

TABLE A.1 Table for Example A.10

(rads/s) |Z| Ohms (rads/s) |Z| Ohms

300 39.339 1700 90.603
400 52.589 1800 81.088
500 71.184 1900 73.588
600 97.665 2000 67.513
700 140.437 2100 62.481
800 222.182 2200 58.240
900 436.056 2300 54.611
1000 1014.938 2400 51.428
1100 469.83 2500 48.717
1200 266.032 2600 46.286
1300 187.052 2700 44.122
1400 145.751 2800 42.182
1500 120.353 2900 40.432
1600 103.111 3000 38.845
Plot the magnitude of the impedance, that is, |Z| versus radian frequency .
Solution:
We cannot type (omega) in the MATLAB Command prompt, so we will use the English letter w
instead.
If a statement, or a row vector is too long to fit in one line, it can be continued to the next line by
typing three or more periods, then pressing <enter> to start a new line, and continue to enter data.

A10 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Appendix B
Introduction to Simulink

T
his appendix is a brief introduction to Simulink. This author feels that we can best introduce
Simulink with a few examples. Some familiarity with MATLAB is essential in understanding
Simulink, and for this purpose, Appendix A is included as an introduction to MATLAB.

B.1 Simulink and its Relation to MATLAB

The MATLAB and Simulink environments are integrated into one entity, and thus we can ana-
lyze, simulate, and revise our models in either environment at any point. We invoke Simulink from
within MATLAB. We will introduce Simulink with a few illustrated examples.

Example B.1
For the circuit of Figure B.1, the initial conditions are i L ( 0 ) = 0 , and v c ( 0 ) = 0.5 V . We will
compute v c ( t ) .
R L
+
1 14 H C
+ vC ( t )
i(t)
43 F

vs ( t ) = u0 ( t )
Figure B.1. Circuit for Example B.1
For this example,
dv
i = i L = i C = C --------C- (B.1)
dt
and by Kirchoffs voltage law (KVL),
di
Ri L + L ------L- + v C = u 0 ( t ) (B.2)
dt
Substitution of (B.1) into (B.2) yields
2
dv d vC
RC --------C- + LC ----------
- + vC = u0 ( t ) (B.3)
dt dt
2

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition B1
Copyright Orchard Publications
Introduction to Simulink

Figure B.12. Model Equation_B_26 complete block diagram

8. The initial conditions i L ( 0 ) = C ( dv C dt )

= 0 , and v c ( 0 ) = 0.5 V are entered by double
t=0
clicking the Integrator blocks and entering the values 0 for Integrator 1 block, and 0.5 for Inte-
grator 2 block. We also need to specify the simulation time. This is done by specifying the simu-
lation time to be 10 seconds on the Configuration Parameters from the Simulation drop
menu. We can also specify the simulation time on the toolbar by overwriting the default value
10.0 to another desired value. We can start the simulation on Start from the Simulation drop

menu or by clicking the icon on the toolbar, and before doing this, we click the Step block
and in the Source Block Parameters dialog box we change the Step time from 1 to 0.
9. To see the output waveform, we double click the Scope block, and then clicking on the

Autoscale icon, we obtain the waveform shown in Figure B.13.

Figure B.13. The waveform for the function v C ( t ) for Example B.1

B12 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Appendix C
A Review of Complex Numbers

T
his appendix is a review of the algebra of complex numbers. The basic operations are defined
and illustrated with several examples. Applications using Eulers identities are presented, and
the exponential and polar forms are discussed and illustrated with examples.

C.1 Definition of a Complex Number

In the language of mathematics, the square root of minus one is denoted as i , that is, i = 1 . In
the electrical engineering field, we denote i as j to avoid confusion with current i . Essentially, j is
an operator that produces a 90degree counterclockwise rotation to any vector to which it is
applied as a multiplying factor. Thus, if it is given that a vector A has the direction along the right
side of the xaxis as shown in Figure C.1, multiplication of this vector by the operator j will result
in a new vector jA whose magnitude remains the same, but it has been rotated counterclockwise
by 90 . Also, another multiplication of the new vector jA by j will produce another 90 counter-
clockwise direction. In this case, the vector A has rotated 180 and its new value now is A .
When this vector is rotated by another 90 for a total of 270 , its value becomes j ( A ) = j A . A
fourth 90 rotation returns the vector to its original position, and thus its value is again A . There-
2 3 4
fore, we conclude that j = 1 , j = j , and j = 1 .

jA y

j ( j A ) = j2 A = A A
x
2
j ( j A ) = j A = A

j ( A ) = j 3 A = j A
Figure C.1. The j operator

Note: In our subsequent discussion, we will designate the xaxis (abscissa) as the real axis, and the
yaxis (ordinate) as the imaginary axis with the understanding that the imaginary axis is just as
real as the real axis. In other words, the imaginary axis is just as important as the real axis.*

* We may think the real axis as the cosine axis and the imaginary axis as the sine axis.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition C1
Copyright Orchard Publications
 Exponential and Polar Forms of Complex Numbers

Figure C.3. Simulink model for Example C.6a

b. The real and imaginary components of this complex number are shown in Figure C.4.

Im
2
5

116.6
63.4
Re
1
Figure C.4. The components of 1 + j 2
Then,
1
j 2-
tan -----
2 2 1 = j116.6
1 + j2 = 1 +2 e 5e = 5 116.6
Check with MATLAB:
y=1+j*2; magy=abs(y); thetay=angle(y)*180/pi; disp(magy); disp(thetay)
2.2361
116.5651
c. The real and imaginary components of this complex number are shown in Figure C.5.
Im
206.6
2 Re
26.6
153.4(Measured
5 Clockwise)
1

Figure C.5. The components of 2 j

Then,

1
j tan
1
------
2 2 2 j206.6 j ( 153.4 )
2 j 1 = 2 +1 e = 5e = 5 206.6 = 5e = 5 153.4

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition C7
Copyright Orchard Publications
Appendix D
Matrices and Determinants

T
his appendix is an introduction to matrices and matrix operations. Determinants, Cramers
rule, and Gausss elimination method are reviewed. Some definitions and examples are not
applicable to the material presented in this text, but are included for subject continuity, and
academic interest. They are discussed in detail in matrix theory textbooks. These are denoted with
a dagger () and may be skipped.

D.1 Matrix Definition

A matrix is a rectangular array of numbers such as those shown below.

1 3 1
2 3 7 or 2 1 5
1 1 5
4 7 6

In general form, a matrix A is denoted as

a 11 a 12 a 13 a 1 n
a 21 a 22 a 23 a 2 n
A = a 31 a 32 a 33 a 3 n (D.1)

a m 1 a m 2 a m 3 a mn

The numbers a ij are the elements of the matrix where the index i indicates the row, and j indi-
cates the column in which each element is positioned. For instance, a 43 indicates the element posi-
tioned in the fourth row and third column.
A matrix of m rows and n columns is said to be of m n order matrix.
If m = n , the matrix is said to be a square matrix of order m (or n ). Thus, if a matrix has five rows
and five columns, it is said to be a square matrix of order 5.
In a square matrix, the elements a 11, a 22, a 33, , a nn are called the main diagonal elements.
Alternately, we say that the matrix elements a 11, a 22, a 33, , a nn , are located on the main diago-
nal.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition D 1
Copyright Orchard Publications
Appendix D Matrices and Determinants

85 R1 C j100
1700 +
IX
V1 V2
R3
VS 100
50
j200 L R2

Figure D.4. Circuit for Example D.18

the current I X can be found from the relation
V1 V2
I X = ------------------- (D.59)
R3

and the voltages V 1 and V 2 can be computed from the nodal equations

V 1 170 0 V 1 V 2 V 1 0
-------------------------------
- + ------------------- + --------------- = 0 (D.60)
85 100 j200
and
V 2 170 0 V 2 V 1 V 2 0
-------------------------------
- + ------------------- + --------------- = 0 (D.61)
j100 100 50

Compute, and express the current I x in both rectangular and polar forms by first simplifying like
terms, collecting, and then writing the above relations in matrix form as YV = I , where
Y = Admit tan ce , V = Voltage , and I = Current
Solution:
The Y matrix elements are the coefficients of V 1 and V 2 . Simplifying and rearranging the nodal
equations of (D.60) and (D.61), we obtain
( 0.0218 j0.005 )V 1 0.01V 2 = 2
(D.62)
0.01 V 1 + ( 0.03 + j0.01 )V 2 = j1.7

Next, we write (D.62) in matrix form as

0.0218 j0.005 0.01 V1 2

0.01 0.03 + j0.01 V2 = j1.7 (D.63)

Y V I

D26 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Solution of Simultaneous Equations with Matrices

where the matrices Y , V , and I are as indicated.

We will use MATLAB to compute the voltages V 1 and V 2 , and to do all other computations. The
script is shown below.
Y=[0.02180.005j 0.01; 0.01 0.03+0.01j]; I=[2; 1.7j]; V=Y\I; % Define Y, I, and find V
fprintf('\n'); % Insert a line
disp('V1 = '); disp(V(1)); disp('V2 = '); disp(V(2)); % Display values of V1 and V2
V1 = 1.0490e+002 + 4.9448e+001i
V2 = 53.4162 + 55.3439i
Next, we find I X from
R3=100; IX=(V(1)V(2))/R3 % Compute the value of IX

IX = 0.5149 - 0.0590i
This is the rectangular form of I X . For the polar form we use the MATLAB script
magIX=abs(IX), thetaIX=angle(IX)*180/pi % Compute the magnitude and the angle in degrees
magIX = 0.5183
thetaIX = -6.5326
Therefore, in polar form
I X = 0.518 6.53

We can also find the current I X using a Simulink model, and to simplify the model we first derive
the Thevenin equivalent in Figure D.4 to that shown in Figure D.5.
By application of Thevenins theorem, the electric circuit of Figure 5.45(a) can be simplified to
that shown in Figure 5.45(b).
X
112 j10.6

100

VTH = 110j6.87 V
Y
Figure D.5. Electric circuit (a) and its Thevenin equivalent (b)
Next, we let V IN = V TH , V OUT = V XY , Z 1 = 112 + j10 , and Z 2 = 100 . Application of the volt-
age division expression yields
Z2
V OUT = ------------------- V IN (D.64)
Z1 + Z2

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition D27
Copyright Orchard Publications
Appendix D Matrices and Determinants

Now, we use the model in Figure D.6 to convert all quantities from the rectangular to the polar
form, perform the addition and multiplication operations, display the output voltage in both polar
and rectangular forms, and show the output voltage on a Scope block in Figure D.7 The Simulink
blocks used for the conversions are discussed in Introduction to Simulink with Engineering Appli-
cationsa, ISBN 9781934404218, Appendix T, Simulink Extras.

Figure D.6. Model for the computation and display of the output voltage for the circuit in Figure 5.45(b)

Figure D.7. Waveform for the output voltage of model in Figure D.6

Spreadsheets have limited capabilities with complex numbers, and thus we cannot use them to com-
pute matrices that include complex numbers in their elements as in Example D.18.

D28 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Appendix E
Window Functions

T
his appendix is an introduction to window functions. The rectangular, triangular, Hanning,
Hamming, Blackman, and Kaiser windows are discussed. An example using each is provided
for illustration of its use.

E.1 Window Function Defined

A window function is a function that is zerovalued outside of some chosen interval. For instance, a
function that is constant inside the interval and zero elsewhere is called a rectangular window, and
describes the shape of its graphical representation. When another function or a signal (data) is mul-
tiplied by a window function, the product is also zerovalued outside the interval: all that is left is
the view through the window. Applications of window functions include spectral analysis, and fil-
ter design.
When selecting an appropriate window function for an application, a comparison graph may be use-
ful. The most important parameter is usually the stop band attenuation close to the main lobe.
All of the window functions that we will discuss are even functions of time when centered at the ori-
gin.

E.2 Common Window Functions

Based on the discussion in the previous section, it appears that a rectangular function would be the
ideal window function to terminate an impulse response with an infinite number of terms. For
instance, let us assume that the impulse response h [ n ] shown in Figure E.1(a) below converges uni-
formly and is represented by a portion of the amplitude response A ( f ) shown in Figure E.1(b).

h[n] A(f)

t f
(a) (b)
Figure E.1. Impulse response with an infinite number of terms

Next, let us assume that the impulse response h [ n ] is terminated abruptly without changing any of
its coefficients, as shown in Figure E.2(a). In this case, the resulting amplitude response A' ( f ) will
be subject to undesired oscillations and poor convergence as shown in Figure E.2(b).

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition E1
Copyright Orchard Publications
Appendix E Window Functions

-50

-100

-150
0 0.2 0.4 0.6 0.8 1

Figure E.10. Normalized frequency domain plot for the triangular window function created with MATLAB
We can also use the MATLAB Window Visualization Tool. The script below generates the plot
shown in Figure E.11.
L=50; wvtool(triang(L))

Time domain Frequency domain

40

1
20

0
0.8

-20
Magnitude (dB)
Amplitude

0.6
-40

0.4 -60

-80
0.2
-100

0 -120
5 10 15 20 25 30 35 40 45 50 0 0.2 0.4 0.6 0.8
Samples Normalized Frequency ( rad/sample)

Figure E.11. Triangular window function generated with the MATLAB Window Visualization Tool

E.2.3 Hanning Window Function

The Hanning* window function is defined as

* This window function is also known as Hann window function or cosinesquared window function.

E6 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Appendix F
Correlation Functions

T
his appendix is an introduction to cross correlation and autocorrelation functions. Cross cor-
relation applies to two different signals and provides of measure of the similarity between one
signal and a timedelayed version of the other signal. Autocorrelation provides a measure of
the similarity between a given signal and its own timedelayed version.

F.1 Cross Correlation

The cross correlation is similar in nature to the convolution* of two functions. Whereas convolution
involves reversing a signal, then shifting it and multiplying by another signal, correlation only
involves shifting it and multiplying. The cross correlation of two signals f 1 ( t ) and f 2 ( t ) denoted as
R 12 ( ) where represents a time delay is defined as

R 12 ( ) =
f 1 ( t )f 2 ( t ) dt (F.1)

where f 2 denotes the complex conjugate of f 2 . The cross correlation R 12 ( ) is an indication of

how well the signals f 1 and f 2 are correlated as the delay varies. If the cross correlation function
peaks for a particular value of the delay , it would indicate a very good correlation, which means
that the two signals match each other very well. Conversely, a very small or zero value of the cross
correlation function would indicate little or no correlation.
In relation (F.1) above the assumes that the time t represents a dummy variable and it is eliminated
in the integration, and thus the cross correlation R 12 ( ) is a function of the delay , whereas in the
convolution integral the delay is the dummy variable. Also, in (F.1) the signal f 2 ( t ) is shifted by
the distance in the positive direction, and if we choose to shift the signal f 1 ( t ) by a distance in
the negative direction, relation (F.1) is expressed as

R 12 ( ) =

f 1 ( t + )f 2 ( t ) dt (F.2)

In probability and statistics, the cross correlation is also referred to as the covariance cov(X, Y)
between two random vectors X and Y.

* Convolution is discussed in Chapter 6.

For a detailed discussion on covariance, please refer to Mathematics for Business, Science, and Technology, ISBN
9781934404010.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition F1
Copyright Orchard Publications
Appendix F Correlation Functions

T
f 1 t + --- T A T 1 7
R 12 ( ) = R 12 --- = Area = A T ---- --- --- = --- A T
2 A 2 2 2 2 8

f2 ( t )

T T/2 T
t
0
(a)

f1 ( t + 0 ) 1
A R 12 ( ) = R 12 ( 0 ) = Area = --- A T
2

f2 ( t )

T T/2 T
t
0
(b)

T
f 1 t ---
2 A
1 A T 1 AT
R 12 ( ) = R 12 --- = Area = ---- --- --- = ------------
2 2 2 2 8
f2 ( t )

T T/2
t
0 T
(c)
Figure F.3. Graphical cross correlation for Example F.2

The cross correlation R 12 ( ) varies with as shown in Figure F.4 where the plot was created with
the MATLAB script below with A = 1 and T = 1 .
x=[-1 -0.5 -0.25 0 0.25 0.5 1]; y=[0 0.25 0.375 0.5 0.375 0.25 0]; p3=polyfit(x,y,5);
x_span=(-1:0.05:1); p3_pol=polyval(p3,x_span); plot(x_span,p3_pol); box off

0.5

R 12 ( )

0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Figure F.4. The cross correlation function for Example F.2

F4 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Appendix F Correlation Functions
F.2 Autocorrelation
Autocorrelation is a timedomain function. It is a measure of how much a signal shape, or wave-
form, resembles a delayed version of itself. It is closely related to the cepstrum.* The value of auto-
correlation can vary between zero and one. A periodic signal, such as a sine wave has an autocorre-
lation which is equal to one at zero time delay, zero at a time delay of onehalf the period of the
wave, and one at a time delay of one period; in other words, it is a sinusoidal wave form itself. Ran-
dom noise has an autocorrelation of one at zero delay, but is essentially zero at all other delays.
Autocorrelation is often used to extract periodic signals from noise.
The autocorrelation of a signal x denoted as R xx ( ) where represents a time delay is a special
case of the cross correlation function defined in relation (T.1) above with y ( t ) = x ( t ) , and it is
defined as
tp
1
R xx ( ) = lim ---
tp t p 0
x ( t )x ( t ) dt (F.5)

and for real signals tp

R xx ( ) = lim --1- x ( t )x ( t ) dt (F.6)
tp t p 0

and this relation states that the autocorrelation is the average of the product of the signal and a
delayed version of the signal as a function of the delay.
Two important properties of the autocorrelation function are:
Property 1:
The autocorrelation is an even functions, i.e.,
Rf ( ) = Rf ( ) when f is a real function (F.7)

R f ( ) = R f ( ) when f is a complex function (F.8)

Property 2:
The maximum value of the autocorrelation function occurs at the origin, i.e.,
[ Rf ( ) ] Rf ( 0 ) for all (F.9)

Example F.4
Determine and sketch the autocorrelation function
at
f( t) = e u0 ( t ) a>0 (F.10)
For the Simulink model use a = 0.5

* Cepstrum is the Fourier transform of the logmagnitude spectrum.

F6 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
 Autocorrelation
Solution:
The autocorrelation function R XX ( ) was created with the MATLAB script
A=100; x=randn(1,A); N=50; alpha=0.9; y=filter([1 zeros(1,N) alpha],1,x);
Rxx=conv(y,fliplr(y)); plot([-A+1:A-1],Rxx); grid
and it is shown in Figure F.8.
150

100

50

-50
-100 -80 -60 -40 -20 0 20 40 60 80 100

Figure F.12. The autocorrelation function of the random noise signal in Figure F.7

Example F.6
For the model in Figure F.13 the Signal Builder block contains the square pulse shown in Figure
F.14, and its autocorrelation function is as shown in Figure 15.

Figure F.13. Model for Example F.6

Figure F.14. Pulse signal in Signal Builder block in Figure F.13

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition F9
Copyright Orchard Publications
Appendix G
Nonlinear Systems

T
his appendix is an introduction to nonlinear systems. Whereas linear systems are expressed in
transfer functions, nonlinear systems are expressed in describing functions defined as the
ratio of the fundamental (first harmonic) frequency of the output response of a nonlinear
device to the amplitude of a sinusoidal input signal to that device. The describing function of a sat-
uration curve is derived.

G.1 Describing Functions

A nonlinear system can be best defined by first defining a linear system.
A linear system satisfies the following condition:
Let c be an arbitrary constant, real or complex. If an input cx1 ( t ) produces an output cy1 ( t ) , and
an input cx 2 ( t ) produces an output cy 2 ( t ) , an input cx 1 ( t ) + cx 2 ( t ) produces an output
cy 1 ( t ) + cy 2 ( t ) . All other systems are considered nonlinear. In fact, all physical systems contain non
linearities and time varying parameters to a certain degree. However, in some cases nonlinear
devices such as relays are introduced in control systems to provide a more effective operation.
Nonlinear systems are in most cases difficult to analyze. The Laplace transformation and the trans-
fer function concept do not apply to nonlinear systems. However, we can apply linearization, a pro-
cess in which a nonlinear state equation can be expanded into a Taylor series about a desired operat-
ing point where only the linear terms of the Taylor series are kept, and all others are discarded.
Nonlinear systems are generally described by describing functions defined as the ratio of the funda-
mental (first harmonic) frequency of the output response of a nonlinear device to the amplitude of
a sinusoidal input signal to that device. Thus, if the fundamental component of the Fourier series of
j 1
the output is C 1 e , and A is the amplitude of the sinusoidal input, the describing function is
j 1
defined as C 1 e A . The procedure is best described with the following example.

The nonlinear characteristics of saturation are shown in Figure G.1

From Figure G.1,
S<A (G.1)
and

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition G1
Copyright Orchard Publications
 Describing Functions
2

b 1 = ----
0
m ( t ) sin t dt (G.7)

and with G.2 we obtain

t1 t2 t3 t4 2

b 1 = ----

0
A sin 2 t dt + ----
t1
S sin t dt + ----
t2
A sin 2 t dt ----
t3
S sin t dt + ----
t4
A sin 2 t dt

but because of the symmetry, the expression for b 1 above can be simplified by integrating over a
quarter of a cycle and multiplying by 4. Thus,
t1 2
4 4
b 1 = -------
0
A sin 2 t dt + -------
t1
S sin t dt (G.8)

Integration of (G.8) yields

2 A
b 1 = --- At 1 ---- sin 2t 1 + 2S cos t 1 (G.9)
2

Substitution of (G.3) into (G.9) and after simplification results in the relation
1 1 1 1
2 S S 2A S S S
b 1 = --- A sin ---- + S cos sin ---- = ------- sin ---- + ---- cos sin ---- (G.10)
A A A A A

and thus the describing function is

1 1
b 2 S S S
-----1 = --- sin ---- + ---- cos sin ---- (G.11)
A A A A

We observe that the describing function for saturation is only a function of the amplitude of the
input and independent of frequency.
The describing function in (G.11) is modeled in Figure G.2.

Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition G3
Copyright Orchard Publications
Appendix G Nonlinear Systems

Note: The XY Graph block does not display

grid lines. They were added manually.

Figure G.2. Model for the describing function in (G.11)

The waveform displayed by the XY Graph block shows the normalized value of the describing
function (Y axis) as a function of the ratio S A (X Axis).

G4 Signals and Systems with MATLAB Computing and Simulink Modeling, Fifth Edition
Copyright Orchard Publications
Index
Symbols Cauer filter - see elliptic filter convolution property of the Fourier
Cayley-Hamilton theorem 5-10 transform - see Fourier transform
% (percent) symbol in MATLAB A-2 characteristic equation 5-18 properties of
cheb1ap MATLAB function 11-35 Cooley and Tukey 10-18
A cheb2ap MATLAB function 11-38 cosine function - Fourier transform of
Chebyshev filters - see filter see Fourier transform of
abs(z) in MATLAB A-22 Chebyshev Type I analog low-pass common functions
active analog filter - see filter filter design - see filter design cosw0t u0(t) function - Fourier transform of
adjoint of a matrix - see matrix Chebyshev Type I filters - see filter see Fourier transform of
admittance Chebyshev Type I low-pass filter common functions
capacitive 4-2 magnitude-square function Cramers rule D-15
inductive 4-2 - see filter cross correlation F-1
complex input 4-11 Chebyshev Type II analog low-pass
algebraic constrain block in Simulink B-18 filter design - see filter design D
aliasing 10-13 Chebyshev Type II filters - see filter
all-pass filter - see filter circuit analysis with Laplace transforms 4-1 d2c MATLAB function 9-43
all-pole approximation 11-20 circuit analysis with state variables 5-21 data points in MATLAB A-13
all-pole low-pass filter circuit transformation from time decade - definition of 11-11
see filter - low-pass to complex frequency 4-1 decimation in frequency
alternate form of the trigonometric clc MATLAB command A-2 see FFT algorithm
Fourier series - see Fourier series clear MATLAB command A-2 decimation in time - see FFT algorithm
alternate method of partial cofactor of a matrix - see matrix deconv MATLAB function A-6
fraction expansion - see partial collect(s) MATLAB symbolic function 3-12 default color in MATLAB A-14
fraction expansion column vector in MATLAB A-17 default line in MATLAB A-14
angle(z) MATLAB function A-22 command screen in MATLAB A-1 default marker in MATLAB A-14
argument 11-2 command window in MATLAB A-1 default values in MATLAB A-11
attenuation rate 11-11 commas in MATLAB A-7 delta (impulse) function
autocorelation F-6 comment line in MATLAB A-2 definition of 1-8
autoscale icon in Simulink B-12 Commonly Used Blocks in Simulink B-7 doublet 1-13
axis MATLAB command A-15 complex admittance - see admittance Fourier transform of - see Fourier
complex conjugate in MATLAB A-4 transform of common functions
B complex conjugate pairs 3-5 higher order 1-13
complex impedance - see impedance nth-order 1-13
band-elimination filter - see filter complex number C-2 sampling property of 1-11
band-elimination filter design complex numbers in MATLAB A-2 sifting property of 1-12
see filter design complex poles 3-5 triplet 1-13
band-limited signal 10-13 Complex to Magnitude-Angle demo in MATLAB A-2
band-pass filter - see filter block in Simulink C-7 DeMoivres theorem 11-14
band-pass filter design computation of the state derivation of the Fourier transform from
see filter design transition matrix 5-10 the Laplace Transform 8-24
band-stop filter - see filter computation of the Z Transform with describing function G-1
Bessel filter - see filter contour integration - see Z transform determinant of a matrix
bilinear MATLAB function 11-59 Configuration Parameters see matrix
bilinear transformation - see in Simulink B-12 DFT - common properties of
transformation methods for mapping congugate of a matrix - see matrix even time function 10-9
analog prototype filters to digital filters conj(A) MATLAB function D-8 even frequency function 10-9
Blackman window function conjugate of a complex number C-3 frequency convolution 10-13
see window functions conjugate time and frequency functions frequency shift 10-12
bode MATLAB function 11-23 of the Fourier transform - see linearity 10-11
box MATLAB command A-12 Fourier transform - properties of odd time function 10-9
buttap MATLAB function 11-17 constant function - Fourier transform of odd frequency function 10-9
buttefly operation 10-23 see Fourier transform - properties of time convolution 10-12
Butterworth analog low-pass Contents Pane in Simulink B-7 time shift 10-11
filter design - see filter design contour integral 9-17 DFT - definition of 10-1
conv MATLAB function A-6 N-point 10-2
C convolution in the complex frequency diagonal elements of a matrix - see matrix
domain - see Laplace transform diagonal matrix - see matrix
c2d MATLAB function 9-42 properties of differentiation in complex frequency
capacitive admittance - see admittance convolution in the discrete-frequency domain property of the Laplace
capacitive impedance - see impedance domain - see Z transform - properties of transform - see Laplace transform
cascade form realization - see digital filter convolution in the discrete-time - properties of
Category I FFT algorithm - see domain - see Z transform - properties of differentiation in time domain property
FFT algorithm convolution in the time domain property of the Laplace transform - see
Category II FFT algorithm - see of the Laplace transform - see Laplace transform - properties of
FFT algorithm Laplace transform - properties of differentiation property of the Fourier
Cauchys residue theorem convolution integral defined 6-7 transform 8-12 - see Fourier
see residue theorem graphical evaluation of 6-8 transform - properties of

IN-1
digital filter 11-1, 11-49, 11-67 exponential form of complex numbers C-5 Finite Impulse Response (FIR) digital filter
Finite Impulse Response (FIR) 11-49 exponential form of the Fourier series - see digital filter
Infinite Impulse Response (IIR) 11-49 see Fourier series FIR - see digital filter
realization of exponential order function - first harmonic - see Fourier series
Direct Form I 11-67 definition of 2-2 harmonics of
Direct Form II 11-68 eye(n) MATLAB function D-6 first-order circuit 5-1
cascade (series) form 11-70 eye(size(A)) MATLAB function D-7 first-order simultaneous
non-recursive 11-49 differential equations 5-1
parallel form 11-70 F Flip Block command in Simulink B-11
recursive 11-49 format in MATLAB A-29
series (cascade) form 11-70 factor(s) MATLAB symbolic function 3-4 fourier MATLAB command 8-33
Digital Filter Design Simulink block 11-70 Fast Fourier Transform (FFT) 10-1, 10-17 Fourier series
digital filter design with Simulink 11-67 FDA Tool Digital Filter Design exponential form of 7-30
dimpulse MATLAB function 9-26 Simulink block 11-78 method used in window functions E-15
Dirac MATLAB function 1-20 FFT algorithm trigonometric form of 7-2
Direct Form I realization - see digital filter Category I 10-19 alternate form of 7-24
realization of Category II 10-19 Fourier series coefficients - evaluation of
Direct Form II realization - see digital filter decimation in frequency 10-19 numerical evaluation using Excel 7-43
realization of decimation in time 10-19 numerical evaluation using MATLAB 7-44
direct term in MATLAB 3-3 double-memory technique 10-19 Fourier series of common waveforms
discontinuous function - definition of 1-2 in-place 10-19 full-wave rectifier 7-20, 7-23
Discrete Fourier Transform (DFT) 10-1 natural input-output 10-19 half-wave rectifier 7-17, 7-20
discrete impulse response 9-36 FFT definition of 10-1, 10-17 square waveform
discrete-time system transfer function 9-36 fft(x) MATLAB function 10-5, 11-65 with even symmetry 7-9, 7-13
discrete unit step function 9-3 figure window in MATLAB A-13 with odd symmetry 7-8, 7-12
discrete-time exponential sequence 9-15 filter - see also digital filter sawtooth 7-9, 7-15
discrete-time systems 9-1 active triangular 7-9, 7-16
discrete-time unit ramp function 9-16 high-pass 4-22, 4-30 Fourier series - harmonics of
discrete-time unit step function 9-14 low-pass 4-22, 4-30 first 7-1, 7-10
Display block in Simulink B-18 all-pass 11-1, 11-89 second 7-1, 7-10
display formats in MATLAB A-29 all-pole 11-20 third 7-1, 7-10
distinct eigenvalues - see eigenvalues band-elimination 11-1, 11-8 Fourier series - symmetry
distinct poles - see poles band-pass 11-1, 11-6 even 7-7
division of complex numbers C-4 band-stop - see band-elimination half-wave 7-7, 7-33
dot operator in MATLAB Bessel 11-90 odd 7-7
division with A-19 Chebyshev 11-10 quarter-wave 7-7 (footnote)
exponentiation with A-19 Inverted 11-36 types of 7-7
multiplication with A-19 magnitude-square function 11-10 Fourier integral - see Fourier transform
double-memory technique prototype 11-10 Fourier transform
see FFT Algorithm Type I 11-24 definition of 8-1
doublet - see delta function Type II 11-36 inverse of 8-1
elliptic 11-37 special forms of 8-2
E high-pass 4-30, 11-1, 11-4 Fourier transform - properties of
low-pass 4-30, 11-1, 11-2 area under f(t) 8-15
Editor Window in MATLAB A-1 low-pass analog filter area under F(w) 8-15
Editor/Debugger in MATLAB A-1 prototypes 11-10 conjugate time and frequency
eig(x) MATLAB function 5-17 maximally flat 11-13 (footnote) functions 8-13
eigenvalues notch (band-elimination) 11-8 constant function 8-18
distinct 15-1 phase shift 11-1, 11-89 frequency convolution 8-14
multiple (repeated) 5-15 RC high-pass 11-4 frequency differentiation 8-12
eigenvector 5-19 RC low-pass 11-2 frequency shifting 8-11
e-jwt u0(t) Fourier transform of - see Fourier RLC band-elimination 4-31, 11-8 imaginary time functions - Fourier
transform of common functions RLC band-pass 4-31, 11-6 transform of 8-6
element-by-element operation in MATLAB filter design - see also digital filter linearity 8-9
division A-19 band-elimination 11-39 Parseval's theorem 8-15
exponentiation A-19 band-pass 11-39 real time functions - Fourier
multiplication A-19 Butterworth analog low-pass 11-13 transform of 8-3
elements of the matrix - see matrix Chebyshev symmetry 8-9
ellip MATLAB function 11-40 Type I 11-24 time convolution 8-14
elliptic filter - see filter Type II 11-36 time differentiation 8-12
elliptic filter design - see filter design elliptic 11-37 time integration 8-13
eps in MATLAB A-21 high-pass 11-39 time scaling 8-10
Eulers identities C-4 low-pass 11-10 time shifting 8-11
even functions 6-4, 7-32 filter MATLAB function 11-60 Fourier Transform derivation from
even symmetry - see Fourier final value theorem in Laplace transform Laplace transform 8-25
series - symmetry see Laplace transform - properties of Fourier transform of common functions
Excel's Analysis ToolPak 10-6 final value theorem in Z transform cosw 0t 8-18
exit MATLAB command A-2 see Z transform - properties of cosw0t u0(t) 8-23
expand MATLAB symbolic function 3-10 find MATLAB function 11-64 delta (d(t) and d(t-a)) 8-16

IN-2
 e-jw0t 8-19 ifourier MATLAB function 8-31 integration in
e-jw0t u0(t) 8-23 IIR - see digital filter complex frequency domain 2-7
signum (sgn(t)) 8-19 ilaplace MATLAB function 3-4 time domain 2-6
sinw0t 8-19 imag(z) MATLAB function A-22 linearity 2-3
sinw0t u0(t) 8-23 imaginary axis - definition of C-2 scaling 2-4
unit step (u0(t)) 8-21 imaginary number - definition of C-1 time periodicity 2-8
Fourier transform of common waveforms imaginary time functions 8-6 time shift 2-3
combined rectangular pulses 8-28 see Fourier transform - properties of Laplace transform of common waveforms
cosine within a rectangular pulse 8-29 impedance full-rectified 2-30
shifted rectangular pulse 8-27 capacitive 4-2 half-rectified 2-25
symmetrical rectangular pulse 8-26 inductive 4-2 linear segment 2-22
periodic time functions 8-29 complex input 4-8, 4-9 pulse 2-21
fourth-order all-pole low-pass filter improper integral - definition of 2-15 rectangular periodic 2-24
see filter, low-pass improper rational function - sawtooth 2-30
fplot in MATLAB A-25 definition of 3-1, 3-12 triangular 2-23
frequency convolution in DFT impulse function - see delta function Laplace transform of common functions
see DFT - common properties of impulse invariant method - see transform of e-at coswt u0(t) 2-20
frequency convolution in Fourier transformation methods for mapping transform of e-at sinwt u0(t) 2-20
transform of - see Fourier transform analog prototype filters to transform of e-at u0(t) 2-18
- properties of digital filters transform of coswt u0(t) 2-19
frequency differentiation in Fourier increments between points transform of d(t) 2-17
transform of - see Fourier transform in MATLAB A-14 transform of d(t-a) 2-17
- properties of inductive admittance - see admittance transform of sinwt u0(t) 2-19
frequency shift in DFT inductive impedance - see impedance transform of tn u0(t) function 2-14
see DFT - common properties of infinite impulse response - see digital filter transform of tn e-at u0(t) 2-18
frequency shift in Fourier transform initial value theorem in Laplace transform transform of u0(t) 2-12
see Fourier transform - properties of see Laplace transform - properties of transform of u1(t) 2-12
frequency shift in Laplace transform initial value theorem in Z transform leakage 10-14
see Laplace transform - properties of see Z transform - properties of left shift in discrete-time domain
freqz MATLAB function 11-54 in-place FFT algorithm 10-20 see Z transform - properties of
full rectification waveform - Laplace see FFT algorithm Leibnitzs rule 2-6
transform of - see Laplace transform integration in frequency in Laplace transform lims= MATLAB command A-25
of common waveforms see Laplace transform - properties of line spectra 7-35
full-wave rectifier - Fourier series of - see integration in time in Laplace transform linear difference equation 9-35
Fourier series of common waveforms see Laplace transform - properties of linearity in DFT
Function Block Parameters in Simulink B-10 Inverse Fourier transform see DFT - common properties of
function files in MATLAB A-25 see Fourier transform linearity in discrete-time domain
fundamental frequency 7-1 Inverse Laplace transform see Z transform - properties of
fzero MATLAB function A-25 see Laplace transform linearity property in Fourier transform
inverse of a matrix - see matrix see Fourier transform - properties of
G Inverse Z transform - see Z transform linearity property in Laplace transform
inversion integral 9-29 see Laplace transform - properties of
Gain block in Simulink B-9, B-18 Inverted Chebyshev filter - see filter linspace MATLAB command A-13
gamma function 2-15 ln (natural log) A-12
Gaussian elimination method D-17 J log in MATLAB A-12
generalized factorial function 2-14 log(x) MATLAB function A-12
Gibbs phenomenon 7-23 j operator C-1 log10(x) MATLAB function A-12
grid MATLAB command A-11 log2(x) MATLAB function A-12
gtext MATLAB command A-13 K loglog MATLAB command A-12
long division of polynomials 9-33
H Kaiser window - see windows lower triangular matrix - see matrix
low-pass analog filter prototypes - see filter
half-wave rectifier - Fourier series of - see L low-pass filter - see filter
Fourier series of common waveforms lp2bp MATLAB function 11-41
half-wave symmetry L Hpitals rule 2-15 lp2bs MATLAB function 11-46
see Fourier series - symmetry laplace MATLAB function 2-26 lp2hp MATLAB function 11-41
Hamming window - see windows Laplace integral - see Laplace transform lp2lp MATLAB function 11-41
Hanning window - see windows Laplace transform
Heaviside MATLAB function 1-19 definition of 2-1 M
help in MATLAB A-2 Inverse of 2-1, 3-1
Hermitian matrix - see matrix Laplace transform - properties of magnitude-squared function 11-10
higher order delta functions - see delta convolution in the complex main diagonal of a matrix - see matrix
function frequency domain 2-11 Math Operations in Simulink B-11
high-pass filter - see filter convolution in the time domain 2-11 MATLAB Demos A-2
high-pass filter design - see filter design differentiation in complex MATLABs Editor/Debugger A-1
frequency domain 2-6 matrix (matrices)
I differentiation in time domain 2-4 adjoint of D-19
final value theorem 2-10 cofactor of D-11
identity matrix - see matrix frequency shift 2-4 conformable for addition D-2
ifft(x) MATLAB function 10-5 initial value theorem 2-9 conformable for subtraction D-2

IN-3
 conformable for multiplication D-4 Nyquist frequency 10-13 recursive realization digital filter
congugate of D-8 see digital filter
definition of D-1 O region of
determinant D-9 convergence 9-3
minor of D-11 octave defined 11-11 divergence 9-3
non-singular D-19 odd functions 6-11, 7-333 relationship between state equations
singular D-19 odd symmetry - see Fourier and Laplace Transform 5-28
diagonal D-1 series - symmetry residue 3-3, 9-37
diagonal elements of D-1 orthogonal functions 7-2 residue MATLAB function 3-3, 3-12
elements of D-1 orthogonal vectors 5-19 residue theorem 9-19
Hermitian D-8 orthonormal basis 5-19 right shift in the discrete-time domain
identity D-6 see Z transform - properties of
inverse of D-220 P RLC band-elimination filter - see filter
left division in MATLAB D-23 RLC band-pass filter - see filter
multiplication in MATLAB A-17 parallel form realization - see digital filter roots of polynomials in MATLAB A-3
power series of 5-9 Parsevals theorem - see roots(p) MATLAB function 3-5, A-3
scalar D-6 Fourier transform - properties of round(n) MATLAB function A-22
size of D-7 partial fraction expansion 3-1 row vector in MATLAB A-3
skew-Hermitian D-9 alternate method of 3-14 Runge-Kutta method 5-1
skew-symmetric D-8 method of clearing the fractions 3-14 running Simulink B-7
square D-1 phase angle 11-2
symmetric D-8 phase shift filter - see filter S
trace of D-2 picket-fence effect 10-14
transpose of D-7 plot MATLAB command A-9 sampling property of the delta function
triangular polar form of complex numbers C-5 see delta function
lower D-6 polar plot in MATLAB A-23 sampling theorem 10-13
upper D-7 polar(theta,r) MATLAB function A-22 sawtooth waveform - see Laplace
zero D-2 poles 3-1 transform of common waveforms
matrix left division in MATLAB - see matrix complex 3-5 sawtooth waveform - Fourier series of
matrix multiplication in MATLAB - see matrix distinct 3-2 see Fourier series of
matrix power series - see matrix multiple (repeated) 3-7 common waveforms
maximally flat filter - see filter poly MATLAB function A-4 scalar matrix - see matrix
mesh(x,y,z) MATLAB function A-15 polyder MATLAB function A-6 scaling property of the Laplace transform
meshgrid(x,y) MATLAB command A-16 polynomial construction from see Laplace transform - properties of
m-file in MATLAB A-1, A-224 known roots in MATLAB A-4 Scope block in Simulink B-12
minor of determinant - see matrix polyval MATLAB function A-5 script file in MATLAB A-2, A-24
MINVERSE Excel function D-25 pre-sampling filter 10-13 second harmonic - see Fourier series
MMULT Excel function D-25 pre-warping 11-52 harmonics of
modulated signals 8-11 proper rational function - semicolons in MATLAB A-7
multiple eigenvalues - see eigenvalues definition of 3-1 semilogx MATLAB command A-12
multiple poles - see poles properties of the DFT semilogy MATLAB command A-12
multiplication by an in discrete-time domain see DFT - common properties of series form realization - see digital filter
see Z transform - properties of properties of the Fourier Transform Shannons sampling theorem
multiplication by e-naT in discrete-time see Fourier transform - properties of see sampling theorem
domain - see Z transform - properties of properties of the Laplace Transform shift of f[n] u0[n] in discrete-time domain
multiplication by n in discrete-time domain see Laplace transform - properties of see Z transform - properties of
see Z transform - properties of properties of the Z Transform sifting property of the delta function
multiplication by n2 indiscrete-time domain see Z transform - properties of see delta function
see Z transform - properties of signal flow graph 10-22
multiplication of complex numbers C-2 Q signals described in math form 1-1
signum function - see Fourier transform
N quarter-wave symmetry - see of common functions
Fourier series - symmetry simout To Workspace block
NaN in MATLAB A-25 quit MATLAB command A-2 in Simulink B-13
natural input-output FFT algorithm simple MATLAB symbolic function 3-6
see FFT algorithm R Simulation drop menu in Simulink B-12
network transformation simulation start icon in Simulink B-12
resistive 4-1 radius of absolute convergence 9-3 Simulink icon B-7
capacitive 4-1 ramp function 1-9 Simulink Library Browser B-8
inductive 4-1 randn MATLAB function 11-65 sine function - Fourier transform of
non-recursive realization digital filter Random Source Simulink block 11-76 see Fourier transform of
see digital filter rationalization of the quotient C-4 common functions
non-singular determinant - see matrix RC high-pass filter - see filter singular determinant - see matrix
nonlinear system G-1 RC low-pass filter - see filter Sinks library in Simulink B-18
normalized cutoff frequency 11-14 real axis C-1 sinw0t u0(t) Fourier transform of - see
notch filter - see filter real number C-2 Fourier transform of common functions
N-point DFT - see DFT - definition of real(z) MATLAB function A-22 size of a matrix - see matrix
nth-order delta function - see delta function rectangular form C-5 skew-Hermitian matrix - see matrix
numerical evaluation of Fourier coefficients rectangular pulse expressed in terms skew-symmetric matrix - see matrix
see Fourier series coefficients of the unit step function 1-4 special forms of the Fourier transform

IN-4
 see Fourier transform Transfer Fcn block in Simulink 4-17 convolution in the discrete
spectrum analyzer 7-35 Transfer Fcn Direct Form II time domain 9-8
square matrix - see matrix Simulink block 11-68 final value theorem 9-9
square waveform with even symmetry - see transfer function of initial value theorem 9-9
Fourier series of common waveforms continuous-time systems 4-13 left shift 9-5
square waveform with odd symmetry - see discrete-time systems 9-35 linearity 9-3
Fourier series of common waveforms transformation between multiplication by an 9-6
ss2tf MATLAB function 5-31 s and z domains 9-20 multiplication by e-naT 9-6
stability 11-13 transformation methods for mapping multiplication by n 9-6
start simulation in Simulink B-12 analog prototype filters to digital filters multiplication by n2 9-6
state equations Impulse Invariant Method 11-50 right shift 9-4
for continuous-time systems 5-1 Step Invariant Method 11-50 shift of f[n] u0[n] 9-3
for discrete-time systems 9-40 Bilinear transformation 11-50 summation 9-7
state transition matrix 5-8 transpose of a matrix - see matrix Z Transform of discrete-time functions
state variables Tree Pane in Simulink B-7 cosine function cosnaT 9-15
for continuous-time systems 5-1 triangular waveform expressed in terms exponential sequence e-naT u0[n] 9-15
for discrete-time systems 9-40 of the unit step function 1-6 geometric sequence an 9-12
State-Space block in Simulink B-13 triplet - see delta function sine function sinnaT 9-15
state-space equations Tukey - see Cooley and Tukey unit ramp function nu0[n] 9-16
for continuous-time systems 5-1 unit step function u0[n] 9-14
for discrete-time systems 9-40 U zero matrix - see matrix
step function - see unit step function zeros 3-1, 3-2
step invariant method - see trans- unit eigenvectors 5-18 zp2tf MATLAB function 11-16
formation methods for mapping analog unit impulse function (d(t)) 1-8
prototype filters to digital filters unit ramp function (u1(t)) 1-9
stop-band filter - see filter unit step function (u0(t)) 1-2
string in MATLAB A-15 upper triangular matrix - see matrix
subplots in MATLAB A-16 using MATLAB for finding the Laplace
summation in the discrete-time Domain transforms of time functions 2-26
see Z transform - properties of using MATLAB for finding the Fourier
symmetric matrix - see matrix transforms of time function 8-31
symmetric rectangular pulse expressed
as sum of unit step functions 1-5 V
symmetric triangular waveform expressed
as sum of unit step functions 1-6 Vandermonde matrix 10-18
symmetry - see Fourier series - symmetry Vector Scope Simulink block 11-78
symmetry property of the Fourier transform
see Fourier transform - properties of W
system function - definition of 8-34
warping 11-52
T window functions
Blackman E-10
Taylor series 5-1 Fourier series method for approximating
text MATLAB command A-13 an FIR amplitude response E-15
tf2ss MATLAB function 5-33 Hamming E-8, E-30
theorems of the DFT 10-10 Hanning E-6, E-26
theorems of the Fourier Transform 8-9 Kaiser E-12, E-33
theorems of the Laplace transform 2-2 other used as MATLAB functions E-14
theorems of the Z Transform 9-3 rectangular E-2
third harmonic - see Fourier triangular E-4, E-22
series - harmonics of Window Visualization Tool in MATLAB E-4
time convolution in DFT
see DFT - common properties of X
time integration property of the Fourier
transform - see Fourier xlabel MATLAB command A-12
transform - properties of
time periodicity property of the Laplace Y
transform - see Laplace
transform - properties of ylabel MATLAB command A-12
time scaling property of the Fourier
transform - see Fourier Z
transform - properties of
time shift in DFT Z transform
see DFT - common properties of computation of with contour
time shift property of the Fourier transform integration 9-17
see Fourier transform - properties of definition of 9-1
time shift property of the Laplace transform Inverse of 9-1, 9-25
see Laplace transform - properties of Z transform - properties of
title(string) MATLAB command A-12 convolution in the discrete
trace of a matrix - see matrix frequency domain 9-9

IN-5