信号与系统陈后金matlab,陈后金信号与系统matlab实验7.pdf

陈后金信号与系统matlab实验7

第八章 离散时间信号与系统的 域分析

z

作者卢未来 日期 年 月 日

: :2012 5 10

8.5利用MATLAB进行离散系统的z域分析

8.5.1部分分式展开式的MATLAB实现

8-27(301 ) MATLAB

例 页 试用 计算

18

F(z)= −1 −2 −3

18+3z −4z − z

的部分分式展开式。

>>%program8_1

num=[18];

den=[183-4-1];

[r,p,k]=residuez(num,den)

r=

0.3600

0.2400

0.4000

p=

0.5000

-0.3333

-0.3333

k=

[]

8.5.2利用MATLAB计算H(z)的零极点与系统特性

例8-28(302页)已知一离散因果LTI系统的系统函数为

−1 −2 −3

z +2z + z

H(z)= −1 −2 −3

1−0.5z −0.005z +0.3z

求该系统的零极点。

>%program8_2

%zerosandpolesofthetransferfunction

b=[121];

a=[1-0.5-0.0050.3];

[r,p,k]=tf2zp(b,a)

r=

-1

-1

p=

0.5198+0.5346i

0.5198-0.5346i

-0.5396

k=

1

例8-29(303页)已知一离散因果LTI系统的系统函数为

2

z +2z+1

H(z)= 3 2 −1

z −0.5z −0.005z +0.3

h[k]

试画出系统的零极点分布图，求系统的单位脉冲响应 和频率响应

jΩ

H(e )，并判断系统是否稳定。

>>%program8_3zerosandpolesofthetransferfunction

b=[121];

a=[1-0.5-0.0050.3];

figure(1);zplane(b,a);

num=[0121];

den=[1-0.5-0.0050.3];

h=impz(num,den);

figure(2);stem(h)

xlabel('k')

title('ImpulseRespone')

[H,w]=freqz(num,den);

figure(3);plot(w/pi,abs(H))

xlabel('Frequency\omega')

title('MagnitudeRespone')

>>