计算动力学第二章辩析.ppt

将式（2)、(3)和（4)代入式(1)，比较等式两边的 同次幂的系数，则可得下列方程组： 由式(5)可解得 把式（8）代入式（6）得： 设初始条件为t＝0时，x(0)=2和 则 得： 为了消除永年项而得到周期解，必须使sin?和cos?的系、 数都等于零，由此可得： 因比方程(9)的解为； 将式(8)、 (10)和(11)代入式(7)，经简化和整理 后，可得， 根据周期性条件，令cos? 的系数等于零，得 把上式代入式(12)，解之，并使它满足初始条件： 可得： 所以，方程(1)的第二次近似解为 其中p=p0+?p1+?2p2=1-?2/16。 方程(1)的稳定振幅为|a|=2，这对应于相平面 上的极限环，读者可自行与上一章几何定性法对比。 [题4—46]试用多尺度法求上题的近似周期解。 [解]仅取两个变量，设 并令 将式(1)和(2)代入范德波方程，并令?0 和?1 的系 数均 等于零，可得： 式(3)的解为， 把上式代入式(4)，得 为了消除永年项，应令： 或 用a乘上式两边，可将上式改写为： 由上式进行分离变量，可得: 设T1=0(即t＝0)时，a=a0，积分上式,可得： 或 由此式可解得： 消除永年项后 ，式(6)的解为： 或范德波方程的一次近似解为： 故范德披方程的零次近似解为： 由式(9)可见，当t时，t→∞时，a（?t）→2，此即为稳 定振幅。这与上题的结果一致。 其中x1(a，?)，x2(a，?)，…是?的以2?为周期的周 期函数。 根据对误差的分析，第一次近似应取为式(4-6)、 (4-7)已如前所述。同理第二次近似则应取为 如此类推。第m次近似在式(4-42)中应取至em-1的项 在式(4-43)的两个式子中均应取至?m的项。确定这 些展开式中的函数在原则上并没有困难，但由于演 算随着近似次数的提高而很快越来越复杂，在实际 上可以进行的只是头二三次近似。这个方法的实用 性不取决于近似次数无限提高时级数的收敛性，而 取决于确定的某一次近似当?→0时的渐近性。 为确定这些展开式中的函数，把式(4-41)对 时间求导数： 利用式(4-43)可以计算下列的量： 把式(4-43，4-47)代入式(4-46)，并且把所得的结 果按小参数?的幂次排列，得 于是，式(4-41)的左端可写为 式(4-41)右端在利用了式(4-42)、(4-48)后可写为 令式(4-49)和式(4-50)的?同幂次系数相等即可得到一系列为求解x1，x2，…，Al，A2，…，P1，P2，…所需的微分方程。而为了使表达式(4-42)以o(?m+1)的精度满足式(4-41)，上述的微分方程须列出m个： 其中 不难看出， fk（a，?）是变量?以2?为周期的周期 函数，又依赖于a0，在求出Aj（a），pj（a）， xj（a），j=1，2…，k之后，即可知道 fk（a，?） 的明显表达式。 首先从式(4-51)第一式解A1（a），p1（a）， x1 （a，? ）。为此研究函数fk（a，?） 和 x1（a，?）的Fourier级数x1 （a， ? ）。为此 研究函数fk（a，?） 和x1（a，?）的Fourier级数。 把本式代人式(4-51)第一式得 令此式相同谐波的系数相等，得 这样，A1(a)和p1(a)以及x1(a，?)展开式中除n＝1 之外的谐波项都已确定. 至于xl展开式中的n＝1的谐波即基本谐波已可 看出根据上面的运算是确定不了的；一般说，基本 谐波的确定带有任意性。因此只要规定x1(a,?) 不 包含基本谐波，就能定出 对x1以至x2，x3，…作这样的规定也使得x的解(4-42) 更加谐调，因为现在基本谐波已全部包含在首项 acos?中，后面的渐近级数 ?x1，?2x2,…内已无需 包含基本谐波了。 于是得 若进行更高一次近似的计算，就把所得的x1， A1，p1代入式(4-52)第二式以得到f1(a，?)的明显 表达式，然后把它也展开成 Fourier级数： 对物x2（a，?）也作不包含基本谐波的规定即 后，即可利用式(4-51)第二式确定A2（a），p2（a） 和x2（a，?）： 如此类推可以解出任意高次近似。以上的运算 是对于普遍情 的表达方式，它显得相当繁琐。对于 具体的 ，运算往往是简单的。作为一个 例子，现计算一个保守系统 的第一和第二近似解。现在 首先计算 式(4-52）第一式： 实际上现在已完成对f0(a，?）的Fourier展开即 从式(4-5