语音信号的低通滤波 和短时分析综合实验.pdfVIP

语音信号处理

实验报告一

实验名称语音信号的低通滤波学院通信工程学院

和短时分析综合实验

教师姓名实验地点6505实验日期2015-03-28

学生姓名学号班级

一、实验内容

1、根据已有语音信号，设计一个低通滤波器，带宽为采样频率的四分之一，求输出信号；

2、辨别原始语音信号与滤波器输出信号有何区别，说明原因；

3、改变滤波器带宽，重复滤波实验，辨别语音信号的变化，说明原因；

4、利用矩形窗和汉明窗对语音信号进行短时傅立叶分析，绘制语谱图并估计基音周期，分析两种

窗函数对基音估计的影响；

5、改变窗口长度，重复上一步，说明窗口长度对基音估计的影响。

二、实验目的

1、读取’.wav’音频文件数据

2、巴特沃斯低通滤波器

3、给语音信号加窗

4、做语音信号短时傅里叶变换

5、做语音信号语谱图

三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况）

计算机一台（装有Matlab）

四、实验试做记录（含程序、数据记录及分析）

结果：

读取’.wav’音频文件数据

[s,Fs,bit]=wavread(‘‘语音信号处理语音信号处理’);%wavread函数读取声音

%soundview(s,Fs);

%sound(s,Fs);

y=s(:,1);

x=(1:length(y))/Fs;%语音信号时间

figure(1);subplot(311);plot(x,y,b)%滤波前信号图形

xlabel(时间);ylabel(幅度);title(原始信号（我到北京去）);

f_N=Fs;%采样率

f_p=f_N/4;f_s=0.4*f_N;R_p=6;R_s=60;%设计要求指标(f_p通带频率、f_s阻带频率、R_p

通带衰减、R_s阻带衰减）

Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);%计算归一化角频率

[n,Wn]=buttord(Wp,Ws,R_p,R_s);%计算阶数和截止频率

[b,a]=butter(n,Wn);%计算H(z)

%freqz(b,a,1000,8000)%作出H(z)的幅频相频图,freqz(b,a,计算点数,采

样率)

%subplot(2,1,1);axis([04000-303]);

ym=filter(b,a,y);%低通滤波器滤波

subplot(312);plot(x,ym,r);%滤波后信号图形

xlabel(时间);ylabel(幅度);title(滤波信号（我到北京去）);

subplot(313);plot(x,y,b);holdon;%将两个图形做在一个坐标中

plot(x,ym,r);legend(原始声音,滤波声音);

xlabel(时间);ylabel(幅度);title(原始与滤波后图形对比);

N=256;

I=s(6001:6000+N,1);%取出其中N个点

w1=window(@rectwin,N);%N个点的矩形窗

w2=window(@hamming,N);%N个点的汉明窗

I1=I.*w1;%对N点信号加矩形窗

I2=I.*w2;%对N