%x:待处理信号
%fs:采样频率
%wp;通带边界频率
%ws:阻带边界频率
%rp:Rp通带最大衰减(dB)
%rs:Rs阻带最小衰减(dB)
function y=BPF(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
07-03
1569
11-30
428
11-30
744