一种数值积分方案用于球函数分析试验.pdfVIP

一种数值积分方案用于球函数分析的试验 郭栋王盘兴房广洋 (南京信息工程大学大气科学系。南京210044) 摘要针对半球月平均位势高度场球函数分析的实际问题。分析比较了吴新元提出的数值积分公式与常用的梯形积分公式、辛卜生积分公 式在标准化缔合勒让德函数模、交角计算的中的误差。结果表明：吴新元数值积分公式明显减小了仨。的交角误差，它对提高月平位势高度场球函 数分析精度有重要意义。 关键词数值积分方案误差分析标准化缔合勒让德函数归—性和正交性 引言 球函数分析方法是分析全(半)球月平均位势高度气候场和异常场集的一种正交函数分解的方法【1-7。。对于经纬均匀格距 矩形格点网．在离散化过程中，经向的球函数正交性不严格成立；它是影响球函数分析精度的重要原因[sl。因此，选择一种 均匀格点域上精度较高的数值积分方案是有意义的工作。 本文用吴新元提出的一种均匀格点域上数值积分方案f9】，与梯形积分方案、辛卜生积分方案【10l用于经向的球函数—— 标准化缔合勒让德函数的归一性和正交性作了比较，论证了吴新元方法有明显的优越性。 1经向数值积分公式 设g(口)为定义在余纬0∈【o，180。】上的函数，为了计算定积分I AO=，rlJ；格点序数J与格点所在余纬的关系为： AO=jIr／J，j=0，J (1) q=j 梯形求积公式(I)【10l： sinoo+qs sin0，)+驴J-In0 fq(O)sinOdO=7一[j-(qo 因0=Oo、q时，sin0=0，简化得： g产n嘭 rq(目)s删p≈了,疗缶J-I 它具有一次代数精确度。 辛普生求积公式(II)： sinoo+qj ∽s…·寺{qo sinOg+4》一n”2J酚12-1sin幺，) 因0=Oo、岛时，sin0=0，简化得： fq(O)sin鲋臼“秀{4J》12一-s嘞一J静12-1咖吼，) 它具有三次代数精确度。 为提高数值积分的精度，吴新元在积分公式中考虑了被积函数的导数，提出了一个新积分方案(8l’这里称吴新元求积公式 (III)。其算式为： 49 fq(O)sin鲋口≈南{7c％s；n岛+乳s；n岛，+·4J蔷12-1％s；n岛，+，6d蔷12％一-s；nB川) sin岛+qo sin以+qjcos已)) +斋{(qo cosoo)一(gj 式中””为导数。因护=岛娩)时，sin护--0、cos秒=1，(-1)，简化得： fq(O)sin…南P和Is血”·6J蔷12q2j-。sin02j_,)+斋¨¨ 它具有五次代数精确度。 2比较方案 1／2=1 酌 rtl,／,i--JI即¨=(印，印)¨2={聒(印)2sinOdO 相同聊不同玎(胛≠∥’)即、印的交角 (9) ‰”．=arccos(印，印)=arccos{r印印sinod0}_90。 离散情况下，数值积分只能给出．磊一、r／m，露矿的近似值： f’=0即11．=|[印，印)‘P] 哟 刁+M∥～叫黼j ㈣ 式中符号“q’既表示数值积分的计算结果，又表示数值积分的计算算符。因为印的计算使用了欧洲中期数值预报中心 (EcwMF)的高精度的算法【101，其计算误差可忽略。故f’。，。、r／’m，”的误差来源于数值积分。 虽然球函数分析对象是定义域是整个球面(口∈【o，180。】，旯E【o，360。】)上的要素场，但是为了分析北、南半球地表的巨大 差异对大