华南理工大学数字信号处理实验报告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这个仿真工具，而且还了解了基音周期提取的自相关函数算法，经过同学间的探讨和老师的教导，理解了自相关函数算法的思想，可谓受益匪浅。