allan方差分析法.docxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
allan方差分析法

《现代导航测试技术》 实验报告 实验名称:光纤陀螺随机误差的Allan方差分析法 班 级:0309103 学 号:030910309 姓 名:许 莹 时 间:2012-12-17 一:实验目的 由于光学陀螺的工作原理和环境干扰等原因,在光学陀螺输出信号中包含很多确定性和随机性的误差项。光学陀螺的随机误差主要包括量化噪声、角速度随机游走、零偏不稳定性】角速率随机游走、速度斜坡和正弦分量,其中前三项误差被认为是其光学性能指标一部分。对于这些随机误差,利用常规的分析方法,例如计算机样本均值和方差并不能揭示潜在的误差源,另一方面,虽然自相关函数和功率谱密度函数分别从时域和频域描述了随机误差的统计特性,但是在实际工作中通过这些函数加以分析将随机误差分离出来是很困难的。 Allan方差法是20世纪60年代由美国国家标准局David Allan提出的,它是一种基于时域的分析方法,不仅可以用来分析光学陀螺的误差特性,而且还可以应用于其他任何精密测量仪器.Allan方差法的主要特点是能非常容易对各种误差源及其对整个噪声统计特性的贡献进行细致的表征和辨识,而且便于计算,易于分离。它提供了一种识别并量化存在于数据中的不同噪声项。 二:实验原理与实验内容 1.Allan方差定义与计算 设以采样时间0对陀螺仪输出角速率进行采样,共采样N个点,把所获得的N个数据分成K,每组包含M个采样点。 K=N/M,M (N-1)/2 如图: 1,2,…,M M+1,M+2,…,2M N-M+1,N-M+2,…,N K=1 k=2 k=K 每一组的持续时间M=M0,称之为相关时间,每一组的平均值为 k(M)= k=1,2,3,…,K Allan方差定义为: 2.Allan方差法最大的优点是可以简便细化分离、辨识光学陀螺的各项误差,同时确定各个误差项对总误差的贡献。 估计方差时,与 实验所用的陀螺仪类型和数据类型获取的环境有关,实验数据中可能存在各种成分的随机噪声。若各项噪声是统计独立的,则Allan方差可以表示成各类误差的平方和,如: 观察之后发现除去正弦噪声含有正弦分量外,其他都可以表示为(n=-2,-1,0,1,2)。因此上式可以简化为: 由于方差一般较小,拟合标准差则可以提高拟合精度,上式可以近似为: 其中=m0,0为采样时间 在最小均方意义上,拟合函数可以求出An,再通过以下计算可以得到量化误差Q,角度随机游走N,零偏不稳定性B,角速率随机游走K和速率斜坡R的估计值,为 R=A23600*sqrt(2); K=A160*sqrt(3); B=A00.6643; N=A-10; Q=A-2qrt(3); 三:实验过程(解算和数据拟合) 选取数据文件:从计算机得到的数据文件中抽取一部分,数据特点是除陀螺仪绕Z轴输出变化以外都比较稳定,可以由Origin软件的图表显示。如图: 选取其中一部分之后得到的z轴角速率输出曲线(50到-50之间变化) 源程序 将取得的数据分为5组,得到不同的总采样时间,再有Allan方差公式得到不同的方差,通过列出五个方程的方程组解出五个系数值。 程序为: %dynamic test about z-rotate %by allan variance analysis clear clc format long e [t,w_x,w_y,w_z,a_x,a_y,a_z,b]=textread(E:\allan\z-rotate.txt); N=size(t,1); sigma_squa=zeros(5,1); w=zeros(5,5);

文档评论(0)

cuotian + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档