Scribd
Upload
92 views4 pages

Lab Report#4.: Control Systems

This lab report discusses three control systems tasks: 1) Proving that tf2zp and zp2tf are inverses, and laplace and illaplace exhibit the same property using a transfer function example. 2) Finding the inverse Laplace transform and solving a differential equation using the transfer function (s^5 +8s^4+23s^3+35s^2+28s+3)/(s^3 +6s^2+8s). 3) Modeling an RLC circuit in Simulink by defining component values, writing the mathematical equation, and creating a Simulink block diagram.

Uploaded by

ShahzebKhurshid
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
92 views4 pages

Lab Report#4.: Control Systems

This lab report discusses three control systems tasks: 1) Proving that tf2zp and zp2tf are inverses, and laplace and illaplace exhibit the same property using a transfer function example. 2) Finding the inverse Laplace transform and solving a differential equation using the transfer function (s^5 +8s^4+23s^3+35s^2+28s+3)/(s^3 +6s^2+8s). 3) Modeling an RLC circuit in Simulink by defining component values, writing the mathematical equation, and creating a Simulink block diagram.

Uploaded by

ShahzebKhurshid
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

Control systems.

Lab report#4.
(Introduction to Matlab and Simulink, Solution of Differential Equations Using
Matlab

and Simulink and Transfer Function Based Modelling )

Group memebers
Muhammad zohaib raza
Rana Imran anwar
Syed arif kazmi

FA13-bet-048
FA13-bet-050
Fa13-bet-026

Submitted to
Sir Ali Ajwad

Lab Tasks:

Task-1:
Considering G(s),
G=s+0.25/[s(s+1)(s+2)(s+3)]
Prove that tf2zp() and zp2tf() are inverse of each other. Also, prove that laplace()
and illaplace() exhibits the same property.

ZP2TF
>> num=[- 0.25]
num =
- 0.2500
>> den=[-1 -2 -3 0]
den =
-1 -2 -3
0
>> g=zpk([-0.25],[0 -1 -2 -3],[1])
Zero/pole/gain:
(s+0.25)
------------------s (s+1) (s+2) (s+3)
>> [num den]=zp2tf([-0.25],[0 -1 -2 -3],[1])
num =
0
0
0 1.0000 0.2500
den =
1
6 11
6
0

TF2ZP
>> num=[1 0.25]
num =
1.0000 0.2500
>> den=[1 6 11 6 0]
den =
1
6 11
6
0
>> g=tf(num,den)
Transfer function:
s + 0.25
-------------------------s^4 + 6 s^3 + 11 s^2 + 6 s
>> [z p k]=tf2zp(num,den)
z=
-0.2500
p=
0
-3.0000
-2.0000
-1.0000

k=
1

Laplace and Ilaplace


>> syms s
>> g=(s+0.25)/(s^4+6*s^3+11*s^2+6*s)
g=
(s + 1/4)/(s^4 + 6*s^3 + 11*s^2 + 6*s)
>> ilaplace(g)
ans =
3/(8*exp(t)) - 7/(8*exp(2*t)) + 11/(24*exp(3*t)) + 1/24
>> laplace(ans)
ans =
3/(8*(s + 1)) - 7/(8*(s + 2)) + 11/(24*(s + 3)) + 1/(24*s)
Task- 2:
Find inverse laplace transform f(t) & hence find inverse from partial fraction form (use MATLAB to compute the
residues). Where

F( s)=

s5 +8 s 4+ 23 s 3+ 35 s2 +28 s+3
s 3 +6 s 2 +8 s

Solve the differential equation to compute x(t)


X + 2X + 10X = t2

Ilaplace
>> syms s
>> g=(s^5+8*s^4+23*s^3+35*s^2+28*s+3)/(s^3+6*s^2+8*s)
g=
(s^5 + 8*s^4 + 23*s^3 + 35*s^2 + 28*s + 3)/(s^3 + 6*s^2 + 8*s)
>> ilaplace(g)
ans =
1/(4*exp(2*t)) + 3/(8*exp(4*t)) + 3*dirac(t) + 2*dirac(t, 1) + dirac(t, 2) + 3/8

Inverse from Partial Fraction


>> num=[1 8 23 35 28 3]
num =
1
8 23 35 28
3
>> den=[0 0 1 6 8 0]
den =
0
0
1
6
8
0
>> [r,p,k]=residue(num,den)
r=
0.3750
0.2500
0.3750
p=

-4
-2
0
k=
1

Lab Task No 3:

Consider the following RLC circuit (Fig 4.3) and take values of your
choice. Make a mathematical equation of this circuit (taking your own
defined values) and then make its Simulink block diagram.

Simulink model.

Result.

You might also like