北京邮电大学MATLAB实验报告

DSP--MATLAB实验报告

实验名称：MATLAB第一次上级实验

学生姓名：

班级：

班内序号：

学号：

日期：

1．实验要求

(1) 常用数字信号序列的产生：

熟悉Matlab 产生数字信号的基本命令，加深对数字信号概念的理解，并能够用Matlab 产生和绘制出一些常用离散信号序列。请用Matlab 画出下列序列的波形（-10

a) δ(n)

b) 单位阶跃序列2 u(n-5)

c) 矩形序列R(n)

d) y(n)＝2sin(0.3πn)+ 0.5cos2(0.6πn)

（2）加、减、尺度（乘除）和移位是数字信号处理中最基本的算术运算，将上述基本序列进行这些基本运算，得到多个序列构成的组合序列。

（3）请用您的计算机声卡采用一段您自己的声音x(n)，长度为45秒，单声道，取样频率44.1kHz，16bit/样值，然后与给定的一段背景音乐y(n) 按下式叠加为一个声音信号z(n)：z(n) = 0.7x(n) + 0.3y(n)

要求在同一个Figure 中,画出采集声音x(n)、背景音乐y(n)和混音z(n) 的时域波形；

2关键算法分析

产生δ(n)，单位阶跃序列 2 u(n-5)，矩形序列R(n)，y(n)＝2sin(0.3πn)+ 0.5cos2(0.6πn)序列，首先在matlab中新建一个m文件，将各种序列的定义写进去，再在命令窗口中写入命令语句，便可得到。

单位取样序列x=impseq(0,-20,120);n=[-20:120]; >> plot(n,x)

单位阶跃序列 2 u(n-5) ：x=stepseq(5,-20,120),n=[-20:120];plot(n,2*x)

矩形序列R(n)：x1=stepseq(-10,-10,60),n1=[-10:60],x2==stepseq(10,-10,60),n2=[-10:60],y=x1-x2,n= [-10:60],plot(n,y)

y(n)＝2sin(0.3πn)+ 0.5cos2(0.6πn)：n=[0:0.1:20];x=2*sin(0.3*pi*n)+0.5*cos(2*0.6*pi*n);plot(n,x)

加法：δ(n)+ 2 u(n-5)：x1=impseq(0,-20,120);n=[-20:120];

x2=2*stepseq(5,-20,120),n=[-20:120];y=x1+x2,plot(n,y)

减法：δ(n)- 2 u(n-5)：x1=impseq(0,-20,120);n=[-20:120]; x2=2*stepseq(5,-20,120),n=[-20:120];y=x1-x2,plot(n,y)

乘法: δ(n)- 2 u(n-5): x1=impseq(0,-20,120);n=[-20:120]; x2=2*stepseq(5,-20,120),n=[-20:120];y=x1.*x2,plot(n,y)

一个figure中显示多个波形：[y,fs,bits]=wavread('我的声音'),t1=(0:length(y)-1)/fs;[x,fs,bits]=wavread('天空之城'),t2=(0:length(x)-1)/fs;subplot(2,2,1),plot(t1,y);subplot(2,2,2),plot(t2,x);

3. 程序运行结果

单位取样序列

单位阶跃序列2 u(n-5)

矩形序列：

y(n)＝2sin(0.3πn)+ 0.5cos2(0.6πn)：

：

加法：δ(n)+ 2 u(n-5)

乘法: u(n)*R5:

移位：

一个figure中显示多个波形:

4. 总结

通过这次实验，我对MATLAB的使用了解了许多，尤其是在对声音和图像的处理方面。MATLAB与其他的编程语言相比，更接近我们现实中的语言，更容易理解，所以学习起来相对容易，但是在编程的过程中，还是因为不了解函数的使用，以及要实现某一功能需要的函数，耗了不少的时间。编程时，我虽然会利用help，但对于每个函数的解释都是英文的，有时读起来有点困难，所以说工科学好英语还是挺重要的，学好MATLAB更是必不可少的。