利用经验模态分解提高极移预报精度（英文）

经验模态分解(Empirical Mode Decomposition,EMD)是一种数据驱动的自适应非线性、非平稳信号分解方法。为提高极移预报精度,将经验模态分解应用于极移预报中。首先利用经验模态分解方法对极移序列进行分解,获得极移的高频分量和低频分量;然后采用最小二乘(Least Squares,LS)外推模型对极移低频分量进行拟合,获得最小二乘拟合残差;其次采用自回归(Autoregressive,AR)模型对极移高频分量和最小二乘拟合残差之和进行建模预报;最后将最小二乘模型和自回归模型外推值相加获得极移预报值。将经验模态分解和LS+AR组合模型预报结果与LS+AR模型预报以及地球定向参数预报比较竞赛(Earth Orientation Parameters Prediction Comparison Campaign,EOP PCC)的预报结果进行比较,结果表明,将经验模态分解应用于极移预报中,可以明显改善极移预报精度。     

一种改正LS+AR模型提高短期极移预报精度的方法

最小二乘(Least Squares,LS)与自回归(Auto Regressive,AR)联合(LS+AR)模型在极移预报(polar motion,PM)中存在以下问题:最小二乘拟合的内部残差值较好而LS外推的残差值较大;LS拟合残差序列是非线性的,故根据预报历元前的残差序列建立的AR模型可能并不适用于待预报的残差序列,存在不匹配预报的情况.针对这两个问题,通过以下方法进行解决:首先对LS拟合数据两端点附加约束条件使其固定到LS拟合曲线上,因此在两端点附近的拟合值与观测值十分接近;然后选取与LS外推残差序列变化趋势接近的内推残差序列作为AR模型的建模对象,进行残差预报.通过实例表明该方法能够有效地提高LS+AR模型在短期极移预报的精度.此外,通过与RLS(Robustified Least Squares)+AR、RLS+ARIMA(Auto Regressive Integrated Moving Average)和LS+ANN(Artificial Neural Network)模型的预报结果对比,证明了该方法在极移预报中的可行性.实例证明了所提出的方法在短期预报中可以取得良好的预报结果,尤其在1–10d超短期的极移预报上可以获得与国际最好预报精度相当的预报结果.     

考虑LS+AR模型基础数据量的极移预报研究

针对目前极移最小二乘(Least Square, LS)+自回归(AutoRegressive, AR)预报模型的单一数据选取方案, 提出分别考虑LS模型数据量和AR残差数据量的组合数据模式, 并对极移预报时单一数据和组合数据预报结果精度进行分析, 探讨模型输入数据量对极移预报精度的影响. 结果表明, 模型输入数据量的变化对极移预报结果影响较大. 采用组合数据预报的方式相比较于单一数据量预报方式精度更高, 特别是针对30--360 d跨度内的中长期预报, 组合数据量的极移预报精度可比单一数据量预报精度有较大改善. 结论证明组合数据在极移预报时具有一定的优势, 可为以后极移预报数据量选取提供一定的借鉴参考意义.     

基于混沌时间序列的Volterra自适应滤波极移预报方法

针对极移复杂的时变特性, 根据混沌相空间坐标延迟重构理论, 提出一种基于Volterra自适应滤波的极移预报方法. 首先, 利用最小二乘拟合算法分离极移序列中的线性趋势项、钱德勒项和周年项, 获得线性极移、钱德勒极移和周年极移的外推值; 其次, 通过C-C关联积分法对最小二乘拟合残差序列进行相空间重构, 并利用小数据量法计算残差序列的最大Lyapunov指数验证其混沌特性, 在此基础上, 构建Volterra自适应滤波器对残差序列进行预测; 最后, 将线性极移、钱德勒极移和周年极移的外推值以及最小二乘拟合残差的预测值相加获得极移最终预报值. 利用国际地球自转服务局(International Earth Rotation and Reference Systems Service, IERS)提供的极移数据进行1--60d跨度预报, 并将预报结果分别与国际地球定向参数预报比较竞赛(Earth Orientation Parameters Prediction Comparison Campaign, EOP PCC)结果和IERS A公报发布的极移预报产品进行对比, 结果表明: 对于1--30d的短期预报, 该方法的预报精度与EOP PCC最优预报方法相当, 当预报跨度超过30d时, 该方法的预报精度低于EOP PCC最优预报方法, 优于参与EOP PCC的其他方法; 与IERS A公报相比, 该方法的短期预报效果较好, 当预报跨度增加时预报精度低于IERS A公报. 预报结果表明该方法更适合于极移短期预报.     

基于经验模式分解和最小二乘支持向量机的卫星钟差预报

