Scribd
Upload
67 views

Lab 06 PDF

This document describes an experiment on linear time-invariant systems and their representation using MATLAB. The objectives are to learn about linear systems modeling and using MATLAB commands to represent systems as transfer functions or pole-zero-gain models. Methods covered include plotting pole-zero maps, simulating responses to impulse, step, and arbitrary inputs using commands like impulse, step, and lsim. Exercises are included to have students obtain and plot impulse/step responses and pole-zero maps for sample systems.

Uploaded by

Abdul Rehman Afzal
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
67 views

Lab 06 PDF

This document describes an experiment on linear time-invariant systems and their representation using MATLAB. The objectives are to learn about linear systems modeling and using MATLAB commands to represent systems as transfer functions or pole-zero-gain models. Methods covered include plotting pole-zero maps, simulating responses to impulse, step, and arbitrary inputs using commands like impulse, step, and lsim. Exercises are included to have students obtain and plot impulse/step responses and pole-zero maps for sample systems.

Uploaded by

Abdul Rehman Afzal
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 7

International Islamic University, Islamabad

Faculty of Engineering and Technology


Iqra College of Technology

Control Technology Lab

LAB # 06: Linear Time Invariant Systems and


Representation using MATLAB

Name of Student: …………………………………………

Registration No.: ……………………………………………

Section: …………………………………………………………

Date of Experiment: ………………………………………

Marks obtained: ………………………………………

Remarks: …………………………………………………

Instructor’s Signature: ……………………………...


Lab 06: Linear Time Invariant Systems and Representation using MATLAB

Lab 06: Linear Time Invariant Systems and Representation


using MATLAB

Objectives: This experiment has following two objectives:


1. Continued with the learning of Mathematical Modeling from previous experiment, we
now start focusing the linear systems. We will learn commands in MATLAB that
would be used to represent such systems in terms of transfer function or pole-zero-
gain representations.
2. We will also learn how to make preliminary analysis of such systems using plots of
poles and zeros locations as well as time response due to impulse, step and arbitrary
inputs.

List of Equipment/Software
Following equipment/software is required:
 MATLAB
 Personal Computer

Mass-Spring System Model


The spring force is assumed to be either linear or can be approximated by a linear
function 𝐹𝑠(𝑥) = 𝐾𝑥, B is the friction coefficient, 𝑥(𝑡) is the displacement and 𝐹𝑎(𝑡) is the
applied force:

The differential equation for the above Mass-Spring system can be derived as follows
𝑑2 𝑥(𝑡) 𝑑𝑥(𝑡)
𝑀 +𝐵 + 𝐾𝑥(𝑡) = 𝐹𝑎 (𝑡)
𝑑𝑡 2 𝑑𝑡

Transfer Function:
Applying the Laplace transformation while assuming the initial conditions are zeros, we get
(𝑀𝑠 2 + 𝐵𝑠 + 𝐾) ∗ 𝑋(𝑠) = 𝐹𝑎 (𝑠)
Then the transfer function representation of the system is given by
𝑂𝑢𝑡𝑝𝑢𝑡 𝐹𝑠 (𝑠) 1
𝑇𝐹 = = = 2
𝐼𝑛𝑝𝑢𝑡 𝑋(𝑠) (𝑀𝑠 + 𝐵𝑠 + 𝐾)

Linear Time-Invariant Systems in MATLAB:


Control System Toolbox in MATLAB offers extensive tools to manipulate and analyze linear
time-invariant (LTI) models. It supports both continuous-time and discrete-time systems.
Systems can be single-input/single-output (SISO) or multiple-input/multiple-output (MIMO).
You can specify LTI models as:

Control Technology Lab (ET-373 L) Page 2


Lab 06: Linear Time Invariant Systems and Representation using MATLAB
Transfer functions (TF), for example,
𝑠+2
𝑃(𝑠) = 𝑠2 +𝑠+10
Note: All LTI models are represented as a ratio of polynomial functions

Examples of Creating LTI Models


Building LTI models with Control System Toolbox is straightforward. The following sections
show simple examples. Note that all LTI models, i.e. TF, ZPK and SS are also MATLAB
objects.
>>num = [1 0];
>>den = [1 2 1];
>>sys = tf(num,den)
Transfer function:
s
-------------
s^2 + 2 s + 1

A useful trick is to create the Laplace variable, s. That way, you can specify polynomials
using s as the polynomial variable.

>>s=tf('s');
>>sys= s/(s^2 + 2*s + 1)

Transfer function:
s
-------------
s^2 + 2 s + 1
This is identical to the previous transfer function.

Example of Creating Zero-Pole-Gain Models


To create zero-pole-gain (ZPK) models, you must specify each of the three components in
vector format. For example,

>>sys = zpk([0],[-1 -1],[1])


Zero/pole/gain:
s
--------
(s+1)^2
produces the same transfer function built in the TF example, but the representation is now
ZPK.

