写一个陷波滤波器程序代码

以下是一个基于MATLAB的陷波滤波器程序代码示例：

function y = notch_filter(x, f0, fs, bw)
% x: 输入信号
% f0: 陷波中心频率
% fs: 采样频率
% bw: 陷波带宽

w0 = 2*pi*f0/fs; % 陷波中心频率的角频率
Q = f0/bw; % 陷波品质因数

% 计算陷波滤波器系数
b0 = 1;
b1 = -2*cos(w0);
b2 = 1;
a0 = 1 + 1/Q*tan(w0/2);
a1 = -2*cos(w0);
a2 = 1 - 1/Q*tan(w0/2);

% 应用陷波滤波器
y = filter([b0, b1, b2], [a0, a1, a2], x);

end

该程序使用了二阶IIR滤波器的实现方式，计算了陷波滤波器的系数，然后使用MATLAB中的`filter`函数进行滤波操作。你可以根据自己的需要进行参数的调整和修改。