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 =
[]
>>
