matlab 2.封装HPF 巴特沃斯高通滤波器

%x：待处理信号
%fs：采样频率
%wp;通带边界频率
%ws:阻带边界频率
%rp:Rp通带最大衰减(dB)
%rs:Rs阻带最小衰减(dB)
function y=HPF(x,fs,wp,ws,rp,rs)
[n,w] = buttord(wp/(fs/2),ws/(fs/2),rp,rs);%确定阶数n和3dB截止频率
[b,a] = butter(n,w,'high');%确定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