球谐函数变换快速计算扰动引力

在球谐函数变换基础上,利用新极下轨道的特殊性,在新坐标下引入Clenshaw求和计算轨道扰动引力。从理论上对比分析了传统方法、球谐函数变换方法和改进方法的计算速度和存储模型需要的物理空间。模拟试验分别采用3种方法计算了一段轨道的扰动引力,试验结果表明,改进的球谐函数变换方法比传统球谐函数变换方法计算速度可提高100倍,数据存储量仅占传统方法的3%。     

广义球谐函数定积分计算方法的改进

运用球谐函数定积分的基本递推公式,推导了在重力场球谐综合与球谐分析中出现的广义球谐函数定积分的计算公式;给出了其适用于超高阶次的改良型递推公式.数值试验表明,该改良公式具有较高的计算精度和计算速度,解决了超高阶次广义球谐函数定积分计算的溢出问题,拓展了这类定积分的计算公式.他们的数值实现为利用位模型计算高分辨率扰动重力场元格网平均值、重力场球谐综合分析等奠定了基础.     

广义球谐函数定积分计算方法的改进

运用球谐函数定积分的基本递推公式，推导了在重力场球谐综合与球谐分析中出现的广义球谐函数定积分的计算公式；给出了其适用于超高阶次的改良型递推公式。数值试验表明，该改良公式具有较高的计算精度和计算速度，解决了超高阶次广义球谐函数定积分计算的溢出问题，拓展了这类定积分的计算公式。他们的数值实现为利用位模型计算高分辨率扰动重力场元格网平均值、重力场球谐综合分析等奠定了基础。     

边界数据稀疏地区似大地水准面模型的建立

研究了斯托克斯移去恢复法和框架约束球谐函数展开模型,分析了数据稀疏、质量不高地区的计算情况,实现了框架约束球谐函数模型精度优于常规斯托克斯移去恢复法。     

一类球谐函数与三角函数乘积积分的计算

本文根据球谐函数的跨次递推公式和三角函数的性质,详细推导了在重力梯度调和分析中出现的一类球谐函数积分的跨次递推公式和递推初始值的计算公式。数值试验表明,球谐函数跨次递推算法具有快速、稳定的优点。该类积分的跨次递推实现,为卫星重力梯度调和分析奠定了算法基础。     

地球引力加速度及其偏导数的计算

本文详细推导了地球引力加速度及其偏导数的计算公式,给出了一组用来计算正规化球谐函数的递推关系。这组递推关系优于一般文献上所给出的用来计算非正规化球谐函数的递推关系,它避免了由于不同阶、级非正规化球谐函数数值相差过大而引起的计算精度的损失。在卫星精密轨道确定中已经证明这种算法是非常有效的。     

基于K-means聚类分析的球冠谐函数拟合高程异常方法

利用球冠谐函数拟合高程异常曲面时,GPS/水准数据的分布和选取,对最后的拟合精度有较大影响。本文提出了一种基于K-Means聚类的球冠谐函数拟合方法,即利用K-Means聚类分析方法对高程异常数据点进行有效分类和选择,再结合球冠谐函数的逼近方法,提高了高程异常曲面的拟合精度。通过与实测高程异常数据进行比较,基于K-Means聚类分析的球冠谐函数拟合方法优于多项式拟合法、移动二次曲面拟合法、多面函数拟合法和球冠谐函数拟合法。     

全球位系数模型构建高精度局部重力场的Clenshaw求和

给出了勒让德函数4种模式的Clenshaw求和公式,在此基础上得到由全球位系数模型计算重力场参数的二阶Clenshaw求和公式,分析了EGM96与WDM94应用于某试验测区局部重力场的精度,为该测区在构建高精度局部重力场时初始全球位系数模型的选取提供了借鉴。     

基于球谐函数模型的电离层预报

据球谐函数模型系数的特点,采用ARMA(p,q)模型对球谐函数模型系数进行预报,由球谐函数模型计算电离层VTEC。提出了针对某一时刻球谐系数进行预报的方法,相比传统按照时刻顺序的预报,预报时间大大延长,预报精度也有所提高。试验结果表明,相比中高纬度地区,低纬度地区预报精度偏低,同时,一天中不同时刻预报结果有所差别,前半天的预报效果明显好于后半天。     

外空扰动引力计算方法的内在联系

利用算子运算的方法,对求定外空扰动引力的Bjerhammar方法、虚拟点质量方法以及虚拟单层密度方法进行了系统地归纳;利用Possion核的再生特性及球谐函数的正交性质,给出了它们之间的换算关系;最后从两个方面分析了实用计算中各种方法的优劣.     

Fourier transform summation of Legendre series and D-functions

