Scribd
Upload
88 views4 pages

Matlab Laplace Transform Techniques

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
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
88 views4 pages

Matlab Laplace Transform Techniques

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
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

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