共查询到19条相似文献，搜索用时 156 毫秒
1.
扩展高阶和超高阶重力场模型的构制与应用的数值稳定性取决于超高阶次缔合勒让德函数的计算方法。文中详细介绍了现有的多种缔合勒让德函数的递推计算方法:标准前向列推法、标准前向行推法、跨阶次递推法和Belikov列推法。从计算速度、计算精度和计算溢出问题3个角度分析比较了阶次高至2 160阶的各种方法的优劣。通过数值试验证明,Belikov列推法和跨阶次递推法是计算超高阶次缔合勒让德函数较优的方法,而其他几种方法不能用于超高阶次缔合勒让德函数的计算。文中结论为超高阶次球谐综合与球谐分析的数值计算提供了可靠的依据。 相似文献
2.
扩展高阶和超高阶重力场模型的构制与应用的数值稳定性取决于超高阶次缔合勒让德函数的计算方法.文中详细介绍了现有的多种缔合勒让德函数的递推计算方法:标准前向列推法、标准前向行推法、跨阶次递推法和Belikov列推法.从计算速度、计算精度和计算溢出问题3个角度分析比较了阶次高至2 160阶的各种方法的优劣.通过数值试验证明,Belikov列推法和跨阶次递推法是计算超高阶次缔合勒让德函数较优的方法,而其他几种方法不能用于超高阶次缔合勒让德函数的计算.文中结论为超高阶次球谐综合与球谐分析的数值计算提供了可靠的依据. 相似文献
3.
4.
首先将级数式中缔合勒让德函数的计算转化为一般勒让德函数的导数计算，并给出相应的递推公式；其次将球谐函数级数式转化为仅含带谐项的级数式，并给出相应系数的计算公式。经实例计算，用机时间节省约90％。这对于应用部门快速计算重力场的有关数据将十分有利。 相似文献
5.
6.
7.
8.
9.
根据经典的球谐函数方法,为满足正交化要求,观测数据需要覆盖整个球面,而对于地表局部测量数据,则无法应用球谐方法解算重力场模型。针对此问题,采用Slepian局部谱分析方法解算中国大陆范围内的实测重力场变化数据,并以GOCE卫星球谐函数解作为已知模型,评估由于实际陆地重力测点的非均匀分布对球谐函数解的误差影响。通过计算多个阶次中国大陆局部范围的Slepian基函数分布;采用GOCE卫星获得重力场模型的前72阶球谐系数作为已知结果,评价实际测点非均匀分布的解算有效性,并针对中国大陆地区采用Slepian基函数进行解算,通过模型对比选择最优截段项数;针对2005—2008年中国大陆地区流动重力测量获得的重力场变化信号进行解算,获得了72阶重力场变化模型。 相似文献
10.
广义球谐函数及其在梯度边值问题的应用 总被引：1，自引：0，他引：1
本文首先指出球谐函数不是球面上唯一的正交函数系,接着证明了体谐函烽的一阶,二阶水平导数也是球面上的正交函数系,最后利用球面上的正交函数系研究了物理大地测量中的二阶,二阶梯度边值问题。 相似文献
11.
Cheinway Hwang 《Journal of Geodesy》1995,70(1-2):110-116
The product of two associated Legendre functions can be represented by a finite series in associated Legendre functions with unique coefficients. In this study a method is proposed to compute the coefficients in this product-sum formula. The method is of recursive nature and is based on the straightforward polynomial form of the associated Legendre function's factor. The method is verified through the computation of integrals of products of two associated Legendre functions over a given interval and the computation of integrals of products of two Legendre polynomials over [0,1]. These coefficients are basically constant and can be used in any future related applications. A table containing the coefficients up to degree 5 is given for ready reference. 相似文献
12.
Toshio Fukushima 《Journal of Geodesy》2012,86(9):745-754
The existing methods to compute the definite integral of associated Legendre function (ALF) with respect to the argument suffer from a loss of significant figures independently of the latitude. This is caused by the subtraction of similar quantities in the additional term of their recurrence formulas, especially the finite difference of their values between two endpoints of the integration interval. In order to resolve the problem, we develop a recursive algorithm to compute their finite difference. Also, we modify the algorithm to evaluate their definite integrals assuming that their values at one endpoint are known. We numerically confirm a significant increase in computing precision of the integral by the new method. When the interval is one arc minute, for example, the gain amounts to 2–4 digits for the degree of harmonics in the range 2 ≤ n ≤ 2,048. This improvement in precision is achieved at a negligible increase in CPU time, say less than 5%. 相似文献
13.
本文所采用的基于输入－输出系统论的谱方法在计算结果的精度上与最小二乘配置方法相当,却很容易用于异性场的计算。用该谱方法对卫星测高及海洋重力资料进行组合求解重力场量（大地水准面差距和重力异常）,其误差估计结果表明各向异性场的计算精度优于各向同性场的精度。 相似文献
14.
万晓云 《武汉大学学报(信息科学版)》2011,36(12):1486-1489
传统的引力场梯度计算公式在两极附近存在奇异性,需要换用其他的非奇异计算公式。从奇异性产生的原因入手,并结合勒让德函数的有关性质,推导了一组新的计算公式。实际计算验证了该公式的正确性和有效性。 相似文献
15.
Improved convergence rates for the truncation error in gravimetric geoid determination 总被引：2，自引：2，他引：0
When Stokes's integral is used over a spherical cap to compute a gravimetric estimate of the geoid, a truncation error results
due to the neglect of gravity data over the remainder of the Earth. Associated with the truncation error is an error kernel
defined over these two complementary regions. An important observation is that the rate of decay of the coefficients of the
series expansion for the truncation error in terms of Legendre polynomials is determined by the smoothness properties of the
error kernel. Previously published deterministic modifications of Stokes's integration kernel involve either a discontinuity
in the error kernel or its first derivative at the spherical cap radius. These kernels are generalised and extended by constructing
error kernels whose derivatives at the spherical cap radius are continuous up to an arbitrary order. This construction is
achieved by smoothly continuing the error kernel function into the spherical cap using a suitable degree polynomial. Accordingly,
an improved rate of convergence of the spectral series representation of the truncation error is obtained.
Received: 21 April 1998 / Accepted: 4 October 1999 相似文献
16.
在局部区域似大地水准面的确定中,应用数学模型来求取高程异常,是工程中常用的方法。对正双曲面、倒双曲面、三次曲面三种不同核函数的多面函数拟合法最佳光滑因子的选取问题进行了比较研究。通过对地形复杂矿区工程实例分析,表明了三种方法在地形较复杂区域似大地水准面确定中的可行性。 相似文献
17.
The integral formulas of the associated Legendre functions 总被引：1，自引：0，他引：1
A new kind of integral formulas for ${\bar{P}_{n,m} (x)}$ is derived from the addition theorem about the Legendre Functions when n ? m is an even number. Based on the newly introduced integral formulas, the fully normalized associated Legendre functions can be directly computed without using any recursion methods that currently are often used in the computations. In addition, some arithmetic examples are computed with the increasing degree recursion and the integral methods introduced in the paper respectively, in order to compare the precisions and run-times of these two methods in computing the fully normalized associated Legendre functions. The results indicate that the precisions of the integral methods are almost consistent for variant x in computing ${\bar{P}_{n,m} (x)}$ , i.e., the precisions are independent of the choice of x on the interval [0,1]. In contrast, the precisions of the increasing degree recursion change with different values on the interval [0,1], particularly, when x tends to 1, the errors of computing ${\bar{P}_{n,m} (x)}$ by the increasing degree recursion become unacceptable when the degree becomes larger and larger. On the other hand, the integral methods cost more run-time than the increasing degree recursion. Hence, it is suggested that combinations of the integral method and the increasing degree recursion can be adopted, that is, the integral methods can be used as a replacement for the recursive initials when the recursion method become divergent. 相似文献
18.
本文推导的椭球谐系数和球谐系数相互之间转换关系的核心思想是在ε~2量级下利用Legendre函数的正交性,从球谐系数求解的积分表示出发,将积分中的椭球坐标变量与球坐标变量相互转换,从而得出椭球谐系数与球谐系数之间的转换关系。本文导出的转换关系有以下优点:①对于第二类Legendre函数的计算采用Laurent级数表示,使计算第二类Legendre函数更为简单;②保留了ε~2量级下,导出的转换关系相比文献[2]的形式更简单,满足物理大地测量边值问题线性化的要求;③顾及了余纬和归化余纬的区别。 相似文献
19.
Toshio Fukushima 《Journal of Geodesy》2018,92(12):1371-1386
We developed an accurate method to compute the gravitational field of a tesseroid. The method numerically integrates a surface integral representation of the gravitational potential of the tesseroid by conditionally splitting its line integration intervals and by using the double exponential quadrature rule. Then, it evaluates the gravitational acceleration vector and the gravity gradient tensor by numerically differentiating the numerically integrated potential. The numerical differentiation is conducted by appropriately switching the central and the single-sided second-order difference formulas with a suitable choice of the test argument displacement. If necessary, the new method is extended to the case of a general tesseroid with the variable density profile, the variable surface height functions, and/or the variable intervals in longitude or in latitude. The new method is capable of computing the gravitational field of the tesseroid independently on the location of the evaluation point, namely whether outside, near the surface of, on the surface of, or inside the tesseroid. The achievable precision is 14–15 digits for the potential, 9–11 digits for the acceleration vector, and 6–8 digits for the gradient tensor in the double precision environment. The correct digits are roughly doubled if employing the quadruple precision computation. The new method provides a reliable procedure to compute the topographic gravitational field, especially that near, on, and below the surface. Also, it could potentially serve as a sure reference to complement and elaborate the existing approaches using the Gauss–Legendre quadrature or other standard methods of numerical integration. 相似文献