陀螺仪动态随机建模概述

  • 格式:pdf
  • 大小:537.16 KB
  • 文档页数:14

《舰船导航》2004年第6期 

陀螺仪动态随机建模概述 

摘要:多年以来,惯性系统界一直在使用各种各样的方法来估算陀螺仪性能参数,包 

括陀螺仪动态漂移和随机漂移特性。自1965年以来,由于数字信号处理技术的进步, 

使得现代建模方法得到了有效的应用,但在系统分析领域外还没有引起广泛的重视。 

为了推广现代方法的应用及其标准化,IEEE hESS陀螺仪和加速度计委员会正在推出 

根据应用最广方法制订的规范化建模方法,即IEEE Std 647的修订版,“单轴激光 

陀螺仪规范格式及测试程序”.因此,下面的讨论仅作为动态随机建模的介缗陛论述。 

其意义在于将其应用于陀螺仪建模、仿真,卡尔曼(Katman)滤波器设计,误差分 

析,性能的规定、预测和评估。模型的通用形式包括:物理性能方程的数学描述; 

由扰动模型和环境敏感模型构成的误差模型:描述随机漂移特性的随机模型;由输 出状态量与附加的测量噪声的缌陛组合构成的测量模型.与之相联系的陀螺仪模型 

通用方程,包括惯性输入的响应,环境敏感项,漂移率和标度因数误差项。本文讨 

论了确定I/0特性的方法.介绍了由时间序列分析建立随机模型的方法,包括频域 

法和时域法。艾伦方差(hl 1an Variance)和PSD的应用是本文的重点。本文介绍 了测试和分析方法,讨论了数据采集、数据整理、数据预处理和结果的评估。最后 列出了按时间顺序排列的相关参考文献,以便读者进一步了解有关内容。 

1 引言 

下文介绍了动态随机建模在陀螺仪建模和性能评估方面的应用。旨在作为概括的评述而 不是全面详尽的讲解。 

在2O世纪4O年代初期,诺伯特·维纳(Norbert Wiener)提出了广义的非线性问题的解 决思路[21:假定待分析系统是一“黑盒”,用已知结构的系统即“白盒”来辨识和表征此待分 

析系统。 

线性问题的解决方法如下:利用各种时域和频域的方法找出可行的黑盒系统的等效体, 此等效体可以由一些“白盒”标准模块组合构建而成。虽然模型的结构与实际的结构可能会 

有所不同,但两者的输入.输出特性是等效的。 

2发展概况 

早在1800年,高斯(Gauss)的最:j ̄--乘估计方法为模型技术奠定了最初的基础【 。。

。目 维普资讯 http://www.cqvip.com 《舰船导航》2004年第6期 

前,各种系统稳态输入.输出特性的确定都是基于这种方法产生的。 

百多年后,也就是1905年,爱因斯坦(Albert Einstein)和吉布斯(Willard Gibbs)创 

建了独立的统计物理学科。到1910年,费希尔(Fisher)将概率密度函数用到最大似然估计 

中【 们。1930年,维纳首次明确使用了频域分析方法,并于1940年创立了现代优化理论【2O]。 

二战期间,出现了博弈论【 和运筹学[67】。这些都涉及到了现代建模技术的最早应用。 

20世纪50年代,随机过程和微分博弈论得到不断发展。1960年,卡尔曼提出了在时域 内进行最佳递归滤波器设计的方法【 们。到1963年,产生了信号识别【2O]和频域时间 序列分析方法【 们。1965年,特克(Tukey)和库里(Coolie)发表了著名的关于快速傅 

里叶变换(FFT)的论文【20】。1966年,大卫.艾伦(David Allan)提出将一种简单的方差分析 

方法用于振荡稳定性研究,此后该方法被用于陀螺仪漂移分析中【 -_ 。到1968年, 

参数辨识方法已被人们熟知【2o】。 

随后的十年间,敏感器的时域和频域特性受到了重视。到1970年,波克斯一詹肯斯 

