如何检测信号中连续上升起始与下降终止？代码识别异常求助

你的代码之所以会把单次活动识别为两个上升沿，核心原因是信号在阈值3000附近出现了震荡波动：当信号短暂低于阈值后又回到阈值以上时，loc = audio_file > 3000生成的二值数组会出现毛刺（比如[1,1,0,1,1]），对其做差分diff后，就会产生两个值为1的上升沿（从0到1的两次跳变），最终导致误识别。

修正方案：先滤除毛刺再检测边沿 #

解决思路是先对二值信号做滤波处理，消除短暂的波动毛刺，再进行边沿检测。以下提供两种适配不同场景的实现方式：

方法一：形态学开运算（需Matlab图像处理工具箱）

形态学开运算可以有效消除二值信号中的小面积噪声毛刺，适合信号震荡持续时间较短的场景：

% 1. 生成初始二值信号
loc = audio_file > 3000;

% 2. 创建结构元素，消除连续少于N个点的毛刺（N根据实际震荡时长调整）
se = strel('line', 5, 0); % 这里设置长度为5，可按需修改
loc_clean = imopen(loc, se); % 开运算：先腐蚀再膨胀，消除小毛刺

% 3. 检测上升/下降沿
d_loc = [false, diff(loc_clean)];
rising_edges_index = find(d_loc == 1); % 连续上升起始点
falling_edges_index = find(d_loc == -1); % 连续下降终止点

方法二：滑动窗口多数投票（无需工具箱）

如果没有Matlab图像处理工具箱，可以用滑动窗口的方式对二值信号做平滑，通过多数投票过滤短暂波动：

% 1. 生成初始二值信号
loc = audio_file > 3000;

% 2. 滑动窗口平滑：窗口大小根据震荡时长调整
window_size = 5; % 示例窗口长度，可按需修改
% 滑动窗口求均值后再阈值化，相当于多数投票
loc_clean = movmean(double(loc), window_size) > 0.5;

% 3. 检测上升/下降沿
d_loc = [false, diff(loc_clean)];
rising_edges_index = find(d_loc == 1);
falling_edges_index = find(d_loc == -1);

注意事项 #

• 窗口大小或结构元素长度需要根据信号的采样率和实际震荡持续时间调整：比如采样率为1kHz，震荡持续10ms，就把窗口大小设为10左右，确保既能滤除毛刺，又不会误删真实的信号边沿。

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