数字信号处理实验报告材料.doc

PAGE 数字信号处理实验讲义 PAGE 1 TOC \o "1-3" \h \z \u HYPERLINK \l "_Toc136100057" 前 言 PAGEREF _Toc136100057 \h 2 HYPERLINK \l "_Toc136100058" 实验一 MATLAB简介 PAGEREF _Toc136100058 \h 3 HYPERLINK \l "_Toc136100059" 实验二 用FFT实现信号的谱分析 PAGEREF _Toc136100059 \h 5 HYPERLINK \l "_Toc136100060" 实验三 IIR数字巴特沃思滤波器的设计 PAGEREF _Toc136100060 \h 8 HYPERLINK \l "_Toc136100061" 实验四 FIR数字滤波器的设计 PAGEREF _Toc136100061 \h 9 PAGE 21 实验一 MATLAB简介 实验目的 1．熟悉MATLAB软件的使用方法； 2．MATLAB的绘图功能； 3．用MATLAB语句实现信号的描述及变换。 实验原理 1．在MATLAB下编辑和运行程序 在MATLAB中，对于简单问题可以在命令窗（command windows）直接输入命令，得到结果；对于比较复杂的问题则可以将多个命令放在一个脚本文件中，这个脚本文件是以m为扩展名的，所以称之为M文件。用M文件进行程序的编辑和运行步骤如下： （1）打开MATLAB，进入其基本界面； （2）在菜单栏的File项中选择新建一个M文件； （3）在M文件编辑窗口编写程序； （4）完成之后，可以在编辑窗口利用Debug工具调试运行程序，在命令窗口查看输出结果；也可以将此文件保存在某个目录中，在MATLAB的基本窗口中的File项中选择Run The Script，然后选择你所要运行的脚本文件及其路径，即可得出结果；也可以将此文件保存在当前目录中，在MATLAB命令窗口，“>>”提示符后直接输入文件名。 2．MATLAB的绘图功能 plot(x,y) 基本绘图函数，绘制 x和y之间的坐标图。 figure(n ) 开设一个图形窗口n subplot(m,n,N) 分割图形窗口的MATLAB函数，用于在一个窗口中显示多个图形，将图形窗口分为m行n列，在第N个窗口内绘制图形。 axis([a0,b0,a1,b1] ) 调整坐标轴状态 title(‘ ’) 给图形加题注 xlabel (‘ ‘) 给x轴加标注 ylabel (‘ ‘) 给y轴加标注 grid 给图形加网格线 3．信号描述及变换 信号描述及变换包括连续时间信号和离散时间信号内容，详细内容请见课本第1章、第2章。 实验内容 1．试用MATLAB绘制出下列信号的波形：(Signal 1.6) ； ； ； （5） 【程序代码】 clear all;close all;clc; syms t; x1=exp(-1.5*t) x2=3*sin(0.5*pi*t) x3=0.5+0.5*sym(('sign(t)')) x4=sym('heaviside(t)')+sym('heaviside(t-1)')-sym('2*heaviside(t-2)') x5=.5*t*(sym('heaviside(t)')-sym('heaviside(t-4)')) subplot(2,3,1); ezplot(x1); axis([-6 3 -500 7000]); title('x1(t)=exp(-1.5t)'); grid on subplot(2,3,2); ezplot(x2); title('x2(t)=3sin(0.5|Dt)'); grid on subplot(2,3,3); fplot('sign(t)/2 + 1/2',[-10 10],1e-8); ezplot(x3,[-10 10]); axis([-10 10 -.2 1.2]); xlabel('t');title('x3(t)=0.5+0.5sgn(t)'); grid on subplot(2,3,4); ezplot(x4,[-1 3]); title('x4(t)=u(t)+u(t-1)-2u(t-2)'); grid on subplot(2,3,5); ezplot(x5,[-2 6]); title('x5(t)=0.5t[u(t)-u(t-4)]'); grid on subplot(2,3,6); axis off 2． 已知连续时间信号(Signal 1.7