The relation between D- and d-functions, spherical harmonic functions and Legendre functions is reviewed. Dmatrices and irreducible representations of the rotation group O(3) and SU(2) group are briefly reviewed. Two new recursive methods for calculations of D-matrices are presented. Legendre functions are evaluated as part of this scheme. Vector spherical harmonics in the form af generalized spherical harmonics are also included as well as derivatives of the spherical harmonics. The special dmatrices evaluated for argument equal to/2 offer a simple method of calculating the Fourier coefficients of Legendre functions, derivatives of Legendre functions and vector spherical harmonics. Summation of a Legendre series or a full synthesis on the unit sphere of a field can then be performed by transforming the spherical harmonic coefficients to Fourier coefficients and making the summation by an inverse FFT (Fast Fourier Transform). The procedure is general and can also be applied to evaluate derivatives of a field and components of vector and tensor fields.     

Methods of harmonic synthesis for global geopotential models and their first-, second- and third-order gradients

Four widely used algorithms for the computation of the Earth’s gravitational potential and its first-, second- and third-order gradients are examined: the traditional increasing degree recursion in associated Legendre functions and its variant based on the Clenshaw summation, plus the methods of Pines and Cunningham–Metris, which are free from the singularities that distinguish the first two methods at the geographic poles. All four methods are reorganized with the lumped coefficients approach, which in the cases of Pines and Cunningham–Metris requires a complete revision of the algorithms. The characteristics of the four methods are studied and described, and numerical tests are performed to assess and compare their precision, accuracy, and efficiency. In general the performance levels of all four codes exhibit large improvements over previously published versions. From the point of view of numerical precision, away from the geographic poles Clenshaw and Legendre offer an overall better quality. Furthermore, Pines and Cunningham–Metris are affected by an intrinsic loss of precision at the equator and suffer from additional deterioration when the gravity gradients components are rotated into the East-North-Up topocentric reference system. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

GPS测站多路径效应建模研究

GPS伪距观测量由于不存在周跳和整周模糊度的问题,可用来实现动态导航定位,但受到多路径效应影响较为严重,定位精度较低。本文以BJFS站为例,计算了伪距多路径效应,首次大量采用单历元数据,以球谐函数的形式建立了伪距多路径效应、高度角和方位角间的函数模型。建模过程中用截断奇异值分解方法解决了单历元数据引起的的病态问题,求得了模型的最小二乘解。在伪距观测数据中加入模型计算的伪距多路径效应改正,用Bernese软件的动态单点定位功能进行伪距单点定位。对BJFS站的伪距动态定位结果进行统计分析,结果表明单点定位精度得到了改善。     

Spatially restricted data distributions on the sphere: the method of orthonormalized functions and applications

 In many geoscientific applications data are irregularly distributed and not globally available, e.g. caps around the poles which are uncovered due to non-polar satellite orbits, or signals being defined solely on bounded regions on the globe. Starting from a sequence of base functions with global support, which in the present case is composed of spherical harmonics being initially non-orthogonal on a bounded subdomain, a set of functions is generated that constitutes an orthonormal basis. Different approaches to realize this transformation are studied and compared with respect to numerical stability and computational effort, and the corresponding effects on the coefficient recovery are investigated. A number of synthetic tests demonstrate the applicability, the benefit, but also the limitations, of this method. Received: 24 March 2000 / Accepted: 9 October 2000     

全张量重力梯度数据的谱表示方法

在文献「１」的基础上,进一步研究全张量重力梯度数据的全局和局部分量的广义球谐谱表示和轨道根数据表示,并给出了广义球谐函数与球谐函数这间的关系,从理论上得到了全张量重力梯度数据的描述方法和由全张量重力梯度网格数据恢复全球重力位谱系数的基本公式,本文对全张量重力梯度数据的谱表示和谱分析所做的工作,对由重力梯度张量的部分分量恢复全球重力位普系数有一定的参考价值。     

The compatibility conditions in altimetry–gravimetry boundary value problems

The role of compatibility conditions for altimetry–gravimetry boundary value problems (AGBVPs) is quite complicated by way of: guaranteeing the existence of the solution; smoothing the data (boundary conditions) at the boundary; and providing a higher level of regularization and smoothness to the solution. Pseudo-differential operators (PDOs) can be applied not only for the reformulation of AGBVPs in a simple way, but also for imposing compatibility conditions at the coastline. It is shown here that spherical PDOs can be combined with the theory of spherical harmonics and spherical wavelets to obtain a numerical solution for AGBVPs with compatibility conditions along the coastline. The role of compatibility conditions, is summarized and emphasized and compatibility conditions are given in an explicit form based on the combined representation of functionals of the disturbing potential in terms of spherical harmonics, spherical PDOs and wavelets. Numerical aspects of the application of the compatibility conditions to a case-study area in Canada are discussed.