自适应滑动窗口法检测p波 M代码

以下是一个简单的MATLAB代码，演示如何使用自适应滑动窗口法检测P波：

function [p_onset, p_offset] = detect_p_wave(ecg_signal, fs)

% 设置自适应滑动窗口的初始长度和增量
win_len = round(0.1 * fs);
win_inc = round(0.05 * fs);

% 设置P波搜索的起始位置和搜索的范围
search_start = round(0.1 * fs);
search_range = round(0.5 * fs);

% 初始化P波起点和终点
p_onset = [];
p_offset = [];

% 初始化自适应滑动窗口的起始位置
win_start = search_start;

% 遍历ECG信号
while win_start + win_len <= length(ecg_signal)

    % 计算当前窗口内的平均值和标准差
    win_mean = mean(ecg_signal(win_start:win_start+win_len-1));
    win_std = std(ecg_signal(win_start:win_start+win_len-1));

    % 根据平均值和标准差计算阈值
    threshold = win_mean + 0.5 * win_std;

    % 在搜索范围内查找第一个超过阈值的样本点
    [p_peak, p_loc] = max(ecg_signal(win_start+search_start:win_start+search_start+search_range));
    p_loc = p_loc + search_start - 1;

    % 如果找到P波峰值
    if p_peak > threshold

        % 在P波峰值左侧查找第一个低于阈值的样本点，作为P波起点
        for i = p_loc-1:-1:win_start+search_start
            if ecg_signal(i) < threshold
                p_onset = [p_onset i];
                break;
            end
        end

        % 在P波峰值右侧查找第一个低于阈值的样本点，作为P波终点
        for i = p_loc+1:win_start+search_start+search_range
            if ecg_signal(i) < threshold
                p_offset = [p_offset i];
                break;
            end
        end

        % 更新自适应滑动窗口的起始位置
        win_start = win_start + win_inc;

    % 如果没有找到P波峰值
    else

        % 将自适应滑动窗口的长度增加一倍
        win_len = win_len * 2;

        % 如果窗口长度超过搜索范围，则跳出循环
        if win_len > search_range
            break;
        end

    end

end

end

该代码使用自适应滑动窗口法来检测ECG信号中的P波。首先，设置自适应滑动窗口的初始长度和增量，以及P波搜索的起始位置和搜索的范围。然后，遍历ECG信号，并在每个窗口内计算平均值和标准差，根据平均值和标准差计算阈值。在搜索范围内查找第一个超过阈值的样本点，作为P波峰值。在P波峰值左侧和右侧查找第一个低于阈值的样本点，作为P波起点和终点。最后，更新自适应滑动窗口的起始位置，并根据需要增加窗口长度。