共查询到17条相似文献，搜索用时 140 毫秒
1.
超高阶缔合勒让德函数的计算在地球重力场研究中是影响计算效率的关键因素之一,其计算耗时会随着截断阶数的增加而呈指数上涨。常用的缔合勒让德函数递推算法有标准前向列推法、标准前向行推法、Belikov递推法和跨阶次递推法。为有效提高缔合勒让德函数递推的计算效率,提出利用基于OpenMP的多核并行算法对上述方法进行并行加速,并通过提高数组维度的方法解决了递推运算并行化过程中的不适用问题。实验表明,所提并行算法在保持相同精度的前提下显著提高了缔合勒让德函数递推的计算效率,加速比在3倍以上,最高可达4.6倍。 相似文献
2.
3.
详细阐述了利用地球重力场模型位系数计算高程异常的数学模型。选取标准向前列递推法和跨阶次递推法作为缔合勒让德函数的计算方法,分析了两种算法的适用范围。以高精度的超高阶EIGEN-6C4重力场模型为例,选择全球范围内5°×5°的格网组成的2701个点和具有代表性的12个点计算高程异常,将计算结果和ICGEM网站计算结果比较验证了算法的可靠性,并绘制了全球范围内5m等高距的高程异常等值线图,实现了全球高程异常的可视化。 相似文献
4.
5.
重力场计算中,经常需要计算以有限阶球谐级数表示的重力场元。常规计算中除需存储(N 1)^2个系数值外,还需迭代计算出(N 2)(N 1)／2个完全正常化勒让德函数值。Clenshaw求和法不需要计算单个球谐函数值而直接计算级数和,因而计算速度上有所提高。总结了球谐函数的零阶导数级数和,并推导了一阶导数级数和。通过数值试验,对于任意点的重力场元,使用C1enshaw求和法计算零阶导数球谐函数和比常规方法节省一半的时间,一阶导数球谐函数之和的计算速度提高幅度不大,并分析了其中的原因。 相似文献
6.
利用地球重力场扰动位系数计算网格剩余垂线偏差 总被引：3，自引：1，他引：3
介绍一种根据重力场扰动位系数计算网络剩余垂线偏差的方法，导出了一组完全正规化连带勒让德函数积分的实用公式，经卫星测高数据处理的实际应用，证明这种方法是有效的。 相似文献
7.
8.
台湾海峡及其邻区岩石圈层下地幔流应力场 总被引：3，自引：0，他引：3
本文应用卫星测到的2～30阶球谐函数系数(GEM8) 计算了台湾海峡及邻近地区的岩石圈层下地幔流应力场。在比较了该区域的构造地质以后,发现应力场的方向,尤其是在板块收敛区域的甚高阶场(26～30阶)的方向与在菲律宾海和南海等区域的板块运动方向一致。 相似文献
9.
全球重力场模型是当今物理大地测量学最为活跃的研究领域之一。本文基于目前国内外最新的重力场模型理论研究成果,提出了利用中国地区细部数据和全球卫星测高2’×2’网格重力异常扩展超高阶位模型的计算方法,详细讨论了数值解算过程中的稳定性和可靠性问题。以EGM96和GPM98CR模型作为参考模型,在全球意义上分别解算得到MOD99a(360阶)、MOD99b(720阶)和MOD99c/d(1800阶),将系列模型MOD99a/b/c/d同中国地区72个GPS水准大地水准面和全球海洋12个地区的卫星测高大地水准面进行了比较,并通过功率谱分析方法检验了4组模型的有效性和可靠性。 相似文献
10.
应用理论推导及数值计算方法，对Stokes随机波的谱特性进行了分析。首先将波面方程，海水质点水平速度用一阶波面分量的非线性组合表示，应用平稳随机高阶短的降阶计算法则，得到了波面方程及海水质点水平速度与一阶波面分量的自相关函数之间的关系，从而确定了Stokes随机波浪的波浪谱密度及海水质点水平速度和加速度谱密度，进而求得有关波浪要素的均方根值。文章还应有数值计算方法，分析了波浪基本参数对均方根值的影响。 相似文献
11.
针对球谐函数定积分计算中Legendre函数递推问题展开研究,分析了标准向前列推法、Belikov法、跨阶次法、X数法以及顾及麦克劳林级数展开式对球谐函数定积分计算的影响。利用Eigen6c-4地球重力场模型计算扰动引力梯度径向分量,分析不同方法之间的差异。实验表明,不考虑麦克劳林级数展开式时4种方法的相对精度在高纬度地区较差,但计算模型扰动引力径向分量的精度一致,结合麦克劳林级数式可提高高纬度地区定积分计算的相对精度,但会降低中低纬度地区定积分计算的精度,并且对高纬度地区扰动引力径向分量的影响极小,但会严重降低低纬度地区扰动引力梯度计算的精度。 相似文献
12.
Note on axisymmetric consolidation of multi-layered poroelastic soils with compressible constituents
In the conventional analysis of Biot’s axisymmetric consolidation, the solid phase and the surrounding fluid were often assumed incompressible for simplicity. Such assumption in soil engineering ignored the effect of compressibility of constituents on the consolidation. In this article, the compressible fluid and solid in soil were taken into account for 3-D consolidation. The pore pressure, the displacements, and the stresses were expressed by two displacement functions, and the Laplace–Hankel transform was applied to set up the stiffness matrix between the generalized displacement and stress. The stiffness matrix consisted of negative exponential functions, ensuring that the computation is efficient and stable. Then the global stiffness matrix is extended by embracing the continuity of the interfaces and boundary conditions of soil base. The relationship between the generalized displacement and stress of the soil base reduces the number of unknowns of the global matrix. Such consideration reduced the number of unknowns of the global matrix and brought in more acceptable boundary conditions where the stiffness and the permeability of semi-infinite soil base can be taken into account. After the inversion of the Laplace–Hankel transform, the real solutions were obtained. The results show that the stiffness and permeability of the soil base can change the development consolidation and that the compressibility of solid particles has a great effect on the settlement in the beginning of consolidation. 相似文献
13.
In order to improve the ocean forecasting in the North Sea and Baltic Sea, an assimilation scheme based on a bottom-topography-dependent anisotropic recursive filter has been used in this study. This scheme can stretch or flatten the shape of a local representative contour surface of the background error covariance function into the form of an ellipse. Furthermore, the computing efficiency has been largely improved due to implicit computation of the background error covariance. A two-month experiment has been used for verifying the impact of assimilating ocean profile observations on ocean forecasting. The results indicate that the use of temperature and salinity profiles can largely improve the oceanic forecasting. The root mean square differences between the forecasts and observations for temperature and salinity have been reduced by 36% and 18% in the experiment period, respectively. Moreover, it is found that the anisotropic recursive filter approach is especially efficient in areas with complex coastlines and sharp fronts, e.g., inner Danish waters. The results also show that the propagation of observation information from an observation position to its neighboring grid points is closely related to currents. 相似文献
14.
A solution is presented for the wave induced drift forces acting on a submerged sphere in a finite water depth based on linearised velocity potential theory. In order to obtain the velocity potential, use has been made of multipole expansions in terms of an infinite series of Legendre functions with unknown coefficients. The series expression for the second order mean forces (drift forces) is provided by integrating the fluid pressure over the body surface. The horizontal drift force is also expressed by a series solution obtained using the far-field method. 相似文献
15.
16.
A model-based approach is proposed to solve the oceanic internal wave signal processing problem that is based on state-space representations of the normal-mode vertical velocity and plane wave horizontal velocity propagation models. It is shown that these representations can be utilized to spatially propagate the modal (depth) vertical velocity functions given the basic parameters (wave numbers, Brunt-Vaisala frequency profile, etc.) developed from the solution of the associated boundary value problem as well as the horizontal velocity components. These models are then generalized to the stochastic case where an approximate Gauss-Markov theory applies. The resulting Gauss-Markov representation, in principle, allows the inclusion of stochastic phenomena such as noise and modeling errors in a consistent manner. Based on this framework, investigations are made of model-based solutions to the signal enhancement problem for internal waves. In particular, a processor is designed that allows in situ recursive estimation of the required velocity functions. Finally, it is shown that the associated residual or so-called innovation sequence that ensues from the recursive nature of this formulation can be employed to monitor the model's fit to the data 相似文献
17.
WANG Yanying QIAN Kun ZHURenchuan Professor School of Naval Architecture Marine Engineering Dalian University of Technology 《中国海洋工程》1998,(4)
In order to determine the extreme sea-state encountered by oceangoing vessels,a method forthe prediction of wave parameters based on the data of long-term distribution sample is developed by useof the fitting approach with the Weibull probability functions and the power functions.As an example aspecial calculation result is given including the analysis of wave parameters for the data of the NorthAtlantic Basin and the computation of motion and load encountered by two container ships CSLR andCSBV in different loading conditions.All computation results are satisfactory compared with corre-sponding design results. 相似文献