应用拟合推估法建立我国大陆地壳水平运动速度场地研究.pdfVIP

应 用 拟 合 推 估 法 建 立 我 国 大 陆 地 壳 水 平 运 动 速 度 场 的 研 究 中 图 分 类 号 = P 2 2 文 献 标 识 码 = B 曾安敏123刘光明23秦显平23唐颖哲23 1.信息工程大学 地理空间信息学院，河南 郑州 450052；2.地理信息工程国家重点实验室， 陕西 西安 710054；3.西安测绘研究所，陕西 西安 710054 摘要：从随机信号协方差函数的建立着手，采用 “地壳运动观测网络工程”直到2008年的观测数据获得的实测速度作为信号协方 差函数的观测数据，建立与实测速度的拟合程度较高的协方差函数，将欧拉矢量模型作为地壳运动的趋势项模型，将局部变化信 息作为信号项模型，采用趋势变化与随机变化相结合的拟合推估法建立中国大陆水平运动模型，提高中国大陆地壳水平运动速度 场的计算精度。 关键词：拟合推估；运动模型；随机信号；协方差函数 12 测绘通报 2012年增刊 即存在互协方差矩阵瓦，=最，。≠0，则拟合推估 4 (10) 解为 嘲锄_【|：】 囊=(A’P￡A)。1AT吒工1 四、协方差函数的拟合 §=磊B’P￡(L—A宕)I (3) §|=-Ys，sz；1§ ＼ 在拟合推估中，关键的问题是在数据处理前已 式中，P￡=(旌。曰7+置)～。 知协方差矩阵最、三驱，，而协方差矩阵中的元素是根 假如事先去除趋向性参数部分，仅保留随机信 据协方差函数计算的。基于地壳连续变形的假设， 号部分，则式(1)改写为 在空间上很靠近的两个站点的速度值必然相关性 三=L—AX=BS+P(4) 强，反之则弱。可以设地壳运动场两点间的协方差 其相应的误差方程也改写为 为距离的连续函数，随机信号协方差函数可采用 V=曰S一乞 (5) 构造如下极值条件 Hirvonen函数为 (11) (6) C(d)=c0／1+k·扩 ．0=矿嚣1y+§■；1§=min 如果同时考虑未测点信号§’，保留雪’与雪相关 协方差函数的Gauss指数函数为 性，则式(6)的拟合推估解为 ‘C(d)=Co·e一脚 (12) V=置(互s+置)叫L， 式中，Co、k为待定常数；C(d)为i、，两点之间的协 方差；d为i√两点之间的距离。 (7) §=三sB’(BYsB7+乏)。1三I 一般来讲，信号的初始方差可按下式计算 §’=乏s．s=；l§ ＼ 通常情况下B=I，如果不考虑观测噪声，即毛 ‘ C(d)=土y弧1 =0，则有 一，一。}．J(13)