数字信号处理实验指导书4要点.docVIP

《》目录 网络资源-基于WEB的数字信号处理实验教学系统 - 1 - 实验一 熟悉MATLAB语言环境 - 2 - 实验二 离散系统分析 - 7 - 实验三 用FFT进行信号的频谱分析 - 9 - 实验四 用窗函数设计FIR数字滤波器 - 12 - 实验五 设计IIR数字滤波器 - 13 - 实验六 随机功率谱估计及MATLAB实现 - 18 - 附录： MATLAB简介 - 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]; x1=1.5*impseq(-1,-4,5)-impseq(3,-4,5); % 列出x1序列 subplot(2,2,1); stem(n1,x1,.); stem(n1,x1,.); title(x(n) 的序列图) ylabel(x1(n)); axis([-5,5,-2,3]);text(5.5,-2,n); 显示的结果如图3所示。 图3 序列x(n)=1.5*δ(n+1)－δ(n－3) 2)打算显示另外一个序列 x(n)=n[u(n)－u(n－8)]－10exp(－0.3(n－10))*[u (n－10)—u(n－16)]，0≤n≤20 新建一个文件stepseq.m,用于生成延迟的单位阶跃序列，其输入参数为序列起始位置ns,序列终止位置nf，及阶跃位置np.将刚才编辑的impseq.m另存为stepseq.m，并改动两条语句： function[x,n]=stepseq(np,ns,nf) x=[(n-np)=0]; 存盘。 下面的操作都在命令窗口(Command Window)中进行。 n2=[0:20];x21 = n2.*(stepseq(0,0,20)-stepseq(8,0,20)); % 列出x21序列 x22 = 10*exp(-0.3*(n2-10)).*(stepseq