matlab 1.封装LPF 巴特沃斯低通滤波器

%x：待处理信号%fs：采样频率%wp;通带边界频率%ws:阻带边界频率%rp:Rp通带最大衰减(dB)%rs:Rs阻带最小衰减(dB)function y=LPF(x,fs,wp,ws,rp,rs)[n,w] = buttord(wp/(fs/2),ws/(fs/2),rp,rs);%确定阶数n和3dB截止频率[b,a] = butter(n,w);%确定H(S)的分子分母系数向量[H,W] = freqz(b,a);%H即频率响应函数，W即数字角频率figure,plot(W*fs/(2*pi),abs(H));%滤波器幅频特性，自变量为模拟频率xlabel('f/Hz'),ylabel('放大倍数（倍）'),title('幅频特性');y=filter(b,a,x);end