再给一个梳状滤波器matlab代码

以下是一个简单的梳状滤波器MATLAB代码示例：

% 设置采样频率 fc 和带宽 bw
fs = 1000;
fc = 50;
bw = 10;

% 计算滤波器参数
delta_f = fc/bw;
Q = fc/delta_f;
beta = 0.5*((1-Q^2)^0.5)/(Q);

% 生成输入信号
t = 0:1/fs:1;
x = sin(2*pi*100*t) + 0.5*sin(2*pi*200*t);

% 进行梳状滤波
y = zeros(size(x));
for i = 1:length(x)
    if i <= Q+1
        y(i) = x(i);
    else
        y(i) = x(i) - beta*y(i-Q);
    end
end

% 绘制原始信号和滤波后的信号
subplot(2,1,1);
plot(t, x);
title('原始信号');
xlabel('时间 (秒)');
ylabel('幅度');

subplot(2,1,2);
plot(t, y);
title('滤波后的信号');
xlabel('时间 (秒)');
ylabel('幅度');

该代码生成了一个包含两个子图的图像，第一个子图显示原始信号，第二个子图显示经过梳状滤波器处理后的信号。