According to OP's comment. It appears that the term 1/s appeared in exponent by mistake. It is actually being multiplied by the denominator. According following method can be used
You can use c2d function to convert it into discrete. Note that the correct way to define the transfer function is like this
s = tf('s')
sys = 1/(s*(s+1)^3)
sys_dis = c2d(sys, 0.1)
Result
>> sys
sys =
1
-----------------------
s^4 + 3 s^3 + 3 s^2 + s
Continuous-time transfer function.
>> sys_dis
sys_dis =
3.925e-06 z^3 + 4.067e-05 z^2 + 3.83e-05 z + 3.278e-06
------------------------------------------------------
z^4 - 3.715 z^3 + 5.171 z^2 - 3.197 z + 0.7408
Sample time: 0.1 seconds
Discrete-time transfer function.
To plot the values
bode(sys_dis)
pzplot(sys_dis)
step(sys_dis)
impulse(sys_dis)
Best Answer