Control systems.
ECN 323-L
Lab codes & notes # 1
Eng. Safa M. Abusaid
Spring 2024
➢ Matlab codes
Please make sure you have Matlab R2015a or newer version.
▪ Laplace and inverse Laplace
Ex_1:
Inserted code Obtained result
syms t s f =
f=3-4*exp(-1*t)+exp(-4*t) exp(-4*t) - 4*exp(-t) + 3
laplace (f)
ans =
1/(s + 4) - 4/(s + 1) + 3/s
2
Ex_2:
Inserted code Obtained result
n=[1 2]; z =
d=[1 4 3 0];
-2
sys=tf(n,d);
[z,p,k]=tf2zp(n,d) p =
pzmap(sys) 0
-3
-1
k =
3
Ex_3:
Inserted code Obtained result
syms s t ft =
Fs=(s+1)/(s^3+s^2+s+0); 1 - exp(-t/2)*(cos((3^(1/2)*t)/2)
- (3^(1/2)*sin((3^(1/2)*t)/2))/3)
ft=ilaplace(Fs)
Ex_4:
Inserted code Obtained result
n=[1 2 8]; sys =
d=[1 2 4 0];
s^2 + 2 s + 8
sys=tf(n,d) -----------------
s^3 + 2 s^2 + 4 s
[z,p]=tf2zp(n,d)
Continuous-time transfer
pzmap(sys) function.
z =
-1.0000 + 2.6458i
-1.0000 - 2.6458i
p =
0.0000 + 0.0000i
-1.0000 + 1.7321i
-1.0000 - 1.7321i