lvboqi.m

用MATLAB实现了该算法

%用fdatool设计滤波器，生成一个.fda文件，不知道怎么导入数据，查找资料说的是用simulink的滤波器模块可以导入.fda生成的滤波器
%所以只用了程序来实现  
c=48;                        %定义滤波器阶数
fc1=8400;fc2=13200;fs=48000;
w1=2*pi*fc1/fs;
w2=2*pi*fc2/fs;               %参数转换，将模拟滤波器的技术指标转换为数字滤波器的技术指标
window=hamming(c+1);        %使用hamming窗函数  根据不同的衰减要求选择窗函数
h=fir1(c,[w1/pi w2/pi],window);  %使用标准响应的加窗设计函数fir1
freqz(h,1,512);                %数字滤波器频率响应
y=wavread('yuanshi.wav');
ynew=filter(h,1,y);        %滤波处理

t=(0:length(ynew)-1)/fs;figure;   
subplot(2,1,1);plot(t,y);xlabel('time(s)');

subplot(2,1,2);plot(t,ynew);xlabel('Time(s)');