大地测量混合边值问题的自然边界元解法武汉理工大学学报交通.PDF

交通科学 第 33 卷 　第 1 期 武 汉 理 工 大 学 学 报 ( 与工程版) Vol . 33 　No . 1 2009 年 2 月 J our nal of Wuhan U niver sit y of Technolo gy Feb . 2009 ( Tran sport ation Science Engineering) 大地测量混合边值问题的自然边界元解法 唐元义1 ,2 ,3) 　许厚泽1) ( 1) ) 中国科学院测量与地球物理研究所 　武汉 　430077 ( 2) ) ( 3) ) 湖北大学数学与计算机科学学院 　武汉 　430062 　 中国科学院研究生院 　北京 　100049 摘要 :在物理大地测量边值问题的球近似下 ,将大地测量混合边值问题转化为 Neumann 外问题 , 用自然边界元法 ,对 Neumann 外问题进行 自然边界归化 ,得到 自然边界积分方程 ,再求 自然边界 积分方程的数值解. 在球界面下 ,利用 自然边界元方法 ,将大地测量混合边值问题的求解 ,转化为 线性方程组的求解 ,而线性方程组的解存在惟一稳定 , 因而大地测量混合边值问题的解存在惟一 稳定. 关键词 :大地测量混合边值问题 ; 自然边界元法 ;数值解 中图法分类号 :P22 DOI :10 . 3963/ j . i ssn .. 2009 . 0 1. 03 1 0 　引　　言 为正常重力; N 为利用卫星测高数据得到的海洋 上大地水准面高; S 1 和 S2 分别为地球平均球面 结合海面上卫星测高数据和陆地上重力观测 上对应于陆地和海洋的 2 个部分, 即: S = S 1 ∪S2 数据 , 以确定地球外部重力场和陆地形状为 目标 , = { r = R} , 且 S 1 ∩S2 = . 而产生了大地测量混合边值问题 ,通常也称为混 T Δ 合边值问题. 在球近似下 ,大地测量混合边值问题 将上述 Sto ke s 边界条件改写为 r = - g 的典型形式表现为下列数学模型[ 1 ] . - 2 T , Sto ke s 边值条件转化为 N eumann 边值 Sto ke sDirichlet 混合边值问题 R ΔT = 0 　　　在球外 Sc 条件, 大地测量混合边值问题改写为 T 2 Sto ke sDirichlet 混合边值问题 Δ r + R T | S 1 = - g ΔT = 0 　　　　在球外 Sc