数字信号处理实验之离散系统的差分方程、冲激响应和卷积分析.docxVIP

数字信号处理实验之离散系统的差分方程、冲激响应和卷积分析 《数字信号处理 A 》实验报告 实验二 实验名称 : 离散系统的差分方程、冲激响应与卷积分析 专业及班级 :电子 131 姓名 :XXX 学号 :XXXXXX 一、实验目的 加深对离散系统的差分方程、冲激响应与卷积分析方法的理解。 二、实验步骤 (附源代码及仿真结果图 ) 1、以下程序中分别使用 conv 与 filter 函数计算 h 与 x 的卷积 y 与 y1,运行程序 ,并分析 y 与 y1 就是否有差别 ,为什么要使用 x[n] 补零后的 x1 来产生 y1;具体分析当 h[n] 有 i 个值 ,x[n] 有 j 个值 ,使用 filter 完成卷积功能 ,需要如何补零？ 编制程序求解下列两个系统的单位冲激响应与阶跃响应 ,并绘出其图形。要求分别用 filter 、 conv、 impz  三种函数完成。 y[n]  0.75 y[ n  1]  0.125 y[ n  2]  x[ n]  x[ n  1] 给出理论计算结果与程序计算结果并讨论。 单位冲激响应 : (1) 用 filter 函数 a1=[1,0 、 75,0、125]; b1=[1,-1]; n=0:20; x1=[1 zeros(1,20)]; y1filter=filter(b1,a1,x1); stem(n,y1filter); title(y1filter); xlabel(x); ylabel(y); 数字信号处理实验之离散系统的差分方程、冲激响应和卷积分析 y1filter 1.5 1 0.5 0 y -0.5 -1 -1.5 -2 0 2 4 6 8 10 12 14 16 18 20 x (2) 用 conv 函数 a1=[1,0 、 75,0、125]; b1=[1,-1]; x1=[1 zeros(1,10)]; [h]=impz(b1,a1,10); y1conv=conv(h,x1); n=0:19; stem(n,y1conv,filled) 1.5 1 0.5 0 -0.5 -1 -1.5 -2 0 2 4 6 8 10 12 14 16 18 20 数字信号处理实验之离散系统的差分方程、冲激响应和卷积分析 (3) 用 impz 函数 a1=[1,0 、 75,0、125]; b1=[1,-1]; impz(b1,a1,21); Impulse Response 1.5 1 0.5 e 0 d u p m -0.5 -1 -1.5 -2 0 2 4 6 8 10 12 14 16 18 20 n (samples) 单位阶跃响应 : 用 filter 函数 a1=[1,0 、 75,0、125]; b1=[1,-1]; n=0:20; x2=ones(1,21); y1filter=filter(b1,a1,x2); stem(n,y1filter); title(y1filter_step); xlabel(x); ylabel(y); y1filter s tep 1 0.8 0.6 0.4 0.2 y 0 -0.2 -0.4 -0.6 -0.8 0 2 4 6 8 10 12 14 16 18 20 x 数字信号处理实验之离散系统的差分方程、冲激响应和卷积分析 (2) 用 conv 函数 a1=[1,0 、 75,0、125]; b1=[1,-1]; x2=ones(1,21); [h]=impz(b1,a1,20); y1=conv(h,x2); y1conv=y1(1:21); n1=0:20; stem(n1,y1conv,filled); title(y1conv); xlabel(n); ylabel(y1[n]); y1conv 1 0.8 0.6 0.4 ] 0.2 n [ 1 y 0 -0.2 -0.4 -0.6 -0.8 0 2 4 6 8 10 12 14 16 18 20 n 2、编制程序求解下列两个系统的单位冲激响应与阶跃响应 ,并绘出其图形。要求分别用 filter 、 conv、 impz 三种函数完成。 y[n] 0.25{ x[ n 1] x[ n 2] x[ n 3] x[ n 4]} 给出理论计算结果与程序计算结果并讨论 a、单位冲激响应 : (1) 用 filter 函数 a1=[1]; b1=[0 、 25,0、 25,0、 25,0、 25]; n=0:20; x1=[1 zeros(1,20)]; y1filter=filter(b1,a1,x1); stem(n,y1filter); title(y1filter); xlabel(x); ylabel(y); 数字信号处理实验之离散系