编写一下matlab为了消除50Hz工频干扰的代码

在MATLAB中消除50Hz工频干扰的基本步骤通常是设计和应用滤波器。这里提供一个简单的例子，使用`equiripple`函数创建一个带阻滤波器来移除50Hz和其附近的干扰。假设我们有一个包含干扰的信号`signal`，频率响应函数(`freqz`)可以帮助我们理解滤波器的效果。下面是代码示例：

% 导入所需的库
clear all; close all; clc;
import mathworks.signalprocessing.*

% 假设你的信号在这里
signal = ...; % 填充实际的信号数据

% 设定滤波器参数
fc = 50; % 要移除的频率 (50Hz)
bw = 5;   % 带宽 (通常选取5Hz左右)
fs = 1000; % 采样率 (假设你的信号是每秒1000样本)

% 创建带阻滤波器
[num, den] = equiripple(butter(2, [fc/bw, fc+bw]/(fs/2), 'bandstop'));

% 应用滤波器到信号
filtered_signal = filtfilt(num, den, signal);

% 检查滤波效果
figure;
plot([signal filtered_signal]);
title('Original Signal vs Filtered Signal');
xlabel('Samples');
ylabel('Amplitude');

% 计算滤波器频率响应
[H, w] = freqz(num, den, fs);
figure;
plot(w, abs(H));
title('Filter Frequency Response');
xlabel('Frequency (Hz)');
ylabel('|H(f)|');

在这个例子中，`butter`函数创建了一个二阶Butterworth滤波器，`equiripple`用于优化滤波器的过渡带平坦度。`filtfilt`函数则进行了无失真线性滤波，`freqz`用于绘制滤波器的幅值频率响应图。