首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 187 毫秒
1.
局部区域模型垂线偏差快速算法   总被引:3,自引:0,他引:3  
针对利用地球重力场模型计算大范围、高分辨率局部区域模型垂线偏差速度慢这一问题,该文在对位系数模型泰勒级数展开的基础上,提出了采用向量运算、混合编程、预先计算缔合勒让德函数Belikov递推方法中与经纬度无关量的方法;最后与严密公式计算结果作对比。实验结果表明,该方法不仅可以满足精度要求,而且可以有效提高计算速度。  相似文献   

2.
远程飞行器和弹道导弹的精确控制与导航会受到地球外部扰动引力场的影响,实践中常利用地球重力场位系数模型计算扰动引力矢量,但其计算效率随着所用模型阶次的升高而显著降低。针对这一问题,文中提出利用OpenMP并行算法快速计算轨迹点扰动引力矢量,在不损失精度的前提下有效提高了计算效率。通过数值实验可知,所提方法可以提高扰动引力矢量的计算效率,加速比达到6倍。这也为重力场元快速计算提供了参考方案。  相似文献   

3.
李伟  刘站科  蒋涛 《测绘科学》2014,39(12):20-23,19
随着地球重力场模型逐渐向超高阶次发展,如何利用超高阶重力场模型快速、精确地计算地球重力场参数成为重要问题.文章给出了利用位系数计算地球外部空间任意点位处重力场参数的数学模型和详细算法,采用比例因子法实现2190阶次缔合勒让德函数的稳定递推计算;并基于VS2008平台开发了重力场参数快速赋值软件.最后利用EGM2008重力场模型和实测GPS水准数据验证了算法、软件的可靠性与效率.  相似文献   

4.
各向异性组合滤波法反演陆地水储量变化   总被引:2,自引:1,他引:1  
超能芳  王正涛  孙健 《测绘学报》2015,44(2):174-182
地球时变重力场模型反演陆地水储量变化已为全球气候变化研究作出巨大贡献,考虑到时变重力场模型球谐系数中存在相关性,其高阶次项具有较大的误差,需采用最优的滤波方法进行空间平滑。本文提出一种新的各向异性组合滤波方法,其基本思想是将改进的高斯滤波法与均方根(root mean square,RMS)滤波法组合,即对球谐系数的低阶次采用改进的高斯滤波法,而高阶次采用RMS滤波法。首先分析了最新的GRACE RL05系列时变重力场模型系数误差特性,基于全球水储量变化反演结果,分析比较了高斯滤波、改进的高斯滤波、RMS滤波和DDK滤波与本文提出的组合滤波法的有效性及精度,并利用模型结果进行了验证,计算结果表明,组合滤波法的中误差最小。研究结果表明,本文提出的组合法相比于先前的滤波方法,可有效地过滤高阶次的噪声,消除南北条带误差,同时减少信号泄漏,提高信噪比,保留更多有效的地球物理信号,进而提高反演精度。  相似文献   

5.
利用最小二乘直接法反演卫星重力场模型的MPI并行算法   总被引:2,自引:0,他引:2  
周浩  罗志才  钟波  陆飚 《测绘学报》2015,44(8):833-839
针对海量卫星重力数据反演高阶次地球重力场模型的密集型计算任务与高内存耗用问题,基于MPI实现了最小二乘直接法恢复高阶次位系数的并行算法。引入并行读写、分块存储与分块计算等方式完成了设计矩阵的构建、法方程的形成与求解等密集型计算任务的并行算法,数值计算结果表明三者的并行相对效率峰值可分别达到95%、68%、63%。利用GOCE轨道跟踪和径向扰动重力梯度数据(共518 400个历元)分别反演了120、240阶次地球重力场模型,计算时间仅为40 min、7 h,内存耗用峰值仅为290 MB、1.57 GB;采用与GOCE同等噪声水平的观测数据恢复的重力场模型精度与GOCE已发布模型的解算精度相一致,联合GRACE和GOCE的解算模型能够实现二者独立信息的频谱互补,表明本文方法可高效稳定地恢复高阶次地球重力场模型。  相似文献   

6.
本文探讨了基于能量守恒方法利用CHAMP卫星精密星历和加速度数据恢复地球重力场模型的原理和方法。给出了地心惯性系下顾及地球自转和非保守力能量损耗的能量守恒方程,并且对日、月摄动位与引潮力附加位的计算方法作了相应的分析,同时介绍了加速度数据的处理方法。基于能量守恒方法,利用2002年1-2月、7-8月和11-12月三个不同时期共180天的CHAMP卫星精密星历和加速度数据恢复了三组50阶次的地球重力场模型GFM01、GFM02和GFM03,并将这些模型与EGM 96重力场模型和GFZ公布的EIGEN-CG01C重力场模型进行比较。结果表明:能量守恒方法恢复的GFM系列模型与EGM 96重力场模型及EIGEN-CG01C重力场模型在低阶位系数上均有较好的一致,但与EIGEN-CG01C模型有更好的一致。这说明了CHAMP卫星对地球中、长波重力场的敏感性,也说明了能量守恒方法恢复低阶地球重力场位系数的有效性。  相似文献   

