多尺度法初识和应用.pdf

《多尺度法初识和应用.pdf》由会员分享，可在线阅读，更多相关《多尺度法初识和应用.pdf（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、多尺度法初识和应用 摘要：简要介绍多重尺度发的中心思想，另外，举例说明多重尺度法在求解方程中的应用。非线性问题的研究 非线性问题的“个性”很强，处理起来十分棘手。历史上曾有过一些解非 线性方程的“精品”，但与大量存在的非线性方程相比，只能算是“凤毛麟角”。因此，长期以来，对非线性问题的研究一直分散在自然科学和技术科学的各个 领域。本世纪六十年代以来，情况发生了变化。人们几乎同时从非线性系统的 两个极端方向取得了突破：一方面从可积系统的一端，即从研究多自由度的非 线性偏微分方程的一端获得重大进展。如在浅水波方程中发现了“孤子”，发 展起一套系统的数学方法，如反散射法，贝克隆变换等，对一些类型的非

2、线性 方程给出了解法；另一方面，从不可积系统的极端，如在天文学、生态学等领 域对一些看起来相当简单的不可积系统的研究，都发现了确定性系统中存在着 对初值极为敏感的复杂运动。促成这种变化的一个重要原因十计算机的出现和 广泛应用。科学家们以计算机为手段，勇敢地探索那些过去不能用解析方法处 理的非线性问题，从中发掘出规律性的认识，并打破了原有的学科界限，从共性、普适性方面来探讨非线性系统的行为。线性与非线性的意义“线性”与“非线性”是两个数学名词。所谓“线性”是指两个量之间所存在 的正比关系。若在直角坐标系上画出来，则是一条直线。由线性函数关系描述的系 统叫线性系统。在线性系统中，部分之和等于整体。

3、描述线性系统的方程遵从叠加原理，即方程的不同解加起来仍然是原方程的解。这是线性系统最本质的特征之一。“非线性”是指两个量之间的关系不是“直线”关系，在直角坐标系中呈一条曲线。最简单的非线性函数是一元二次方程即抛物线方程。简单地说，一切不是一次的函 数关系，如一切高于一次方的多项式函数关系，都是非线性的。由非线性函数关系 描述的系统称为非线性系统。多尺度法的基本思想 多尺度法首先是由 Sturrock(1957)、Cole(1963)、Nayfeh(1965)等提出的,此后得到进一步的发展。上面介绍该法的基本思想与方法。我们考虑形式为 的方程所控制的系统，设方程的解为 将原点移至中心位置 是合适

4、的。于是有 0qfq 22100 xxqxqq0qq 此时第一式可写成 假设 f 可以展为泰勒级数，则上式可写为 其中 而 f(n)表示关于自变量的 n 阶导数，对于中心，而 我们可以把方程的解看成是多个自变量的函数，而不是一个自变量的函数。也就是们可以把x看成是 t 和 ,的函数。多尺度法的基本思想是，将表示响应的展开式考虑成为多个自变量(或多个尺度)的函数。即 因此关于t的导数变成了关于的偏导数的展开式，即 然后代入方程进行求解，求出 。这时，方程的 解可写成：然后按照小参数法(摄动法)建立 的 各阶方程，进而 求出 多重尺度法的应用 一、求解自治系统 例 1.4.1 求 Duffing

5、方程（1.1.4）0qqx00qxfx Nnnnxax10 01qfnann！00qf 0qfnt,2,1,0ntTnntTtTtT2210101100DDddddddTtTTtTt20212102022D2DDDD2Dddt,321xxx,321xxx30(1)xxx 自由振动的二次近似解（用多尺度法）解：求二次近似解可选三个变量，设 2001210122012(,)(,)(,)xx T T Tx T T Tx T T T 代入原方程，并用到式（1.4.3），可得到下列方程组 200200 xxT （1.4.4a）22311020012xxxTT T (1.4.4b)222232220122

