matlab设计一个fir带阻滤波器

下面是MATLAB代码，设计一个FIR带阻滤波器，具有以下规格：

- 带阻频带范围：2000 Hz - 3000 Hz
- 阻带频带范围：1500 Hz - 3500 Hz
- 采样率：8000 Hz
- 阻带衰减：60 dB
- 通带最大纹波：0.5 dB

fs = 8000; % 采样率
fpass = [2000 3000]/(fs/2); % 通带频率范围
fstop = [1500 3500]/(fs/2); % 阻带频率范围
apass = 0.5; % 通带最大纹波
astop = 60; % 阻带衰减

% 使用firpm函数设计FIR带阻滤波器
n = ceil(6.6/(fstop(2)-fstop(1))); % 滤波器阶数
b = firpm(n, [0 fstop(1) fstop(2) 1], [1 0 1], [apass astop]); % 设计滤波器

% 绘制滤波器的幅频响应曲线
freqz(b, 1, 1024, fs);

% 将滤波器应用于信号
% x为输入信号，y为输出信号
y = filter(b, 1, x);

请注意，上述代码中使用了firpm函数来设计FIR滤波器。该函数可以根据指定的通带和阻带频率范围、通带最大纹波和阻带衰减等参数，自动计算滤波器的阶数和系数。设计出的滤波器通常比手动设计的滤波器更优秀。