信号与系统课程设计.ppt

《信号与系统课程设计.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、课程设计班级：成员：众所周知，语音信号的处理分析已变得非常流行，基于语音处理分析技术的产品也开始流入市场，遍布人们的生活。这也是本小组致力于该方面研究的原因。为了研究不同类型的声音信号性质，我们以数学知识为基础，通过快速傅里叶变换及其逆变换等一系列技术手段，从时域图，频域图入手，经对语音信号素材的处理，分析，对比，类比，对各类语音信号性质有了一定的了解，得出了不少有意义的结论。1、对语音信号进行频域分析，找到语音信号的主要频谱成分所在的带宽，验证为何电话可以对语音信号采用8KHz的采样速率。分析：由于要验证电话可以对语音信号用8KHz的采样速率，根据采样定理，就要验证语音信号主要频带范围内的

2、最高频率不超过4KHz，因此需要画出原语音信号的频谱图，并对其频谱图进行分析。相关图像设计思路：首先用44100KHz的采样速率采集语音信号，由于采样信号周期非常小，所以采样后得到的信号几乎和原信号一致，所以可大致用采样信号大致代替原信号，然后通过对采样信号做傅里叶变换，得到信号的频谱图。通过对频谱图分析可知，该信号主要频带宽度在0~4000KHz左右，因此根据采样定理，电话用8KHz的采样速率去采集语音信号可以从采样信号中恢复出原信号。2、分析男声和女声的差别。我们知道男声和女声在频域上是有些差别的，一般大家都会认为女声有更多高频的成分，验证这种差别。同时，提出一种方法，能够对一段音频信号

3、是男声信号、还是女声信号进行自动的判断。分析：首先，选择了一段男生声音素材和女生声音素材，用matlab分别绘出两段声音信号的频谱图,从图中可以看出，男生的声音频率主要分布在100Hz~800Hz,女生的声音频率主要分布在200Hz~1800Hz,女生的声音高频成分较多。这是因为，声波是由物体振动产生的机械波，男人声带宽而厚，振动频率低；女人声带窄而薄，振动频率高。而我们平时所感受得男生声音低沉，女生声音尖细，则是由于发声时男女声带的振动频率的高低不同，所以男女音调的高低不同。相关图像对于判断男女生声音，搜索资料发现男声声音大至在60—550Hz，女声声音大至在150—1200Hz。观察男女

4、声频谱图可知男女声音信号在幅度值相对较高的频率值上男生的频率较低，女声的频率较高，多次分析验证知此结果较正确，可以利用这一点粗略判断，即实现一般情况下男女声的自动判断。找到幅度最大点的频率值，将界限值设定在150Hz，编写程序实现自动判断。程序代码：y1=wavread('C:Q2.wav');Fs=44100;%采样频率%yt1=fft(y1);%傅里叶变换%L=length(yt1);df=Fs/L;Fx=df*(0:L-1);%将横轴变为频率轴%figure(1)plot(Fx,abs(yt1));axis([01000005000]);xlabel('频率/Hz');[a,b]=m

5、ax(abs(yt1))f=b/L*Fsiff<150title('判断结果：男生声音')elsetitle('判断结果：女生声音')end经验证，此方法正确率较高，实现了男女声音的自动正确判断。3、采样后的序列，验证是否会对信号的质量产生影响。降采样的方法很简单，例如命令y=wavread(‘SpecialEnglish.wav’);将语音文件读入后保存在向量y中，这时对应的采样频率为44.1KHz。使用y1=y(1:2:length(y))命令，就可以将原序列y每隔1个采样后放入序列y1中，这时y1序列对应的采样频率即为22KHz。.wav文件的采样速率为44.1KHz，仍然远远高于我们

6、通常说的语音信号需要的频谱宽度，例如在电话对语音信号的采样中，我们仅仅使用8KHz的采样速率。对读入的音频数据进行不同速率的降采样，使用wavplay()命令播放。x=wavread(‘C:Q2.wav');%将声音信号采样称字符串%fs=44100x1=x(1:2:length(x));%对原声音信号进行1/2降采样%x2=x(1:5:length(x));%对原声音信号进行1/5降采样%x3=x(1:10:length(x));%对原声音信号进行1/10降采样%subplot(411);wavplay(x,fs);%播放原声音信号%y=fft(x);df=fs/length(y);fx

7、=df*(0:length(y)-1);plot(fx,abs(y));%绘制原声音信号的频谱分析图%axis([0800002000]);title('原声音频谱分析图');subplot(412);wavplay(x1,fs/2);%播放1/2降采样声音信号%y1=fft(x1);df=fs/length(y1);fx=df*(0:length(y1)-1);plot(fx,abs(y1));%绘制1/2