你目前给出的技术指标不完全,还需要采样频率、通带波纹、阻带衰减。给全了倒是可以帮你编一个。
那就假设几个吧:
采样频率:Fs=80KHz
经计算得相应数字滤波器的指标:
通带截止频率:fp=0.6*Fs/2
通带波纹:rp=0.01(假设的,单位dB,下同);
阻带截止频率:fs=0.7*Fs/2
阻带衰减:rs=0.1(假设的,同上)
程序如下,给出了幅频特性曲线。
clear;
Fs=80;
fp=0.6*Fs/2;
rp=0.01;
fs=0.7*Fs/2;
rs=0.1;
f= [fp fs];
A=[1 0];
dev=[rp rs];
[n,f0,m0,W]=remezord(f,A,dev,Fs);
b=remez(n,f0,m0,W);
[h,W]=freqz(b,1,256,1);
h=abs(h);
h=20*log10(h);
plot(w,h);grid;
xlabel('频率(归一化)');
ylabel('幅度(dB)');