电子科大2012信号与系统语音信号分析课程设计报告.pdfVIP

电子科大2012信号与系统语音

信号分析课程设计报告

LT

绘图命令：plot、stem等

求频谱、频率响应：fft、freqz等

卷积、滤波：conv、deconv、filter等

声音文件读写：wavread、wavwrite、wavplay等

六、设计内容

1、对语音信号进行频域分析，找到语音信号的主要频谱成分所在的带宽，说明

为何电话可以对语音信号采用8KHz的采样速率。

程序代码:

[x,fsx,bits]=wavread(FriendsCut.wav);

subplot(211);

plot(x);

title(时间域特性);

subplot(212);

y=fft(x,fsx);

df=fsx/length(y);

fx=df*(0:length(y)-1);

plot(fx,abs(y));

axis([080000500]);

title(频率域特性);

说明：

找到一个音频文件，通过截取其中某一段并进行格式转换，得到wav格式

的音频文件（FriendsCut.wav）。对声音信号的频谱图进行分析，使用Matlab

绘制该语音信号的频谱图，观察频谱图，读出声音信号的频率范围，由采样定

理可知，如果需要重建声音信号，需产生一个周期冲激串，其冲激幅度就是采

样得到的样本值，给冲激串加上一个增益即可恢复原声音信号，此为采样定理。

使用Matlab中的快速傅里叶变换(fft)，绘制出声音文件的时域波形和频

域波形并进行分析，观察可得声音信号的主要频率范围为0~2500Hz，而根据采

样定理可知，该声音信号采样频率应不小于5000Hz，故电话使用8000Hz的采样

频率可以保证声音在失真较少的条件下进行信号的采样与恢复。

输出结果：

2、分析男声和女声的差别。我们知道男声和女声在频域上是有些差别的，一般

大家都会认为女声有更多高频的成分，验证这种差别。同时，提出一种方法，

能够对一段音频信号是男声信号、还是女声信号进行自动的判断。

程序代码：

y1=wavread(Man.wav);

Fs=44100;

yt1=fft(y1);

df=Fs/length(yt1);

Fx=df*(0:length(yt1)-1);

figure(1)

subplot(221);

plot(y1);title(男声时域特性);

subplot(222);

plot(Fx,abs(yt1));

axis([010000010000]);

title(男声频域特性);

xlabel(频率/Hz);

y2=wavread(Wom.wav);

Fs=44100;

yt2=fft(y2);

df=Fs/length(yt2);

Fx=df*(0:length(yt2)-1);

figure(1)

subplot(223);

plot(y2);title(女声时域特性);

subplot(224);

plot(Fx,abs(yt2));

axis([01000005000]);

title(女声频域特性);

xlabel(频率/Hz);

s1=0;

forj1=2700:4000;

s1=s1+abs(y1(j1));

end

s0=0;

fori1=200:4000;

s0=s0+abs(y1(i1));

end

p1=s1/s0

s2=0;

forj=2700:4000;

s2=s2+abs(y2(j));

end

s=0;

fori=200:4000;

s=s+abs(y2(i));

end

p2=s2/s

说明：

首先，通过截取一段男女对话的音频并进行格式转化，得到wav格式的男声

语音文件(Man.wav)和女声语音文件（Wom.wav）。通过wavread命令和快速傅里

叶变换fft命令对语音信号进行相应的处理，并绘制出相应的时间域与频率域特

性图。通过观察得到的图形，我们可以看出男声与女声的声音大部分集中于

200-4000Hz之间，而男声较女声更多的声音集中在较为低的频率段。所以我们

通过计算男声和女声在高频部分（我们选择了2700-4000这一频段）的声音比例，

来较为粗略地判断男女声的差异，并通过这样的思路来形成一个简单的男女声

判断机制。

输出结果：

男声高频部分所占比例：

p1=0.2432