上述效果用matlab实现并给出代码

以下是使用MATLAB实现notch滤波器来消除心电信号中的50Hz工频干扰的代码：

% 导入心电信号数据
ecg_data = load('ecg_data.mat');
ecg_signal = ecg_data.ecg_signal;

% 设置notch滤波器参数
Fs = 1000; % 采样率
Fc = 50; % 中心频率
BW = 10; % 带宽
N = 4; % 阶数

% 设计notch滤波器传递函数
[b, a] = iirnotch(Fc/(Fs/2), BW/(Fs/2), N);

% 应用notch滤波器
ecg_filtered = filtfilt(b, a, ecg_signal);

% 绘制滤波前后的心电信号
t = 1:length(ecg_signal);
subplot(2,1,1)
plot(t, ecg_signal)
title('Original ECG Signal')
xlabel('Time (samples)')
ylabel('Amplitude')
subplot(2,1,2)
plot(t, ecg_filtered)
title('Filtered ECG Signal')
xlabel('Time (samples)')
ylabel('Amplitude')

上述代码中，我们使用了IIR数字滤波器（`iirnotch`函数）来设计notch滤波器的传递函数，并使用`filtfilt`函数来应用滤波器。最后，我们绘制了滤波前后的心电信号，以便比较滤波效果。