MATLAB含噪信号超阈值峰值统计：需消除相邻峰值干扰

我来帮你搞定这个噪声干扰导致峰值统计过多的问题～你的原代码思路是检测“从低于阈值跳到高于阈值”的点，但强噪声会让信号在阈值附近反复波动，导致连续多个相邻点都被误判为独立峰值，这就是问题所在。下面给你两种实用的解决方案：

方法1：用MATLAB自带的findpeaks函数（最简便） #

findpeaks是MATLAB专门为峰值检测设计的函数，内置了很多过滤噪声的参数，非常适合你的场景。核心是通过设置最小峰高（确保超过阈值）和最小峰距（过滤相邻的噪声小峰值）来精准统计有效峰值：

thresh = 3; % 你的阈值
min_peak_distance = 5; % 可根据信号采样率/噪声强度调整，单位是采样点数
x = VarName1;

% 调用findpeaks，设置关键参数
[peaks_values, peak_locations] = findpeaks(x, ...
    'MinPeakHeight', thresh, ... % 只保留高于阈值的峰值
    'MinPeakDistance', min_peak_distance); % 两个峰值至少间隔指定采样点数

% 输出结果
fprintf('有效峰值数量：%d\n', length(peaks_values));
disp('峰值位置索引：');
disp(peak_locations);

你可以根据信号的实际情况调整min_peak_distance——如果噪声是高频密集的，就把这个值设大一点，比如10或20，直到过滤掉那些干扰性的相邻小峰值。

方法2：手动处理连续超阈值区间（自定义性更强） #

如果你不想依赖内置函数，可以手动把连续超过阈值的信号片段归为一个“峰值块”，每个块只保留一个真正的峰值（比如块内的最大值位置）：

thresh = 3;
x = VarName1;

% 第一步：标记所有超过阈值的点
above_thresh = x >= thresh;

% 第二步：找到所有连续超阈值区间的起始和结束索引
diff_mask = diff([0; above_thresh; 0]); % 用差分找区间边界
start_indices = find(diff_mask == 1);
end_indices = find(diff_mask == -1) - 1;

% 第三步：遍历每个区间，提取该区间内的峰值位置（最大值所在点）
valid_peak_locs = [];
for i = 1:length(start_indices)
    current_segment = x(start_indices(i):end_indices(i));
    [~, max_pos_in_segment] = max(current_segment); % 找到区间内最大值的位置
    valid_peak_locs = [valid_peak_locs, start_indices(i) + max_pos_in_segment - 1];
end

% 输出结果
fprintf('有效峰值数量：%d\n', length(valid_peak_locs));
disp('峰值位置索引：');
disp(valid_peak_locs);

这个方法的好处是你可以自定义峰值的选择逻辑——比如不想用最大值，也可以用区间的中点或者其他规则，灵活性更高。

为什么你的原代码会统计过多峰值？ #

你的原代码idx(yest)会把每一次“从低于阈值到高于阈值”的跳变都算作一个峰值，但强噪声会让信号在阈值附近反复横跳（比如[2,3,4,3,2,3,4]这种序列），原代码会把第2、6个点都标记为峰值，但这其实是同一个噪声波动带来的连续超阈值段，我们需要的是每个这样的段只保留一个有效峰值。

内容的提问来源于stack exchange，提问作者Simon