共查询到15条相似文献，搜索用时 171 毫秒
1.
2.
3.
针对球谐函数定积分计算中Legendre函数递推问题展开研究,分析了标准向前列推法、Belikov法、跨阶次法、X数法以及顾及麦克劳林级数展开式对球谐函数定积分计算的影响。利用Eigen6c-4地球重力场模型计算扰动引力梯度径向分量,分析不同方法之间的差异。实验表明,不考虑麦克劳林级数展开式时4种方法的相对精度在高纬度地区较差,但计算模型扰动引力径向分量的精度一致,结合麦克劳林级数式可提高高纬度地区定积分计算的相对精度,但会降低中低纬度地区定积分计算的精度,并且对高纬度地区扰动引力径向分量的影响极小,但会严重降低低纬度地区扰动引力梯度计算的精度。 相似文献
4.
详细阐述了利用地球重力场模型位系数计算高程异常的数学模型。选取标准向前列递推法和跨阶次递推法作为缔合勒让德函数的计算方法,分析了两种算法的适用范围。以高精度的超高阶EIGEN-6C4重力场模型为例,选择全球范围内5°×5°的格网组成的2701个点和具有代表性的12个点计算高程异常,将计算结果和ICGEM网站计算结果比较验证了算法的可靠性,并绘制了全球范围内5m等高距的高程异常等值线图,实现了全球高程异常的可视化。 相似文献
5.
重力场计算中,经常需要计算以有限阶球谐级数表示的重力场元。常规计算中除需存储(N 1)^2个系数值外,还需迭代计算出(N 2)(N 1)／2个完全正常化勒让德函数值。Clenshaw求和法不需要计算单个球谐函数值而直接计算级数和,因而计算速度上有所提高。总结了球谐函数的零阶导数级数和,并推导了一阶导数级数和。通过数值试验,对于任意点的重力场元,使用C1enshaw求和法计算零阶导数球谐函数和比常规方法节省一半的时间,一阶导数球谐函数之和的计算速度提高幅度不大,并分析了其中的原因。 相似文献
6.
超高阶缔合勒让德函数的计算在地球重力场研究中是影响计算效率的关键因素之一,其计算耗时会随着截断阶数的增加而呈指数上涨。常用的缔合勒让德函数递推算法有标准前向列推法、标准前向行推法、Belikov递推法和跨阶次递推法。为有效提高缔合勒让德函数递推的计算效率,提出利用基于OpenMP的多核并行算法对上述方法进行并行加速,并通过提高数组维度的方法解决了递推运算并行化过程中的不适用问题。实验表明,所提并行算法在保持相同精度的前提下显著提高了缔合勒让德函数递推的计算效率,加速比在3倍以上,最高可达4.6倍。 相似文献
7.
基于势流理论和时域高阶边界元方法,建立了三维完全非线性数值波浪水槽模型.利用源造波法产生入射波浪,应用五阶斯托克斯波理论给定波浪速度;采用混合欧拉-拉格朗日方法追踪流体瞬时水面,将二阶泰勒级数展开法应用于更新下一时间步的波面和速度势;通过加速势的方法准确计算自由水面速度的法向导数和物面速度势的时间导数.对完全非线性波浪进行了模拟,得到了稳定的波形.当波浪非线性较小时,与四阶Runge- Kutta法(RK4)计算结果和五阶斯托克斯波理论解均吻合良好;随着波浪非线性的增大,计算结果误差逐渐增大.通过数值试验分析,在满足精度要求的基础上,本方法计算时间略大于四阶Runge- Kutta法的四分之一,大大减少了计算量. 相似文献
8.
9.
10.
11.
The Seasat altimeter data has been completely adjusted by a crossing arc technique to reduce the crossover discrepancies to
approximately ±30 cm in five regional adjustments. This data was then used to create sea surface heights at 1° intersections
in the ocean areas with respect to the GRS80 ellipsoid. These heights excluded the direct tidal effects but included the induced
permanent deformation.
A geoid corresponding to these sea surface heights was computed, based on the potential coefficients of the GEML2 gravity
field up to degree 6, augmented by Rapp's coefficients up to degree 180. The differences between sea surface heights and the
geoid were computed to give approximate estimates of sea surface topography. These estimates are dominated by errors in both
sea surface heights and geoid undulations. To optimally determine sea surface topography a spherical harmonic analysis of
raw estimates was carried out and the series was further truncated at degree 6, giving estimates with minimum wavelengths
on the order of 6000 km.
The direction of current flow can be computed on a global basis using the spherical harmonic expansion of the sea surface
topography. Ths has been done, not only for Seasat/GEML2 estimates, but also using the recent dynamic topography estimates
of Levitus. The results of the two solutions are very similar and agree well with the major circulation features of the oceans. 相似文献
12.
The constant and harmonic parts of the global ocean tide are modeled by up to nine major tidal constituents, namely, S2, M2,
N2, K1, P1, O1, Mf, Mm, and Ssa. Our computations start with the Fourier sine and cosine series expansion for the tidal constituents,
including the constant Mean Sea Level (MSL). Although the frequencies of the tidal constituents are considered known, the
coefficients of the sine and cosine functions are assumed to be unknown. Subsequently, the coefficients of the sine and cosine
functions, as well as the constant part of the Fourier expansion, were expanded into spherical harmonics up to degree and
order n, where n corresponds to the number of linearly independent spherical harmonic base functions needed to model the tidal constituents,
determined via independent columns of the Gram matrix. The unknown coefficients of the spherical harmonic expansions are computed
using sea level observations within cycles #1–#350 of the TOPEX/Poseidon satellite altimetry over 11 years of its mission.
A set of orthonormal base functions was generated for the marine areas covered by TOPEX/Poseidon observations from the spherical
harmonics using a Gram-Schmidt orthogonalization process. These were used for modeling the dominant tidal constituents. The
computed models based on orthonormal base functions for the nine tidal constituents and the constant part of the Fourier expansion,
were tested numerically for their validity and accuracy, proving centimeter accuracy. 相似文献
13.
The approach presented is directed toward a specific adaptation of the least‐squares collocation with noise, yielding smooth predictions of geophysical quantities. The smoothing corresponds here to a truncated gravity field equivalent to an (n’, n') spherical‐harmonic expansion. This is reflected in the truncation, at the degree n‘, of the pertinent covariance and cross‐covariance functions in most (but not all) instances. The smooth predictions of geophysical quantities, made in an equilateral grid corresponding to the truncation degree n‘, serve in constructing contour maps after having been densified for the needs of a contour routine. Such a densification is carried out efficiently via errorless collocation with the degree truncation n‘ throughout. Consistent with this procedure, “residuals” at observation points (i.e., discrepancies between the contour map and the data) are computed using the same algorithm. The complete collocation approach is utilized for a 2° resolution of the earth's gravity field with emphasis on the oceanic geoid, based on the residuals from a global spherical‐harmonic adjustment of SEASAT altimetry. The presented results include contour maps of geoid undulations and gravity anomalies. They are compared to the results of a point‐mass adjustment, another technique based on the spherical‐harmonic adjustment. The agreement between these two techniques is found to be excellent. 相似文献
14.
15.
基于106个月的Grace（gravity recovery and climate experiment）月重力场模型（120阶次）,消除了月重力场的月、季度及年度变化,得到了稳态的地球重力场模型（Grace_sta）。在2～120阶次之间,Grace_sta与已有高阶重力场EGM2008及EGM96三个模型的阶方差是一致的。在2～100阶次之间,Grace_sta模型误差阶方差要小于EGM2008与EGM96误差阶方差。在全球范围内,Grace_sta重力场的大地水准面与EGM2008相应阶次的大地水准面标准差约为3cm,与EGM96模型大地水准面差异则高达52cm。结果表明,Grace_sta足可以取代EGM2008重力场模型2～100阶次的低阶部分,新得到的稳态重力场模型可为海面地形分析提供了可靠的参考场。 相似文献