matlab代码写50Hz带阻滤波

在MATLAB中，编写一个50Hz带阻滤波器通常需要使用数字信号处理工具箱（Digital Signal Processing Toolbox），特别是`fir1`函数用于设计线性相位有限 impulse response (FIR) 滤波器。50Hz带阻滤波器通常是为了抑制某些频段内的干扰，比如在电力系统分析中，可能希望移除工频噪声。

下面是一个简单的例子，假设我们想要设计一个低通主响应，然后在其截止频率附近添加两个带阻区域，包括50Hz。这里我们选择一个高斯窗函数作为滤波器的设计方法：

% 设定滤波器参数
fs = 1000;       % 采样率，单位Hz
f_c = 50;        % 中心频率，这里是50Hz
bw = [10 10];     % 带宽，上下各10Hz
order = 40;      % 滤波器阶数

% 计算截止频率和带阻区范围
cutoff_freq = f_c / 2;
low_edge = cutoff_freq - bw(1)/2;
high_edge = cutoff_freq + bw(2)/2;

% 设计低通滤波器
lpf_design = fir1(order, cutoff_freq / fs);
% 添加带阻特性
[b,a] = butter(order, [low_edge high_edge]/fs, 'bandstop');
h_bp = impz(b,a);

% 结合低通滤波器和带阻滤波器
h_total = lpf_design .* h_bp;

% 显示滤波器频率响应
freqz(h_total, fs)
xlabel('Frequency (Hz)')
ylabel('Magnitude')
title('50Hz Band-stop Filter')

% 这里省略了实际应用中的滤波操作部分