' (A) ' S Z N T
' (A) ' S Z N T
syms s z n T
Gs=(s+4)/((s+2)*(s+5));
pretty(Gs)
gt=ilaplace(Gs);
gnT=compose(gt,n*T);
Gz=ztrans(gnT);
Gz=simplify(Gz);
Gz=subs(Gz,T,0.5);
Gz=vpa(simplify(Gz),6);
Gz=vpa(factor(Gz),6);
'G(z) evaluated for T=0.5'
pretty(Gz)
'(b)'
Gs=(s+1)*(s+2)/(s*(s+3)*(s+4));
'G(s)'
pretty(Gs)
gt=ilaplace(Gs);
gnT=compose(gt,n*T);
Gz=ztrans(gnT);
Gz=simplify(Gz);
Gz=subs(Gz,T,0.5);
Gz=vpa(simplify(Gz),6);
Gz=vpa(factor(Gz),6);
'G(z) evaluated for T=0.5'
pretty(Gz)
'(c)'
Gs=20/((s+3)*(s^2+6*s+25));
pretty(Gs)
gt=ilaplace(Gs);
gnT=compose(gt,n*T);
Gz=ztrans(gnT);
Gz=simplify(Gz);
Gz=subs(Gz,T,0.5);
Gz=vpa(simplify(Gz),6);
Gz=vpa(factor(Gz),6);
'G(z) evaluated for T=0.5'
pretty(Gz)
'(d)'
Gs=15/(s*(s+1)*(s^2+10*s+81));
'G(s)'
pretty(Gs)
gt=ilaplace(Gs);
gnT=compose(gt,n*T);
Gz=ztrans(gnT);
Gz=simplify(Gz);
Gz=subs(Gz,T,0.5);
Gz=vpa(simplify(Gz),6);
Gz=vpa(factor(Gz),6);
'G(z) evaluated for T=0.5'
pretty(Gz)
ans =
(a)
s+4
--------------(s + 2) (s + 5)
ans =
G(z) evaluated for T=0.5
z (z - 0.177350)
1.00001 -----------------------------(z - 0.0820851) (z - 0.367881)
ans =
(b)
ans =
G(s)
(s + 1) (s + 2)
----------------s (s + 3) (s + 4)
ans =
G(z) evaluated for T=0.5
z (z - 0.296742) (z - 0.840812)
1.00000 -------------------------------------(z - 0.135335) (z - 0.223130) (z - 1.)
ans =
(c)
20
----------------------2
(s + 3) (s + 6 s + 25)
ans =
G(z) evaluated for T=0.5
(z + 0.223130) z
0.394980 -------------------------------------------2
(0.0497871 + 0.185700 z + z ) (z - 0.223130)
ans =
(d)
ans =
G(s)
15
-------------------------2
s (s + 1) (s + 10 s + 81)
ans =
G(z) evaluated for T=0.5
(z + 0.289173) (z + 0.116365) z
0.0578297 -----------------------------------------------------2
(z - 0.606531) (z - 1.) (0.00673795 + 0.135489 z + z )
>>