用matlab设计一个梳状滤波器，滤波器能够让 1KHz、2KHz、3KHz、4KHz、 5KHz、6KHz、7KHz 和 8KHz 的信号通过，其他信号一律滤除，可以使用一个 1KHz 的方波信号作为激励，利用梳状滤波器滤出各次谐波分量。

在MATLAB中设计一个梳状滤波器，通常用于音频处理和信号分析中，它允许特定频率范围内的信号通过而抑制其他频率。对于这个特定的需求，我们可以使用数字IIR（无限 impulse response）滤波器，例如巴特沃斯滤波器或者椭圆滤波器，来构建一个带通滤波器，中心频率设置为1KHz到8KHz，并且有良好的阻带衰减特性。

以下是基本步骤：

1. **确定滤波器参数**：
- 中心频率：1KHz到8KHz，选择适当的数量（比如8个）对应于你需要滤过的各次谐波。
- 带宽：为了包含所有的频率，需要选择合适的衰减点（如-3dB截止点），并计算间隔（比如每增加1KHz）。
- 滤波器类型：一般建议使用巴特沃斯滤波器，因为它对过渡带的要求不高，适合带通应用。

2. **设计滤波器**：

   fs = 96000; % 假设采样率 (Hz)
   f_cents = [1:8]*1000; % 中心频率列表 (Hz)
   bw = f_cents(2) - f_cents(1); % 宽度

   % 设计巴特沃斯滤波器
   b, a = butter(length(f_cents), 2*bw/(fs/2), 'bandpass');

3. **创建梳状滤波器**：
将多个滤波器链起来，每个滤波器只针对其中一个中心频率：

   bandpassFilters = cellfun(@(fc) fir1(100, fc*(fs/2)/bw, 'bandpass'), f_cents);
   
   % 将单个滤波器应用于激励信号
   inputSignal = rectpuls(fs, 1, 1000); % 1kHz方波信号
   filteredSignals = cellfun(@conv2, bandpassFilters, {inputSignal}, 'UniformOutput', false);

4. **结果分析**：
现在`filteredSignals`包含了每个频率的单独信号。你可以查看它们以确认是否达到了预期的效果。