论文部分内容阅读
重力测量的目的在于研究地球的重力场,重力场的研究帮助人们了解地球形状、内部物质分布、内部构造等,对导航、军事、地震研究和资源勘探等方面具有重要意义,获取高精度的绝对重力值是重力测量系统的关键。目前能在野外精确测量重力加速度的是落体式激光干涉绝对重力仪。其测量方法一般是通过真空腔内的自由落体下落产生激光干涉条纹,被光电转换器转化成电信号后再被高速数字化仪转化成数字信号,通过处理可以得到落体的重组轨迹-时间距离坐标,最后通过算法拟合得到高精度g值,采用合适的拟合算法对保证获取高精度重力场信息至关重要。最小二乘法是曲线拟合和参数估计的常用方法之一,由于最小二乘算法简单,计算量少,收敛性能好,被广泛应用在绝对重力仪数据拟合中,但针对实际测量中不断出现的新问题,仍需要在常规最小二乘的基础上对算法进行探究和改进。本文考虑到绝对重力干涉信号的频率变化特性,提出采用局部最小二乘算法进行处理,但是针对绝对重力仪数据段中存在异常点等方面,常规最小二乘和局部最小二乘处理得到的估计值将受到严重的歪曲甚至不能接受,因此提出采用抗差最小二乘处理,然而抗差最小二乘算法涉及到了高维矩阵乘积以及逆矩阵的求解,为了进一步提高算法的可靠性,从参数估计求解矩阵的形式和针对绝对重力仪数据的特点出发,提出了加权递推最小二乘算法,它将复杂庞大的计算过程转化为简单过程的多次重复,避免高维逆矩阵乘积以及逆矩阵求解,代码实现容易且编程量少。本文采用MATLAB软件产生模拟数据,分别采用常规最小二乘、局部最小二乘、抗差最小二乘和加权递推最小二乘对模拟数据进行参数估计,通过估计值与设定值的偏差、标准差及相应的处理时间等方面对它们的估计结果进行比较分析,分析结果表明加权递推算法在绝对重力仪数据拟合中具有更好的效果,并且稳定性高,拟合精度好,在异常数据的抑制上也具有独特的优势。然后分别使用上述四种算法对实际绝对重力仪数据进行处理后,通过每组平均值、标准差以及每组处理时间等方面进行对比分析,分析结果表明,针对实际数据出现的不确定性因素,加权递推算法有着其他三种算法不具备的抗干扰能力,这与模拟数据计算分析的结果基本一致。加权递推最小二乘应用到自由落体式绝对重力仪的数据处理中,极大地抑制了异常点对重力加速度计算结果的影响,提高了重力值的准确度和稳定性,相对于常规最小二乘、局部最小二乘和抗差最小二乘,加权递推最小二乘反映出更好的处理结果,鲁棒性能良好,更适合于绝对重力仪数据处理中。