6、001021223xxxx xTT TT TT (1.4.4c)设式（1.4.4a）的解为 01200(,)exp()exp()xA T TiTAiT 其中 A 是未知复函数，A是 A 的共轭。用复数形式表示是为了运算方便。把 0 x代入式（1.4.4b）223110020123exp()exp(3)xAxiA AiTAi TccTT 其中 cc 表示前面各项的共轭。为使 x1,不出现永年项，必须 21230AiA AT （1.4.4d）又求得 3101exp(3)8xAi Tcc 把01,xx代入（1.4.4c）,并利用条件（1.4.4d），有 2324522000202152132exp(

7、)exp(3)exp(5)888xAxiA AiTA Ai TAiTccTT 消除永年项 32215208AiA AT （1.4.4e）解2x为 45200211exp(3)exp(5)6464xA Ai TAiTcc 利用式（1.4.4d）,（1.4.4e）求 A（T1,T2）如下：由（1.4.4d）2132AiAAT 由（1.4.4c）2321516AiA AT 利用式（1.4.3a）并注意到00AT，就得到 223315216dAiA AiA Adt 令1exp()2Aai，其中,a是 t 的实函数，将之代入上式，实、虚部展开，有 0a 2243158256aa 积分得 0aa 2240

8、3158256aat 00,a为积分常数，所以 224001315exp()28256Aaiaa ti 于是，原方程二阶近似解为 322500001211cos(1)cos3cos532321024xaaaa 其中 2240315(1)8256aa t 二、无限传输方程的近似解 （一）稳定性分析 对于系统 ()()()()()x tx tx tx tf x t （2.1.1）对于方程（2.1.1）的根0 x,如果对0 x的任一邻域 U,存在0 x的一个属于 U 的邻域1U，使系统（2.1.1）的解()x t，若有01xU,则对一切0t，有()x tU，就称0 x是稳定的，否则就称为不稳定的。如

9、果0 x稳定，并且有0()limtx tx，就称0 x是渐近稳定的。定义：若（2.1.1）的零解对 都是渐近稳定的。则称（2.1.1）为全时滞稳定的。或叫无条件稳定或绝对稳定。可求(2.1)的特征方程：将txce代人到方程（2.1.1）中则有，()tx tce()()tx tce ()()tx tc e 所以有：()()0tttc ec ece 即有：0ee （2.1.2）1ee 若0时，则1为其特征根。如果其特征根位于左半平面，而当由 0 增至时，不越过虚轴，则系统（2.1.2）的更全具有负实部，这样系统（2.1）的零解为全时滞稳定的。因此，要使（2.1.1）为全时滞稳定，首先要使（2.1.

10、2）的根具有负实部。只有当(2.1.1)的特征根为纯虚数时，方程的解才有近似周期解。用i代人（2.1.1）中，有 0iiiiee 即(cossin)cossin0iiii 所以有 cossin0sincos0 令22()(1cos)cosf 当1cos0时，在区间上0,2上，22()2(1cos)sinsin0f 函数 f单调 当 0时，2()(0)0ff 当2时，22()()024ff 函数与 X 轴有交点，方程有解，即 特征方程（2.1.2）有纯虚根。（二）近似周期解 在3x的非线性扰动的情况下，可求系统的一次近似周期解（利用多尺度法）设2001210122012()(,)(,)(,)x

11、tx T T Tx T T Tx T T T （2.2.1）其中2012,nnTt Tt TtTt 应用微分算子，记00DT，11DT，知：2201010()0()dDDdtTT （2.2.2）由001101()(,)(,)x txTTx TT20()，知 001101()(,)(,)x tx TTx TT20()(2.2.3)根据二元函数的泰勒展开：00(,)f xh yk 0000(,)()(,)f xyhkf xyxy 令0010(,0,)Tx hTyk 知 100100110(,)(,)(0)Txx TTx TTxTT 0011(,)xxTTT 1011 1(,)x TTDx （2.2

12、.4）110110110(,)(,)(0)Txx TTx TTxTT 1011(,)xxTTT 1011 1(,)x TTD x （2.2.5）将（2.2.4），（2.2.5）代人(2.2.3)中 得到时滞项：2001101()(,)(,)0()x tx TTx TT 10111(,)x TTD x+210111(,)x TTD x+20()20011011001(,)(,)(,)0()x TTx TTD x TT（2.2.6）3301()()x txx 32001001101(,)3(,)(,)xT TxT T x T T 22330011011013(,)(,)(,)x T Tx T Tx

