Ollllaaaaaaa

A=[0 1 0 0;0 -.1818 2.6754 0;0 0 0 1;0 -.454 31.213 0]

A =

0 1.0000 0 0
0 -0.1818 2.6754 0
0 0 0 1.00000 -0.4540 31.2130 0

>> B=[0;1.8181;0;4.5454]

B =

0
1.8181
0
4.5454

>> C=[0 0 1 0]

C =

0 0 1 0

>>D=[0]

D =

0

>> [num,dem]=ss2tf(A,B,C,D)

num =

0 -0.0000 4.5454 0.0009 0

dem =

1.0000 0.1818 -31.2130 -4.4599 0

>>H=tf([0 -0.0000 4.5454 0.0009 0],[1.0000 0.1818 -31.2130 -4.4599 0])

Transfer function:
4.545 s^2 + 0.0009 s
-------------------------------------
s^4 + 0.1818 s^3 - 31.21 s^2 - 4.46 s>> s=tf('s')

Transfer function:
s

>> H1=H/s

Transfer function:
4.545 s^2 + 0.0009 s
---------------------------------------
s^5 + 0.1818 s^4 - 31.21 s^3 - 4.46s^2

>> H2=c2d(H1,.01
??? H2=c2d(H1,.01
|
Error: Expression or statement is incorrect--possibly unbalanced (, {, or [.

>> H2=c2d(H1,0.01)

Transfer function:7.573e-007 z^3 + 2.271e-006 z^2 - 2.272e-006 z - 7.567e-007
-----------------------------------------------------------
z^4 - 4.001 z^3 + 6.001 z^2 - 3.998 z + 0.9982

Sampling time:0.01
>> Z=tf('z')

Transfer function:
z

Sampling time: unspecified
>> GZ=((Z-1/Z)*H2)

Transfer function:
7.573e-007 z^5 + 2.271e-006 z^4 - 3.029e-006 z^3 - 3.028e-006 z^2 +2.272e-006 z + 7.567e-007
---------------------------------------------------------------------------------------------
z^5 - 4.001 z^4 + 6.001 z^3 - 3.998 z^2 + 0.9982 zSampling time: 0.01
>> a=sym('a')

a =

a

>> b=sym('b')

b =

b

>> Kd=sym('Kd')

Kd =

Kd

>> conv([1 0.23],conv([1 0.099],conv([1 0.0215],conv([1...