如何准确计算timespec结构时间间隔的移动平均值

问题根源 #

你现在的方案有两个核心错误，整数截断反而是次要影响因素：

1. 逻辑错误：不能独立平均tv_sec和tv_nsec字段
   tv_nsec的取值范围永远在0~999999999之间，两个字段是进位关系，不是独立的统计维度。举个最简单的例子：两个间隔分别是1.9秒（tv_sec=1, tv_nsec=900000000）和2.1秒（tv_sec=2, tv_nsec=100000000），真实平均间隔是2秒，但独立平均两个字段会得到tv_sec=1、tv_nsec=500000000，也就是1.5秒，误差直接到0.5秒，和整数截断没有关系。
2. 计算顺序错误放大截断误差
   你用的平均公式拆成了两次除法：old_avg*(n-1)/n和new_value/n，每一步除法都向下取整丢精度，累积下来误差会越来越大。

最简可行方案 #

不要拆分时间字段做计算，所有间隔统一转成单值的64位纳秒整数计算，全程不需要碰浮点数，精度足够，代码也简单：

1. 你要统计的是消息的接收间隔，不是时间戳本身，所以每次收到新消息时，先和上一条消息的时间戳做差，把差值转成int64_t类型的纳秒数
2. 调整平均公式的计算顺序，把两次除法合并成一次，所有分子算完再做除法，把截断误差压到最小
3. 只有在需要输出timespec格式结果的时候，再把纳秒值拆成秒和纳秒字段即可

参考实现代码：

#include <stdint.h>
#include <time.h>

// 全局状态，按需改成结构体封装也行
static struct timespec last_ts = {0};
static int64_t avg_interval_ns = 0;
static uint64_t msg_cnt = 0;

/**
 * 每收到一条目标类型消息，传入当前时间戳调用即可
 * 计算完成后全局变量avg_interval_ns就是平均间隔（单位：纳秒）
 */
void update_msg_interval_avg(struct timespec now) {
    msg_cnt++;
    // 第一条消息没有前置间隔，仅记录时间戳
    if (msg_cnt == 1) {
        last_ts = now;
        return;
    }

    // 计算本次消息和上一条的间隔，统一转纳秒
    int64_t sec_delta = (int64_t)(now.tv_sec - last_ts.tv_sec);
    int64_t nsec_delta = (int64_t)(now.tv_nsec - last_ts.tv_nsec);
    int64_t interval_ns = sec_delta * 1000000000LL + nsec_delta;

    // 改进的移动平均公式，仅做一次除法，截断误差不超过1ns
    avg_interval_ns = (avg_interval_ns * (msg_cnt - 1) + interval_ns) / msg_cnt;

    // 更新上一条消息时间戳
    last_ts = now;
}

/**
 * 需要取timespec格式的平均间隔时调用这个函数转换
 */
struct timespec get_avg_interval_timespec() {
    struct timespec res;
    res.tv_sec = avg_interval_ns / 1000000000LL;
    res.tv_nsec = avg_interval_ns % 1000000000LL;
    // 处理负余数（正常时间单调递增的场景不会触发，加个兜底）
    if (res.tv_nsec < 0) {
        res.tv_sec -= 1;
        res.tv_nsec += 1000000000LL;
    }
    return res;
}

补充说明 #

• 这套实现全程用整数运算，没有浮点开销，哪怕是内核态、嵌入式无浮点单元的场景也能直接用
• 单次计算的截断误差最大不超过1纳秒，没有跨字段的逻辑误差，完全满足消息间隔统计的精度要求
• 如果你的场景允许用浮点数，直接把avg_interval_ns改成double类型，连1纳秒的截断误差都能消掉，代码逻辑不变
• 64位纳秒值可以表示正负292年的时间长度，哪怕msg_cnt到百万级，计算时的中间值也不会溢出，绝大多数场景不需要考虑溢出问题。如果确实要处理千万级以上的样本量，可以给平均值预留32位的小数位做定点运算，精度还能再提。

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