MATLAB BASICS ECEN 605 Linear Control Systems Instructor

MATLAB BASICS ECEN 605 Linear Control Systems Instructor: S. P. Bhattacharyya

STARTING MATLAB l You can start MATLAB by double-clicking on the MATLAB icon or invoking the application from the Start menu of Windows. The main MATLAB window will then pop-up.

BASIC OPERATION MATLAB is based on matrix and vector algebra; even scalars are treated as 1 x 1 matrices. l e. g. v = [1 3 5 7]; t = 0: . 1: 10; k = 0: 10; M = [1 2 4; 3 6 8]; l

APPROACH l l l TRY EXPLORING ALL THE COMMANDS INITIALLY TO GET FAMILIARIZED. TAKE UP EXAMPLES AND SOLVE THEM TRY THE ASSIGNMENT QUESTIONS

SPECIAL MATRICES There a number of special matrices that can be defined: null matrix: M = [ ]; nxm matrix of zeros: M = zeros(n, m); nxm matrix of ones: M = ones(n, m); nxn identity matrix: M = eye(n); l

GENERAL INFORMATION l l l Matlab is case sensitive so "a" and "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 help followed by a particular function name. 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 vector x and size(x) returns the dimension of the matrix x.

CREATING FUNCTIONS l l l 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;

USING THE FUNCTIONS l l The following commands typed from within MATLAB demonstrate how this M-file is used: x = 2; y = 3; z = yplusx(y, x)

PLOTTING l l l l plot xlabel ylabel title grid axis stem subplot

EXAMPLE OF A PLOT l l l Z = [0 : pi/50 : 10*pi]; X = exp(-. 2. *Z). *cos(Z); Y = exp(-. 2. *Z). *sin(Z); plot 3(X, Y, Z); grid on; xlabel('x-axis'); ylabel('y-axis'); zlabel('z-axis');

Continuous Time System Analysis l l l Transfer Function Representation Time Simulations Frequency Response Plots Control Design State Space Representation

Transfer Function Representation l l l Tf 2 zp Zp 2 tf Cloop Feedback Parallel series

Cont… Transfer functions are defined in MATLAB by storing the coefficients of the numerator and the denominator in vectors. Given a continuous-time transfer function B(s) H(s) = ---- A(s) l

Cont… Where B(s) = b. Ms. M+b. M-1 s. M-1+…+b 0 and A(s) = a. Ns. N+a. N-1 s. N-1+…+a 0 Store the coefficients of B(s) and A(s) in the vectors num = [b. M-1 … b 0] den = [a. N-1 … a 0] l

Example 5 s+6 H(s) = -------- s 3+10 s 2+5 l num = [5 6]; l den = [1 10 0 5]; l all coefficients must be included in the vector, even zero coefficients

Cont… l To find the zeros, poles and gain of a transfer function from the vectors num and den which contain the coefficients of the numerator and denominator polynomials: [z, p, k] = tf 2 zp(num, den)

Example l num = [5 6]; den = [1 10 0 5]; l [z, p, k] = tf 2 zp(num, den) l z = -1. 2000 p = -10. 0495 0. 0248+0. 7049 i 0. 0248 -0. 7049 i k = 5 l l l

Example Cont… (s-z 1)(s-z 2). . . (s-zn) H(s) = K ------------- (s-p 1)(s-p 2). . . (s-pn) 5*(s+1. 2) ----------------------------- (s+10. 0495)(s-{0. 0248+0. 7049 i})(s-{0. 0248 -0. 7049 i})

Cont… To find the numerator and denominator polynomials from z, p, and k: l [num, den] = zp 2 tf(z, p, k) l To reduce the general feedback system to a single transfer function: T(s) = G(s)/(1+G(s)H(s)) [num. T, den. T]=feedback(num. G, den. G, num. H, den. H);

Cont… For a unity feedback system: [num. T, den. T] = cloop(num. G, den. G, -1); l To reduce the series system to a single transfer function, T(s) = G(s)H(s) [num. T, den. T] = series(num. G, den. G, num. H, den. H); l To reduce the parallel system to a single transfer function, T(s) = G(s) + H(s) [num. T, den. T] = parallel(num. G, den. G, num. H, den. H); l

Example l [num, den] = zp 2 tf(z, p, k) l num = 0 0 5 6 den = 1. 0000 10. 0000 5. 0000 l

Time Simulations l l residue Step Impulse lsim

Cont… l l [R, P, K] = residue (B, A) finds the residues, poles and direct term of a partial fraction expansion of the ratio of two polynomials B(s)/A(s) The residues are stored in r, the corresponding poles are stored in p, and the gain is stored in k.

Example If the ratio of two polynomials is expressed as b(s) 5 s 3+3 s 2 -2 s+7 ------- = --------- a(s) -4 s 3+8 s+3 l l l b = [ 5 3 -2 7] a = [-4 0 8 3]

Example Cont… l l l r = -1. 4167 -0. 6653 1. 3320 p = 1. 5737 -1. 1644 -0. 4093 k = -1. 2500 B(s) R(1) R(2) R(n) ------ = ---- + ----- +. . . + ----- + K(s) A(s) s - P(1) s - P(2) s - P(n)

Find the response of a system to a particular input First store the numerator and denominator of the transfer function in num and den, respectively. l To plot the step response: step(num, den) l To plot the impulse response: impulse(num, den) l

Cont… For the response to an arbitrary input, use the command lsim (linear simulation) l Create a vector t which contains the time values in seconds t = a: b: c; l Define the input x as a function of time, for example, a ramp is defined as x = t lsim(num, den, x, t); l

Frequency Response Plots l l l Freqs Bode Logspace Log 10 Semilogx

Cont… l l To compute the frequency response H of a transfer function, store the numerator and denominator of the transfer function in the vectors num and den. Define a vector w that contains the frequencies for which H) is to be computed, for example w = a: b: c where a is the lowest frequency, c is the highest frequency and b is the increment in frequency. H = freqs(num, den, w)

Cont… To draw a Bode plot of a transfer function which has been stored in the vectors num and den: bode(num, den) l

Cont… To customize the plot, first define the vector w which contains the frequencies at which the Bode plot will be calculated. l Since w should be defined on a log scale, the command logspace is used. l For example, to make a Bode plot ranging in frequencies from 0. 1 to 100, define w by w = logspace(-1, 2); l The magnitude and phase information for the Bode plot can then be found by: [mag, phase] = bode(num, den, w); l

Cont… To plot the magnitude in decibels, convert mag using the following command: magdb = 20*log 10(mag); l To plot the results on a semilog scale where the yaxis is linear and the x-axis is logarithmic: semilogx(w, magdb) l For the log-magnitude plot : semilogx(w, phase) l

Control Design l Rlocus Consider a feedback loop where G(s)H(s) = KP(s) and K is a gain and P(s) contains the poles and zeros of the controller and of the plant. l Suppose that the numerator and denominator coefficients of P(s) are stored in the vectors num and den. rlocus(num, den) l

State Space Representation l l l Ss Step Lsim Ss 2 tf Tf 2 ss ss 2 ss