综合评价线性多步积分公式轨道积分性能的两项新指标.pdfVIP

维普资讯 中 国 科 学 院 上 海 天 文 台 年 刊 2O02年第23期 ANNAI．S OF SHANGHAI OBSERVATORY ACADEMIA SINICA No．23，2O02 综合评价线性多步积分公式 轨道积分性能的两项新指标 徐 继 鸿 (中国科学院国家天文台，北京 100012) 提 要 概述 了制约线性多步积分公式轨道积分状态的多种因素；提出了综合评价线性多步积分公式 积分性能的两项新指标。建议在对数值计算有较高精度要求 的科研项 目中，应将构造并选择适合 研究项 目的线性多步积分公式以及高效的积分方式列为课题前期工作的重要部分。 主题词 ：线性多步积分公式 一 轨道积分 ～ 技术指标 分类号 ：P135 ． 1 引 言 考察一个线性多步积分公式的优劣 ，通常采用以下四项技术指标 ： (1)局部截断误差 ：简言之 ，就是把微分方程的理论解代人线性多步积分公式后公式的两 端之差，误差常数是其表征。采用线性多步积分方法积分天体轨道时的积分误差 ，通常按轨道 的切向与法向分解为沿迹误差和垂迹误差两个分量。相对而言，沿迹误差 的数值总是要比垂 迹误差大得多(这是因为沿迹误差 A2~．AnT，在 △ 中，轨道平均角速度 rt的误差 An将被时 间 放大)，因而 ，通常讨论沿迹误差。线性多步积分方法的沿迹误差又都源于两个方面，它 们分别是截断误差和舍人误差。在积分过程中由数值存储时的四舍五人会造成舍人误差，舍 人误差除与数值积分公式、积分方式有关外，还依赖于计算机的字长。积分公式的中途截断会 带来截断误差。截断误差是用以判断该线性多步积分公式精确度的一个重要指标。假设前面 各步无截断误差 ，局部截断误差 + = +，()(p+1(t)+D( )，其中， 为积分步长 ， P为阶数 ，cD+是与积分系数有关的常数(即误差常数)，D是余项略写。假设前面各步也有 截断误差 ，全局截断误差 e+ = (t+)一 + ，其 中， (t+)是理论解 ； +是数值解。全 局截断误差 的上限为 Ie+ I≤P {A +GyM(t一口)}exp[P 凹 (t一口)]。其中，A，日， L，P ，G等是与积分系数有关 的常数 ，t一口=nh，yg I “(t)I在 [口，b]区间 内的最 大值。 2o01年5月 22日收到 ． 维普资讯 第 23期 综合评价线性多步积分公式轨道积分性能的两项新指标 35 (2)数值稳定性 ：它取决于该积分公式在逐步的积分过程中误差传递与累积的方式，积分 间隔越长，其表现将越充分。既然均方差数值的大小 ，在客观上反映了一个积分器的误差累积 与扩散的状况 ，因而 ，均方差随时间的变化趋势也就表征了该积分器的数值稳定性。 (3)计算的复杂性 ：表示一种积分公式逐步积分过程的复杂程度 ，也即每一步积分所需计 算量的大小 ，用于衡量该线性多步积分公式的工作效率。 (4)方法的普适性 ：反映该线性多步积分公式对不同类型轨道的适应能力。 然而 ，为了对一个线性多步积分公式在轨道积分时的综合性能作出确切地评价 ，仅凭这四 条是不够的，必须补充考虑新的技术指标。 在本文的讨论 中，坐标系、单位制、比对标准、编译系统等的选取均同文献 [5](勿再赘述)。 初值如下 ： i0 =0．444317792473； 0=0．77； 口0= 1．3021148741； 0= 3．56388560259； =0．053903982449；0=一0．00264108192654； (e0 0．054；c￡J0—0．049)。 2 反映线性多步积分公式轨道积分性能的两个特征 任何一种线性多步积分方法的计算精度都必然要受到积分步长、积分间隔、被积轨道 的类 型 (有无耗散力、偏心率大小)