数字语音处理上机实验.doc

文档介绍：数字语音处理
上
机
实
验
指导老师:
姓名:
学号:
班序号:
实验一:MATLAB试验
实验目的
了解信号谱分析特点,掌握各种窗函数功率谱的异同及其对短时信号谱分析的作用。
明确短时傅立叶变换与声谱图的关系,理解其对语音信号的时频分析作用。
学习音频文件的基本操作方法。
学习用MATLAB实现基本的语音信号处理,提高自学和动手能力,培养学习兴趣。
实验内容
输入并运行以下MATLAB代码:
程序如下:
%(1)信号波形和对数功率谱示例
% 波形:100样点,抽样频率Fs=8kHz,信号频率
%f=2kHz
x=sin(2*pi*(0:99)*2000/8000);
% 频谱:是波形的快速傅立叶变换
X=fft(x);
% 对数功率谱:功率谱取10lg或幅度谱取20lg
P=20*log10(abs(X));
% 绘制波形和对数功率谱(注意观察25和75处的尖峰)
figure(1);
plot(x);
title('信号的波形');
figure(2);
plot(P);
title('信号的对数功率谱');
%(2)窗函数的对数幅度谱示例
% 避免对0取对数
floor=;
% 图1:矩形窗
W_rec=10*log10(abs(fft(boxcar(512),1024))+floor);
figure(1);
plot(W_rec(1:end/2));
title('矩形窗的对数幅度谱');
% 图2:海明窗
W_hamming=10*log10(abs(fft(hamming(512),1024))+floor);
figure(2);
plot(W_hamming(1:end/2));
title('海明窗的对数幅度谱');
% 图3:汉宁窗
W_hanning=10*log10(abs(fft(hanning(512),1024))+floor);
figure(3);
plot(W_hanning(1:end/2));
title('汉宁窗的对数幅度谱');
(2)窗函数的对数幅度谱示例
% 避免对0取对数
floor=;
% 图1:矩形窗
W_rec=10*log10(abs(fft(boxcar(512),1024))+floor);
figure(1);
plot(W_rec(1:end/2));
title('矩形窗的对数幅度谱');
% 图2:海明窗
W_hamming=10*log10(abs(fft(hamming(512),1024))+floor);
figure(2);
plot(W_hamming(1:end/2));
title('海明窗的对数幅度谱');
% 图3:汉宁窗
W_hanning=10*log10(abs(fft(hanning(512),1024))+floor);
figure(3);
plot(W_hanning(1:end/2));
title('汉宁窗的对数幅度谱');
2、编程绘制以下窗序列的对数幅度谱(选取N=512):
(1)指数窗
W1(n) = a 0nN-1
= 0 otherwise
【编程如