Matlab中Filter Builder的使用(以滤除音乐中的人声为例)

Matlab中Filter Builder的使用(以滤除音乐中的人声为例)

以matlab2017a为例
先设计一个带通滤波器,修改如下参数为

Input sample rate 为信号采样率,可以通过右上角按钮看滤波器通带
确定之后就会生成一个Hbs滤波函数

这样既可以调用该函数
代码如下
[y,Fs] = audioread('01_1.wav');    %[y,Fs] = audioread(filename) y为音频数据,Fs为采样率
t=0:1/Fs:(size(y)-1)/Fs;           
f=Fs*(0:4095)/4096;
y1=fft(y,4096);
fz=abs(y1(1:4096));     %取y1前2048个构成一维数组
figure(1);
plot(f,fz);
title('滤波前的频谱')
xlabel('Hz');
ylabel('fuzhi');
axis([0,5000,0,90]);
pks=findpeaks(fz);       
[b,c]=sort(pks,'descend');%对fz值进行排列 
wz1=find(fz==b(1));
wz2=find(fz==b(2));
wz3=find(fz==b(3));
F1=wz1/4096*Fs; %找出幅度最大时对应的频率的值,滤除该值
F2=wz2/4096*Fs;
F3=wz3/4096*Fs;
fil_1=filter(Hbs,y);%Hbs滤波函数 y为滤波输入
y2=fft(fil_1,4096);
fz_2=abs(y2(1:4096));
figure(2);
plot(f,fz_2);
axis([0,5000,0,90]);
audiowrite('300_1k.wav',fil_1,Fs)
如果对你有用就表扬一下我哦~ 害羞


评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

文锦渡

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值