重力场数学拟合模型解的影响因素探究     

鹿荻  陈锴 《测绘科学》2016,41(3):28-33

为了提高地球重力场模型不适定方程求解的精度,该文采用谱分析方法从级数展开阶数、数据采样率及数据缺失量3个方面探索影响数学拟合效果的根本因素:从常用的三角级数及勒让德级数模型出发,引出重力场拟合模型球谐函数模型,观察在改变级数展开阶数、数据采样率及数据缺失量等情况下所对应设计矩阵谱结构的变化,并从微观上研究影响误差分配的有关因素及最小奇异值对误差的决定性作用,为探求重力场模型解不准的原因及实现更高精度的全球重力场模型的建立提供参考。  相似文献   

Moritz解析延拓解与Bjerhammar虚拟球面解的等价性          下载免费PDF全文

程芦颖  许厚泽 《海洋测绘》2008,28(3):6-7

利用调和函数的性质以及球谐函数展开理论,并根据实数域上函数的幂级数展开式证明了Moritz解析延拓解与Bjerhammar虚拟球面解的等价性,同时分析了两种解的内在区别。  相似文献   

球面距离计算方法及精度比较     

樊东卫  何勃亮  李长华  韩军  许允飞  崔辰州 《天文研究与技术》2019,(1):69-76

球面距离(角间距)计算是天文或地理学中极常用的计算之一,也是目标查找、锥形检索、交叉证认等方法的基础。数学上,通过球面几何可以直接计算出两点的距离,前人已经推导出了多个复杂程度不一的计算方法。但是由于计算机的精度有限,在进行数值计算时有舍入误差,导致公式计算结果出现偏差。对几个常用的球面距离计算公式进行了考察,测试并对比它们在不同计算环境下的精度与优缺点。此外还展示并比较了几种常用天文软件包、数据库的球面距离计算方法,以期有助于天文工作者选择适合自己当前需要的计算方法。  相似文献   

Higher harmonics induced by a submerged horizontal plate and a submerged rectangular step in a wave flume     

J. Brossard  G. Perret  L. BlonceA. Diedhiou 《Coastal Engineering》2009

The decomposition of a monochromatic wave over a submerged plate is investigated experimentally in a wave flume. Bound and free higher harmonic modes propagating upstream and downstream the structure are discriminated by means of moving resistive probes. The first-order analysis shows a resonant behaviour linked to the ratio of the plate's width and the fundamental mode wavelength over the plate. The second-order analysis shows an energy transfer from the fundamental mode towards free harmonics propagating downstream the structure. This transfer is linked to the ratio of the width of the plate and the bound harmonic wavelength over the plate. We also performed experiments with a submerged step to compare the efficiency of both structures. The submerged plate is shown to be a more efficient breakwater than the step, at the first as well as the second-order.  相似文献   

Triangulated spherical splines for geopotential reconstruction     

M. J. Lai  C. K. Shum  V. Baramidze  P. Wenston 《Journal of Geodesy》2009,83(8):695-708

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.  相似文献   

A study on the reference frame consistency in recent Earth gravitational models   总被引：1，自引：1，他引：0  

Christopher Kotsakis 《Journal of Geodesy》2009,83(1):31-50

All gravity field functionals obtained from an Earth gravitational model (EGM) depend on the underlying terrestrial reference frame (TRF), with respect to which the EGM’s spherical harmonic coefficients refer to. In order to maintain a coherent framework for the comparison of current and future EGMs, it is thus important to investigate the consistency of their inherent TRFs, especially when their use is intended for high precision studies. Following the methodology described in an earlier paper by Kleusberg (1980), the similarity transformation parameters between the associated reference frames for several EGMs (including the most recent CHAMP/GRACE models at the time of writing this paper) are estimated in the present study. Specifically, the differences between the spherical harmonic coefficients for various pairs of EGMs are parameterized through a 3D-similarity spatial transformation model that relates their underlying TRFs. From the least-squares adjustment of such a parametric model, the origin, orientation and scale stability between the EGMs’ reference frames can be identified by estimating their corresponding translation, rotation and scale factor parameters. Various aspects of the estimation procedure and its results are highlighted in the paper, including data weighting schemes, the sensitivity of the results with respect to the selected harmonic spectral band, the correlation structure and precision level of the estimated transformation parameters, the effect of the estimated differences of the EGMs’ reference frames on their height anomaly signal, and the overall feasibility of Kleusberg’s formulae for the assessment of TRF inconsistencies among global geopotential models.  相似文献   