针对卫星钟差(satellite clock bias,SCB)呈现非线性、非平稳变化的特性,提出结合经验模式分解(empirical mode decomposition,EMD)和最小二乘支持向量机(least squares support vector machines,LSSVM)的钟差预报方法.首先对钟差相邻历元间作一次差,并利用经验模式分解将差分序列分解成若干不同频率的平稳分量,分解后的分量突出了差分序列不同的局部特征;然后根据各个分量的变化规律,选择合适的核函数和相关参数构造不同的最小二乘支持向量机模型分别预报;最后将各分量预报值叠加得到一次差预报值,再将其还原得到钟差预报值.实验结果表明,所提方法与常用的二次多项式(quadratic polynomial,QP)模型、灰色系统(grey model,GM)模型和单一的最小二乘支持向量机模型相比,具有较高的预报精度和较强的泛化能力.     

基于极限学习机的极移中长期预报

传统的极移预报多是基于最小二乘外推和自回归等线性模型,但极移包含了复杂的非线性成分,线性模型的预报效果往往不甚理想。将一种新型神经网络极限学习机(Extreme Learning Machine, ELM)用于极移中长期预报。首先利用最小二乘外推模型对极移序列进行拟合,获得趋势项外推值,然后采用极限学习机对最小二乘拟合残差进行预报,最终的极移预报值为趋势项外推值与残差预报值之和。将极限学习机的预报结果同反向传播(Back Propagation, BP)神经网络与地球定向参数预报比较活动(Earth Orientation Parameters Prediction Comparison Campaign, EOP PCC)的预报结果进行对比,结果表明,极限学习机用于极移中长期预报是高效可行的。     

经验模态分解方法在脉冲星计时残差分析及预报中的应用

高精度脉冲星计时分析在脉冲星诸多应用领域中占重要地位,而多个未能消除的误差源的存在,限制了脉冲星计时精度的进一步提高.分析其观测特征,这种影响将导致脉冲星计时残差波动的多样性.利用从北美计时阵(NANOGrav)观测的4颗毫秒脉冲星数据,采用经验模态分解(EMD)方法对计时残差的波动特征进行了研究.通过分析发现:计时残差的振荡可用几个内在本征模态函数和一个剩余趋势揭示其多时间尺度的振荡结构和非平稳性.进一步比对脉冲星间多振荡结构对计时残差波动贡献的差异,发现经验模态分解非常有利于寻找脉冲星计时残差数据间的共性及个性.同时,采用自回归滑动平均模型进行残差预报,发现应用经验模态分解方法将残差序列分解后再预报,在短期预报时具有明显的优势.上述结果有利于拓展脉冲星计时残差在时域的分析.     

基于有效角动量函数的多参数极移预报

