数字信号处理--实验一 熟悉MATLAB环境

1、学生实验报告开课学院及实验室： 年 月 日学院机械与电气工程学院年级、专业、班姓名学号实验课程名称数字信号处理实验成绩实验项目名称实验一 熟悉MATLAB环境指导老师一、实验目的 1熟悉MATLAB的主要操作命令。2学会用MATLAB创建时域离散信号。3学会创建MATLAB函数。二、实验原理参阅附录MATLAB基本操作及常用命令。三、使用仪器、材料1、硬件：计算机 2、软件：Matlab四、实验步骤完成以下操作。1数组的加、减、乘、除运算。 输入A=1 2 3 4；B=3 4 5 6；计算：C=A+B；D=A-B；E=A.*B；F=A./B；G=A.B；并用stem语句画出A、B、C、D、E、

2、F、G。2用MATLAB实现以下序列(1)单位抽样序列(2)单位阶跃序列(3)矩形序列(4)正弦序列x(n)=5sin(0.5n+ /4)(5)指数序列x(n)=exp(-0.5n)3用MATLAB生成以下两个序列：并作以下运算，并绘制运算后序列的波形。（1）（2）（3）（4）（5）4利用MATLAB读取一个WAV文件，并画出其波形图。将此WAV文件的信号幅度衰减一半后再存为另一个WAV文件，并画出其波形图。五、实验过程原始记录（数据、图表、计算等） 1.程序如下：A=1 2 3 4;B=3 4 5 6;C=A+B;D=A-B;E=A.*B;F=A./B;G=A.B;stem(A, .);ti

3、tle(A);stem(B, .);title(B);stem(C, .);title(C);stem(D, .);title(D);stem(E, .);title(E);stem(F, .);title(F);stem(G, .);title(G);运行结果：2.（1）程序如下：%取n0=0 （2）程序如下：%取n0=1n0=0;n1=-10;n2=10; n0=1;n1=-10;n2=10;n=n1:n2; n=n1:n2;delta=(n-n0)=0; u=(n-n0)=0;stem(n,delta,.) %产生单位抽样序列 stem(n,u,.) %产生单位阶跃序列xlabel(n)

4、,ylabel(delta(n) xlabel(n),ylabel(u(n)运行结果： 运行结果：（3）程序如下: %取N=6 （4）程序如下:N=6;n1=-10;n2=10; n=-10:0.2:10;n=n1:n2; x=5*sin(0.5*pi*n+pi/4);R=(N-1)=n)&(n=0); stem(n,x,.) %产生正弦序列stem(n,R,.) %产生矩形序列 xlabel(n)xlabel(n),ylabel(R(n) ylabel(x(n)运行结果： 运行结果：（5）程序如下：n=-10:10;x=exp(-0.5*n);stem(n,x,.) %产生指数序列xlabe

5、l(n)ylabel(x(n)运行结果：3.生成x(n)、h(n)序列的程序：n=-10:10;x=zeros(1,10),1,2,3,4,5,zeros(1,6);subplot(1,2,1);stem(n,x,.); title(x(n); %产生x(n)序列n=-10:10;h=zeros(1,10),1,2,1,2,zeros(1,7);subplot(1,2,2);stem(n,h,.); title(h(n); %产生h(n)序列运行结果：产生(1)-(5)的程序：n=-10:10;x=zeros(1,10),1,2,3,4,5,zeros(1,6);h=zeros(1,10),1

6、,2,1,2,zeros(1,7);x1 =zeros(1,15),1,2,3,4,5,zeros(1, 1);subplot(3,2,1);stem(n,x1,.); title(x(n-5); %产生x(n-5)序列x2=zeros(1, 5),1,2,3,4,5,zeros(1,11);subplot(3,2,2);stem(n,x2,.); title(x(n+5); %产生x(n+5)序列x3= fliplr(x);subplot(3,2,3);stem(n,x3,.); title(x(-n); %产生x(-n)序列x4=x+h;subplot(3,2,4);stem(n,x4,.

7、); title(x(n)+h(n); %产生x(n)+h(n)序列x5=3*x;subplot(3,2,5);stem(n,x5,.); title(3x(n); %产生3x(n)序列x6=x.*h;subplot(3,2,6);stem(n,x6,.); title(x(n)h(n); %产生x(n)h(n)序列运行结果：4. 将WAV文件放在E盘，路径为“E:玻璃破碎.wav”。程序如下：y,fs,bits=wavread(E:玻璃破碎.wav); %读取文件subplot(1,2,1);plot(y);title(未处理); s=1/2*y;subplot(1,2,2);plot(s);title(处理后);wavwrite(s,E:玻璃破碎副本.wav) %形成新的声音文件运行结果：六、实验结果及分析 本次实验中，通过MATLAB软件编程，熟悉了一些基本的操作，包括数组的产生和运算、一些常见序列的产生和运算及基本函数的使用。通过本次实验，我对于MATLAB用于信号处理的操作有了一定的了解，之前也曾在信号与系统实