1 / 4 数 字 信 号 处 理 实 验 报 告实验项目名称:离散系统的 MATLAB 实现实验日期:2012 年 10 月 11日实验成绩:实验评定标准:1)实验结果是否正确A()B()C()2)实验结果分析A()B()C()3)实验报告是否按照规定格式A()B()C()一、 实验目的离散系统的模型时域表示的MATLAB 实现函数有 filter 和 impz函数两种,本实验通过以上函数实现离散系统的时域模型。二、 实验器材PC 机一台MATLAB 软件三、 实验内容1.filter 函数:利用递归滤波器或非递归滤波器对数据滤波,将一个离散系统可以看成一个滤波器, 系统的输出就是输入经过滤波器的结果。●y=filer(b,a,x); 说明: y 表示由向量 b 和 a 组成的系统对输入信号 x 进行滤波,系统的输出为y。2 / 4 2.impz 函数:直接给出系统的单位冲击响应。●impz(b,a);实验要求:(1)当系统的输入输出差分方程为:0.9(n)-y(n-1)-0.8y(n-2)=0.3x(n)+0.7x(n-1)+0.7x(n-2),并且系统的输入为 N 点单位冲激序列的时候,分别用filter 和 impz 函数得到系统单位冲激相应的曲线,并对其进行比较。四、 实验结果系统的输入方程为: y(n)-0.8y(n-1)-0.5y(n-2)=0.7x(n)+0.3x(n-1)程序: clear all; Pulse=[1,zeros(1,63)]; B=[0.7 0.3]; A=[1 -0.8 - 0.5]; H1=filre(b,a,pulse); H2=impz(b,a,4); Suplot(2,1,1),stem(h1),title( ‘filer function’); Suplot(2,1,2),stem(h2),title( ‘impz function’); 系统的输入方程为: 0.9y(n)-y(n-1)-0.8y(n-2)=0.3x(n)3 / 4 +0.7x(n-1)+0.7x(n-2) :clear all; Pulse=[1,zeros(1,63)]; B=[0.3 0.7 0.7]; A=[0.9 -1 - 0.8]; H1=filre(b,a,pulse); H2=impz(b,a,4); Suplot(4,1,1),stem(h1),title( ‘filer function’); Suplot(4,1,2),stem(h2),title( ‘impz function’); 系统的输入方程为:y(n)-0.9y(n-1)-0.3y(n-2)=0.6x(n)+0.2x(n-1)clear all; Pulse=[1,sin(1*63)]; B=[0.6 0.2]; A=[ 1 -0.9 -0.3]; H1=filre(b,a,pulse); H2=impz(b,a,4); Suplot(2,1,1),stem(h1),title( ‘filer function’); 4 / 4 五、 实验结果分析y(n)-0.8y(n-1)-0.5y(n-2)=0.7x(n)+0.3x(n-1)0.9y(n)-y(n-1)-0.8y(n-2)=0.3x(n)+0.7x(n-1)+0.7x(n-2) y(n)-0.9y(n-1)-0.3y(n-2)=0.6x(n)+0.2x(n-1)