杭电数字信号处理实验2.pdf

信号、系统与信号处理实验Ⅱ

实验报告

姓名：王健

学号：

班级：

上课时间：周五-六七八

实验名称：离散系统频率响应和零极点分布

一、实验目的

通过使用MATLAB函数仿真简单的离散时间系统，研究其时频域特性，以加深对离

散时间系统的冲击响应，频率响应分析和零极点分布概念的理解。

二、实验原理与要求

对于一个线性时不变离散时间系统，其输入输出关系可用以下常系数线性差

分方程描述：

故此系统的系统函数可用下式表示：

其频率响应为：

将系统函数的分子，分母分别作因式分解，可得到LTI系统的零极点增益表

达式为:

通过系统的零极点增益表达式，可以判断一个LTI离散时间系统的稳定性。对

一个因果的离散时间系统，若所有的极点都位于单位圆内，则系统是稳定的，同

理，由零极点分布图可以大致估计出系统的频率响应：

(1)单位圆附近的零点对幅度响应的谷点的位置与深度有明显影响，当零点位于

单位圆上时，谷点为零，零点可在单位圆外。

(2)单位圆附近的极点对幅度响应的峰点位置和高度有明显影响。

要求

一个LTI离散时间系统的输入输出差分方程为：

y(n)-1.6y(n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)

（1）编程求此系统的单位冲击响应，并画出波形

（2）若输入:编程求

此系统的输出，并分别画出波形

（3）编程得到系统频响的幅度响应和相位响应，并画图。

（4）编程得到系统的零极点图，并分析系统的因果性和稳定性。

（5）若系统为y(n)=0.45x(n)+0.5x(n-1)+0.45x(n-2)+0.53y(n-1)-0.46y(n-2)

在其上述解果。

三、实验程序与结果

（1）y(n)-1.6y(n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)系统的冲击响应、输

出信号、幅频曲线、相频曲线、零极点图

clear;clc;closeall;

b=[0.50.1];

a=[1-1.61.28];

h=impz(b,a,50);

stem(h);title(冲击响应)

x=[12345zeros(1,50)];

y=conv(x,h);

figure

subplot(2,1,1)

stem(y);title(conv方法输出信号)

y1=filter(b,a,x);

subplot(2,1,2)

stem(y1);title(filter方法输出信号)

[h,f]=freqz(b,a,512,2000);

mag=abs(h);

ph=angle(h)*180/pi;

figure

subplot(2,1,1),plot(f,mag),xlabel(频率),ylabel(幅度);title(幅频曲线)

subplot(2,1,2),plot(f,ph);xlabel(频率),ylabel(相位);title(相频曲线)

[zpk]=tf2zp(b,a);

figure

zplane(z,p,k);title(零极点图);

结果:

（2）y(n)=0.45x(n)+0.5x(n-1)+0.