(Box—Jenkins)时间序列分析方法[14,20,331连同系统识别和卡尔曼自适应滤波技术【2O】一 

起得到了发展。2O世纪8O年代,非线性多输入多输出(MIMO)随机最优控制/估计得到人 

们的关注。到2O世纪9O年代,人工智能专家系统思想在建模技术中得到了应用。 

3单元模型 

“单元”是指 

能够实现某一特 

定功能的可运行 的实体,它可以是 

敏感器、系统,或 其他装置。单元模 

型可分解为几部 

分,如图1所示, 第一部分,实体, 

或称为物理模型, 

用微分方程或代 

数方进行描述,表 

征了其运行的物 理原理。这是应用 图1系统模型 

动态建模技术研究实体的关键性部分。第二部分,误差模型,包括扰动模型和环境模型,扰 动模型涵概了对实体方程参数变动的灵敏性,环境模型包含了对环界干扰的灵敏性第三部 分,是随机模型,包含了在其它运行都良好的条件下观测到的随机漂移。测量模型则由系统 

2 维普资讯 http://www.cqvip.com 《舰船导航》2004年第6期 

状态量和输出噪声的线性组合构成。 

其想法如下:如果可以对受测单元建立模型,就可以充分利用模型来估计其性能,甚至 

可以将性能提高到更高的水平。也就是说,可以通过最优滤波将系统误差降低到最小残留白 

噪声的范围内。该模型也可以用来根据规范进行性能预测评估。 

4陀螺仪模型方程 种适用于多种敏感器的通用模型方程如图2所示。它包括惯性输入因素项、(含不对准 

在内)环境敏感因素项和随机(含量化噪声在内)因素项。 

图2通用模型方程 

为了更好地组织各种模型组分,在RLG标准中提出了这种陀螺仪模型方程的分解方法

对于一些陀螺仪,如动力调谐陀螺仪(DTG:IEEE Std 813),完整的模型方程含有很多项, 详见IEEE陀螺仪标准。 

5建模 建模技术主要应用于仿真研究 性能估计和卡尔曼滤波器的设计。动态建模与随机建模 

最基本的差别如下: 

在动态建模中,给定输入向量(一个或多个输入)和预期的输出向量(一个或多个输出), 从两者的时间序列来确定输入/输出关系。应用包括在输出端叠加随机噪声的情况

而在随机建模中,可能无法直接得到输入。可以假定这样一种(输入白化)模型,由白 

3 维普资讯 http://www.cqvip.com 《舰船导航》2004年第6期 

噪声激厉JJ、与被测单元具有相同的输出特性。这类模型并不罕见,故可以选出几种规范的形 

式。例如,美国标准技术局(NIST)的大卫.艾伦(David Allan)利用幂级数表示功率谱密度 

(PSD)并利用时域相应的方差分析进行振荡稳定性分析【8】。在IEEE Std 647附件C(Annex C of the IEEE Std 647)中介绍了这种方差分析。其想法为,由一个或多个假定的、强度为Ni 

的白噪声源激励一规范化的传递函数,得到的响应与实际研究对象(黑盒模型)具有相同的 

统计特性与谱特性。这也是陀螺仪漂移分析的目标。 

6动态建模 除去环境灵敏性外,光学陀螺仪的动态建模仅涉及到标度因数、偏值平"不对准。在确定 

性输入的驱动下,利用回归法可以得出陀螺仪的动态模型。 动态建模最适合于是自旋式转子(Spinning whee1)陀螺仪,如两自由度陀螺仪或表现出 

很强的两自由度效应的单自由度陀螺仪【 。后一类陀螺仪其径向刚度系数很小,或者角动 量相对于刚度系数米说很大,因此角动量对陀螺仪响应具有很大的影响作用。得到的结果是 与动力调谐陀螺仪(DTG)模型方程类似的两自由度陀螺仪模型。该动态模型的特征响应为 

四阶,四个特征值取决于陀螺仪的角动量。其频域传递函数可田下述方法进行分析。 僧伽(Sinha)和库斯坦(Kuszta)[331给出了多种估算传递函数的方法。例如,频率响应 

法就是较早使用的经典方法之--I 。在FFT出现之前,将输入信号进行离散频率阶跃变化, 

同时可测出由输入到输出的相对幅值与相位。目前使用的方法,白噪声常常是与求和点成比 

例的(力矩)速率插入。利用数字处理技术计算出开环传递函数和闭环传递函数。该方法在 

出现附加的不相关噪声时仍可进行估算。 时间响应法是利用阶跃响应和脉冲响应的另一种估算传递函数的经典方法。通过对采样 

数据进行卷积积分(离散形式)可以由输入、输出确定其脉冲响应,此方法被称作反卷积p引。 

利用一矩阵方程可以将各时间点的输出表示为输入的函数。该方程的解即是被测对象黑盒模 

型的脉冲响应。 另一种利用脉冲(或阶跃)响应法的途径是将离散传递函数(Z变换形式)表示成为两幂 级数的比值式【3 。得出两个矩阵方程,其系数以输出脉冲响应表示。根据脉冲响应数据可以 

求解出分母系数,进而根据分母系数和输出数据可以求解出分子系数。 有儿种方法可用与离散时间传递函数相对应的微分方程来模拟输出【3 。输出可以由其过 

去值(即自回归部分)和输入的过去值及当前值(滑动平均部分)来表示。一矩阵方程把构 

成传递系数的参数向量和一系列输入、输出数据(表示为矩阵)相联系,此输出向量又构成 

了另外一系列的连续输出数据。可以利用LSE(最小平方)或MLE(最大似然估计)方法将该 

问题转化为噪声数据来处理。推导出的递归公式用于实时估计。 

另一种方法是先引入白噪声输入,再分析输出相对输入的关系。该时域方法称为相关法【 ’ 乃J计算输入输出的互相关值,由此推导出特性曲线。该方法仅限于由线性时不变系统产生的 

4 维普资讯 http://www.cqvip.com 《舰船导航》2004年第6期 

平稳时间序列。此方法的主要问题在于:当模型变得比较复杂时,很难确定其互相关函数。 尤其困难的是,当相关项分布于很大的动态范围内或不同频率成分相互交迭时,很难把它们 

进行分解。如果模型的结构形式也是未知的,那么识别工作的难度将很大。这是Box—Jenkins 

时域时间序列分析法以及大多数动态模型法存在的问题之一。 频域方法【l 1(也称为谱分析)是一种较好的模型研究工具。相关函数的傅立叶变换即为 

功率谱密度(PSD)。广泛采用的方法是利用FFT从原始数据计算PSD。求总体均值的实时处 

理大大提高了估算传递函数的能力。利用频域方法,可以用XPSD(输入输出的互功率谱密度) 除以PSD(输入的功率谱密度)来估算传递函数[91,同时还给出了传递函数的幅值与相位。 

7随机模型 

对于随机模型来说,很重要的一点就是引用白噪声及由此构造传递函数的思想。之所以 

这样做是因为:如果输入是白噪声,就可以直接由输出的PSD估算出最小相位线性时不变系 

统的传递函数。此时仅由输出的PSD就可以推算出传递函数,而不需要再求出输入输出的 

XPSD。相位信息则可以由幅值响应来确定。 

因此,对于线性时不变系统,假设输入为白噪声,只要输出已知,就可以表述未知的模 

型。除了输入不可观测以外,随机建模的许多方法与一些动态建模的方法是相似的。利用PSD 

估算传递函数的频域方法属于直接估计类型。即使在特定的失常状况下,比如表现为l/f过程 

的偏值不稳状态(颠振速率噪声)、角度量化噪声(在特征上不同于连续的角度白噪声),也 

可以利用精确的分析技术进行辨识。 

如同动态建模的情况一样,在时域进行随机建模也有几种方法可行。与PSD方法对应的 

是相关函数法,两者的联系为傅立叶变换对关系。类似于相应的动态建模方法,可以由输出 

数据计算出自协方差序列进而推导出等效的脉冲响应。有一种方法是将协方差函数拟合成为 

指数项与阻尼正弦项之和,利用最小二乘估计法来求取模型参数。这类似于用部分分式展开 

式来表示频率响应函数。 

另一相关法表示为自回归滑动平均(ARMA)过程,它把自协方差序列与差分方程的系