极移的变化与多种激发息息相关,这些激发包括大气表面压力和大气风、海底压力和洋流、陆地水分布以及气候变暖导致的海平面变化,并且可以通过有效角动量函数来估计.在极移预报中,通过刘维尔方程融合有效角动量函数,并利用最小二乘与自回归组合的方法进行拟合及外推,同时,对自回归模型的可调节参数设置更多的选择,在不同的极移预报阶段,对于不同分量的预报匹配更优的参数,有效地提高了极移的预报精度.在441次1–90 d的极移预报实验中,短中期的预报改善更为明显,在1–6 d和7–30 d的极移X预报结果中,分别有56.9%和53.5%优于国际地球自转服务(International Earth Rotation Service, IERS)的预报;在1–6 d和7–30 d的极移Y预报结果中,分别有66.5%和59.7%优于IERS的预报.整体上,极移Y的预报精度比极移X的预报精度有更多的提升,以IERS的地球定向参数(Earth Orientation Parameters, EOP)产品EOP 14C04 (IAU2000A)为参考,极移X预报在第1 d、第5 d的MAE (Mean Absolute...     

基于小波分析和神经网络的卫星钟差预报性能分析

为了有效地进行卫星钟差预报和更好地反映卫星钟差特性,提出了一种基于小波分析和神经网络的4阶段混合模型来实现卫星钟差的预报,并给出了基于小波分析和径向基函数(Radial Basis Function,RBF)神经网络进行卫星钟差预报的基本思想、预报模型和实施步骤.采用"滑动窗"划分数据,利用神经网络预测小波分解和去噪后的钟差序列各层系数,更精确地把握钟差序列复杂细致的变化规律,从而更好地逼近钟差序列.为验证该混合预报模型的可行性和有效性,利用GPS卫星钟差数据进行钟差预报精度分析,并将其与灰色系统模型和神经网络模型进行比较分析.仿真结果显示,该模型具有较好的预报精度,可为实时GPS动态精密单点定位提供较高精度的卫星钟差.     

基于Takagi-Sugeno模糊神经网络模型的卫星钟差预报方法

卫星钟差预报精度的不断提升是精密导航的关键问题.为了进一步提高钟差的预报精度和更好地反映钟差的变化特性,提出一种基于Takagi-Sugeno模糊神经网络(Fuzzy Neural Network,FNN)的钟差预报方法.该方法首先根据钟差数据的特点对钟差进行预处理,然后以预处理后的数据建立一种高精度预报钟差的Takagi-Sugeno模糊神经网络算法.采用IGS(International Global Navigation Satellite System Service)不同采样间隔的精密钟差数据进行了短期预报试验,并与ARIMA(Auto-Regressive Integrated Moving Average)模型、GM(1,1)模型及QP(Quadratic Polynomial)模型进行了对比试验,分析结果表明:对不同类型原子钟,该方法用于钟差短期预报是可行的、有效的,其获得的卫星钟差预报结果明显优于常规方法.     

Acceleration of Numerical Integration of the Equations of Motion of Asteroids

Solar System Research - Finding and studying possible collisions of asteroids approaching the Earth requires a significant amount of computation. This paper describes the R0 program created to...     

Equation of state and anisotropy of pressure of magnetars

Magnetars are the neutron stars with the highest magnetic fields up to 10¹⁵–10¹⁶ G. It has been proposed that they are also responsible for a variety of extra-galactic phenomena, ranging from giant flares in nearby galaxies to fast radio bursts. Utilizing a relativistic mean field model and a variable magnetic field configuration, we investigate the effects of strong magnetic fields on the equation of state and anisotropy of pressure of magnetars. It is found that the mass and radius of low-mass magnetars are weakly enhanced under the action of the strong magnetic field, and the anisotropy of pressure can be ignored. Unlike other previous investigations, the magnetic field is unable to violate the mass limit of the neutron stars.     

Observations of the intensity of the penumbra of sunspots

Observations of the penumbral intensity of sunspots in 13 wavelength regions are presented. In 4 wavelength regions 54 sunspots are measured. In the other wavelength regions the number of sunspots considered ranges from 3–19.The penumbral intensity alters with position within the spot. This intensity variation is found to be comparable with the change in intensity from one spot to another. The penumbral intensity is found to be independent of spot size in the sample considered.The penumbra model of Kjeldseth Moe and Maltby (1969) with = 0.055 is supported by the measurements.     

A survey of orbits of co-orbitals of Mars

Many asteroids with a semimajor axis close to that of Mars have been discovered in the last several years. Potentially some of these could be in 1:1 resonance with Mars, much as are the classic Trojan asteroids with Jupiter, and its lesser-known horseshoe companions with Earth. In the 1990s, two Trojan companions of Mars, 5261 Eureka and 1998 VF₃₁, were discovered, librating about the L₅ Lagrange point, 60° behind Mars in its orbit. Although several other potential Mars Trojans have been identified, our orbital calculations show only one other known asteroid, 1999 UJ₇, to be a Trojan, associated with the L₄ Lagrange point, 60° ahead of Mars in its orbit. We further find that asteroid 36017 (1999 ND₄₃) is a horseshoe librator, alternating with periods of Trojan motion. This asteroid makes repeated close approaches to Earth and has a chaotic orbit whose behavior can be confidently predicted for less than 3000 years. We identify two objects, 2001 HW₁₅ and 2000 TG₂, within the resonant region capable of undergoing what we designate “circulation transition”, in which objects can pass between circulation outside the orbit of Mars and circulation inside it, or vice versa. The eccentricity of the orbit of Mars appears to play an important role in circulation transition and in horseshoe motion. Based on the orbits and on spectroscopic data, the Trojan asteroids of Mars may be primordial bodies, while some co-orbital bodies may be in a temporary state of motion.     

Graphical representation of periodic motion of the line of nodes of the Moon

In the text-books of astronomy, sections generally related to the Moon deal with the orbital elements of the Earth-Moon system such asa, e, i, , and the time of perigee passage. While the MEAN of the first of the three elements do not vary, mean longitude of the ascending node-mean longitude of the lunar perigee and the time of perigee passage undergoes secular as well as periodic changes due predominantly to the action of the Sun's gravitational attraction. While to a certain degree, explanations related to the calculation of the lunar orbit parameters are given, not a single graphical representation of these short- or long-periodic changes are presented. We allow the number of data related to these periodic changes must cover a large span of time; and if regression of the line of nodes or advances of the line of apses are to be graphically seen, data covering 18.61 and 8.85 yr, respectively, are needed. In this work we particularly aim at the graphical representation of the periodic changes of the line of nodes.