butfilter.m

来自「随机减量方法」· M 代码 · 共 4 行

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);