使用LP法和同伦摄动法求解van der pol方程作业.docxVIP

求解van der pol方程 （1）LP法 对比可知 令，则 原方程化简为： 取 将(2)、(3)式代入方程(1)，整理得： 又因为 的同次幂项系数相等，可得： 初始条件， 解方程（4）中第一式并代入初始条件得： 代入(4)的第二式，得： 消去久期项得： 于是有 ， 初始条件，解得： 于是方程(4)的第三式可整理得： 消去久期项，令上式中和的系数为零，求得： ， 于是有 整理可得最终解： 其中 （2）同伦分析法 令 则 （1） 即 物理上，动力系统极限环可用周期函数描述，即 （2） 选取为初始猜测解。 根据（2）取辅助线性算子： 定义非线性算子： 构造零阶形变方程： 满足条件： 其中 当时， 当时， 假定和的选取合适，使得 于是 假定和的选取合适，上述级数收敛，于是 将零阶形变方程 对q求n次导数，再除以，最后令q=0。得高阶形变方程。 当n=1时 消去久期项得 为简便，取得 边界条件：，得 当n=2时 代入 代入 代入 代入 化简得 消去久期项得 ... 我们仅需保证级数： 收敛。 当时，若则上述级数收敛。 当时， 所以 结果和L-P法解得的结果一样。 （3）作图 使用Matlab编程得到相轨线图如下： 由图1,2,3可以看出，LP法与RK法之间的误差比较小，随着小参数的增长，二者之间的偏差越来越大。 图1，的相轨迹图，实线为LP法，虚线是RK法 图2，的相轨迹图，实线为LP法，虚线是RK法 图3，的相轨迹图，实线为LP法，虚线是RK法