共查询到20条相似文献，搜索用时 15 毫秒
1.
为了更充分利用低轨重力卫星的高精度观测数据,根据卫星轨道的扰动理论,导出了应用卫星轨道与星间距离观测值联合反演地球重力场模型的算法.该算法的实质是将牛顿运动方程在卫星轨道处进行展开,转化为第二类Volterra积分方程,并采用基于移动窗口的9次多项式内插公式进行数值求解.给出了该算法的观测方程,用QR分解法消去局部参数矩阵,最后采用预条件共轭梯度法求解法方程.利用GRACE卫星2008-01-01~2008-08-01时间段内的轨道及星间距离观测数据,解算了120阶次的地球重力场模型SWJTU-GRACE01S,该模型在120阶处的阶方差为1.58×10-8,大地水准面差距累计误差为22.29 cm,与美国GPS水准网比较的标准差为0.793 m,结果表明:SWJTU-GRACE01S模型精度介于EIGEN-GRACE01S与EIGEN-GRACE02S模型之间,从而验证了该算法的有效性. 相似文献
2.
连续傍轴近似公式及其多种离散形式 总被引：5，自引：2，他引：3
基于单向波动法向透射的一般运动学模型和局部多项式插值函数,推导建立了连续形和2阶离散形法向透射傍轴近似公式。指出连续形傍轴近似公式具有多种离散形式,本文给出的标准形式、已有的多次透射公式和平面波多项式外推公式都是连续形傍轴近似公式的不同离散形式,因此具有基本相当的精度、稳定性和计算效率以及相同的普适性和局限性。 相似文献
3.
变异函数是地统计学中区域化变量空间结构分析和空间局部插值的主要分析工具.理论变异函数模型的获取是地质统计学中的基础性工作,它是了解区域化变量的变异特征、进一步对地质统计学计算的必要环节.针对现有的理论变异函数的拟合方法,如人工拟合法、线性规划拟合法、加权多项式拟合法、目标规划拟合法等的不足之处,充分利用粒子群优化算法在求解非线性优化问题时具有的全局寻优的特点,提出基于粒子群优化的理论变异函数拟合方法.在实例应用中,分别利用粒子群优化算法和加权多项式拟合方法进行理论变异函数拟合,交叉验证结果表明粒子群优化算法预测精度较高,具有较强的稳健性. 相似文献
4.
高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率,本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程,在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数,在求解大型稀疏线性方程组时利用直接求解器进行快速求解,从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法,能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应. 相似文献
5.
利用CHAMP卫星轨道和加速度计数据推求地球重力场模型的一种有效方法是能量守恒方法.本文基于能量守恒方程,推导了整体求解尺度、偏差、积分常量和位系数的计算公式,提出了整体求解位系数、积分常量、尺度和偏差参数的重力场恢复方法-改进的能量法.该方法摆脱了目前使用能量法时,必须预先通过其它外部手段标定加速度计数据以及积分常量通过近似计算获取的现状.基于德国慕尼黑技术大学(TUM)提供的约300天的CHAMP卫星几何法轨道和GFZ提供的加速度计数据,计算出了60×60阶地球重力场模型XISM-CHAMP1S,并与EIGEN-CG03C、EIGEN-CHAMP03S、EIGEN2、ENIGN1S、EGM96模型进行了比较.结果表明:XISM-CHAMP1S模型精度明显优于相同阶次EGM96、EIGEN1S和EIGEN2,与同阶次的EIGEN-CHAMP03S模型精度最为接近. 相似文献
6.
地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法（如有限元法、谱元法、高阶有限差分法等）,其计算速度较慢;计算速度较快的算法（如低阶有限差分法、付氏伪谱法等）计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法（基于Forsyte广义正交多项式的褶积微分算子法）,该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性. 相似文献
7.
针对目前建立区域地磁场模型截断阶数受限,地磁图精度偏低等缺陷,研究了运用Taylor多项式拟合法建立区域地磁场模型所存在的边界效应问题.从理论上深入分析了产生边界效应误差的根源,提出了采用边界插值约束的Taylor多项式拟合法,在改善区域地磁模型边界效应的同时,提高了模型拟合的截断阶数,进而绘制出较高精度的区域地磁场图,同时对该技术和现有补充国际地磁参考场值方法进行了比较分析,最后进行了仿真试验.试验证明,该技术具有优于补充国际地磁参考值的方法,能克服现有区域地磁场模型的边界效应问题,提高了区域地磁场模型的拟合精度. 相似文献
8.
利用CHAMP卫星矢量和标量地磁测量数据,通过三维Taylor多项式模型建立了2010.0年中国及邻近地区在300km高度附近的X、Y、Z和F分量的三维磁场模型.为了比较验证所建模型,分析了其与Taylor多项式模型(二维模型)不同截断阶数所对应的均方偏差(RMSE)、残差及分布等.结果表明,三维模型每一阶的RMSE和残差绝对值的平均值均要比二维模型的小约45%.由于采用了系数完全展开的建模方式,三维模型的系数数量约为二维模型的2倍,三维模型较低的截断阶数可以反映更多的地磁信息.本研究中的5阶三维Taylor模型基本可达到8阶Taylor模型的精度.两种模型绘制的地磁场及残差分布有较好的一致性. 相似文献
9.
10.
数字强震记录“精确”插值、积分和微分的权函数算法——I.时域离散序列数值分析方法与应用研究 总被引：4，自引：0，他引：4
现有强震记录常规处理程度中采用的插值、积分、微分运算方法计算简单、省时,但在某些情况下,不能确保所要求的计算精度,针对这一问题,本文假定强震记录是理想的有限长记录并满足采定理,在此基础上提出了对强震记录进行“精确”插值,积分和微分运算以及权函数的概念,并推导建立了相应的计算公式,理论和算例分析表明：（１）本文算法中使用的插值,积分和微分运算权函数均具有不变和快速收敛的性质;（２）只需要少量的权系数,本文算法就能获得足够的计算精度;（３）选取足够多的权系数,本文算法能以任意精度副近真实的精确解序列,（４）本文算法可以成为强震数据处理和相关结构动力分析的基础。 相似文献
11.
本文基于新型能量插值法,利用美国喷气推进实验室(JPL)公布的2008年的GRACE-Level-1B实测数据,反演了120阶GRACE地球重力场.首先,由于GPS轨道测量精度相对较低,通过将K波段测距仪高精度的星间距离观测量插值引入双星动能差中,进而建立了新型能量插值卫星观测方程.其次,详细对比分析了2点、4点、6点和8点能量插值观测方程对地球重力场反演精度的影响.研究结果表明:基于最优的信噪比,6点能量插值公式有利于提高120阶GRACE地球重力场的反演精度.最后,基于美国、欧洲和澳大利亚的GPS/水准观测数据检验了本文新建立的WHIGG-GEGM03S地球重力场模型的正确性和有效性. 相似文献
12.
基于双二次插值的探地雷达有限元数值模拟 总被引：3，自引：0，他引：3
从探地雷达(GPR)满足的波动方程出发,详细介绍了二维GPR模型单元剖分、二次插值、数值积分和有限元刚度矩阵总体合成的GPR有限元求解过程.为解决数值模拟时截断边界处的超强反射,采用Clay Bout透射边界条件对雷达波进行衰减,进而压制了来自截断边界处的反射波.在满足时间步长与空间网格差分稳定性前提下,采用中心差分法对GPR有限元方程进行离散,并用不完全LU分解预处理的BICGSTAB算法求解系数方程组,然后编制了基于双二次插值的GPR有限元正演模拟matlab程序.运用该程序分别对矩形和"V"字形两个典型地电模型进行正演计算,得到了正演剖面图,将该正演剖面图与基于线性插值的FEM算法的正演剖面图做了对比分析.结果表明基于双二次插值FEM算法相比基于双线性插值FEM算法异常响应更明显,具有更高的模拟精度,更有利于指导雷达剖面的数据解译. 相似文献
13.
本文提出了一种新的基于正交多项式变换来压制动校正后CMP道集上随机噪声的方法,将地震资料经过正交多项式变换,建立正交多项式系数谱,该谱描述了地震资料在正交多项式域的能量分布,可以较好地分离有效波和随机噪声,相对于固定阶次的基于正交多项式变换压制噪声的方法,所提出的方法根据能量随阶次变化的规律,自适应确定表达有效信号的正交多项式阶次,既提高了信号和噪声的分离效果,又有效地保护了地震信号中的AVO信息,对人工合成数据和实际资料的处理结果表明了所提出方法的有效性. 相似文献
14.
《地球物理学报》2020,(9)
本文在法方程层面融合GOCE卫星的V_(xx)、V_(yy)、V_(zz)和V_(xz)重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用IIR带通滤波器在5～41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明:引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型. 相似文献
15.
16.
Spectral analysis for recovering the Earth's gravity potential by satellite gravity gradients 
下载免费PDF全文
下载免费PDF全文 根据仪器功率谱密度和重力位系数阶方差的定义,本文建立了卫星重力梯度测量噪声功率谱密度与重力场模型的误差阶方差的直接对应关系,并基于此讨论了重力梯度测量精度、卫星轨道高度以及运行时间对地球重力场恢复精度的影响.相比于传统的基于最小二乘法评估卫星载荷噪声对地球重力场恢复精度的影响而言,本文提出的方法简单、直接,有助于快速设计和确定卫星重力测量计划的有关参数. 相似文献
17.
《地震工程与工程振动》2017,(3)
本文简要综述了时域连续Galerkin(TCG),改进TCG,以及时域间断Galerkin(TDG)以及高阶单步法等几种典型的高阶时域积分算法。若采用p阶多项式对基本未知量进行插值近似,采用单场格式的常规TCG算法、改进TCG算法、两场格式的TDG算法求得的位移和速度分别具有p、2p、2p+1阶收敛精度,在计算成本相当的情况下,改进TCG算法具有最高收敛精度。高阶单步法多为后处理的算法,大多具有4阶精度。相比最常用的二阶算法,高阶单步法操作简单,易于实施,保留二阶算法无条件稳定、计算效率高等优点,可对二阶算法进行误差估计,并可根据精度需求,指导时间步长自适应调整。 相似文献
18.
本文基于星间加速度法开展了插值公式、相关系数和采样间隔对GRACE Follow-On星间加速度精度影响的研究. 模拟结果表明:1)适当增加数值微分公式的插值点数可有效提高插值精度. 基于9点Newton插值公式,星间加速度的插值误差为4.401×10-13 m·s-2,分别基于7点、5点和3点插值公式,插值误差增加了1.192倍、6.912倍和274.029倍. 2)适当增大相关系数可有效降低星间加速度的误差. 基于相关系数0.99,星间加速度方差为3.777×10-24 m2·s-4,分别基于相关系数0.90、0.70、0.50和0.00,方差增加了9.780倍、22.404倍、26.217倍和26.820倍. 3)随着采样间隔增大,星间加速度方差逐渐降低,但卫星观测值的空间分辨率也同时降低,因此合理选取采样间隔有利于地球重力场精度的提高. 4)基于9点Newton插值公式、相关系数(K波段测量系统星间距离和星间速度0.85、GPS轨道位置和轨道速度0.95、星载加速度计非保守力0.90)和采样间隔10 s,利用预处理共轭梯度迭代法,精确和快速反演了120阶GRACE Follow-On地球重力场,在120阶处累计大地水准面精度为4.602×10-4 m. 相似文献
19.
本文在法方程层面融合GOCE卫星的Vxx、Vyy、Vzz和Vxz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型. 相似文献
20.