常微分方程值问题的数值解法上机报告.docxVIP

常微分方程初值问题的数值解法上机报告（一）问题：考虑著名的Lorenz方程其中为变化区域有一定限制的实参数。（1）对取定的参数值，选取不同的初值，观察计算的结果有什么特点？解的曲线是否有界，是不是周期的或者趋于某个固定的点？（2）在问题允许的范围内适当改变其中的参数值，再选取不同的初值，观察并记录计算的结果有什么特点？是否发现什么不同的现象？（二）解决问题的算法考虑一阶常微分方程组初值问题其中可以用四级四阶古典Runge-Kutta方法求解：（三）使用的软件IDL（四）数值结果及结果分析取h=0.001s显然，当初值取为(0,0,0)T时，x,y,z恒为0.(1)取定参量值仅对初值的x分量作一个小扰动，取为(0.0001,0,0)T，此时求解该常微分方程组，可以得到如下的轨迹：解的曲线有界，但是曲线不是周期的，而且不会趋于某个固定的点，而是在两个值附近振荡。仅对初值的y分量作一个小扰动，取为(0,0.0001,0)T，可以得到如下轨迹：可以得到和上面相同的结论。仅改变初值的z分量，取为(0,0,20)T，可以得到如下轨迹：可以看到只改变z方向的初值，结果为一条直线。这是由于根据方程可知，当x=0,y=0时，x，y不随时间变化，仅z随时间变化。下面改变方程的初值来求解。当初值取为(-5,-5,-5)T时，可以得到如下轨迹：当初值取为(200,200,200)T时，可以得到如下轨迹：将