《数字信号处理》实验指导书
150010005000012345678910150010005000012345678910
2、用直接Ⅰ型结构实现滤波器: 已知某滤波器的差分方程(见下),试用直接Ⅰ型结构实现该滤波器,求冲激响应。
y(n)=x(n)+1.4x(n-1)+0.3x(n-2)-0.3y(n-1)+0.2y(n-2),
具体操作步骤:
(4) 打开IIR02.mdl文件,窗口中已有所需模块,请按下图连接各模块;
x(n)¥?μ?DòáD1z1.41zy(n)y1??ê3á?1?±£′??±μ??á?D1z-0.31z0.30.2
(5) 用鼠标点击仿真窗口工具条中的?图标开始仿真;
(6)回到Matlab命令窗口中,输入以下程序,观察冲激响应序列; 程序代码
n=0:10; stem(n,y1);
程序说明
作出冲激序列
- 15 -
《数字信号处理》实验指导书
1500100050000123456789101.510.50-0.5012345678910
3、用直接Ⅱ型结构实现滤波器:
试用直接Ⅱ型结构实现内容2中的滤波器,求冲激响应。具体操作步骤: (1) 打开IIR03.mdl文件,窗口中已有所需模块,请按下图连接各模块;
x(n)¥?μ?DòáD1z-0.31zy(n)y2??ê3á?1?±£′??±μ??á?D1.40.20.3
(2)用鼠标点击仿真窗口工具条中的?图标开始仿真;
(3)回到Matlab命令窗口中,输入以下程序,观察冲激响应序列(比较结果); 程序代码
subplot(2,1,1); stem(n,y1); subplot(2,1,2); stem(n,y2);
- 16 -
程序说明
直接Ⅰ型冲激序列
直接Ⅱ型冲激序列
《数字信号处理》实验指导书
1.510.50-0.50123456789101.510.50-0.5012345678910
4、用级联型结构实现滤波器(选作):
已知某滤波器的系统函数H(z)可进行如下分解,试用级联型结构实现该滤波器
H(z)?4(1?z(1?0.5z?4??1)(1?1.4z?1?z?2)?2?1)(1?0.9z??1?0.8z?1)?2?21?z?1?11?1.4z1?0.9z?z1?0.5z?1?0.8z
?4?H1(z)?H2(z)(1) 打开Filter04.mdl文件,点击?图标进行仿真;
(2) 回到Matlab命令窗口中,输入以下程序,观察冲激响应序列; 程序代码
stem(n,y3(1,1:11));
程序说明
级联型冲激序列
- 17 -
《数字信号处理》实验指导书
1.510.50-0.50123456789101050-5012345678910
(3) 打开Filter05.mdl文件(交换子系统的顺序),点击?图标进行仿真; (4) 回到Matlab命令窗口中,输入以下程序,观察冲激响应序列(比较结果); 程序代码 subplot(2,1,1); stem(n,y3(1,1:11)); subplot(2,1,2); stem(n,y4(1,1:11));
程序说明
- 18 -
《数字信号处理》实验指导书
1050-50123456789101050-5012345678910
七、实验报告要求
1、实验步骤按实验内容指导进行;
2、对于实验内容1、2和3的数据必须给出的仿真模块图和冲激响应图; 3、实验内容4仅作为练习,不用给出相关数据;; 4、实验报告按规定格式填写,要求如下: (1)实验步骤根据自己实际操作填写;
(2)各小组实验数据不能完全相同,否则以缺席论处;
5、实验结束,实验数据交指导教师检查,得到允许后可以离开,否则以缺席论处;
八、实验注意事项
1、Matlab编程、文件名、存盘目录均不能使用中文。 2、实验结果要保存,打印后贴在实验报告上。 3、仿真模块图的拷贝
(1)在仿真界面下,先把整个窗口缩放到合适的大小; (2)点击菜单Edit / Copy model to clipboard; (3)在Word界面下,进行粘贴操作;
九、思考题
- 19 -
《数字信号处理》实验指导书
1、要得到一个方便调节零极点的IIR滤波器应用什么结构,为什么? 注释 无
主要参考文献:
[1]《MATLAB符号运算及应用》,黄忠霖编,国防工业出版社,2004
- 20 -