资源描述:
《数字信号处理实验报告》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数字信号处理实验讲义前言2实验一MATLAB简介3实验二用FFT实现信号的谱分析5实验三IIR数字巴特沃思滤波器的设计8实验四FIR数字滤波器的设计911实验一MATLAB简介实验目的1.熟悉MATLAB软件的使用方法;2.MATLAB的绘图功能;3.用MATLAB语句实现信号的描述及变换。实验原理1.在MATLAB下编辑和运行程序在MATLAB中,对于简单问题可以在命令窗(commandwindows)直接输入命令,得到结果;对于比较复杂的问题则可以将多个命令放在一个脚本文件中,这个脚本文件是以m为扩展名的,所以称之为M文件。用M文件进
2、行程序的编辑和运行步骤如下:(1)打开MATLAB,进入其基本界面;(2)在菜单栏的File项中选择新建一个M文件;(3)在M文件编辑窗口编写程序;(4)完成之后,可以在编辑窗口利用Debug工具调试运行程序,在命令窗口查看输出结果;也可以将此文件保存在某个目录中,在MATLAB的基本窗口中的File项中选择RunTheScript,然后选择你所要运行的脚本文件及其路径,即可得出结果;也可以将此文件保存在当前目录中,在MATLAB命令窗口,“>>”提示符后直接输入文件名。2.MATLAB的绘图功能plot(x,y)基本绘图函数,绘制x和y
3、之间的坐标图。figure(n)开设一个图形窗口nsubplot(m,n,N)分割图形窗口的MATLAB函数,用于在一个窗口中显示多个图形,将图形窗口分为m行n列,在第N个窗口内绘制图形。axis([a0,b0,a1,b1])调整坐标轴状态title(‘’)给图形加题注xlabel(‘‘)给x轴加标注ylabel(‘‘)给y轴加标注grid给图形加网格线3.信号描述及变换信号描述及变换包括连续时间信号和离散时间信号内容,详细内容请见课本第1章、第2章。实验内容1.试用MATLAB绘制出下列信号的波形:(Signal1.6)(1);(2)(
4、3);21(1);(5)【程序代码】clearall;closeall;clc;symst;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([-63-5007000]);ti
5、tle('x1(t)=exp(-1.5t)');gridonsubplot(2,3,2);ezplot(x2);title('x2(t)=3sin(0.5¦Ðt)');gridonsubplot(2,3,3);fplot('sign(t)/2+1/2',[-1010],1e-8);ezplot(x3,[-1010]);axis([-1010-.21.2]);xlabel('t');title('x3(t)=0.5+0.5sgn(t)');gridonsubplot(2,3,4);ezplot(x4,[-13]);title('x4(t)=
6、u(t)+u(t-1)-2u(t-2)');gridonsubplot(2,3,5);ezplot(x5,[-26]);title('x5(t)=0.5t[u(t)-u(t-4)]');gridonsubplot(2,3,6);axisoff2.已知连续时间信号(Signal1.7),,试用MATLAB绘制出下列信号的波形:(1);(2);(3);(4);(5)。21【程序代码】clearall;closeall;clc;figure(2)symst;x1=(4-t)*(sym('heaviside(t)')-sym('heaviside
7、(t-4)'))x2=exp(-2*t)*sym('heaviside(t)')x3=sin(2*pi*t)x4=subs(x1,t,t/2)x5=subs(x4,t,t-2)x6=subs(x2,t,-t)x7=x2+x6x8=x7*x3subplot(2,3,1)text(0,0.9,'x1(t)=(4-t)[u(t)-u(t-4)]');text(0,0.7,'x2(t)=exp(-2t)u(t)');text(0,0.5,'x3(t)=sin(2¦Ðt)');axisoff;boxoff;subplot(2,3,2);ezplot
8、(x5,[-3,10]);gridtitle('x4(t)=x1(t/2)');subplot(2,3,3);ezplot(x5,[-1,12]);gridtitle('x5(t)=x4(t-