7.
在利用卫星重力梯度数据求解地球重力场模型位系数时,为进一步提高调和分析法的计算效率,推导出了利用5种组合{Tzz}、{Txx+Tyy}、{Tzz-Txx-Tyy}、{Txz,Tyz}、{Txx,Tyy,Txy}构建重力场模型位系数的向量化调和分析法和向量化调和分析的FFT算法.实验表明,利用30′×30′格网Tzz数据最快可在1 s内反演300阶次重力场模型.与传统FFT算法相比,向量化调和分析的FFT算法计算效率最高可提高16倍.  相似文献   

8.
扰动位的综合确定   总被引:2,自引:1,他引:2  
利用地球重力场任意一种有关信息都可以描述地球重力场一定的情况。根据卫星轨道摄动观测求定的引力位球谐系数只能表示地球重力场的长波分量。大地水准面起伏是地球扰动引力场越来越丰富的有用信息,但目前用其计算的引力位系数也只是在低阶较准确。重力异常、垂线偏差、单层密度和纯重力异常都利于求定高阶位系数,其中与大地水准面起伏有关的量,如纯重力异常和单层密度,用它们计算位系数等于联合应用大地水准面和重力异常,故用其计算的位系数在低阶次精度也较好。重力异常垂直梯度是描述扰动引力场细部最有利的信息。本文给出利用各种类型观测资料计算位系数精度估计式,提出综合利用各种资料求定位系数依资料类型的谱特性赋权的方法。  相似文献   

9.
基于能量守恒原理,建立了SST-ll星间距离变率观测噪声谱与重力场误差谱的关系,以GRACE相关指标模拟分析了卫星间距、卫星高度和距离变率精度对恢复地球重力场的影响.结果表明,增大卫星间距可提高恢复低阶次位系数的精度,卫星间距超过500 km对提高恢复重力场精度的作用已不明显;降低轨道高度可提高恢复高阶次位系数的精度,卫星高度每降低100 km,恢复位系数的有效阶次提高20阶以上;提高星间距离变率精度可大幅度提高恢复重力场的精度,距离变率精度每提高一个量级,恢复位系数的有效阶次提高约28阶.将模拟结果与GGM02S和EIGEN-GRACE02S模型进行比较,初步验证了本文方法的可行性.  相似文献   

10.
用GRACE星间速度恢复地球重力场   总被引:9,自引:2,他引:7  
本文首先给出了用星间速度恢复地球重力场的数学模型,然后用GRACE卫星30天的星间速度观测值计算了一个100阶的地球重力场模型DQM2006S2。为了对这一模型的精度进行评述,将它与EGM 96,EIGEN-CHAMP03S和GGM01S3个地球重力场模型作了比较,并用这一模型计算了高程异常与GPS/水准实际观测值进行了比较,结果表明:DQM2006S2模型精度优于EGM 96和EIGEN-CHAMP03S模型精度,但是不及GGM01S模型精度。精度不及GGM01S的原因是GGM01S模型使用了111天的星间速度数据,其数据量约为DQM2006S2模型使用数据量的4倍。  相似文献   

11.
在评估重力场模型计算空间扰动引力精度时,对模型截断误差常采用阶方差方法。文中将6种经典的重力异常阶方差模型与现有超高阶重力场模型的阶方差进行比较,TSD模型与重力场模型的差值最小。根据重力异常阶方差模型TSD,文中分析不同高度、不同阶次利用重力场模型计算空中扰动引力时截断误差的影响。实验结果表明:36阶模型截断误差最大径向和水平方向分别为26.455 1mGal、25.946 3mGal;360阶模型截断误差最大径向和水平方向分别为9.969 0mGal、9.960 9 mGal;2160阶模型截断误差最大径向和水平方向分别为2.538 5 mGal、2.538 1mGal;2160阶模型计算空中扰动引力时,即使在低空附近,截断误差在2.5mGal以内,计算高度超过5km,截断误差可以忽略;超过400km的高度,都可以用36阶模型计算,截断误差在1mGal以内。  相似文献   

12.
四种改进积分法的低空扰动引力计算   总被引:1,自引:0,他引:1  
针对Stokes积分方法计算扰动引力中计算点从空中趋近地面时存在积分奇异和不连续的问题,该文提出了去中央奇异点法、奇异点积分值修正法、中央格网加密算法和改进积分式法4种改进Stokes积分的计算公式,并进行了实验计算。计算结果表明:近地空间范围内,4种改进算法都能在一定程度上改进原始积分的奇异性问题;相同条件下,奇异点积分值修正法和改进积分式法计算精度最高,适宜于低空计算;改进积分式法通过理论推导,得到了从球外部到球面统一、连续且无奇异的改进Stokes积分公式,理论严谨。  相似文献   

