基于主成分分析的相对重力观测零漂和固体潮提取.pdfVIP

第 35卷 第 5期 大 地 测 量 与 地 球 动 力 学 VoI．35No．5 2015年 10月 JournalofGeodesyandGeodynamics Oct．。20l5 文章编号 2015)05—090605 基于主成分分析的相对重力观测 零漂和固体潮提取 于红娟 郭金运 。 李九龙 穆大鹏。 孔巧丽 。 1 山东科技大学测绘科学与工程学院，青岛市前湾港路579号 ，266590 2 海岛(礁)测绘技术国家测绘地理信息局重点实验室，青岛市前湾港路 579号，266590 3 中国科学院测量与地球物理研究所 ，武汉市徐东大街 340号 ，430077 摘 要 ：利用主成分分析 PCA提取 2014-08—15～23CG-5相对重力仪静态相对观测的零 点漂移和 固体潮 。 结果显示 ，PCA提取 的零点漂移值与最小二乘线性拟合结果非常一致 ，其差仅在 10 ~Gal／d量级 ；PCA估 算 的固体潮值与 CG-5重力仪 内部软件提供 的固体潮模型计算值统计结果基本一致 ，两者之差均小于 8 uGal，均方根均 小于 5uGal。 关键词 ：PCA；零点漂移 ；固体潮 ；相对重力测量 ；C 5重力仪 中图分类号 ：P223 文献标识码 ：A 本文通过实验 ，利用主成分分析 (principal 阵计算得到的线性变换矩阵。假设协方差矩阵为 componentanalysis，PCA)提取CG-5重力仪 口静 ，则它是一个 阶对称 阵。通过求解特征根和 态相对观测的零点漂移和 固体潮。将经过固体潮 特征 向量 ，矩阵 邑 可 以分解为 ： 改正的9d观测数据作为原统计量 ，使用 PCA提 邑 一 EDE (2) 取Lt零漂，并与最小二乘线性拟合值作对 比分析 ， 式中，D是对角阵，对角线上的元素 ≥ 。…≥ 用含有固体潮和零漂影 响的观测数据 ，在进行零 是矩阵 c的特征根 ，同时等于对应重力观测数据 点漂移改正的基础上使用 PCA提取固体潮 ，并与 主成分 的方差 ，E 是 由特征 向量构成 的正交矩 阵_3]。因此，可 由式 (1)得到 ： CG-5重力仪内部软件提供 的固体潮模型计算值 G — PCA(；E (3) 作对 比。通过两组实验对比结果来验证本文方法 前 k个主成分的综合能力以及所包含的重力 的可靠性和适用性 。 值信息量可用累积贡献率 a表示 ： 1 主成分分析方法 P a 一 ∑ ／∑ ×100G (4) i= 】 i 重力观测的主成分分析_l2 是基于协方差矩 式 中，表示第 i个特征值 ，P表示特征值总个数 ， 阵将重力数据作为原统计变量 ，通过线性变换得 k表示选取的个数 。 到新变量 ，新变量之间统计不相关 。PCA 的主要