This tutorial is available as a supplement to the textbook Fundamentals of Signals and Systems Using rd the Web and Matlab, 3 edition by Edward Kamen and Bonnie Heck, published by Prentice Hall. The tutorial covers basic MATLAB commands that are used in introductory signals and systems analysis. It is meant to serve as a quick way to learn MATLAB and a quick reference to the commands that are used in this textbook. For more detailed information, the reader should consult the official MATLAB documentation. An easy way to learn MATLAB is to sit down at a computer and follow along with the examples given in this tutorial and the examples given in the textbook.
The tutorial is designed for students using either the professional version of MATLAB (ver. 7.01) with the Control Systems Toolbox (ver. 6.1), the Signal Processing Toolbox (ver. 6.2.1), and the Symbolic Math Toolbox (ver. 3.1) or using the Student Edition of MATLAB (ver. 7.0).
The topics covered in this tutorial are:
1. MATLAB Basics 2 A. Definition of Variables 2 B. Definition of Matrices 4 C. General Information 6 D. M-files 6 2. Fourier Analysis 9 3. Continuous Time System Analysis 11 A. Transfer Function Representation 11 B. Time Simulations 15 C. Frequency Response Plots 17 D. Analog Filter Design 18 E. Control Design 19 F. State Space Representation 20 4. Discrete-Time System Analysis 22 A. Convolution 22 B. Transfer Function Representation 22 C. Time ...
MATLAB Tutorial This tutorial is available as a supplement to the textbookFundamentals of Signals and Systems Using the Web and Matlab, 3rd by Edward Kamen and Bonnie Heck, published by Prentice Hall. edition The tutorial covers basic MATLAB commands that are used in introductory signals and systems analysis. It is meant to serve as a quick way to learn MATLAB and a quick reference to the commands that are used in this textbook. For more detailed information, the reader should consult the official MATLAB documentation. An easy way to learn MATLAB is to sit down at a computer and follow along with the examples given in this tutorial and the examples given in the textbook. The tutorial is designed for students using either the professional version of MATLAB (ver. 7.01) with the Control Systems Toolbox (ver. 6.1), the Signal Processing Toolbox (ver. 6.2.1), and the Symbolic Math Toolbox (ver. 3.1) or using the Student Edition of MATLAB (ver. 7.0). The topics covered in this tutorial are: 1. MATLAB Basics A. Definition of Variables B. Definition of Matrices C. General Information D. M-files 2. Fourier Analysis 3. Continuous Time System Analysis A. Transfer Function Representation B. Time Simulations C. Frequency Response Plots D. Analog Filter Design E. Control Design F. State Space Representation 4. Discrete-Time System Analysis A. Convolution B. Transfer Function Representation C. Time Simulations D. Frequency Response Plots E. Digital Filter Design F. Digital Control Design G. State Space Representation 5. Plotting 6. Loading and Saving Data
1. MATLAB BasicsMATLAB is started by clicking the mouse on the appropriate icon and is ended by typingexitor by using the menu option. After each MATLAB command, the "return" or "enter" key must be depressed. A. Definition of VariablesVariables are assigned numerical values by typing the expression directly, for example, typing a = 1+2yields: a = 3The answer will not be displayed when a semicolon is put at the end of an expression, for example typea = 1+2;. MATLAB utilizes the following arithmetic operators: + addition - subtraction * multiplication / division ^ power operator ' transpose A variable can be assigned using a formula that utilizes these operators and either numbers or previously defined variables. For example, sinceawas defined previously, the following expression is valid b = 2*a;To determine the value of a previously defined quantity, type the quantity by itself: byields: = b 6If your expression does not fit on one line, use an ellipsis (three or more periods at the end of the line) and continue on the next line. c = 1+2+3+... 5 6+7;+
2
There are several predefined variables which can be used at any time, in the same manner as user-defined variables: i sqrt(-1) j sqrt(-1) pi 3.1416... For example, y= 2*(1+4*j)yields: y= 2.0000 + 8.0000iThere are also a number of predefined functions that can be used when defining a variable. Some common functions that are used in this text are: abs of a number (absolute value for real numbers) magnitude angle of a complex number, in radians angle cos function, assumes argument is in radians cosine sin function, assumes argument is in radians sine exp exponential function For example, withydefined as above, c = abs(y)yields: = c 8.2462 c = angle(y)yields: c = 1.3258Witha=3 as defined previously, c = cos(a)yields: = c 0.9900- c = exp(a)yields: c = 20.0855
3
Note thatexp For example, withcan be used on complex numbers.y = 2+8ias defined above, c = exp(y)yields: = c -1.0751 + 7.3104iwhich can be verified by using Euler's formula: c = e2cos(8) + je2sin(8) B. Definition of Matrices MATLAB is based on matrix and vector algebra; even scalars are treated as 1x1 matrices. Therefore, vector and matrix operations are as simple as common calculator operations. Vectors can be defined in two ways. The first method is used for arbitrary elements: v = [1 3 5 7];creates a 1x4 vector with elements 1, 3, 5 and 7. Note that commas could have been used in place of spaces to separate the elements. Additional elements can be added to the vector: v(5) = 8;yields the vector v = [1 3 5 7 8] defined vectors can be used to define a new. Previously vector. For example, withvdefined above a = [9 10]; b = [v a];creates the vectorb = [1 3 5 7 8 9 10]. The second method is used for creating vectors with equally spaced elements: t = 0:.1:10;creates a 1x101 vector with the elements 0, .1, .2, .3,...,10. Note that the middle number defines the increment. If only two numbers are given, then the increment is set to a default of 1: k = 0:10;creates a 1x11 vector with the elements 0, 1, 2, ..., 10.
4
Matrices are defined by entering the elements row by row: M = [1 2 4; 3 6 8];creates the matrix M=é1 2 4ùëê3 6 8ûú There are a number of special matrices that can be defined: null matrix: M = []; nxm matrix of zeros: = zeros(n,m); M nxm matrix of ones: = ones(n,m); M nxn identity matrix: M = eye(n);A particular element of a matrix can be assigned: M(1,2) = 5;places the number 5 in the first row, second column. In this text, matrices are used only in Chapter 12; however, vectors are used throughout the text. Operations and functions that were defined for scalars in the previous section can also be used on vectors and matrices. For example, a = [1 2 3]; b = [4 5 6]; c a + b= yields: = c 5 7 9Functions are applied element by element. For example, t = 0:10; x = cos(2*t);creates a vectorxwith elements equal to cos(2t) for t = 0, 1, 2, ..., 10.
5
Operations that need to be performed element-by-element can be accomplished by preceding the operation by a ".". For example, to obtain a vector xthat contains the elements of x(t) = tcos(t) at specific points in time, you cannot simply multiply the vectortwith the vectorcos(t). Instead you multiply their elements together: t = 0:10; x t.*cos(t);= C. General Information " MATLAB is case sensitive so "aand "A" are two different names. Comment statements are preceded by a "%". On-line help for MATLAB can be reached by typing help for the full menu or typing helpfollowed by a particular function name or M-file name. For example, help cosgives help on the cosine function. The number of digits displayed is not related to the accuracy. Tochange the format of the display, typeformat short efor scientific notation with 5 decimal places,format long efor scientific notation with 15 significant decimal places and bank formatfor placing two significant digits to the right of the decimal. The commands who and whos give the names of the variables that have been defined in the workspace. The command length(x)returns the length of a vectorxandsize(x)returns the dimension of the matrixx. D. M-filesM-files are macros of MATLAB commands that are stored as ordinary text files with the extension "m", that isfilename M-file can be either .m. An lista function with input and output variables or it can be a of commands. All of the MATLAB examples in this textbook are contained in M-files that are available from the textbook website. The following describes the use of M-files on a PC version of MATLAB. MATLAB requires that the M-file must be stored in a directory that is specified in the MATLAB path list. For example, consider using MATLAB on a PC with a user-defined M-file stored in a directory called "\MATLAB\MFILES". ThentoaccessthatM-file,addthatdirectorytothepathbyclicking“Filethenclicking“SetPathand following directions.
6
The M-files associated with this textbook should be downloaded from www.ece.gatech.edu/users/192/book/M-files.html and copied to a user-specified directory. That directory should be added to the path. The M-files that come with MATLAB are already in the path. To create an M-file, click on “File then click on “New then click on “M-file and an editor screen will appear. Type in your MATLAB commands and then save the file to a directory that is in the path. The M-file is run by typing the name of the M-file from the MATLAB command window (without the .m extension). For example, suppose an M-file named example.m is located in the path. Then typing example from the command prompt runs that M-file. As example of an M-file that defines a function, create a file in your working directory named yplusx.m that contains the following commands: function z = yplusx(y,x) z = y + x;The following commands typed from within MATLAB demonstrate how this M-file is used: x = 2; y = 3; z = yplusx(y,x)MATLAB M-files are most efficient when written in a way that utilizes matrix or vector operations. Loops and if statements are available, but should be used sparingly since they are computationally inefficient. An example of the use of the commandforis for k=1:10, x(k) = cos(k); endThis creates a 1x10 vector x containing the cosine of the positive integers from 1 to 10. This operation is performed more efficiently with the commands k = 1:10; x = cos(k);which utilizes a function of a vector instead of a for loop. An if be used to definestatement can conditional statements. An example is if(a <= 2), b = 1; elseif(a > 4) = b = 2; else b = 3;
7
endThe allowable comparisons between expressions are >=, <=, <, >, ==, and ~=. Several of the M-files written for this textbook employ a user-defined variable which is defined with the command input example, suppose that you want to run an M-file with different values of a. For variableT. The following command line within the M-file defines the value: T = input('Input the value of T: ') Whatever comment is between the quotation marks is displayed to the screen when the M-file is running, and the user must enter an appropriate value.
8
2. Fourier AnalysisCommands covered: dft idft fft ifft contfftThe dftcommand uses a straightforward method to compute the discrete Fourier transform. Define a vectorxand compute the DFT using the command X = dft(x)The first element in X function The to the value of X(0). corresponds dft available from the is MathWorks ftp site and is defined in Figure C.2 of the textbook. The command idft a straightforward method to compute the inverse discrete Fourier uses transform. Define a vectorXand compute the IDFT using the command x = idft(X)The first element of the resulting vector xis x[0]. The functionidftis available at the MathWorks ftp site and is defined in Figure C.3 of the textbook. For a more efficient but less obvious program, the discrete Fourier transform can be computed using the command fftwhich performs a Fast Fourier Transform of a sequence of numbers. To compute the FFT of a sequence x[n] which is stored in the vectorx, use the command X = fft(x)Used in this way, the command fft interchangeable with the command is dft more. For computational efficiency, the length of the vector xshould be equal to an exponent of 2, that is 64, 128, 512, 1024, 2048, etc. The vector xcan be padded with zeros to make it have an appropriate length. MATLAB does this automatically by using the following command where N is defined to be an exponent of 2: X = fft(x,N);The longer the length ofx Due to a wrap around effect, only the, the finer the grid will be for the FFT. first N/2 points of the FFT have any meaning. Theifftcommand computes the inverse Fourier transform:
9
x = ifft(X);The FFT can be used to approximate the Fourier transform of a continuous-time signal as shown in Section 6.6 of the textbook. A continuous-time signal x(t) is sampled with a period of T seconds, then the DFT is computed for the sampled signal. The resulting amplitude must be scaled and the corresponding frequency determined. An M-file that approximates the Fourier Transform of a sampled continuous-time signal is available from the ftp site and is given below: function [X,w] = contfft(x,T); [n,m] = size(x); if n<m, x = x'; end Xk = fft(x); N = length(x); n = 0:N-1; n(1) = eps; X = (1-exp(-j*2*pi*n/N))./(j*2*pi*n/N/T). Xk.'; * w = 2*pi*n/N/T;The input is the sampled continuous-time signal x and the sampling time T. The outputs are the Fourier transform stored in the vectorXand the corresponding frequency vectorw.
10
3. Continuous Time System AnalysisA. Transfer Function Representation Commands covered: tf zpk tf2zp zp2tf feedback parallel seriesTransfer functions are created in MATLAB from the vectors that contain the coefficients of the numerator and the denominator. Given a continuous-time transfer function B(s) H (s) = A(s) where B(s) = bMsM+bM-1sM-1+...+b0and A(s) = sN+aN-1sN-1+...+a0 the coefficients of B(s) and. Store A(s) in the vectors num = [bM bM-1 ... b0] and = [1 a denN-1 a ...0]. In this text, the names of the vectors are generally chosen to benumandden, but any other name could be used. Then, type sys = tf(num,den). Wheresysis the desired name of the system For example, is defined by num = [2 3]; den = [1 4 0 5]; H = tf(num,den);Note that all coefficients must be included in the vector, even zero coefficients. A transfer function may also be defined in terms of its zeros, poles and gain: