数值分析实验报告模板精选 .doc

数值分析实验报告模板 常微分方程数据值解:用龙格—库塔法分析Lorenz方程的特性 一、 问题的提出 考虑著名的Lorenz方程 ?dx?dt?s ??dy??rx?y?xz ?dt ?dz?dt?xy?bz? 其中s，r，b为变化区域内有一定限制的实参数，该方程形式简单，表面上看并无惊人之处，但由该方程揭示出的许多现象，促使“混沌”成为数学研究的崭新领域，在实际应用中产生了巨大的影响。 二、问题分析 Lorenz方程实际上是一个四元一阶常微分方程，用解析法精确求解是不可能的，只能用数值计算，最主要的有欧拉法、亚当法和龙格- 库塔法等。为了得到较高精度的数值解,本实验要求采用经典四阶龙格—库塔方法求解该问题。 三、实验程序及注释 算法程序 function [T]=Runge_Kutta ％定义算法，其中f为待解方程组， x0是初始自变量，y0是初始函数 值，h是步长，n为步数 if nargin n=100; ％如果输入参数个数小于5，则步数 n=100 end r=size;r=r; ％返回初始输出矩阵的行列数，并将 值赋给r s=size;s=s; ％返回初始输入矩阵的行列数，并 将值赋给s r=r+s; T=zeros; T=[y0;x0]; for t=2:n+1 ％以下是具体的求解过程 k1=feval); k2=feval+T;x0+h/2]); k3=feval+T;x0+h/2]); k4=feval;x0+h]); x0=x0+h; T=[T+*;x0]; end 主程序 function dy=fun ％定义函数 s=; ％给参数s，r，b赋值 r=; b=/3; dy=s*-x); ％Lorenz方程表达式 dy=-x*x-x); dy=x*x-b*x; dy=dy; 运行程序 T=Runge_Kutta; ％调用前面的算法程序 plot3,T,T); ％显示三分量的关系图 axis ％定义坐标轴长度 view ％设定观察角度 四、实验数据结果及分析 各初始变量相同时的图像分析 各初始变量取相同的值[10，10，10]，运行上述程序后，得到如下图像： 从图中可以看出，各初始变量相同时，曲线总是被吸引回奇怪吸引子附近作来回跳跃。初始变量值取为[-10，-10，-10] ，[20，20，20]时，依然如此。图像如下： [-10，-10，-10] [20，20， 20] 初始值的每个分量变化对图像的影响 y分量： [0，2，0] [0，5， 0] [0，15，0] [0，20， 0] 从上面可以看出，随着初始y值的增大，奇怪吸引子中曲线在其附近来回跳跃的两个位置中的一个吸引力变弱，另一个吸引力变强。初始y继续增大到某一特定值，情况又会变回来。这说明在空间存在一些区域，当初始位置位于这些区域外时解将出现奇怪吸引子的性质，而在这些区域以内解将呈现普通吸引子的性质。 z分量：[0，0， 20] 从上图可以看出解的曲线为一直线，这可以从方程的角度来解释。当x=0,y=0时在方程中dx/dt=0,dy/dt=0，x，y 方向的值不发生变化，仅z方向的值变化，因此解为一直线。 调整参数r、s、b对图像的影响 为便于分析，我们只调整r、s、b三个参数中的任意一个。当只调整b且将初始变量取为[0，eps，0]。具体情况如下： s=，r=，b=/3 s=，r=，b=/3 s=，r=，b=/3 s=，r=， b=/3 增大b值时，Lorenz曲线在其附近来回跳跃的两个位置会一个加强，一个减弱。当b达到某一值时，个位置丧失吸引力，另一位置则将曲线完全吸引过来变成普通吸引子。改变s和r的值也有类似的现象。 五、实验结论 本实验利用龙格—库塔法对Lorenz方程进行了分析，从实验中我们得出，Lorenz方程的解对初始变量和参数r、s、b具有很强的敏感性。 山西大学计算机与信息技术学院 湖北民族学院理学院《数值分析》课程实验报告 湖北民族学院理学院《数值分析》课程实验报告 xn?) 湘潭大学数学与计算科学学院 学生实验报