2018-EE-183
LAB 4: MODELING AND SIMULATION OF SECOND
ORDER SISO DYNAMIC SYSTEM
TABLE OF CONTENTS
Define the Transfer Function:
Poles And Zeroes:
Finding all the Parameters:
Conversion of Transfer Function to State Space
Conversion of State Space to Transfer Function
System in Time domain
Step System Time Domain:
Finding Poles and Zeros from "tf2zp" Command
Bode Plot
PZMAP of the System
Step Response of System
DEFINE THE TRANSFER FUNCTION:
m = 2
m = 2
b =0.5
b = 0.5000
k = 0.9
k = 0.9000
s = tf('s')
s =
Continuous-time transfer function.
H = s/(m*s^2+b*s+k)
H =
-------------------
2 s^2 + 0.5 s + 0.9
Continuous-time transfer function.
POLES AND ZEROES:
pole(H)
ans = 2×1 complex
-0.1250 + 0.6591i
-0.1250 - 0.6591i
zero(H)
ans = 0
FINDING ALL THE PARAMETERS:
[Z,P,K] = zpkdata(H)
Z = 1×1 cell array
{[0]}
P = 1×1 cell array
{2×1 double}
K = 0.5000
pole(H)
ans = 2×1 complex
-0.1250 + 0.6591i
-0.1250 - 0.6591i
zero(H)
ans = 0
CONVERSION OF TRANSFER FUNCTION TO STATE SPACE
[A,B,C,D] = tf2ss([0 1 0],[2 0.5 0.9])
A = 2×2
-0.2500 -0.4500
1.0000 0
B = 2×1
C = 1×2
0.5000 0
D = 0
CONVERSION OF STATE SPACE TO TRANSFER FUNCTION
A = [0,1;-0.45,-0.25]
A = 2×2
0 1.0000
-0.4500 -0.2500
B = [0;0.5]
B = 2×1
0.5000
C = [0 1]
C = 1×2
0 1
D = 0
D = 0
[M,P] = ss2tf(A,B,C,D)
M = 1×3
0 0.5000 0
P = 1×3
1.0000 0.2500 0.4500
SYSTEM IN TIME DOMAIN
syms w
Y = ilaplace(w/(2*w^2+0.5*w+0.9))
Y =
vpa(Y,4)
ans =
STEP SYSTEM TIME DOMAIN:
syms w
Y = ilaplace(w/(2*w^3+0.5*w^2+0.9*w))
Y =
vpa(Y,4)
ans =
FINDING POLES AND ZEROS FROM "TF2ZP" COMMAND
[z,p] = tf2zp([0 1 0],[2 0.5 0.9])
z = 0
p = 2×1 complex
-0.1250 + 0.6591i
-0.1250 - 0.6591i
[z,p,k] = tf2zpk([0 1 0],[2 0.5 0.9])
z = 0
p = 2×1 complex
-0.1250 + 0.6591i
-0.1250 - 0.6591i
k = 0.5000
BODE PLOT
bode(H)
PZMAP OF THE SYSTEM
pzmap(H)
STEP RESPONSE OF SYSTEM
stepplot(H)