基于双差分值和RR间隔处理的心电信号R峰检测算法（MATLAB R2018A）

心电信号中的R峰是确定心率和节律、以及检测其它波形特征点（图1A）的基础。R峰的准确检测是心率变异性分析、心拍分割和心律失常识别重要的处理步骤。

现有的心电信号R峰检测方法主要为基于规则的决策法和基于深度学习的检测方法。基于规则的决策法通常对心电信号进行时域或频域变换，增强R波分量，抑制其他波形和各种噪声分量，在增强R波分量的变换信号中以阈值法进行点检测，再通过启发式的后处理校正R峰位置。然而，规则难以涵盖心电信号的各种复杂情况，如QRS波群具有多种形态（图1B）、可穿戴式心电信号存在复杂干扰（图1C）、以及不同病人的心率差异大等；另外，规则的归纳依赖于丰富的领域知识，难以在各种心律失常事件中精确地检测出R峰，漏检率和误检率高。

鉴于此，提出一种基于双差分值和RR间隔处理的心电信号R峰检测算法，运行环境为MATLAB R2018A。

for i = 1:size(qrs_regions, 1)qrs_start = qrs_regions(i, 1);qrs_end = qrs_regions(i, 2);% Ensure QRS window does not extend beyond the signalqrs_start = max(qrs_start, 1);qrs_end = min(qrs_end, length(y));% Extract the QRS windowqrs_window = y(qrs_start:qrs_end);% Calculate the maximum and minimum amplitudesmax_amplitude = max(qrs_window);min_amplitude = min(qrs_window);% Calculate the mean of maximum and minimum amplitudesmean_amplitude = (max_amplitude + min_amplitude) / 2;% Calculate relative magnitudes by subtracting the mean from all data pointsrelative_magnitudes = qrs_window - mean_amplitude;% Find the position of the maximum relative magnitude within the window[~, max_rel_magnitude_index] = max(relative_magnitudes);% Convert the relative index to the absolute indexr_peak_location = qrs_start - 1 + max_rel_magnitude_index;% Add the R-peak location to the listr_peak_locations = [r_peak_locations, r_peak_location];完整代码：https://mbd.pub/o/bread/mbd-ZZiUl55v
end