MATLAB实现心电信号R波检测

要在MATLAB中实现心电信号R波检测，可以使用基于波峰检测的方法。以下是一个简单的示例代码：

% 读取心电信号数据
ecgData = load('ecgData.mat');

% 滤波
fs = 1000;  % 采样率为1000Hz
fc = 50;  % 截止频率为50Hz
[b, a] = butter(2, fc / (fs/2), 'high');  % 高通滤波器
ecgFilt = filtfilt(b, a, ecgData.ecg);

% 求导
ecgDiff = diff(ecgFilt);

% 平方
ecgSqr = ecgDiff .^ 2;

% 移动平均
winSize = round(0.15 * fs);  % 窗口大小为150ms
ecgMovAvg = movmean(ecgSqr, winSize);

% 波峰检测
[pks, locs] = findpeaks(ecgMovAvg, 'MinPeakHeight', mean(ecgMovAvg), 'MinPeakDistance', round(0.5 * fs));

% 绘制心电信号和R波位置
figure;
plot(ecgData.ecg);
hold on;
scatter(locs, pks);
xlabel('Sample Number');
ylabel('Amplitude');
title('R Peak Detection');

该代码通过高通滤波、求导、平方、移动平均和波峰检测等步骤，实现了心电信号R波检测。当然，这只是一个简单的示例，实际应用中可能需要更复杂的方法来提高检测准确率。