广州大学数字信号处理实验一

开课学院及实验室： 年 月 日 学 院 机 电 年级、专业、班 06信工（1） 姓名 业 学号 成绩 指导 教师 0615050081 实验课程名称 实验项目名称 一、 实验目的

数字信号处理 实验一 熟悉MATLAB环境 1．熟悉MATLAB的主要操作命令。 2．学会用MATLAB创建时域离散信号。 3．学会创建MATLAB函数。

二、 实验原理

参阅附录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、F、G。 2．用MATLAB实现以下序列

(1)单位抽样序列

(n-n0)1,nn00,nn0 (2)单位阶跃序列

1,nn0u(n-n0)0,nn0(3)矩形序列

1,(0nN1)RN(n)0,(n0,nN)

(4)正弦序列

X(n)=5sin(0.5πn+ π/4) (5)指数序列

X(n)=exp(-0.5n)

3．用MATLAB生成以下两个序列：

x(n)(n)2(n1)3(n2)4(n3)5(n4) h(n)(n)2(n1)(n2)2(n3)

并作以下运算，并绘制运算后序列的波形。 （1）x(n5),（2）x(n) （3）x(n)h(n) （4）3x(n) （5）x(n)h(n)

x(n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

subplot(3,3,1) stem(A,'.') subplot(3,3,2) stem(B,'.') subplot(3,3,3) stem(C,'.')

subplot(3,3,4) stem(D,'.') subplot(3,3,5) stem(E,'.') subplot(3,3,6) stem(F,'.') subplot(3,3,7) stem(G,'.')

实验内容2 结果与程序如下(已改正):

离散时间信号的产生 (1)单位抽样序列

1,nn0

(n-n0)0,nn0

先定义delta函数,并保存. function[x,n]=delta(n0,n1,n2) n=[n1:n2]; x=[(n-n0)==0];

然后执行下面程序.(以δ(n-3)为例) [x,n]=delta(3,-1,10); stem(n,x);

(2)单位阶跃序列

1,nn0u(n-n0)0,nn0

先定义stepseq函数,并保存.

function[x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)>=0];

然后执行下面程序.(以u(n-3)为例)

[x,n]=stepseq (3,-1,10); stem(n,x);

(3)矩形序列

1,(0nN1)RN(n)0,(n0,nN)

先定义aaa函数,并保存.

function[x,n]=aaa(N,n1,n2) n=[n1:n2];

x=[(N>n)&(n>=0)];

然后执行下面程序.(以R3)为例)

[x,n]=aaa (3,-1,10); stem(n,x);

(4)单位斜坡序列

x(n)n,n0

0,n0

先定义ramp函数,并保存.

function[x,n]=ramp(n1,n2) n=[n1:n2]; x=n;

然后执行下面程序 [x,n]=ramp (0,10); Stem(n,x);

(5)正弦序列

例: x=5*sin(0.5*pi*n+ pi/4) n=-pi:0.1:pi;

x=5*sin(0.5*pi*n+pi/4); stem(n,x)

(6)指数序列

例:x=5*exp(-0.5*n) n=-1:0.1:1;

x=5*exp(-0.5*n); stem(n,x)

(7)任意序列 例:

x=[1,5,-4,2,5,-1,5]; n=1:length(x); stem(n,x)

实验内容3 结果与程序如下:

n=[-10:10]

h=[(n)==0]+2*[(n-1)==0]+[(n-2)==0]+2*[(n-3)==0] subplot(3,3,1);stem(n,h,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('h(n)')

x=[(n)==0]+2*[(n-1)==0]+3*[(n-2)==0]+4*[(n-3)==0]+5*[(n-4)==0] subplot(3,3,2);stem(n,x,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('x(n)') subplot(3,3,3);stem(n+5,x,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('x(n-5)') subplot(3,3,4);stem(n-5,x,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('x(n+5)') subplot(3,3,5);stem(-n,x,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('x(-n)') subplot(3,3,6);stem(n,x+h,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('x(n)+h(n)') subplot(3,3,7);stem(n,3*x,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('3x(n)') subplot(3,3,8);stem(n,x.*h,'.')

axis([-10,10,0,5]);xlabel('n');ylabel('x(n)h(n)')

实验内容3波形与程序如下:

[y,fs,bits]=wavread('C:\\Documents and Settings\\Administrator\\桌面\\数字信号处理课后答案\\爱的代价.wav') t=(0:length(y)-1)/fs subplot(2,1,1); plot(t,y);

axis([0,13,-0.4,0.6]) xlabel('时间(s)'); ylabel('幅度'); subplot(2,1,2) plot(t,y/2)

axis([0,13,-0.4,0.6])

xlabel('时间(s)'); ylabel('幅度')

wavwrite(y/2,fs,bits, 'C:\\Documents and Settings\\Administrator\\桌面\\数字信号处理课后答案\\爱的代价2.wav')