语音滤波 matlab
需积分: 0 161 浏览量
2009-06-06
22:56:04
上传
评论 2
收藏 25KB DOC 举报
[x,fs,bite]=wavread('c:lsndmgr.wav',[1000 20000]);
sound(x,fs,bite);
n=1024;
subplot(2,1,1);
y=plot(x(50:n/4));
grid on ;
title('时域信号')
X=fft(x,256);
subplot(2,1,2);
plot(abs(fft(X)));
grid on ;
title('频域信号');
%采样频率10Hz，通带截止频率fp=3Hz，阻带截止频率fs=4Hz
%通带衰减小于1dB，阻带衰减大于20dB
%使用双线性变换法由模拟滤波器原型设计数字滤波器
T=0.1; FS=1/T;
fp=3;fs=4;
wp=fp/FS*2*pi;
ws=fs/FS*2*pi;
Rp = 1; % 通带衰减
As = 20; % 阻带衰减
% 频率预畸
OmegaP = (2/T)*tan(wp/2); % Prewarp Prototype Passband freq
OmegaS = (2/T)*tan(ws/2); % Prewarp Prototype Stopband freq
%设计butterworth低通滤波器原型
N =
ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS
)));
OmegaC = OmegaP/((10^(Rp/10)-1)^(1/(2*N)));
[z,p,k] = buttap(N); %获取零极点参数
p = p*OmegaC;
k = k*OmegaC^N;
B = real(poly(z));
b0 = k;
cs = k*B;
ds = real(poly(p));
% 双线性变换
czy0001
- 粉丝: 0
- 资源: 1
最新资源
- mysql语句大全及用法
- Online Maps v3 3.7.15版本
- OpenGL渲染引擎-设计与实践
- meta-llama-3-8b-instruct 的 model-00001-of-00004.safetensors 的1/3
- C# yolov8 TensorRT +ByteTrack Demo.rar
- WX小程序源码无后台同城拼车(完整带PHP后台)
- 基于Swift的空气质量监测应用设计源码 - BadAir
- java开发的基于kafka、xlog的web日志实时分析storm topology
- k8s1.20.4-高可用集群部署-新增项目-kubernetes安装包和详细文档笔记整理
- 基于UC3842-S开关电源PDF原理图+反激开关电源全套学习教程资料.zip
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
评论0