实验二离散系统的差分方程、冲激响应和卷积分析.docx

广西科技大学实验报告 数字信号处理实验 Guangxi University of Science and Technology实验报告实验课程： 数字信号处理 实验内容：离散系统的差分方程、冲激响应和卷积分析 院 （系）： 计算机科学与通信工程学院 专 业： 通信工程 班 级： 学生姓名： 学 号： 指导教师： 2016年 5 月 25 日 实验二、离散系统的差分方程、冲激响应和卷积分析一、实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。二、实验原理：离散系统其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号，。记系统单位冲激响应，则系统响应为如下的卷积计算式： 当时，hn是有限长度的（n：0，M），称系统为FIR系统；反之，称系统为IIR系统。 在MATLAB中，可以用函数y=Filter(p,d,x) 求解差分方程，也可以用函数 y=Conv(x,h)计算卷积。三、实验内容及步骤：1求yn+0.75yn-1+0.125yn-2=xn-xn-1的系统响应以及零极点分布图。程序如下：% 用impz()函数求冲激响应：M=60;num=1 -1 0;den=1 0.75 0.125;y=impz(num,den,M);subplot(2,2,1);stem(y);xlabel(时间);ylabel(振幅);title(单位冲激响应：impz()函数方法)n=0:60;x=cos(2*pi*0.47*n);num=1 -1 0;den=1 0.75 0.125;y=filter(num,den,x); %注意：“”后加一个空格再写指令，这样就可以直subplot(2,2,2);stem(n,y)xlabel(时间);ylabel(振幅);title(余弦输入响应: filter()函数方法)% 系统零极点分布图：num=1 -1 0;den=1 0.75 0.125;subplot(2,2,3);zplane(num,den);grid图: 2-12.求yn=0.25xn-1+xn-2+xn-3+xn-4系统零极点分布图。程序如下：%用impz()函数求冲激响应：M=60;num=0.25 0.25 0.25 0.25;den=1 0 0 0;y=impz(num,den,M);subplot(2,2,1);stem(y);xlabel(时间);ylabel(振幅);title(单位冲激响应：impz()函数方法)n=0:60;x=cos(2*pi*0.47*n);num=0.25 0.25 0.25 0.25;den=1 0 0 0;y=filter(num,den,x); %注意：“”后加一个空格再写指令，这样就可以直subplot(2,2,2);stem(n,y)xlabel(时间);ylabel(振幅);title(余弦输入响应: filter()函数方法) % 系统零极点分布图：num=0.25 0.25 0.25 0.25;den=1 0 0 0;subplot(2,2,3);zplane(num,den);grid实验结果：图: 2-2四、实验总结与分析：通过本次实验，基本学会了用MATLAB软件编程求离散