数字信号处理实验指导书2016试题.doc

PAGE - 2 - PAGE - 27 - 《数字信号处理》 实验指导书 王宏志 吕洪武 张李梅 孙慧然 计算机科学与工程学院 2016年3月 目录  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc382845360 网络资源-基于WEB的数字信号处理实验教学系统  PAGEREF _Toc382845360 \h - 1 -  HYPERLINK \l _Toc382845361 实验一 熟悉MATLAB语言环境  PAGEREF _Toc382845361 \h - 2 -  HYPERLINK \l _Toc382845362 实验二 离散系统分析  PAGEREF _Toc382845362 \h - 7 -  HYPERLINK \l _Toc382845363 实验三 用FFT进行信号的频谱分析  PAGEREF _Toc382845363 \h - 9 -  HYPERLINK \l _Toc382845364 实验四 用窗函数设计FIR数字滤波器  PAGEREF _Toc382845364 \h - 12 -  HYPERLINK \l _Toc382845365 实验五 设计IIR数字滤波器  PAGEREF _Toc382845365 \h - 13 -  HYPERLINK \l _Toc382845366 实验六 随机功率谱估计及MATLAB实现  PAGEREF _Toc382845366 \h - 18 -  HYPERLINK \l _Toc382845367 附录： MATLAB简介  PAGEREF _Toc382845367 \h - 19 - 网络资源 基于WEB的数字信号处理实验教学系统 / 单击”软件下载”即可下载实验系统,如下图所示: 软件安装界面如下所示: 软件运行界面如下所示: 实验一 熟悉MATLAB语言环境 一、实验目的 熟悉Matlab的基本使用方法，重点掌握常用于数字信号处理的相关指令。 利用Matlab实现序列的显示，运算等，加深对信号处理原理课程所学内容的理解。 二、实验内容及步骤 1. 双击Matlab 图标，进入主窗口，如图1所示。其中右侧的是命令窗口(Command Window)。Current Directory显示的是当前的目录，如果要在命令窗口中直接调用一个M文件，则必须保证其在这个目录下。不同的机器上，Current Directory可能不同。 图1 进入Matlab时的画面 2. 下面的操作都在命令窗口(Command Window)中进行。 键入help subplot 并回车，将会看到命令subplot的使用方法的帮助。以后遇到不会的命令可以用键入“ help 某命令”的方式来查看其使用方法。 subplot的用法 ：subplot(m,n,p)是将一个窗口分成m*n个小窗口，p是小窗口的编号，方向是从左至右。 键入help stem 并回车，将会看到命令stem的使用方法的帮助。stem用来画离散序列的柄状图。 3. 下面我们新建立一个M文件。首先点击图1所示的主菜单File下的“新建”图标，将弹出图2 所示的窗口。接着在这个弹出的窗口点击“保存”图标，建立一个文件impseq.m，将其保存在Current Directory下（当前的目录，见图1）。 图2 新建一个文件impseq.m 4. 我们在这个impseq.m文件中输入如下语句并保存，注意以 %开头来写注释。 function[x,n]=impseq(np,ns,nf) % 生成x(n)=delta(n-np);ns≤n≤nf，即单位冲激序列δ(n) % np代表脉冲的位置, ns为序列起始位置，nf为序列终止位置。 % 调用方式[x,n]=impseq(np,ns,nf) if nsnp|nsnf|npnf error(输入位置参数不满足ns=n=nf) else n=[ns:nf]; x=[(n-np)==0];%if判断语句到这里结束 %上面的这条语句是关键语句，作用是如果满足条件(n-np)==0，则x=1。 %在n=ns:nf的一串值中，只有一个值会满足这个逻辑式，因而只在这个n=np %处,x为1，其余的n值 处，x均为0.这样就构成了延迟np的单位冲激序列 end 5.下面我们再回到图1所示的窗口 1)打算显示一个序列x(n)=1.5*δ(n+1)－δ(n－3),下面的操作都在命令窗口(Command Window)中进行。注意：如果一个语句后面有分号，就不显示中间结果。 n1=[-4:5