位场球谐分析的基本理论.docVIP

PAGE 40 2007.5 2007.5 PAGE 41  第二章 位场球谐分析的基本理论 球谐分析是卫星重力和磁场数据分析解释主要工具。地球外部的重力场可以表示成球谐级数形式，同时，通过不同阶次的球谐展开，可以对地球重力场进行分析，以达到显示地球重力场特征，进而研究地球重力异常各种成因的目的。本章主要介绍论文所涉及的球谐分析的一些基本理论。 §2.1 位场拉普拉斯方程的解 可以证明，地球外部引力场是调和的，其满足拉普拉斯方程。在不同的坐标系中，拉普拉斯方程有不同的形式，为了便于讨论，本节分别对直角坐标系和球坐标系中拉普拉斯方程的求解进行讨论。 直角坐标系中拉普拉斯方程的解 假设V表示地球引力位，其为直角坐标系中空间点(x, y, z)的调和函数，则有?2V=0，即 (2-1) 用分离变量法解方程，令 代入方程（2-1）则有 式中X?, Y?, Z?分别为X, Y, Z对x, y, z的二次导数。 解方程，得V(x, y, z)的一般表达式 其中为待定常数，根据边界条件来确定。 球坐标系中拉普拉斯方程的解 在球坐标系中，引力位V可表示为空间点(r, ?, ? )的函数，即V(r, ?, ? )，其中r为点的坐标径向距离，? 为余纬度，? 为经度，如图2-1所示。引力位V的拉普拉斯方程可表示成 (2-2) 图2-1 球坐标系 同理，可以用分离变量法解方程(2-2)，即令 代入方程(2-2)，则有 (2-2a) 上式中只第三项与有关，则第三项是一常数。 令 (2-3) 其中m是整数，可解得 代入方程(2-2)，则有 (2-4) 从上式可以看出，左边第一项只与有关，第二项只与有关，则两项都是常数。 令 (2-5) 这里n为整数。 解方程(2-5)，可得 当r→ ? 时，V有限，可知C=0. 将(2-4)和(2-5)代入方程(2-2a) (2-6) 方程(2-6)解的形式为 其中为连带勒让德函数，其中m是整数，是常参数。 当m, n取不同的整数值时，方程(2-3)和(2-5)的解是特解，考虑到位场的叠加性质，将所有这些特解累加起来，的一般表达式可写成 (2-7) 其中和是常参数。 可以证明，地球外部引力位的球谐表达式可以写成 (2-8) 其中和是通常所称的球谐系数。 §2.2 勒让德函数与连带勒让德函数 勒让德与连带勒让德函数的一般形式及其递推公式 勒让德函数球谐函数的核心组成部分。根据上一节方程(2-6)，把解代入，可得 当m＝0时，把x = cos? 代入，可得到勒让德方程，即 （2-9） 解方程(2-9)，可得勒让德函数表达式，即 显然，勒让德函数的级数表达形式为 (2-10) 其中为中最大的整数值。 如图2-2所示，勒让德函数具有以下性质： (1) 当是偶数时是偶函数，是奇数时是奇函数; (2) 在[–1，+1]的区域内有个零点。 图2-2 (a) 函数P0(x)~P6(x)在区间[-1,1]上的图形；(b) 函数P6(x)在球面上的图形，其中x=cos? ，灰色和白色分别代表正负区域。 当m ?0时，把x = cos? 代入，方程(2-9)变成连带勒让德方程，即 解方程可得连带勒让德函数微分表达式为 连带勒让德函数级数展开形式可写成 (2-11) 利用(2-10)及(2-11)式可以计算勒让德函数。然而，计算高阶勒让德多项式时不方便的，需要导出其递推公式。 由 (2-12) (2-13) (2-14)