豆搜网 文档下载 文档下载导航
设为首页 | 加入收藏
搜索 请输入内容:  
 导航当前位置: 文档下载 > 所有分类 > 工程科技 > 信息与通信 > 基于MATLAB的语音信号录制采集和窗函数法滤波器设计
侵权投诉

基于MATLAB的语音信号录制采集和窗函数法滤波器设计

[y,fs,bits]=wavread('G:\你好.wav');
sound(y,fs) % 回放语音信号
n=length(y) %选取变换的点数
y_p=fft(y,n); %对n点进行傅里叶变换到频域
f=fs*(0:n/2-1)/n; % 对应点的频率
figure(1)
subplot(2,1,1);
plot(y); %语音信号的时域波形图
title('原始语音信号采样后时域波形');
xlabel('时间轴')
ylabel('幅值 A')
subplot(2,1,2);
plot(f,abs(y_p(1:n/2))); %语音信号的频谱图
axis([0 22500 0 0.6]);
title('原始语音信号采样后频谱图');
xlabel('频率Hz');
ylabel('频率幅值');
wp=2*pi*1000/fs;
ws=2*pi*1200/fs;
Rp=1;
Rs=53;
wdelta=ws-wp;
N=ceil(8*pi/wdelta);%取整
wn=(wp+ws)/2;
[b,a]=fir1(N,wn/pi,hamming(N+1));%选择窗函数,并归一化截止频率
figure(2)
freqz(b,a,512);title('FIR低通滤波器');
f2=filter(b,a,y);
figure(3)
subplot(2,1,1);plot(y);title('FIR低通滤波器滤波前的时域波形');
subplot(2,1,2);plot(f2);title('FIR低通滤波器滤波后的时域波形');
sound(f2,44100);%播放滤波后的语音信号
F0=fft(f2,1024);
f=fs*(0:511)/1024;
figure(4)
y2=fft(y,1024);
subplot(2,1,1);
plot(f,abs(y2(1:512)));
title('FIR低通滤波器滤波前的频谱')
xlabel('频率/Hz');
ylabel('幅值');
subplot(2,1,2)
F2=plot(f,abs(F0(1:512)));
title('FIR低通滤波器滤波后的频谱')
xlabel('频率/Hz');
ylabel('幅值');

第1页

热门文档

相关文档

站点地图 | 文档上传 | 侵权投诉 | 手机版
新浪认证  诚信网站  绿色网站  可信网站   非经营性网站备案
本站所有资源均来自互联网,本站只负责收集和整理,均不承担任何法律责任,如有侵权等其它行为请联系我们.
文档下载 Copyright 2013 doc.docsou.com All Rights Reserved.  闽ICP备15022310号-9  闽公网安备 35021102001881号  email
返回顶部