13.
根据地球重力场参数的计算模型,借助全球地球重力场模型EGM2008,利用Microsoft Visual Studio 2010面向对象的功能,设计和开发了可以计算任意点重力异常、高程异常、垂线偏差分量的计算程序,以德国地学中心(GFZ)研制的相关软件对程序的正确性进行了验证,以吉林省各城市的经纬度数据计算了对应的重力异常、高程异常、垂线偏差分量。  相似文献   

14.
马志伟  陆洋  涂弋  朱传东  郗慧 《测绘学报》2016,45(9):1019-1027
多种类型高分辨率重力场数据的不断增加,使得在局部范围内精化重力场模型成为了可能。本文采用Abel-Poisson核将重力场量表示成有限个径向基函数线性求和的形式,对局部区域的多种重力场数据进行联合建模。为了提高运算速度,运用了基于自适应精化格网算法的最小均方根误差准则(RMS)来求解径向基函数平均带宽。以南海核心地区为例,联合两种不同类型、不同分辨率的重力场资料(大地水准面起伏6'×6'、重力异常2'×2'),构建了局部区域高分辨率的重力场模型。所建模型表示的重力场参量达到了2'×2'的分辨率,对原始的重力异常数据(2'×2')拟合的符合程度达到±0.8×10-5m/s2。结果表明,利用径向基函数方法进行局部重力场建模,避免了球谐函数建模收敛慢的问题,有效提高了模型表示重力场的分辨率。  相似文献   

15.
利用EGM2008地球重力场模型可以将地球重力场解析化,由此我们可以十分方便快速地表示和计算重力异常、高程异常、垂线偏差等其它重力场参数。本文运用Surfer软件将EGM2008重力场模型可视化,获取了全球和局部地区重力场相关参数的等值线图、线框图、3D表面图、影像图。  相似文献   

16.
地球重力场模型正演计算软件系统的设计   总被引:1,自引:0,他引:1  
本文介绍了自主开发的地球重力场模型正演(计算高程异常、重力异常、垂线偏差等)计算软件—EGMCAL,给出了该软件的整体设计框架和功能。利用GFZ(德国地学中心)研制的相关软件验证了该软件的正确性,并利用EGMCAL对我国东部某大桥高程控制网数据进行了处理,得出三种精度最高的模型在此区域的计算结果。  相似文献   

17.
由CHAMP星载GPS相位双差数据解算地球引力场模型   总被引:1,自引:0,他引:1  
利用7d的CHAMP星载GPS相位观测数据和48个IGS跟踪站的观测数据,构造星地双差相位观测量,进行GPS数据预处理;利用Cowell Ⅱ数值法进行轨道积分和分块Bayes最小二乘参数估计,解算了地球引力场位系数。该模型与EGM96相比(70阶次),大地水准面起伏差异最大为2.872m,差弄精度为0.522m,平均差异为-0.003m,这说明本文解算的地球重力场模型与EGM96没有系统性差异。  相似文献   

18.
最小二乘配置最初是在组合各种资料来研究地球形状与重力场的一种数学方法,目前最小二乘配置已经在测绘数据处理中得到广泛应用。本文首先分析了目前采用的最小二乘配置法解算方法,在讨论了矩阵的奇异值分解(Singular Value Decomposition,SVD)方法的基础上,推导得出了矩阵SVD分解与广义逆矩阵的关系,得出了可以直接利用SVD分解求解矩阵的Moore-Penrose广义逆,并推导了应用SVD分解求解最小二乘配置的估值计算公式和精度估算公式,最后通过重力异常实例进行了计算,得出矩阵的SVD分解用于最小二乘配置解算的正确性和可行性,为最小二乘配置的求解提供了一种新方法。  相似文献   

19.
2009年GOCE卫星升空以后,卫星重力梯度数据参与解算的GOCE系列重力场模型已有多家研究机构相继公布。本文分别采用青藏地区的GPS/水准和重力异常实测数据对GOCE重力场模型进行了外部测试,并在重力异常验证过程中引入了一种新的滤波方法,验证结果表明在青藏地区GOCE重力场模型相比其它系列模型的优势在于中波段。同时,探讨了GOCE重力场模型与其他系列模型在青藏地区主要差异值的空间分布以及首次利用统计分析方法找出模型之间主要差异值的阶次分布,得出如下结论:模型之间的较大差异值在空间水平方向上主要分布在喜马拉雅山脉、天山等地形起伏较大的区域,在垂直方向上主要集中在岩石圈。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号