EX 13.1 

%---------------------------------------------------------------------------------

delta = 0.1;  %采样时间为0.1s

s= tf('s');

sysc= 1/(s*(s+1));  %连续系统传递函数

sysd = c2d(sysc, delta); % 得到离散系统传递函数，采样时间为0.1s

[num, den ] = tfdata(sysd, 'v');%从传函中提取分子分母多项式系数；v功能，表示返回分子分母多项式系数；

[z,p,k] = tf2zp(num, den)  % 得到零极点

sysd2 = zpk(z,p,k,delta) %得到零极点形式，0.1s 为采样时间

%计算结果：

                 0.0048374 (z+0.9672)

sysd2 =     --------------------

                    (z-1) (z-0.9048)

%---------------------------------------------------------------------------------

%式（13.3.8）

delta = 2;  %采样时间

s= tf('s');

sysc= (s+4)/(s+1)^3;  %连续系统传递函数

sysd = c2d(sysc, delta); % 得到离散系统传递函数，采样时间为delta

[num, den ] = tfdata(sysd, 'v');%从传函中提取分子分母多项式系数；v功能，表示返回分子分母多项式系数；

[z,p,k] = tf2zp(num, den)  % 得到零极点

sysd2 = zpk(z,p,k,delta) %得到零极点形式，delta为采样时间

%计算结果

              1.564 (z+0.6082) (z+0.02809)

sysd2 = ----------------------------

                      (z-0.1353)^3

%---------------------------------------------------------------------------------

%式（13.3.8）

delta = 0.5;  %采样时间

s= tf('s');

sysc= (s+4)/(s+1)^3;  %连续系统传递函数

sysd = c2d(sysc, delta); % 得到离散系统传递函数，采样时间为delta

[num, den ] = tfdata(sysd, 'v');%从传函中提取分子分母多项式系数；v功能，表示返回分子分母多项式系数；

[z,p,k] = tf2zp(num, den)  % 得到零极点

sysd2 = zpk(z,p,k,delta) %得到零极点形式，delta为采样时间

%计算结果

              0.13337 (z+1.097) (z-0.1286)

sysd2 = ----------------------------

                      (z-0.6065)^3