华南理工大学数字信号处理实验报告1(曹老师)

一、实验目的

1. 加深对常用离散信号的理解。

2. 熟悉使用Matlab 在时域中产生一些基本的离散时间信号。

二、实验内容

1、代码：

N = input('Input N = ');

x = [1 zeros(1,N-1)];

n = 0:N-1;

stem([[-(N-1):-1],n],[zeros(1,N-1),x]);

axis([-(N-1),N-1,0,1]);

运行结果：

2、代码：

N = input('Input N = ');

x = ones(1,N);

n = 0:N-1;

stem([[-(N-1):-1],n],[zeros(1,N-1),x])

axis([-(N-1),N-1,0,1]);

运行结果：

3、代码：

N = input('N = ');

A = input('A = '); f = input('f = ');

Fs = input(

'Fs = '); fai = input('fai = '); n = 0:N-1;

x = A*sin(2*pi*f.*n/Fs+fai); stem(n,x);

运行结果：

0123456789

4、代码：

N = input('N = '); r = input('r = '); w = input('w = '); n = 0:N;

x = r*exp(j*w*n); subplot(1,2,1); stem(n,abs(x)); title('·ù?è'); subplot(1,2,2); stem(n,angle(x)); title('?à??');

运行结果：

5

10

幅

度

5

10

相位

5、代码：

N = input('N = '); a = input('a = '); n = 0:N; x = a.^n; stem(n,x)

运行结果：

012345678910

基音周期的提取：自相关函数法

代码如下：

function pitch

x=wavread('d:\HuangHaoyu\a0.wav');

figure(1);

stem(x,'.');

n=160;

for m=1:length(x)/n;

for k=1:n;

Rm(k)=0;

for i=(k+1):n;

Rm(k)=Rm(k)+x(i+(m-1)*n)*x(i-k+(m-1)*n);

end

end

p=Rm(10:n);

[Rmax,N(m)]=max(p);

end

N=N+10;

T=N/8;

figure(2);stem(T,'.');axis([0 length(T) 0 10]);

xlabel('帧数(n)');ylabel('周期(ms)');title('各帧基音周期'); T1= medfilt1(T,5);

figure(3);stem(T1,'.');axis([0 length(T1) 0 10]); xlabel('帧数(n)');ylabel('周期(ms)');title('各帧基音周期');

运行结果如下：

2

4

6

8

101214

16

18

20

012345678910帧数(n)

周期(m s )

各帧基音周期

三、实验总结

通过这次实验，我学会了如何使用matlab产生各种基本的离散序列，并进一步熟悉了matlab这个仿真工具，而且还了解了基音周期提取的自相关函数算法，经过同学间的探讨和老师的教导，理解了自相关函数算法的思想，可谓受益匪浅。