电子设计工程
Electronic Design Engineering
第 25卷
Vol.25
第 22期
No.22
2017年 11月
Nov. 2017
收稿日期:2016-09-21 稿件编号:201609188
作者简介:罗 棋(1990—),男,湖北天门人,硕士。研究方向:绝对重力仪落体控制。
绝对重力测量在计量学,地球物理和大地测量
有着重要的作用,其精度的提高对地震活动、地质勘
测都有不可估量的作用
[1]
。一般测量绝对重力加速
度 g 的方法是基于落体每下落半个波长距离产生一
个干涉条纹的规律,获得物体下落的位移信息,对比
铷原子钟给出的时间,得到时间距离对,拟合得到重
力加速度值
[2]
。
目前普遍采用最小二乘和非线性拟合时间距离
对
[3]
,与非线性相比,最小二乘原理简单,易于理解和
掌握,且在编程中易于实现,但是最小二乘算法适用
于近似服从正态分布的测量数据,参数拟合结果对
非正态分布的观测粗差极为敏感
[4]
,没有考虑观测误
差的不同对算法解的影响
[5]
。用双采样过零检测,时
刻数据将越来越密集,时间间隔越来越短,测量的误
差影 响将 会放 大,因此 需要 对最 小二 乘进行改进。
文献[4]提出了抗差估计,但是涉及到求解复杂多维
的逆矩阵。基于上述问题,为了进一步提供算法的
可行性和精度要求,提出加权递推形式,避免了复杂
逆矩阵的求解。文中基于文献[4]中的迭代算法,采
用文献[5]的加权阵,对拟合的时间距离对进行加权
处理,拟合误差方差越大的赋予的权重越小,误差方
差小的赋予的权重大,这也充分体现了加权递推最
小二乘算法的优越性
[5]
。
1 加权递推最小二乘算法
1.1 加权递推最小二乘
绝对重力仪研究的对象是做自由落体运动的
落体,针对落体运动,绝对重力仪数据符合二次多项
式曲线(包含重力梯度和光速有限性请参考文献[6]
中的算法,在拟合公式中增加正弦项和余弦项以仰
制激光器调制频率对输出的影响参考文献[4],将其
整理后依旧可以带入最小二乘拟合参数,并不影响
算 法 的 通 用 性),通 常 假 设 时 间 距 离 对 数 据 为
(t
i
,s
i
)(i = 1,2,⋯m)
,最小二乘的基本思想即选择一个
n 次多项式
s(t)
用一种使均方误差极小的方式来逼近
加权递推最小二乘在绝对重力仪数据处理中的应用
罗 棋,庞 聪,李查玮
(中国地震局地震研究所(地震大地测量重点实验室),湖北 武汉 430071)
摘要:针对绝对重力仪高精度的需求,在抗差参数估计的基础上提出了加权递推算法及其改进算
法,加权递推和改进算法改良抗差权重因子并且避免求复杂多维的逆矩阵。通过仿真表明,该算
法对含有白噪声或者异常数据都起到很好的仰制作用并提高了重力测量精度。
关键词:绝对重力仪;抗差估计;加权递推;最小二乘
中图分类号:P315 文献标识码:A 文章编号:1674-6236(2017)22-0065-04
Application of weighted recursive least square to data processing of absolute
gravimeter
LUO Qi,PANG Cong,LI Cha⁃wei
(Institute of Seismology,CEA,Wuhan 430071,China)
Abstract: In view of the high precision demand of absolute gravimeter,the weighted recursive algorithm
and its improved algorithm are proposed based on the robust parameter estimation,weighted recursive
algorithm and improved robust weighting factor and inverse matrix for complex multidimensional.
Simulation results show that the proposed algorithm has a good effect on the white noise or abnormal data
and improves the accuracy of gravity measurement.
Key words: absolute gravimeter;robust estimation;weighted recursive;least square
-- 65