Matlab中findpeaks识别平峰区域起始峰值的问题求助

针对你遇到的平台型峰区（红色块）中定位上升沿结束后首个峰值的需求，以下是几种精准的自定义实现方案：

方案1：基于连续高值区间的定位 #

先识别所有高于阈值的连续峰区，再提取每个区间内上升沿结束的第一个点：

threshold = 4.3; % 匹配你之前用的最小峰值高度
highRegions = data >= threshold;

% 定位连续高值区间的起始、结束索引
regionStarts = strfind([false, highRegions], [false, true]);
regionEnds = strfind([highRegions, false], [true, false]);

targetPeaks = [];
for i = 1:length(regionStarts)
    startIdx = regionStarts(i);
    endIdx = regionEnds(i);
    regionData = data(startIdx:endIdx);
    
    % 计算差分，找到上升趋势结束的位置（差分从正转平/负）
    diffData = diff(regionData);
    % 可根据数据波动调整这个差分阈值
    flatPos = find(diffData <= 0.01, 1, 'first');
    
    if isempty(flatPos)
        % 若区间全程上升，取区间最后一个点作为峰值
        peakInRegion = length(regionData);
    else
        % 上升沿结束后的第一个点，即平台起始点
        peakInRegion = flatPos + 1;
    end
    
    targetPeaks = [targetPeaks, startIdx + peakInRegion - 1];
end
peakHeights = data(targetPeaks);

方案2：直接遍历捕捉上升沿转峰值的节点 #

这种方法更直观，直接遍历数据，捕捉“从低于阈值上升到高值后，第一个停止上升的点”：

threshold = 4.3;
targetPeaks = [];
n = length(data);
isRising = false;

for i = 2:n
    % 触发上升沿检测：当前点达标，前一个点未达标
    if data(i) >= threshold && data(i-1) < threshold
        isRising = true;
    end
    
    % 当处于上升状态，且下一个点不再上升（或已到数据末尾），记录当前点
    if isRising && (i == n || data(i+1) <= data(i))
        targetPeaks = [targetPeaks, i];
        isRising = false; % 重置状态，寻找下一个峰区
    end
end
peakHeights = data(targetPeaks);

方案3：结合局部最大值验证（可选） #

如果需要确保目标点是严格的局部最大值，可以在上述方案基础上加入islocalmax验证：

validPeaks = [];
validHeights = [];
% 窗口大小可根据峰区宽度调整
isLocalMax = islocalmax(data, 'WindowSize', 5, 'FlatSelection', 'first');

for idx = targetPeaks
    if isLocalMax(idx)
        validPeaks = [validPeaks, idx];
        validHeights = [validHeights, data(idx)];
    end
end

关键调整点 #

• threshold：必须匹配你的数据峰值下限，和之前用的minPeakHeight保持一致
• 差分阈值（方案1中的0.01）：根据数据的波动幅度调整，若数据噪声大，可适当增大
• WindowSize（方案3）：根据峰区的宽度设置，确保覆盖峰区的局部范围

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