编写一段matlab的代码，找出功率高于预选功率阈值的所有局部最大值，计算这些峰值之间的成对频率距离，并仅保留成对距离在有效人类心率范围内的对，找出成对距离大致相等并且形成连续路径的峰值

以下是一个可能的实现，其中假设有一个名为power的向量包含功率值，预选功率阈值为threshold，有效人类心率范围为[40, 200] bpm：

% 找出功率高于预选功率阈值的局部最大值
[pks, locs] = findpeaks(power, 'MinPeakHeight', threshold);

% 计算成对频率距离
fs = 1000; % 采样频率为1000 Hz
freqs = locs / length(power) * fs; % 将位置转化为频率
dists = diff(freqs);

% 仅保留成对距离在有效人类心率范围内的对
valid_pairs = dists > 60 / 200 * fs & dists < 60 / 40 * fs;

% 找出成对距离大致相等并且形成连续路径的峰值
pairs = [locs(1:end-1), locs(2:end)];
valid_pairs_locs = pairs(valid_pairs,:);
diffs = diff(valid_pairs_locs,[],2);
diffs_thresh = mean(diffs) * 0.5; % 假设相邻峰值间隔的阈值为平均间隔的一半
valid_pairs_mask = diffs > diffs_thresh;
valid_pairs_mask = [false; valid_pairs_mask];
valid_pairs_locs = valid_pairs_locs(valid_pairs_mask,:);
valid_pairs_locs = unique(valid_pairs_locs(:)); % 去除重复的峰值位置

% valid_pairs_locs 就是符合条件的峰值位置