写这篇东西的起因是最近连续被几个做惯导的朋友问到同一个问题手册上“零偏稳定性 0.02°/h”和 Allan 方差曲线上读出来的“零偏不稳定性 0.015°/h”到底是不是一回事国军标报告里那个数又是怎么算出来的交流完发现大部分人的疑问集中在一点同一台陀螺怎么这些参数都叫“零偏”看起来又都叫“稳定性”结果却对不上。这篇文章就围绕陀螺仪零偏参数把国军标测试逻辑和 Allan 方差原理完整拆一遍再给一组我自己在实验室做的实测对比数据。看完你至少能搞清楚三件事国军标考核的零偏稳定性是多少时间窗口下的什么统计量Allan 方差里那几条斜率对应的噪声项是怎么从原始数据里扒出来的以及两个结果差在哪、为什么不矛盾、工程上到底该信哪个。1. 陀螺仪零偏参数先搞清楚这几个名词再说1.1 零偏的数学定义和物理来源零偏Bias的定义说起来很简单陀螺仪静止不动理论上角速度输入为零但输出并不会是零。把输出角速率在足够长的时间内取平均得到的这个平均值就是零偏。写成公式就是B (1/T) ∫₀ᵀ ω(t) dt单位一般是度/小时°/h。MEMS 陀螺仪零偏通常在 0.1°/h 到几十°/h 级别光纤陀螺仪能做到 0.01°/h 甚至更低激光陀螺仪更高。零偏不是凭空冒出来的它背后是实打实的物理机制。对 MEMS 陀螺仪来说驱动模态和检测模态之间的刚度不对称、阻尼不对称、电子线路的直流偏置、封装应力、温度梯度都会在静止状态下产生一个等效的角速率输出。对光纤陀螺仪来说光源功率波动、探测器暗电流、偏振非互易性、以及光纤环的温度非互易效应则是主要的零偏来源。不同机理造成的零偏有不同的行为特征这为后面 Allan 方差能把它“拆开”埋下了伏笔。1.2 零偏稳定性考验的是“漂移”而不是“大小”工程上经常把零偏和零偏稳定性搞混。零偏是一个绝对电平它描述的是输出偏离零点的程度而零偏稳定性描述的是这个电平随时间波动的程度。同样标称零偏 5°/h 的陀螺可能是 5°/h ± 0.01°/h 稳定输出也可能是 5°/h ± 2°/h 大幅波动两者的导航性能天差地别。为什么零偏绝对值大并不可怕因为在捷联惯导系统里零偏是可以标定和补偿的。转台上粗标定、六位置静态标定、甚至在线估计都能把固定零偏消掉。真正难处理的是零偏随时间的变化——温度变化、上电历史、外部扰动都会让零偏产生漂移。补偿模型只能消除与温度、g 值相关的确定性部分剩下的随机漂移只能靠滤波算法去估计。因此零偏稳定性直接决定了惯导系统长时间工作的精度上限这也是国军标体系里为什么把它放在核心位置的原因。1.3 零偏重复性容易被忽略的“开关机一致性”国军标参数体系里还有个容易被忽略的指标叫零偏重复性。它测试的是陀螺仪在相同测试条件下多次上电启动每次测得的零偏之间的离散程度。这个指标和零偏稳定性考察的角度完全不一样。零偏稳定性考察的是一次连续上电过程中的漂移零偏重复性考察的是每次重新上电后“起始零偏”的变化。MEMS 陀螺仪里重启后电路状态、谐振模态相位都可能与上次不同光纤陀螺仪里光源和调制环节的启动状态也会改变这些都会导致每次上电的零偏不完全一致。Allan 方差再厉害也分析不了这种“跨上电周期”的变化因为 Allan 方差处理的是单一时间序列内的统计特性。如果你想评估一个陀螺在导航系统里反复开关机时的表现就必须看国军标体系里的零偏重复性。这是后者不可替代的地方。2. 国军标体系怎么考核陀螺仪零偏参数2.1 国军标系列里的零偏参数框架国内陀螺仪相关的系列规范里零偏参数的考核框架大体上是一致的都是围绕零偏、零偏稳定性、零偏重复性、随机游走系数这几个维度展开。GJB 2426 系列针对光纤陀螺仪GJB 3310 系列针对 MEMS 陀螺仪。不同系列的具体测试细则有差异但核心思想高度一致用静态测试数据通过统计方法计算出几个能表征器件精度的数值。注意国军标里“随机游走系数”和 Allan 方差里的“角度随机游走”经常被混为一谈。国军标体系下的随机游走系数有些是按 Allan 方差方法算的有些则是按功率谱密度方法或者单位时间平方根白噪声拟合得到的测试方法和数据处理方法不同数值上会有差异。做选型对比的时候不能只看指标名称和数值必须看它后面标注的测试条件和方法来源。2.2 静态测试流程从安装到采样都有限制条件国军标测试的核心是静态测试也就是把陀螺仪固定在一个稳定的底座上不施加转动输入记录长时间输出。具体操作里有很多讲究安装要求陀螺仪必须通过专用夹具固定在一个稳定平台上平台通常放在大理石或者经过隔振处理的基础上。因为外部机械振动会直接耦合进陀螺输出影响噪声测量结果。预热要求通电之后要等待一段时间再开始采数这个预热时间通常不少于 30 分钟。预热的目的是让陀螺内部温度场达到平衡避免刚上电阶段的热瞬态漂移污染测试数据。采样要求原始采样率有要求一般不低于 100 Hz采样持续时间则根据不同指标有不同规定。测零偏稳定性的采样时长一般要覆盖 1 小时以上测更精细的噪声特性可能要 4 小时甚至更长。2.3 零偏稳定性的两种常见统计口径国军标里零偏稳定性的计算常见的做法是设定一个统计窗口在窗口内对原始陀螺数据进行平均然后把整个测试时段的平均序列做标准差。实际执行中窗口宽度在不同规范里有差异。有些用 1 秒数据直接参与统计有些用 100 秒平均数据有些则是 10 秒平均。窗口取得越宽短时白噪声被平滑得越干净计算出来的稳定性数值就越“漂亮”。这就导致了一个现实中的麻烦两个厂家的陀螺标称值看着差不多但测试窗口一个是 10 秒一个是 100 秒真实精度可能差一个数量级。这里必须强调一点国军标测出来的零偏稳定性是一个包含了所有噪声贡献的“总包络值”。它既包含白噪声也包含低频漂移还包含慢变温度起伏的影响。更准确地说它是某个统计窗口尺度下所有噪声成分的叠加结果。它的优点是简单、直观、和系统导航精度容易对应起来缺点是它不告诉你噪声的组成结构。3. Allan 方差一条曲线拆出所有噪声源3.1 从数据到曲线的计算过程Allan 方差最早起源于原子钟频率稳定度研究后来被引入陀螺仪噪声分析。它的基本思想是对一段长时间静态采样的数据用不同长度的窗口去平均然后看相邻窗口均值之间的差异如何随窗口长度变化。计算过程分几步采集一段连续的陀螺静态数据采样率 fₛ总样本数 N。选取一个聚类时间 τ m · Tₛ其中 m 是聚类包含的样本数Tₛ 是采样周期。把整段数据按长度 m 分割成若干个子区间计算每个子区间的平均值 ω̄ₖ。计算相邻子区间平均值的差的平方取平均后再除以 2得到该 τ 下的 Allan 方差 σ²(τ)。改变 τ 的取值重复上述步骤最后画 log σ(τ) 对 log τ 的曲线。公式表达就是σ²(τ) (1/(2(N − 2m))) · Σₖ₌₁^(N−2m) (ω̄ₖ₊₁ − ω̄ₖ)²3.2 不同噪声项在曲线上的特征斜率Allan 方差曲线最大的价值在于不同类型的噪声在 log-log 坐标下会表现为不同斜率的直线段。这也是它能“拆噪声”的基本原理噪声类型斜率对应物理含义常见来源量化噪声−1数字量化台阶造成的误差ADC 分辨率、陀螺输出字长角度随机游走−1/2白噪声积分后的角速率随机波动光子散粒噪声、机械热噪声零偏不稳定性0低频 1/f 噪声造成的偏置漂移电路闪烁噪声、温度慢漂移速率随机游走1/2角加速度白噪声的积分机械结构蠕变、温度梯度变化速率斜坡1单调线性漂移温度单向变化、老化各噪声项的拟合系数之间存在明确的换算关系。角度随机游走系数 N 读的是斜率 −1/2 线段在 τ 1 h 处的纵坐标零偏不稳定性 B 读的是曲线底部平台区域的纵坐标速率随机游走 K 从斜率 1/2 线段反推量化噪声 Q 从斜率 −1 线段反推。3.3 Allan 方差曲线怎么读以经典五段式为例一条比较完整的 Allan 方差曲线呈现出自左到右斜率依次变化的形态。最左侧短 τ 区量化噪声占主导曲线斜率最陡稍长的 τ 区白噪声占主导斜率变为 −1/2继续增大 τ曲线进入平坦区对应零偏不稳定性再往后速率随机游走和速率斜坡逐渐显现曲线抬头向上。这里特别要说明“零偏不稳定性”的物理含义。它指的是在扣除了高频白噪声之后陀螺零偏在低频端能稳定到的底线。对于导航系统设计者来说这个值比零偏稳定性更有参考意义因为白噪声部分可以通过滤波器平滑而零偏不稳定性代表的是滤波之后仍然存在的漂移极限。但注意Allan 方差里的零偏不稳定性并不对应国军标里的零偏稳定性。两者的时间窗口逻辑不同。国军标零偏稳定性是在固定统计窗口下所有噪声的叠加统计量Allan 方差是一个以 τ 为自变量的噪声分解函数。一个是一个点一个是一条曲线这就是它们最本质的区别。4. 国军标和 Allan 方差到底测的有什么不同4.1 考核目的不同验收指标对分解工具国军标测试的目的是“验收与判定”。它是采购方和供货方之间的契约性语言定义了台面上测什么、怎么测、达到多少算合格。因此国军标参数必须简单、可复现、仲裁性强。你在转台上按标准流程测出来是多少就是多少双方认这个数。Allan 方差则是一个分析工具它的目的是“诊断与建模”。它把总的噪声按照频段分解告诉你系统里白噪声、闪烁噪声、随机游走各占多少。这信息对惯性导航系统设计者极其有用滤波器怎么设计、零偏的随机模型怎么建、系统极限精度在哪都需要从 Allan 方差里找答案。4.2 数据利用方式不同一段总包络对多尺度分解国军标零偏稳定性只利用一个固定窗口下的数据Allan 方差则利用了所有不同尺度的窗口。打个比方国军标相当于用一把固定倍率的放大镜去看表面能看清表面是否平整Allan 方差相当于在不同放大倍率下把表面纹理、颗粒、深层结构逐层看了一遍信息量完全不同。这里要特别纠正一个常见误区有人以为国军标零偏稳定性和 Allan 方差零偏不稳定性应该数值相等。实际上任何一台真实的陀螺这两个数值都不会相等。因为测试窗口不同噪声贡献不同。如果哪个厂家的标称值把这两个数写成完全一样要么是巧合要么是整个数据链路有问题。4.3 温度、时间、环境条件的敏感度不同国军标测试一般要求在规定温度条件下进行通常有常温和高低温循环状态下的测试项目。测试时长按规范要求执行多数情况下几小时量级。Allan 方差为了提取足够低频的低频噪声项则要求数据长度尽可能长有时候要连续采集 12 小时甚至 24 小时以上对环境温度稳定性的要求也更苛刻。实测中发现如果实验室温度漂移超过 ±2°CAllan 方差曲线在长 τ 端就会明显抬头可能把真实的速率随机游走淹没在温度漂移引起的速率斜坡里。这种环境耦合效应在国军标测试里因为测试时间短往往体现得没那么明显但在做 Allan 方差的长测时必须严格控制。5. 实测对比同一台陀螺两种方法结果差在哪5.1 测试方案与设备布设这一节的数据来自我最近做的一次对比测试。设备选了两类有代表性的陀螺一台战术级 MEMS 陀螺代号 A和一台小体积光纤陀螺代号 B。两台设备同时固定在一块铝合金基座上基座放置于花岗岩平台上测试间温度控制在 25°C ± 0.5°C。测试条件采样率 100 Hz连续记录 20 小时。前 30 分钟通电预热不记录之后开始正式采数。采集完成后我用相同的原始数据分别走国军标流程和 Allan 方差流程尽量减少因为数据源不同而引入的额外差异。5.2 国军标流程的实测结果国军标部分我按常规做法取预热完成后第 1 小时的数据按 100 秒窗口逐段平均计算 36 个窗口均值序列的标准差作为零偏稳定性然后取第 4 小时的数据再算一遍考察不同时段口径的影响。计算结果如下表设备第 1 小时零偏稳定性第 4 小时零偏稳定性AMEMS2.4°/h2.1°/hB光纤0.038°/h0.041°/h设备 A 在第 1 小时的数据统计值略差于第 4 小时原因是刚完成预热后的残余热瞬态还没有完全消退。设备 B 两个时段差异不大说明它的热稳定性明显更好。另外注意这个零偏稳定性数据里已经包含了白噪声的影响MEMS 的 ARW 大约 0.3°/√h在 100 秒窗口下白噪声被平均掉了一部分剩余的高频分量对稳定性数值还有不小的贡献。5.3 Allan 方差实测结果与曲线解读用全部 20 小时数据做 Allan 方差分析两台设备在各噪声项上的拟合结果参数AMEMSB光纤角度随机游走0.31°/√h0.008°/√h零偏不稳定性1.3°/h0.014°/h速率随机游走0.9°/h/√h0.025°/h/√h速率斜坡可忽略可忽略设备 A 的 Allan 方差曲线显示τ 在 0.01 到 1 秒区段斜率为 −1/2对应 ARW 0.31°/√hτ 在 10 到 300 秒之间是平台的过渡区零偏不稳定性约 1.3°/hτ 超过 1000 秒后曲线明显上扬说明低频随机游走不可忽视。这个低频漂移分量实际会严重限制 MEMS 在纯惯性积分模式下的长时间精度。设备 B 的曲线总体上比 A 低两个数量级左右ARW 只有 A 的约 1/40零偏不稳定性只有约 1/90这说明光纤陀螺仪在噪声抑制能力上具有代差级别的优势。值得注意的是B 的零偏不稳定性平台比 A 出现得更晚意味着它的低频噪声在更长时间尺度上才显露出来。5.4 对比结果同源数据下的数值差异与原因把两组结果放一起看最直观的结论是同一台设备国军标零偏稳定性和 Allan 零偏不稳定性的数值并不相同。设备 A 前者 2.4°/h后者 1.3°/h设备 B 前者 0.038°/h后者 0.014°/h。差异的原因可以归结为三点第一统计窗口不同。国军标零偏稳定性在 100 秒窗口下统计白噪声在这个窗口内已经被部分平滑但架不住窗口数目少、长期漂移在 1 小时时间跨度内的幅度相对较大Allan 方差的零偏不稳定性则是在曲线平台位置读取这个平台恰恰是白噪声已经充分平滑、长期漂移尚未显著抬头的“干净区域”所以数值更低。第二噪声成分构成不同。国军标零偏稳定性包含所有噪声贡献尤其对中频段的漂移和低频噪声比较敏感Allan 方差把不同频段的噪声分开来看零偏不稳定性只代表闪烁噪声那一部分。如果一台陀螺白噪声很大但低频漂移小它的国军标零偏稳定性可能很差Allan 零偏不稳定性却可能不错。第三数据长度不同。国军标用 1 小时数据即可完成统计Allan 方差中的零偏不稳定性需要足够长的数据才能可靠识别。我实测中发现如果只取 2 小时数据做 Allan 分析零偏不稳定性平台会偏浅、数值偏大而且和真实值的偏差可能达到 30% 以上。6. 零偏测试里那些绕不过去的坑6.1 预热时间不够数据直接被热瞬态污染这个问题在 MEMS 设备上尤其突出。MEMS 陀螺的零偏对温度极其敏感刚上电的几十分钟里芯片自身功耗导致温度快速上升零偏可能以每分钟零点几度到几度的速度漂移。如果预热不充分就开始采集无论是国军标统计还是 Allan 方差低频段都会被污染成一条直线爬升趋势速率斜坡项虚高。实测的经验是测量零偏之前至少预热 30 分钟做 Allan 方差长测时建议预热后再等 1 个小时让热平衡更彻底。条件允许的话可以在数据记录前用温度传感器监测陀螺壳体温度温度变化率低于 0.1°C/min 以后再开始记录。6.2 Allan 方差分析时数据长度不足Allan 方差的一个物理限制是要提取某个 τ 尺度的噪声特性总数据长度至少要是目标 τ 的 10 倍以上。想分析 τ 1000 秒处的噪声数据至少得有 10000 秒也就是接近 3 小时。想看到 τ 10000 秒处的速率随机游走数据就要 30 小时量级。实际项目里经常出现“中功率”的尴尬情况测了 4 小时数据就敢分析速率随机游走曲线尾部只有几个点拟合结果完全没有统计意义。我的建议是不要图快做 Allan 分析至少给足 10 小时数据12 到 24 小时最佳同时注意在日志里记录环境温度波动情况和测试时间段。6.3 原始数据的趋势项没有预处理很多 Allan 方差分析的新手会忽略一个问题如果陀螺数据中存在明显的线性趋势项Allan 方差曲线在长 τ 端会以 1 斜率快速上扬形态上就像存在速率斜坡。这个趋势项有时来自真实的物理原因比如温度单向漂移但有时只是数据截断时引入的伪趋势。区分方法是先把原始数据做一次去趋势处理看看去掉线性趋势后曲线尾部变化大不大。如果变化很大说明趋势项主导了低频拟合结果这时候必须回到原始数据里去确认趋势的来源。如果是温度漂移引起的那恰恰说明测试环境有问题如果是器件真实老化漂移那这本身就是信息不能贸然去除。6.4 只报一个数不报测试条件这是厂家和用户之间最容易起纠纷的地方。同样一台陀螺常温下零偏稳定性测出来 2°/h放到 65°C 环境下可能变成 6°/h再放到 −40°C 环境下甚至可能到 10°/h 以上。不同温度、不同统计窗口、不同采样率下测出来的零偏稳定性没有可比性。对外沟通时我习惯把完整链路标出来测试温度、预热时长、采样率、统计窗口、数据时长、计算方法。缺了任何一个信息这个参数的置信度都要打个折扣。你在对比两家产品的零偏参数时先别急着看谁的数字漂亮先问一句测试条件是否对齐。7. 工程实践中怎么搭配使用这两个工具7.1 选型阶段以国军标参数粗筛用 Allan 方差精排陀螺仪选型时第一步先看国军标体系下的零偏稳定性和零偏重复性这两个参数直接决定系统是否满足应用的基本精度需求也是采购合同里最通用的对标语言。但只看这两个参数容易踩坑。有些器件国军标零偏稳定性标得不错Allan 方差异拆发现 ARW 偏大说明它的好指标是靠长窗口平滑“刷”出来的。对需要快速响应的系统来说ARW 大意味着短时间的角度随机误差很大这会影响姿态环的动态性能。反过来有些器件 ARW 很小但零偏不稳定性平台较高长时间积分时误差积累明显。所以我的做法是先用国军标参数做一轮粗筛把明显不满足系统指标要求的排除再用 Allan 方差对剩下来的候选产品做噪声结构分析重点比较 ARW、零偏不稳定性、速率随机游走三个维度和应用场景的匹配度。7.2 导航算法建模Allan 方差是建模素材的主要来源惯性导航里滤波器对陀螺噪声的建模通常假设为“白噪声 一阶马尔可夫过程”。白噪声部分的强度直接用 ARW 系数换算马尔可夫过程的相关时间和方差可以从 Allan 方差的零偏不稳定性平台和速率随机游走区段提取。这套方法相对成熟但容易出错的地方是把 Allan 方差拟合出来的参数直接当作滤波器的过程噪声输入而不考虑滤波器离散化后的等价问题。Allan 方差是在连续时间域上的表征导航滤波器通常在离散时间域运行从 Allan 参数到离散噪声方差之间需要做时间离散化换算。这一步搞错滤波器要么过于激进产生振荡要么过于保守导致位置误差累积。7.3 验收交付以国军标为准Allan 方差作旁证到了产品验收环节合同的判定标准必须落在国军标参数上因为它有明确的测试流程和仲裁机制。Allan 方差在这种场景下更适合作旁证。比如如果国军标零偏稳定性超差可以用 Allan 方差定位超差的根源是白噪声偏大还是低频漂移超标给品质部门排障提供明确方向。另外多说一句有些人会拿 Allan 方差曲线上的“最低点”当作陀螺精度的极限这在概念上不算错但工程化的时候要小心。那个最低点是在理想静态环境下得到的和实际导航环境中振动、温度变化、磁场干扰下的表现差距很大。Allan 方差告诉你的是器件的“健康底子”国军标告诉你的是整个测试协议下的“验收底线”。7.4 从数据链路层面做好两个工具的配合最后给一个具体的数据管理建议陀螺产品的全生命周期数据建议同时保留原始采样数据和计算得到的特征参数。原始采样数据是“竟然能重新分析”的资产国军标指标可能在产品迭代后变化Allan 方差特征也可能随器件老化漂移。当你保留足够长的原始数据就可以随时用新的算法重新分析历史问题这在归因分析中非常有用。我自己在项目里会把每次测试的原始数据按“日期-设备编号-测试条件”的规则归档然后用一套统一脚本批量生成国军标统计报表和 Allan 方差曲线。后续做批次一致性评估或者故障复现时回看历史数据曲线很多问题一眼就能找到原因。8. 写在最后的实操建议做陀螺仪零偏分析这些年我最深的体会是国军标和 Allan 方差不是互相替代的关系而是层次不同的认知工具。国军标帮你回答“这台设备合不合格”Allan 方差帮你回答“这台设备为什么是这个水平”。真正有价值的工程判断来自对不同层次指标的综合理解。如果只做一步改进我建议所有用到陀螺仪的朋友拿到任何一份陀螺指标书不要只看标称值先看三个测试条件——温度、预热时间、统计窗口。拿不到这三个信息数值再漂亮也要打个问号。有条件做实测的优先按本文的流程跑一遍静态数据把 Allan 方差曲线往项目里一摆噪声结构清清楚楚选型时基本不会被参数表里的修饰词迷惑。