更新时间:2024-05-21 14:04
摄动问题可分为正则摄动和奇异摄动两类形式。如果令 ε=0,Γε的表达式可化为Γ0,而且是一致有效的,就称这个摄动问题是正则摄动问题。如果在Sε中令ε=0会导致问题无解或多解,或者虽然当ε=0时Sε能化为s0并有解Γ0,但表达式Γε不一致有效,则称这个摄动问题为奇异摄动问题。
摄动方法(perturbation method)是把系统视为理想模型的参数或结构作了微小扰动的结果来研究其运动过程的数学方法。这种方法最早应用于天体力学,用来计算小天体对大天体运动的影响,后来广泛应用于物理学和力学的理论研究。摄动方法作为一般的数学方法,也是控制理论研究中的一种工具。摄动方法的基本思路是:如果一个系统Sε中包含有一个难以精确确定或作缓慢变化的参数ε,就可以令 ε=0,使系统Sε退化为s0,而把Sε看作是s0受到(由于ε≠0而引起的)摄动而形成的受扰系统。问题因而简化成为在求解S0的基础上来找出系统Sε的运动表达式。这样做往往能达到简化数学处理的目的。摄动方法所提供的系统Sε的运动Γε的形式是s的幂级数(可能包含负幂次项),级数的各项系数是有关变量(时间、状态变量等)的函数。如果在这些变量的容许变化范围内,当ε趋于零时,Γε的表达式一致地(均匀地)趋于S0的运动表达式Γ0,就称表达式Γε为一致有效的。
奇异摄动问题则复杂得多,当ε 趋于0时系统Sε的行为或结构往往发生本质的或剧烈的改变,出现各种复杂的现象。奇异摄动问题的研究已发展为控制理论的一个重要分支。其中常用的方法有伸缩坐标法、匹配渐近展开法、复合展开法、参数变易法、平均法、多重尺度法等。
对于弱非线性系统,若把非线性部分看作是对线性部分的摄动,常能用摄动方法(这种情况常称为小参数法)得到相当好的结果。奇异摄动理论与分岔理论、突变论等也有比较密切的关系。
研究天体在真实轨道上的坐标和在中间轨道上的坐标之差,这个差值称为坐标摄动。在经典方法中,常把坐标摄动表示为某个小参量(例如摄动行星的质量)的幂级数,然后逐项进行计算。由于计算技术的发展,微分方程近似解法中皮卡迭代法正逐步代替原来的小参量幂级数展开方法。它的主要优点是有统一的迭代过程,使计算过程能高度自动化。
这是1858年恩克在研究彗星的运动时提出的,它讨论坐标摄动在直角坐标系中的表示式,经常用于计算短周期彗星和月球火箭的轨道。这种方法的优点是:摄动方程的推导简单,形式对称,可以直接得到坐标,便于计算天体的历表。它的缺点是:以直角坐标表示的摄动量难于显示出摄动的几何特性和力学含义;随着时间跨度的增长,直接坐标的三个摄动量往往同时变大,以致不能把它们所服从的方程作线性化处理,否则就要多次更换零点。
自然天体一般总是围绕着某个主天体运动,例如行星绕着太阳运动,卫星绕着行星运动。因此,球坐标或极坐标的摄动就有较明显的几何意义。克莱洛和拉普拉斯在研究彗星的运动和大行星运动理论时最早提出了球坐标摄动方法。后来,纽康对拉普拉斯方法作了改进,特别是在展开摄动函数时运用了算符运算,使展开过程不仅有简洁的数学表示式,而且有规则的处理过程,便于以后在电子计算机上进行计算。纽康成功地运用这个方法研究了水星、金星、地球、火星四颗内行星以及天王星、海王星的运动,据此编成的内行星的历表,一直是二十世纪以来编算天文年历的基础。希尔提出了一种以真近点角为引数的球坐标摄动法,它曾被成功地用于计算第一号小行星──谷神星的摄动。
1963年穆森提出了另一种计算坐标摄动的方法,用于计算天体坐标在向径、速度和角动量三个方向上的摄动量。尽管这样的分解不正交,但由于它有不少优点,如有较明显的力学意义,推导方便,积分直接、运用算符运算、各阶摄动方程具有统一而紧凑的形式,并便于计算自动化,现正用于建立新的大行星运动理论。
这是轨道要素作为基本变量的摄动方法。如果行星只受太阳的吸引,正如开普勒定律所描述的,它将沿着一个固定的椭圆运动,决定椭圆运动的六个轨道要素应是常数。若考虑到其他因素的影响,行星将偏离原来的椭圆,六个轨道要素就不再是常数,它们将遵循由常数变易法导出的规律而变化。在这种情况下,可得到一族椭圆,它们逐个地与真实轨道相切,在相切点,二者不仅有相同的坐标,而且有相同的速度;只是加速度彼此不同,一个是真实加速度,另一个是椭圆加速度,二者之差正是摄动力引起的摄动加速度。由于种摄动加速度的作用,天体在下一时刻将离开这个椭圆,走上邻近的一个瞬时椭圆;相反,一旦摄动作用消失,天体将沿着消失点的瞬时椭圆一直运动下去。天体在太阳辐射压摄动下的运动正是这样:当辐射压起作用时,天体的瞬时椭圆不断变化;但当天体进入一个阳光照不到的阴影区时,辐射压消失,天体就沿着入影点的瞬时椭圆运动下去,直到跑出这个影子为止。 天体的真实轨道就是瞬时椭圆族的包络线。与坐标摄动相比,椭圆轨道要素的变化一般要缓慢得多,因而便于处理。瞬时椭圆法最早是欧拉在十八世纪中叶研究木星与土星的相互摄动时提出的,后由拉格朗日加以改进。他根据常数变易法,利用拉格朗日括号,严格地导出了描述椭圆轨道要素变化的摄动方程──拉格朗日方程。这种方法的应用十分广泛,特别是被勒威耶成功地用来研究大行星的运动。
关于正则摄动问题有下面几种常用的方法:
① 林斯泰特-庞加莱方法 它原用于天体力学。但这个方法对于消除一些问题中高阶近似的强奇性是无效的。1949年,M.J.莱特希尔作了重要推广,引进了自变量的非线性变换,求得了一系列物理问题的一致有效渐近解。
② 克雷洛夫-博戈柳博夫方法 由于取了函数的平均值,此法又称平均法。
③ 调和均衡法
常微分方程奇异摄动问题主要是大参数问题。下面列举两种类型问题。
第二类,在微分方程的系数中具有转向点的奇点问题。
近十多年来,研究奇异摄动问题的数值方法也有了发展,根据奇异摄动问题的特点构造特殊的数值方法(有限差分方法和有限元方法) ,以便取得良好的数值结果。
中国力学工作者在对摄动方法的发展有开创性贡献。钱伟长在1948年解圆板大挠度问题时,即提出称为合成展开法的方法,郭永怀在1953年把由庞加莱和莱特希尔发展起来的方法推广应用于边界层效应的粘性流问题;钱学森1956年又深入阐述了这个方法的重要性,并称之为PLK方法。
这是一种以分析力学为基础的方法。其基本思想是:对变量进行一系列适当的正则变换,以求降低运动方程的阶次,使新的方程具有较简单的形式,例如得出一个描述等速直线运动或简谐振动的方程,从而使问题得解。
十九世纪,德洛内从这个观点出发建立了著名的德洛内月球运动理论。他首先将月球的摄动函数展开成四百多个三角项,然后进行一系列的正则变换,使每次变换都能消去其中的一项。他花了差不多二十年的时间,总共进行了上千次变换,找到了三个合适的角速度,将月球的轨道要素都表示成时间的三角多项式,而不包含任何长期项。德洛内的工作为天体力学中的变换理论奠定了基础。这种方法是由一系列形式统一的循环过程组成的,因此非常便于用电子计算机进行计算。
德洛内之所以要进行那样多的变换,是为了对摄动函数中的每一项都给以严格的数学处理。这在实用上是没有必要的,某些高阶项尽可以略去。以这种想法为指导,蔡佩尔在二十世纪初建立了蔡佩尔变换。他先把摄动函数中的角变量按它们变化快慢排队,然后在一定精度范围内寻找适当的变换,以便一次消去所有含快变量的项,得出一组平均化的方程,进而对新的方程重复类似的过程,直至消去全部角变量为止。与德洛内方法相比,这种方法的工作量小得多,因此,它一出现就被成功地用来研究小行星的运动。人造卫星上天后,它得到了更广泛的应用。但是,蔡佩尔变换也有一些缺点,其中最突出的是:决定新旧变量转换关系的母函数是混合型的,同时含有新旧两种变量,使用颇为不便。为了克服这一缺点,堀源一郎在二十世纪六十年代提出了一种以李变换为基础的理论──堀源-李变换。其优点是:不仅新旧变量之间的变换具有显函数的形式,同时其结果在正则变换之下保持不变,因此它与用哪一组正则变量进行计算无关,而具有通用性。 电子计算机的创制和发展不仅大大提高数值计算的精度和速度,而且代替人们完成大量机械的重复的推导,今天已广泛用于摄动理论研究。德普里特、亨拉德、罗姆利用电子计算机编制了一个分析月球历表。单就计算太阳主要摄动项而言,摄动函数就有近3,000项,并通过李变换,得到了近50,000项月球坐标表示式。其规模之大,远非德洛内理论所能相比。
影响天体运动的摄动因素多种多样:有万有引力引起的保守力,有介质阻尼引起的耗散力,有连续作用的力,也有诸如辐射压引起的间断力等。影响大行星动的主要摄动因素是行星间的相互吸引;地球大气的阻尼使卫星陨落于地面;太阳辐射压决定着彗尾的形状;潮汐摩擦则是卫星轨道演化的主要动力。只有准确地掌握各种摄动因素,才能准确无误地计算天体的运动,解释各种壮丽的天象。反之,通过精密的观测和准确掌握天体的运动规律,就可以根据摄动理论的分析,弄清天体周围的力学环境,如测定摄动天体的质量、主天体的学扁率和弹性模量、大气密度和各种引力场参数等等,至还能预告一些未知天体的存在与行迹。因此,摄动理论不仅有丰富的理论内容,也有较高的实用价值。