Scribd
Upload
53 views

1/22/19 11:44 AM MATLAB Command Window 1 of 1

The document contains the work of a student submitting solutions to problems on Laplace and inverse Laplace transforms. It includes: 1) Taking the Laplace transform of two functions. 2) Taking the inverse Laplace transform of two functions. 3) Analyzing transfer functions derived from two sets of polynomials. 4) Finding the residues and poles of a transfer function.

Uploaded by

Nico
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
53 views

1/22/19 11:44 AM MATLAB Command Window 1 of 1

The document contains the work of a student submitting solutions to problems on Laplace and inverse Laplace transforms. It includes: 1) Taking the Laplace transform of two functions. 2) Taking the inverse Laplace transform of two functions. 3) Analyzing transfer functions derived from two sets of polynomials. 4) Finding the residues and poles of a transfer function.

Uploaded by

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

1/22/19 11:44 AM MATLAB Command Window 1 of 1

>> %Midterm Practical Exam in Feedback Lab (January 22, 2019)


>> %Submitted to: Engr. Joey P. Sarmiento
>> %Submitted by: Monsalud, Robert Russell N.
>> %Section: E41
>>
>> %Number 1 letter a
>> syms t
>> laplace (8*(t^2)*cos((3*t)+(pi/4)))

ans =

4*2^(1/2)*(6/(s^2 + 9)^2 - (24*s^2)/(s^2 + 9)^3) - 4*2^(1/2)*((6*s)/(s^2 + 9)^2 -


(8*s^3)/(s^2 + 9)^3)

>> %Number 1 letter b


>> syms t
>> laplace (3*t*exp(-2*t)*sin((4*t)+(pi/3)))

ans =

(6*(2*s + 4))/((s + 2)^2 + 16)^2 - (3*3^(1/2)*(1/((s + 2)^2 + 16) - ((2*s + 4)*(s +


2))/((s + 2)^2 + 16)^2))/2

>>
MATLAB Command Window Page 1

>> %Number 2 letter a


>> syms s
>> ilaplace (((s^2+3*s+10)*(s+5))/((s+3)*(s+4)*(s^2+2*s+100)))

ans =

(20*exp(-3*t))/103 - (7*exp(-4*t))/54 + (5203*exp(-t)*(cos(3*11^(1/2)*t) - (11^(1/2)


*sin(3*11^(1/2)*t))/57233))/5562

>> %Number 2 letter b


>> syms s
>> ilaplace ((s^3+4*s^2+2*s+6)/((s+8)*(s^2+8*s+3)*(s^2+5*s+7)))

ans =

(1199*exp(-4*t)*(cosh(13^(1/2)*t) - (4262*13^(1/2)*sinh(13^(1/2)*t))/15587))/417 -
(65*exp(-(5*t)/2)*(cos((3^(1/2)*t)/2) + (131*3^(1/2)*sin((3^(1/2)*t)/2))/15))/4309 -
(266*exp(-8*t))/93

>>
MATLAB Command Window Page 1

>> %Number 3AA


>> nume = [45 2925 51390 147240 133200 53280];
>> deno = [1 115 4499 70700 553692 5201463 3483390 2749500];
>> tf(nume,deno);
>> g = tf(nume,deno);
>> [z1,p1,k1]=zpkdata(g);
>> zpk(z1,p1,k1)

ans =

45 (s+34.88) (s+26.83) (s+2.122) (s^2 + 1.17s + 0.5964)


-----------------------------------------------------------------
(s+47) (s+39) (s+26.34) (s^2 + 0.6618s + 0.5695) (s^2 + 2s + 100)

Continuous-time zero/pole/gain model.

>> %Number 3AB


>> nume = [45 2925 51390 147240 133200 53280];
>> deno = [1 115 4499 70700 553692 5201463 3483390 2749500];
>> tf(nume,deno)

ans =

45 s^5 + 2925 s^4 + 51390 s^3 + 147240 s^2 + 133200 s + 53280

-----------------------------------------------------------------------

s^7 + 115 s^6 + 4499 s^5 + 70700 s^4 + 553692 s^3 + 5.201e06 s^2

+ 3.483e06 s + 2.75e06

Continuous-time transfer function.

>>
MATLAB Command Window Page 1

>> %Number 3BA


>> nume = [56 3528 41888 51576 41552];
>> deno = [1 225 16778 427711 1093333 1188715 753676 165165];
>> tf(nume,deno);
>> g = tf(nume,deno);
>> [z2,p2,k2]=zpkdata(g);
>> zpk(z2,p2,k2)

ans =

56 (s+47.72) (s+14) (s^2 + 1.276s + 1.111)

---------------------------------------------------------------------------

(s+87.62) (s+80.06) (s+54.59) (s+1.411) (s+0.3766) (s^2 + 0.9391s + 0.8119)

Continuous-time zero/pole/gain model.

>> %Number 3BB


>> nume = [56 3528 41888 51576 41552];
>> deno = [1 225 16778 427711 1093333 1188715 753676 165165];
>> tf(nume,deno)

ans =

56 s^4 + 3528 s^3 + 41888 s^2 + 51576 s + 41552

-----------------------------------------------------------------------

s^7 + 225 s^6 + 16778 s^5 + 427711 s^4 + 1.093e06 s^3 + 1.189e06 s^2

+ 753676 s + 165165

Continuous-time transfer function.

>>
MATLAB Command Window Page 1

>> %Number 3CA


>> nume = [5 565 16710 140400];
>> deno = [1 143 6849 123717 788690 3468760 4804800 0];
>> tf(nume,deno);
>> g = tf(nume,deno);
>> [z3,p3,k3]=zpkdata(g);
>> zpk(z3,p3,k3)

ans =

5 (s+72) (s+26) (s+15)


---------------------------------------------------
s (s+56) (s+55) (s+24.91) (s+2.087) (s^2 + 5s + 30)

Continuous-time zero/pole/gain model.

>> %Number 3CB


>> nume = [5 565 16710 140400];
>> deno = [1 143 6849 123717 788690 3468760 4804800 0];
>> tf(nume,deno)

ans =

5 s^3 + 565 s^2 + 16710 s + 140400

-----------------------------------------------------------------------

s^7 + 143 s^6 + 6849 s^5 + 123717 s^4 + 788690 s^3 + 3.469e06 s^2

+ 4.805e06 s

Continuous-time transfer function.

>>
MATLAB Command Window Page 1

>> %Number 4
>> nume = [10000 750000 350000];
>> deno = [1 113 4022 58200 754275 4324375 2586375 0];
>> [r4,p4,k4]=residue(nume,deno)

r4 =

-0.0023 + 0.0000i
0.0088 + 0.0000i
0.0066 + 0.0537i
0.0066 - 0.0537i
-0.2210 + 0.0000i
0.0660 + 0.0000i
0.1353 + 0.0000i

p4 =

-55.2886 + 0.0000i
-44.3110 + 0.0000i
-2.2730 +13.6568i
-2.2730 -13.6568i
-8.1812 + 0.0000i
-0.6732 + 0.0000i
0.0000 + 0.0000i

k4 =

[]

>>

You might also like