数字信号处理实验指导书(修订版)要点.docVIP

目录 实验一 卷积实验…………………………………………… 1 实验二DFT和FFT实验…………………………………… 5 实验三 用双线性变换法设计IIR数字滤波器……………11 实验四 用窗函数法设计FIR数字滤波器…………………16 实验一 卷积实验 一.实验目的 1.熟悉离散信号和系统的MATLAB 表示和产生方法； 2.熟悉线性卷积的MATLAB编程方法； 3 利用卷积方法观察分析系统的时域特性。 二 实验原理 (一) 离散时间信号的MATLAB表示 序列x(n)={2,6,1,2,0,3,4,5,6},用MATLAB语句表示为： n=[-4 -3 -2 -1 0 1 2 3 4]; x=[ 2 6 1 2 0 3 4 5 6]; 其中：n表示取样时间点，x表示对应时间点的值。 如果不需要采样位置信息或这个信息是多余的时候(例如序列从n=0开始)，用户可以只使用x向量来表示序列。 1.单位样值序列产生函数 IMPSEQ.M function[x,n]=impseq(n0,n1,n2) %产生x(n)=delta(n-n0);n1=n=n2; n=[n1:n2]; x=[(n-n0)==0]; 在命令行窗口里执行： x= impseq(20,0,100); stem(x) 有： 2.单位阶跃序列产生函数 STEPSEQ.M function[x,n]=stepseq(n0,n1,n2) %产生x(n)=u(n-n0);n1=n=n2; n=[n1:n2]; x=[(n-n0)=0]; 在命令行窗口里执行： x=stepseq(20,0,100); stem(x) 有： 3正余弦信号的产生： % x(n)=1.5cos(0.02*pi*n) -150=n=150 n=-150:150 x=1.5*cos(0.02*pi*n); stem(n,x) 执行结果: 4.产生均值为0，方差为1的高斯随机噪声序列。 n=-150:150; Noise=randn(1,301); stem(n, Noise) 执行结果: (二) 序列的卷积 一个线性系统一般应满足下式： y(n)=T[x(n)] 由数字信号处理理论可知，线性系统的脉冲响应即为h(n)，上式一般称为线性卷积，一般可以表示为： y(n)=x(n)*h(n) 如果任意序列是无限长度的，就不能用MATLAB来直接计算卷积。MATLAB提供了一个内部函数conv( )来计算两个有限长序列的卷积。 conv( )函数假定两个序列都从n=0开始，调用格式如下所示： y=conv(x,h) 例：设x(n)=[2,3,6,5,7,6,4],h(n)=[1,2,1,2,0,1],则 y(n)=x(n)*h(n) 可由以下语句实现： x=[2,3,6,5,7,6,4]; h=[1,2,1,2,0,1]; y=conv(x,h) y = Columns 1 through 8 2.00 7.00 14.00 24.00 29.00 39.00 36.00 34.00 Columns 9 through 12 21.00 15.00 6.00 4.00 三 实验内容 分别实现单位样值序列产生函数, 单位阶跃序列产生函数, 正余弦信号和高斯随机噪声序列; 2.计算线性卷积 任意产生有限长序列x1(n)和x2(n)，用conv（）函数计算它们的线性卷积，会出结果的图形。 3.已知某系统的输入序列为x(n)=0.9nu(n),系统的单位样值响应为h(n)=[1,0.3,1,0.8,0,1]，要求： 编程产生上述序列; 计算序列x(n)输入系统后的响应,绘出0≤n≤100时的输出序列。 四 实验报告要求 (1) 简述实验目的及实验原理。 (2) 按实验内容要求,给出Matlab 程序,附上实验过程中的信号序列、系统响应序列的图形。 (3) 总结实验中的主要结论。 实验二 DFT和FFT实验 一.实验目的 1加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法， 所以FFT的运算结果必然满足DFT的基本性质)。 2掌握DFT、IDFT的MATLAB编程方法； 3