13、 T T（2.2.7）223000111010012()xxxxxxx tTTTTTT （2.2.8）将（2.2.1）（2.2.2）(2.2.3)（2.2.4）（2.2.5）（2.2.7）（2.2.8）代人原方程得 ()()()x tx tx t 20 0011 0010 0011 101(,)(,)(,)(,)D x T TDx T TD x T TDx T T 200011 0010 1011 101(,)(,)(,)(,)D x T TD x T TD x T TD x T T0011011 001(,)(,)(,)x TTx TTD x TT 32 23 24 300100110100

14、1101101(,)3(,)(,)3(,)(,)(,)x T Tx T T x T Tx T Tx T Tx T T 这样根据多项式的性质，可知，指数012,的系数在等式两边相等。这样就有，00 0010 001001:(,)(,)(,)0D x T TD x T Tx TT （2.2.9）则，当(,)a bD时，系统可形如（2.1.1），这样0i是特征方程的根。易见方程（2.2.9）有如下形式的谐波解：000011(,)()T ixTTA Tecc 其中 cc 表示前面各项的共轭，000000111(,)()()T iT ixTTA TeA Te 0 00 00 032332001111(,

15、)()3()()T iT iT ixT TAT eAT eA T e 000000223231113()()()T iT iT iA TeATeATe 0 00 00 00 0333223111111()3()()3()()()T iT iT iT iA T eA T A T eAT A T eA T e 11 0010 1011 0010 101:(,)(,)(,)(,)Dx T TD x T TDx T TD x T T 1011001(,)(,)x TTD xTT 3001(,)xTT 又有，000001001111(,)T iT ixAAD xTTeeTTT 这样，0101010110

16、1(,)(,)(,)D xTTD xTTxTT 3100110011001001(,)(,)(,)(,)D x T TD x T TD x TTx T T 0 000 000 000 033111111()T iiT iiT iiT iAAAAAeeeeeeA T eTTTTT=0 000 000 0011111T iiT iiT iiAAAAAeeeeeeTTTTT 0 00 00 00 0333223111111()3()()3()()()T iT iT iT iA T eA T A TeA T A T eA T e 0002111113()()iT iAAAeATA TeTTT 0 00

17、 00 0332331111113()()()()T iT iT iAAA T A TeA T eA T eTT 而对齐次方程0 1010 101101(,)(,)(,)0D x T TD x T Tx TT的特征方程有：00000wiwiw iw iee 得，000()wiwiewi 为此，我们可以设 1()111()()2ib TA Ta Te 可令11,ababDDTT 这样，11()()111111()22ibTibTAabea TeiTTT 11()()1122ib Tib TabD eaD ie 1()1()2ibTabeDaDi 由于所求的为方程的近似周期解，所以其永年项为 0.

18、则，02111113()()0iAAAeATA TTTT 即，01()2111()13()()02iib TabeDaD ieA T A T 而，11()()111111()()()()22ib Tib TA T A Ta Tea Te214a 这样有，011()()2111()1 30242iib Tib TabeDaDieaae 即，033()104iabDaD iea 22000aaaaaDwiDD wiDDwi 22000bbbbbaDwaDi a wDaDi awD 33033044aw ia 分离实部和虚部 得 2230002230000304304aabbbaaabbDDaD w

19、a w Daw DaDwDwDwaDaDaw 根据克拉默法则解方程组，得 320320220220003434aaaaawaaDaawwwaa 32320000220220003434baawwwwDaawwwaa 把回代，因此有：0 1010 101101(,)(,)(,)D x T TD x T Tx TT0 033T iA ecc（注：其余那些项为永年项为零）因此，1x应有形如0033TicAe形式的周期解。将0033101(,)T ixTTcA e回代，则有 003330033iiciciceA e 0330033iAciie 所以方程有如下形式的周期解：0033101300(,)33iiA ex T Tiie 结论：当0时，当1cos0时，在区间上0,2上，22()2(1cos)sinsin0f 函数 f单调 当 0时，2()(0)0ff 当2时，22()()024ff 函数与 X 轴有交点，方程有解，即 特征方程（2.1.2）有纯虚根。参考文献 1陈予恕，唐云.非线性动力学中的现在分析方法 第 1 版 北京：科学出版 1992 2张锦炎，冯贝叶.常微分方程几何理论与分支问题 第二次修订本 北京大学出版社，1997 3魏俊杰，中立型微分方程零解的稳定性征与全局 Hopf 分支 J.数学学报，2002，45（1）：93-104.4 多尺度法论文 郑国金