数字信号处理实析.docVIP

数字信号处理基础实验教程 徐振宇 编 浙江树人大学信息科技学院 2008年5月实验 一 离散时间信号的时域分析 一、实验目的 1、熟悉MATLAB软件的使用 2、掌握离散时间信号的表示方法 二、实验原理 2.1基础知识回顾 R1.1离散时间信号用一个称为样本的数字序列来表示。典型的离散时间信号或序列的样本值用表示，其中自变量的取值范围是从到之间的整数。为了表示方便，序列通常不用大括号表示。 R1.2离散时间信号可以是一个有限长的序列，也可以是一个无限长序列。有限长（也称为有限时宽或有限范围）序列仅定义在有限的时间间隔中： （1.1） 其中而，。有限长序列的长度或时宽为 （1.2） R1.3单位样本序列（通常称为离散时间冲激或单位冲激）用表示，其定义为 （1.3） 单位阶跃序列用表示，其定义为 （1.4） 2.2用到的MATLAB命令 本实验中将用到的MATLAB命令如下： 运算符和特殊符号 ： . + - * / ； % 基本矩阵和矩阵控制 i ones zeros 二维图形 axis stem title xlabel ylabel 通用图形函数 clf 2.3单位样本和单位阶跃序列的实现 产生的两个基本离散时间序列分别是式（1.3）和式（1.4）给出的单位样本序列和单位阶跃。使用MATLAB命令，可产生长度为N的单位样本序列u[n] u=[1 zeros(1,N-1)]； 延时M个样本且长度为N的单位样本序列ud[n]，可用如下MATLAB命令产生，其中MN： ud = [zeros(1,M) 1 zeros(1,N-M-1)]; 同样，使用下面的MATLAB命令可产生长度为N的单位阶跃序列s[n] s=[ones(1,n)]; 产生延时单位阶跃序列的方法，类似于产生延时单位样本序列的方法。 三、实验内容 绘制一个单位样本序列： % Generation of a Unit Sample Sequence clf; % Generate a vector from -10 to 20 n = -10:20; % Generate the unit sample sequence u = [zeros(1,10) 1 zeros(1,20)]; % Plot the unit sample sequence stem(n,u); xlabel(Time index n);ylabel(Amplitude); title(Unit Sample Sequence); axis([-10 20 0 1.2]);  实验二 用窗口法设计FIR滤波器 实验目的 1、掌握MATLAB软件的使用 2、了解一个实际滤波器设计过程，加深掌握用窗口法设计FIR滤波器的原理和方法 3、了解窗函数对滤波器性能的影响 基本原理 设所希望得到的滤波器的理想频率响应为。那么FIR滤波器的设计就在于寻找一个传递函数去逼近。在这种逼近中最直接的一种方法是从单位取样响应序列h(n)着手，使h(n)逼近理想的单位取样响应。又因为 但是一般来说，这样得到的单位取样响应往往都是无限长序列；而且是非因果的。我们以一个截止频率为的线形相位理想低通为例说明。设低通滤波器的时延为，即 这是一个以为中心的偶对称的无限长非因果序列。这样一个无限长的序列怎样用一个有限长序列去近似呢？最简单的办法就是直接截取它的一来代替它。例如把n=0到n=N-1的一段截取来作为h(n)，但是为要保证所得到的是线形相位滤波器。必须满足h(n)的对称性，所以时延应该取h(n)长度的一半，即=（N-1）/2 这种直接截取的办法可以形象地想象为，h(n)好比是通过一个“窗口”所看到的一段。h(n)可表达为和一个“窗口函数”的乘积。在这里，窗口函数就是矩形脉冲函数，即： 但是一般来说，窗口函数并不一定是矩形函数，可以在矩形以内还对作一定的加权处理，因此，一般可以表示为 这里w(n)就是窗口函数。这种对理想单位取样响应加窗的处理对频率响应会产生以下三点影响： 1使理想特性不连续的边沿加宽，形成一过渡带，过渡带的宽度取决于窗口频谱的主瓣宽度。 2在过渡带两旁产生肩峰和余振，它们取决于窗口频谱的旁瓣；旁瓣越多，余振也越多；旁瓣相对值越大，肩峰则越强。 3增加截取长度N，只能缩小窗口频谱的主瓣宽度而不能改变旁瓣的相对值；旁瓣与主