butfilter.m

随机减量方法

function y=butfilter(x,wp,ws)%实现带通滤波，输入为需滤波的信号x，频率带宽[wp ws],输出为带通滤波后的信号。
wp=wp/(20/2);ws=ws/(20/2);
[b,a]=butter(4,[wp ws],'bandpass');
y=filter(b,a,x);