Least squares and integral methods for the spherical harmonic analysis of theSq-field     

Koji Kawasaki  S. Matsushita  J. C. Cain 《Pure and Applied Geophysics》1989,131(3):357-370

Spherical harmonic coefficients (SHCs) for the daily magnetic variation fields (solar and lunar) and the main field of the earth are usually estimated by the method of least squares applied to a truncated spherical harmonic series. In this paper, an integral method for computing the SHCs for the solar quiet daily magnetic variation fieldSq is described and applied toSq data for May and June 1965. TheSq SHCs thus derived are then compared with the results obtained using both unweighted and weighted versions of the least squares method. The weighting used tends to orthogonalize the least squares terms. The integral and weighted least squares results agree closely for terms up to order 4 and degree 30, but both disagree considerably for the higher degree terms with the results of the unweighted least squares. Errors introduced by the numerical integration can be shown to be small, hence the disagreement between integral and unweighted least squares coefficient sets arises from improper weighting. Also, it is concluded that discrepancies between the geomagnetic northward and eastward component-derived coefficient sets arise from either time-dependent external sources that produce non-local-time, based fields or nonpotential sources and not from truncation of the spherical harmonic series as has previously been suggested.Deceased.  相似文献   

Modelling of mantle postglacial relaxation in axisymmetric geometry with a composite rheology and a glacial load interpolated by adjusted spherical harmonics analysis     

G. Dal Forno  P. Gasperini 《Geophysical Journal International》2007,169(3):1301-1314

Although studies on glacial isostatic adjustment usually assume a purely linear rheology, we have previously shown that mantle relaxation after the melting of Laurentide ice sheet is better described by a composite rheology including a non-linear term. This modelling is, however, based on axially symmetric geometry and glacial forcing derived from ICE-3G and suffers from a certain amount of arbitrariness in the definition of the ice load. In this work we apply adjusted spherical harmonics analysis to interpolate the ice thicknesses of ICE-3G and ICE-1 glaciological models. This filters out the non-axisymmetric components of the ice load by considering only the zonal terms in the spherical harmonics expansion. The resulting load function is used in finite-element simulation of postglacial rebound to compare composite versus purely linear rheology. Our results confirm that composite rheology can explain relative sea level (RSL) data in North America significantly better than a purely linear rheology. The performance of composite rheology suggests that in future investigations, it may be better to use this more physically realistic creep law for modelling mantle deformation induced by glacial forcing.  相似文献   

外部质量产生的扰动位的球谐展开式     

于锦海 《测绘学院学报》1996,(3)

本文讨论了参考椭球面之外，地球质量产生的扰动位的球谐展开式。分别在地球的平均球面上以及地球表面之外，给出了该扰动位的球谐展开式，为进一步研究外部质量对扰动位的影响建立了理论基础。  相似文献   

实际预报可预报性的时空依赖性分析   总被引：4，自引：1，他引：4  

李志锦  纪立人 《大气科学》1996,20(3):290-297

本文利用国家气象中心1990年1月1日至1992年2月29日的1～5 d、500 hPa高度场业务预报结果研究了不同球谐函数谱分量和经验正交函数(EOF)谱分量的可预报性。对球谐函数谱分量的研究表明可预报性并不是随着空间尺度的增加呈现单调的递减关系，主要表现为在总波数n等于5附近具有最大可预报性。可预报性主要依赖于总波数n，经卷大圆上的波数和纬向圆周上的波数对可预报性具有完全相同的重要意义。对EOF谱分量的分析表明，随着EOF指数(即序号数)的增大，可预报性依次减小。从EOF和球谐函数的关系发现前几个EOF分量正是对应着最可预报的球谐函数分量。进一步分析表明，最可预报的分量对应大气运动的慢变流型。这些结果对如何使用数值预报产品以及如何进行延伸预报具有重要意义。  相似文献