北邮matlab实验一实验报告

数字信号处理Matlab实验数字信号的产生和基本运算

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

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

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

c) 矩形序列 R(n)

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

a)

利用.m文件

function[x,n]=a(n0,n1,n2)

if((n0<n1)||(n1>n2)||(n0>n2))

error('cuowu')

end

n=[n1:n2];

x=[(n-n0)==0];

图像如下图所示

数字信号处理Matlab实验数字信号的产生和基本运算

b)

建立.m文件建立函数u

命令如下：

% step_function.m

% for example ss2_6

function[x,n]=u(n0,n1,n2)

if((n0<n1)||(n0>n2)||(n1>n2))

error('cuowu')

end

n=n1:n2;

x=[(n-n0)>=0];

实现2u(n-5)的命令是

Y=2*u(5,-10,10);

N=-10:10;

Stem(n,y)

图像如下：

c) 命令如下：

数字信号处理Matlab实验数字信号的产生和基本运算

function[x,n]=R(n0,n1,n2,n3)

if((n0<n2)||(n0>n1)||(n1>n3)||(n2>n3))

error('cuowu')

end

n=n2:n3;

x=(n>=n0)-(n>=n1);

>> y=R(1,3,-10,10);

>> n=-10:10;

>> stem(n,y)

>>

图像如下：

d)

利用.m文件

命令如下：

n=-10:10;

y=2*sin(0.3*n.*n)+0.5*cos(0.6*n.*n);

stem(n,y)

函数图像如下：

数字信号处理Matlab实验数字信号的产生和基本运算

2 以上数列相加

命令如下：

>> n=-10:10;

y1=a(0,-10,10);

y2=2*u(5,-10,10);

y3=R(1,3,-10,10);

y4=2*sin(0.3*pi*n)+0.5*cos(0.6*pi*n);

>> subplot(3,2,1);stem(n,y1+y2);

>> subplot(3,2,2);stem(n,y1+y3);

>> subplot(3,2,3);stem(n,y1+y4);

>> subplot(3,2,4);stem(n,y2+y3);

>> subplot(3,2,5);stem(n,y2+y4);

>> subplot(3,2,6);stem(n,y3+y4);

图像：

数字信号处理Matlab实验数字信号的产生和基本运算

音频处理：

[record,fs,bits]=wavread('record.wav');

[background,fs,bits]=wavread('background.wav');

>> x=record((fs*0+1:fs*45),:);

>> y=background((fs*0+1:fs*45),:);

>> z=0.7*x+0.3*y;

>> wavwrite(z,fs,bits,'music.wav');

>> subplot(3,1,1),plot(x),title('record');

>> subplot(3,1,2),plot(y),title('background');

>> subplot(3,1,3),plot(z),title('music');

//此时的录音是用audition录的，也可以用matlab自带的wavrecord函数来录音。

数字信号处理Matlab实验数字信号的产生和基本运算

实验心得：

做音频处理的时候由于操作失误，导致背景音乐声音较小，需要改进。