详细写一下MATLAB对下载了MIT-BIH心电信号.mat文件预处理的代码

MIT-BIH心电信号是一个经典的心电信号数据集，包含48个记录的心电信号文件和相应的注释文件。如果已经下载了MIT-BIH心电信号.mat文件，可以使用MATLAB对其进行预处理。

1. 加载数据

load('mitdb_ecg_signals.mat'); % 加载心电信号数据
load('mitdb_ecg_annotations.mat'); % 加载心电信号注释数据

2. 选择一个记录进行预处理

record = 1; % 选择第1个记录
ecg_signal = ecg_signals{record};
annotations = ecg_annotations{record};

3. 滤波

fs = 360; % 采样率
f1 = 0.5; % 截止频率1
f2 = 30; % 截止频率2
Wn = [f1 f2]/(fs/2);
[b,a] = butter(1,Wn); % 一阶巴特沃斯滤波器
filtered_signal = filtfilt(b,a,ecg_signal); % 无相位延迟滤波

4. 去基线漂移

baseline = medfilt1(filtered_signal,round(fs*0.2)); % 中值滤波
baseline_corrected_signal = filtered_signal - baseline; % 去基线漂移

5. 去除噪声

noise_std = std(baseline_corrected_signal); % 计算信号噪声标准差
threshold = 3*noise_std; % 设置阈值
noisy_samples = abs(baseline_corrected_signal) > threshold; % 标记噪声样本
baseline_corrected_signal(noisy_samples) = NaN; % 去除噪声样本

6. 绘制心电信号和注释

t = (1:length(ecg_signal))/fs; % 时间轴
figure;
plot(t,ecg_signal,'k');
hold on;
plot(t,annotations*max(ecg_signal),'r');
xlabel('Time (s)');
ylabel('Amplitude (mV)');
title('Original ECG Signal');
legend('ECG Signal','Annotations');

figure;
plot(t,baseline_corrected_signal,'k');
hold on;
plot(t,annotations*max(baseline_corrected_signal),'r');
xlabel('Time (s)');
ylabel('Amplitude (mV)');
title('Preprocessed ECG Signal');
legend('ECG Signal','Annotations');

以上代码示例仅供参考，具体的预处理过程还需根据实际情况进行调整和优化。