Scribd
Upload
38 views

Lab Report

The document describes 7 tasks related to Laplace transforms: 1) Finding the Laplace transform of t*exp(t) in Matlab. 2) Finding the Laplace transform of 2t^3*u(t-4) in Matlab. 3) Verifying the shifting property of Laplace transforms in Matlab. 4) Verifying the scaling property for cos(2t) in Matlab. 5) Verifying the linearity property for cos(2t) and sin(2t) in Matlab. 6) Finding the inverse Laplace transform of (3s^3 + 4s^2 + s + 11)/(s^2 + 41s + 12) in Matlab. 7) Calculating the

Uploaded by

HarisArshad
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
38 views

Lab Report

The document describes 7 tasks related to Laplace transforms: 1) Finding the Laplace transform of t*exp(t) in Matlab. 2) Finding the Laplace transform of 2t^3*u(t-4) in Matlab. 3) Verifying the shifting property of Laplace transforms in Matlab. 4) Verifying the scaling property for cos(2t) in Matlab. 5) Verifying the linearity property for cos(2t) and sin(2t) in Matlab. 6) Finding the inverse Laplace transform of (3s^3 + 4s^2 + s + 11)/(s^2 + 41s + 12) in Matlab. 7) Calculating the

Uploaded by

HarisArshad
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

LAB REPORT 10

LAB TASKS

1. Find the Laplace transform of tet in Matlab.


clc
syms f t;
f=t*(exp(t));
laplace(f)

ans =

1/(s - 1)^2

0, t4
2. Find the Laplace transform of the signal in Matlab: f (t )   3 or f (t)  2t3u(t  4).
2t , t  4

clc
syms f t;
f=2*t^3*(heaviside(t-4));
laplace(f)
ans =

(128*exp(-4*s))/s + (96*exp(-4*s))/s^2 + (48*exp(-4*s))/s^3 + (12*exp(-


4*s))/s^4

3. Verify the shifting property on above signal f(t) in Matlab as:


 
L [ f (t  a)u (t  a)]   f (t  a)u (t  a)e  st dt   f (t  a)e  st dt
0 a

Let   t  a, then
 
L [ f (t  a)u (t  a)]   f ( )e  s (  a ) d  e  sa  f ( )e  s d  e  sa F ( s)
0 0
 
F ( s  a)   f (t )e ( s a )t dt   [e at f (t )]e  st dt  L [e at f (t )]
0 0

clc
syms f t ;
f1=2*(t-3)^3*(heaviside(t-7));

laplace(f1)

ans =

2*exp(-3*s)*((exp(-4*s)*(64*s^3 + 48*s^2 + 24*s + 6) - 6)/s^4 + 6/s^4)


4. Verify the scaling property on the following signal f(t)=cos(at), a=2, (a=scaling factor) in
Matlab as:

L [ f (at )]   f (at )e  st dt
0

Let   at , then


s
 s 1    1 s
L [ f (at )]   f ( )e a
d   f ( )e a d  F ( )
0 a a 0 a a
clc
syms f t ;
f1 = cos (2*t);
laplace(f1)

ans =

s/(s^2 + 4)

5. Verify the linearity property for two appropriate signals in Matlab as:
  
ℒ [af (t) + bg(t)]   [af (t )  bg (t )]e  st dt  a  f (t )e  st dt  b  g (t )e  st dt  aF ( s )  bG ( s )
0 0 0

clc
syms f t c d e;
f1 = cos (2*t);
f2 = sin (2*t);
f = (a*f1) + (b*f2);
laplace(f)
c = laplace (a*f1);
d = laplace (b*f2);
e = c+d

ans =

[ (5*s)/(s^2 + 4) + 18/(s^2 + 4), (2*s)/(s^2 + 4) + 8/(s^2 + 4), (45*s)/(s^2 + 4) + 10/(s^2 + 4),


s/(s^2 + 4) + 6/(s^2 + 4), s/(s^2 + 4) + 4/(s^2 + 4)]

e=

[ (5*s)/(s^2 + 4) + 18/(s^2 + 4), (2*s)/(s^2 + 4) + 8/(s^2 + 4), (45*s)/(s^2 + 4) + 10/(s^2 + 4),


s/(s^2 + 4) + 6/(s^2 + 4), s/(s^2 + 4) + 4/(s^2 + 4)]

6. Find Inverse Laplace Transform for the following Laplace transform:

3𝑠 3 + 4𝑠 2 + 𝑠 + 11
𝐹(𝑠) =
𝑠 2 + 41𝑠 + 12
clc
syms s;
f = ((3*s^3) + 4*s^2 + s +11)/(s^2+ 41*s + 12)
ilaplace (f)

f=

(3*s^3 + 4*s^2 + s + 11)/(s^2 + 41*s + 12)

ans =

3*dirac(1, t) - 119*dirac(t) + 4844*exp(-(41*t)/2)*(cosh((1633^(1/2)*t)/2) -


(97863*1633^(1/2)*sinh((1633^(1/2)*t)/2))/3955126)

7. Calculate the partial fraction form for the following Laplace transform:

3𝑠 3 + 4𝑠 2 + 𝑠 + 11
𝐹(𝑠) =
𝑠 2 + 41𝑠 + 12
clc
b = [3 4 1 11];
a = [0 1 41 12];
[r p k] = residue (b,a)
r1 =17.25;
r2= 0.0003;
p1 =-40.7052;
p2= -0.2948;
k1=3;
k2=-119;
k=k1+k2;
f = r1/(s-p1) + r2/(s-p2) + k

r=

1.0e+03 *

4.8437
0.0003

p=

-40.7052
-0.2948

k=
3 -119

f=

1.0e+02 *

Columns 1 through 3

-1.1558 - 0.0001i -1.1600 + 0.0000i -1.1600 + 0.0000i

Column 4

-1.1600 + 0.0000i

You might also like