如何检测方差突变时刻并从此刻开始计算平均距离？

核心需求 #

通过Arduino读取超声波传感器数据并传输到Python，需实现：

• 精准识别传感器数据方差骤升的时刻（对应物体进入检测范围的节点）
• 从该时刻开始，计算后续距离数据的平均值

现有问题：直接计算方差平均值的方法无法有效识别骤升节点。

可行解决方案 #

1. 滑动窗口方差+阈值触发 #

采用滑动窗口实时计算数据方差，当窗口内方差超过预设阈值（基于正常状态方差的倍数或标准差偏移）时，判定为骤升时刻。

代码示例

import numpy as np
from collections import deque

# 替换为你的Arduino串口读取逻辑
sensor_data = [20.1, 20.3, 20.2, 20.0, 15.2, 18.7, 16.5, 14.3, 19.1, 20.2]
window_size = 5  # 根据采样频率调整窗口大小
# 用前N个数据计算正常状态方差基准
normal_var = np.var(sensor_data[:window_size])
threshold = normal_var * 3  # 阈值设为正常方差的3倍，可根据实际数据调整

window = deque(maxlen=window_size)
triggered = False
valid_data = []

for distance in sensor_data:
    window.append(distance)
    if len(window) < window_size:
        continue
    
    current_var = np.var(window)
    if not triggered and current_var > threshold:
        triggered = True
        print(f"方差骤升触发，当前距离：{distance}")
    
    if triggered:
        valid_data.append(distance)

if valid_data:
    avg_distance = np.mean(valid_data)
    print(f"触发后的平均距离：{avg_distance:.2f} cm")
else:
    print("未检测到方差骤升事件")

2. 方差差分法 #

计算连续窗口方差的差值，当差值超过设定阈值时判定为骤升时刻，适合方差上升趋势明显的场景：

import numpy as np

sensor_data = [20.1, 20.3, 20.2, 20.0, 15.2, 18.7, 16.5, 14.3, 19.1, 20.2]
window_size = 5
var_diff_threshold = 2  # 根据实际数据调整差值阈值

var_list = []
# 预计算所有滑动窗口的方差
for i in range(len(sensor_data) - window_size + 1):
    window_data = sensor_data[i:i+window_size]
    var_list.append(np.var(window_data))

# 定位方差骤升点
trigger_idx = -1
for i in range(1, len(var_list)):
    diff = var_list[i] - var_list[i-1]
    if diff > var_diff_threshold:
        # 映射到原始数据的索引
        trigger_idx = i + window_size - 1
        break

if trigger_idx != -1:
    valid_data = sensor_data[trigger_idx:]
    avg_distance = np.mean(valid_data)
    print(f"方差骤升时刻索引：{trigger_idx}，触发后的平均距离：{avg_distance:.2f} cm")
else:
    print("未检测到方差骤升事件")

3. 优化建议 #

• 阈值校准：先采集一段无物体的正常数据，计算其方差的均值和标准差，将阈值设为正常方差均值 + 2~3倍标准差，减少误触发
• 防抖处理：触发后持续监测2-3个窗口，确认方差保持高值再正式启动平均计算，避免短暂干扰导致误判
• 数据预处理：对原始距离数据做滑动平均平滑，过滤传感器随机噪声，提升方差检测的准确性

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