共查询到20条相似文献，搜索用时 31 毫秒
1.
Computations of Fourier coefficients and related integrals of the associated Legendre functions with a new method along with their application to spherical harmonics analysis and synthesis are presented. The method incorporates a stable three-step recursion equation that can be processed separately for each colatitudinal Fourier wavenumber. Recursion equations for the zonal and sectorial modes are derived in explicit single-term formulas to provide accurate initial condition. Stable computations of the Fourier coefficients as well as the integrals needed for the projection of Legendre functions are demonstrated for the ultra-high degree of 10,800 corresponding to the resolution of one arcmin. Fourier coefficients, computed in double precision, are found to be accurate to 15 significant digits, indicating that the normalized error is close to the machine round-off error. The orthonormality, evaluated with Fourier coefficients and related integrals, is shown to be accurate to O(10?15) for degrees and orders up to 10,800. The Legendre function of degree 10,800 and order 5,000, synthesized from Fourier coefficients, is accurate to the machine round-off error. Further extension of the method to even higher degrees seems to be realizable without significant deterioration of accuracy. The Fourier series is applied to the projection of Legendre functions to the high-resolution global relief data of the National Geophysical Data Center of the National Oceanic and Atmospheric Administration, and the spherical harmonic degree variance (power spectrum) of global relief data is discussed. 相似文献
2.
M. K. Paul 《Journal of Geodesy》1978,52(3):177-190
Recurrence relations for the evaluation of the integrals of associated Legendre functions over an arbitrary interval within
(0°, 90°) have been derived which yield sufficiently accurate results throughout the entire range of their possible applications.
These recurrence relations have been used to compute integrals up to degree 100 and similar computations can be carried out
without any difficulty up to a degree as high as the memory in a computer permits. The computed values have been tested with
independent check formulae, also derived in this work; the corresponding relative errors never exceed 10−23 in magnitude.
Contribution from the Earth Physics Branch No. 719 相似文献
3.
Fourier transform summation of Legendre series and D-functions 总被引：4，自引：1，他引：3
T. Risbo 《Journal of Geodesy》1996,70(7):383-396
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. 相似文献
4.
Non-Singular Expressions for the Gravity Gradients in the Local North-Oriented and Orbital Reference Frames 总被引：4，自引：3，他引：1
The conventional expansions of the gravity gradients in the local north-oriented reference frame have a complicated form, depending on the first- and second-order derivatives of the associated Legendre functions of the colatitude and containing factors which tend to infinity when approaching the poles. In the present paper, the general term of each of these series is transformed to a product of a geopotential coefficient and a sum of several adjacent Legendre functions of the colatitude multiplied by a function of the longitude. These transformations are performed on the basis of relations between the Legendre functions and their derivatives published by Ilk (1983). The second-order geopotential derivatives corresponding to the local orbital reference frame are presented as linear functions of the north-oriented gravity gradients. The new expansions for the latter are substituted into these functions. As a result, the orbital derivatives are also presented as series depending on the geopotential coefficients multiplied by sums of the Legendre functions whose coefficients depend on the longitude and the satellite track azimuth at an observation point. The derived expansions of the observables can be applied for constructing a geopotential model from the GOCE mission data by the time-wise and space-wise approaches. The numerical experiments demonstrate the correctness of the analytical formulas.An erratum to this article can be found at 相似文献
5.
Spherical harmonic series, commonly used to represent the Earth’s gravitational field, are now routinely expanded to ultra-high
degree (> 2,000), where the computations of the associated Legendre functions exhibit extremely large ranges (thousands of
orders) of magnitudes with varying latitude. We show that in the degree-and-order domain, (ℓ,m), of these functions (with full ortho-normalization), their rather stable oscillatory behavior is distinctly separated from
a region of very strong attenuation by a simple linear relationship: , where θ is the polar angle. Derivatives and integrals of associated Legendre functions have these same characteristics.
This leads to an operational approach to the computation of spherical harmonic series, including derivatives and integrals
of such series, that neglects the numerically insignificant functions on the basis of the above empirical relationship and
obviates any concern about their broad range of magnitudes in the recursion formulas that are used to compute them. Tests
with a simulated gravitational field show that the errors in so doing can be made less than the data noise at all latitudes
and up to expansion degree of at least 10,800. Neglecting numerically insignificant terms in the spherical harmonic series
also offers a computational savings of at least one third. 相似文献
6.
7.
Toshio Fukushima 《Journal of Geodesy》2018,92(2):123-130
In order to accelerate the spherical harmonic synthesis and/or analysis of arbitrary function on the unit sphere, we developed a pair of procedures to transform between a truncated spherical harmonic expansion and the corresponding two-dimensional Fourier series. First, we obtained an analytic expression of the sine/cosine series coefficient of the \(4 \pi \) fully normalized associated Legendre function in terms of the rectangle values of the Wigner d function. Then, we elaborated the existing method to transform the coefficients of the surface spherical harmonic expansion to those of the double Fourier series so as to be capable with arbitrary high degree and order. Next, we created a new method to transform inversely a given double Fourier series to the corresponding surface spherical harmonic expansion. The key of the new method is a couple of new recurrence formulas to compute the inverse transformation coefficients: a decreasing-order, fixed-degree, and fixed-wavenumber three-term formula for general terms, and an increasing-degree-and-order and fixed-wavenumber two-term formula for diagonal terms. Meanwhile, the two seed values are analytically prepared. Both of the forward and inverse transformation procedures are confirmed to be sufficiently accurate and applicable to an extremely high degree/order/wavenumber as \(2^{30}\,{\approx }\,10^9\). The developed procedures will be useful not only in the synthesis and analysis of the spherical harmonic expansion of arbitrary high degree and order, but also in the evaluation of the derivatives and integrals of the spherical harmonic expansion. 相似文献
8.
S.J. Claessens 《Journal of Geodesy》2005,79(6-7):398-406
Several new relations among associated Legendre functions (ALFs) are derived, most of which relate a product of an ALF with trigonometric functions to a weighted summation over ALFs, where the weights only depend on the degree and order of the ALF. These relations are, for example, useful in applications such as the computation of geopotential coefficients and computation of ellipsoidal corrections in geoid modelling. The main relations are presented in both their unnormalised and fully normalised (4π-normalised) form. Several approaches to compute the weights involved are discussed, and it is shown that the relations can also be applied in the case of first- and second-order derivatives of ALFs, which may be of use in analysis of satellite gradiometry data. Finally, the derived relations are combined to provide new identities among ALFs, which contain no dependency on the colatitudinal coordinate other than that in the ALFs themselves. 相似文献
9.
球冠谐分析中非整阶Legendre函数的性质及其计算 总被引：9，自引：4，他引：5
局部重力场的谱方法是当前重力学的研究方向,该方法的核心问题是如何构造合适的谱函数以及如何对谱函数实施快速、有效的计算。当所研究的区域近似一个球冠时,迂冠谐函数是该区域对应的谱函数,它由非整阶勒让德（Ｌｅｇｅｎｄｒｅ）函数和三角函数组成,显然非整阶勒让德函数的构造和计算是研究球冠谐函数的关键。本文研究了非高小阶勒让德函数的性质和实用计算方法,包括如何对非整阶勒让德函数实施规格化处理。 相似文献
10.
Methods of harmonic synthesis for global geopotential models and their first-, second- and third-order gradients 总被引：2，自引：2，他引：0
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. 相似文献
11.
Efficient numerical computation of integrals defined on closed surfaces in ℝ3 with non-integrable point singularities that arise in physical geodesy is discussed. The method is based on the use of polar
coordinates and the definition of integrals with non-integrable point singularities as Hadamard finite part integrals. First
the behavior of singular integrals under smooth parameter transformations is studied, and then it is shown how they can be
reduced to absolutely integrable functions over domains in ℝ2. The correction terms that usually arise if the substitution rule is formally applied, in contrast to absolutely integrable
functions, are calculated. It is shown how to compute the regularized integrals efficiently, and, numerical efforts for various
orders of singularity are compared. Finally, efficient numerical integration methods are discussed for integrals of functions
that are defined as singular integrals, a task that typically arises in Galerkin boundary element methods.
Received: 15 April 1997 / Accepted: 7 May 1998 相似文献
12.
A function having some properties of a wavelet and being harmonic around a given point in R
3 is defined, and three models showing the local relationships between the disturbing density, the disturbing potential and
the disturbing gravity are established by using the function as the kernel function of the integrals in the models. The local
relationship has two meanings. One is that we can evaluate with a high accuracy the integrals in the models by using mainly
high-accuracy and high-resolution data in a local area. The other is that we can obtain a stable solution with high resolution
when we invert the integrals in the models because of the rapid decrease of the kernel function of the integrals. As a result,
with these models we evaluate one quantity with high resolution, in a band limited by the maximum degree of a set of geopotential
coefficients or by the resolution (spacing) of the local data, from another quantity (or quantities) in a local area, and
the resulting solution is stable.
Received: 6 April 1998 / Accepted: 16 June 1999 相似文献
13.
An iterative method is presented which performs inversion of integrals defined over the sphere. The method is based on one-dimensional fast Fourier transform (1-D FFT) inversion and is implemented with the projected Landweber technique, which is used to solve constrained least-squares problems reducing the associated 1-D cyclic-convolution error. The results obtained are as precise as the direct matrix inversion approach, but with better computational efficiency. A case study uses the inversion of Hotine’s integral to obtain gravity disturbances from geoid undulations. Numerical convergence is also analyzed and comparisons with respect to the direct matrix inversion method using conjugate gradient (CG) iteration are presented. Like the CG method, the number of iterations needed to get the optimum (i.e., small) error decreases as the measurement noise increases. Nevertheless, for discrete data given over a whole parallel band, the method can be applied directly without implementing the projected Landweber method, since no cyclic convolution error exists. 相似文献
14.
本文推导的椭球谐系数和球谐系数相互之间转换关系的核心思想是在ε~2量级下利用Legendre函数的正交性,从球谐系数求解的积分表示出发,将积分中的椭球坐标变量与球坐标变量相互转换,从而得出椭球谐系数与球谐系数之间的转换关系。本文导出的转换关系有以下优点:①对于第二类Legendre函数的计算采用Laurent级数表示,使计算第二类Legendre函数更为简单;②保留了ε~2量级下,导出的转换关系相比文献[2]的形式更简单,满足物理大地测量边值问题线性化的要求;③顾及了余纬和归化余纬的区别。 相似文献
15.
Due to the fact that the spectrum of a convolution is the product of the spectra of the two convolved functions, the convolution
integrals of physical geodesy can be evaluated very efficiently by the use of the fast Fourier transform (FFT) provided that gravity and/or terrain data are available on a regular grid. All Fourier transform-based methods usually treat
the gridded data as point values despite the fact that these discrete values may have been obtained by averaging and they
represent mean values over the whole area of a grid element. In the frequency domain, this fact can be taken into account
very easily, because the spectra of mean and point data are related via a two-dimensional (2D) sinc function. The paper shows explicitly this relationship using the convolution integrals for the evaluation of geoid
undulations, deflections of the vertical, and gravity and gradiometry terrain effects. Numerical tests are presented, indicating
that the differences in the two approaches may become significant when highly accurate results are wanted. The application
of the2D sinc function in the evaluation, update, and inversion of other convolution integrals is briefly discussed as well. 相似文献
16.
通过求解负荷勒夫数满足的微分方程的渐近解研究负荷勒夫数的渐近表达式。得出的结论比Farrell的经典结论精确一个数量级。 相似文献
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.
19.
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%. 相似文献
20.
M. K. Paul 《Journal of Geodesy》1973,47(4):413-425
Neglecting distant zones in the computation of geoidal height using Stokes' formula gives rise to some truncation error. This
truncation error is expressible as a weighted summation of the zonal harmonic components of the gravity anomaly. Making use
of the well-known properties of Legendre polynomials, a compact method of computing these theoretical coefficients has been
developed in this paper. 相似文献