This example shows a more complicated ZPK model.


>>sys=zpk([1 0], [-1 -3 -.28],[.776])
Zero/pole/gain:

0.776 s (s-1)
--------------------
(s+1) (s+3) (s+0.28)

Control Technology Lab (ET-373 L) Page 3


Lab 06: Linear Time Invariant Systems and Representation using MATLAB

Plotting poles and zeros of a system:


Pzmap
Compute pole-zero map of LTI models

pzmap(sys)
pzmap(sys1,sys2,...,sysN)
[p,z] = pzmap(sys)

Description:
pzmap(sys) plots the pole-
zero map of the continuous- or
discrete-time LTI model sys. For
SISO systems, pzmap plots the
transfer function poles and zeros. The
poles are plotted as x's and the zeros
are plotted as o's.
pzmap(sys1,sys2,...,sysN)
plots the pole-zero map of several LTI
models on a single figure. The LTI
models can have different numbers of
inputs and outputs. When invoked
with left-hand arguments, [p, z] =
pzmap(sys) returns the system poles
and zeros in the column vectors p and
z. No plot is drawn on the screen. You
can use the functions sgrid or zgrid to
plot lines of constant damping ratio
and natural frequency in the s-plane or
z-plane.

Example:

Plot the poles and zeros of the continuous-time system.

2𝑠 2 + 5𝑠 + 1
𝐻(𝑠) = 2
𝑠 + 2𝑠 + 3

>>H = tf([2 5 1],[1 2 3]); sgrid


>>pzmap(H)

Control Technology Lab (ET-373 L) Page 4


Lab 06: Linear Time Invariant Systems and Representation using MATLAB
Simulation of Linear systems to different inputs
impulse, step and lsim

You can simulate the LTI systems to inputs


like impulse, step and other standard inputs
and see the plot of the response in the
figure window. MATLAB command
‘impulse’ calculates the unit impulse
response of the system, ‘step’ calculates the
unit step response of the system and ‘lsim’
simulates the (time) response of continuous
or discrete linear systems to arbitrary
inputs. When invoked without left-hand
arguments, all three commands plots the
response on the screen.

For example:

To obtain an impulse response


>> H = tf([2 5 1],[1 2 3]);
>>impulse(H)

To obtain a step response type


>>step(H)

Time-interval specification:
To contain the response of the system you
can also specify the time interval to
simulate the system too.
For example,
>> t = 0:0.01:10
>> impulse(H,t)

Or

>> step(H,t)

Simulation to Arbitrary Inputs:


To simulates the (time) response of
continuous or discrete linear systems to
arbitrary inputs use ‘lsim’. When invoked
without left-hand arguments, ‘lsim’ plots
the response on the screen.
lsim(sys,u,t) produces a plot of the time response of the LTI model sys to the input time
history ‘t’,’u’. The vector ‘t’ specifies the time samples for the simulation and consists of
regularly spaced time samples.

Control Technology Lab (ET-373 L) Page 5


Lab 06: Linear Time Invariant Systems and Representation using MATLAB
T = 0:dt:Tfinal

The matrix u must have as many rows as


time samples (length(t)) and as many
columns as system inputs. Each row
𝑢(𝑖), specifies the input value(s) at the
time sample 𝑡(𝑖).

Simulate and plot the response of the


system to a square wave with period of
four seconds.
2𝑠 2 + 5𝑠 + 1
𝑠 2 + 2𝑠 + 3

First generate the square wave with gensig. Sample every 0.1 second during 10 seconds:

>>[u,t] = gensig(‘square’,4,10,0.1);

Then simulate with lsim.

>> H = tf([2 5 1],[1 2 3])

Transfer function:
2 s^2 + 5 s + 1
-----------------
s^2 + 2 s + 3

>> lsim(H,u,t)

Simulink Representation of LTI System:

Control Technology Lab (ET-373 L) Page 6


Lab 06: Linear Time Invariant Systems and Representation using MATLAB
Exercise 1:

Consider the transfer function


6𝑠 2 + 1
𝐺(𝑠) =
𝑠 3 + 3𝑠 2 + 3𝑠 + 7
Using MATLAB plot the pole zero map of the above system

Exercise 2:

a. Obtain the unit impulse response for the following system


𝐶(𝑠) 1
= 2
𝑅(𝑠) 𝑠 + 0.2𝑠 + 1
b. Obtain the unit step response for the following system
𝐶(𝑠) 𝑠
= 2
𝑅(𝑠) 𝑠 + 0.2𝑠 + 1
c. Explain why the results in a. and b. are same?

Exercise 3:

A system has a transfer function


15
𝐶(𝑠) ( ⁄𝑧)(𝑠 + 𝑧)
= 2
𝑅(𝑠) 𝑠 + 3𝑠 + 15

Plot the response of the system when R(s) is a unit impulse and unit step for the
parameter z=3, 6 and 12.

Control Technology Lab (ET-373 L) Page 7

You might also like