福建农林大学常微分课程论文.docVIP

福建农林大学计算机与信息学院 （数学类课程） 课程实习报告 课程名称： 常微分方程课程实习 实习题目： 常微分方程数值求解问题的实习 姓 名： XXX 系： 应用数学 专 业： 数学与应用数学 年 级： 2010 学 号： 102260002XXX 指导教师： 陈永雪 职 称： 讲师 2011 年 12 月 1 日  福建农林大学计算机与信息学院数学类 课程实习报告结果评定 评语： 成绩： 指导教师签字： 评定日期： 目 录 1. 实习的目的和任务 1 2. 实习要求 1 3. 实习地点 1 4. 主要仪器设备 1 5. 实习内容 1 5.1数值求解Lorenz方程，模拟混沌现象……………………………………1 5.1.1数值求解方程 2 5.1.2数值解对初值的敏感性 3 5.1.3模拟混沌现象 5 5.2用不同格式对同一个初值问题的数值求解及其分析 7 5.2.1数值求解 7 5.2.2 用欧拉法求解………………………………………………………...9 5.2.3 用改进欧拉法求解………………………………………………….10 5.2.4用4阶龙格—库塔求解…………………………………………….11 5.2.5问题讨论与分析…………………………………………………….13 5.3一个算法不同步长求初值问题及其分析………………………………..15 5.3.1取步长……………………………………………………………….15 5.3.2分析和结论………………………………………………………….19 6. 结束语 19 参考文献 19 常微分方程课程实习 1.实习的目的和任务 目的：通过课程实习能够应用MATLAB软来计算微分方程（组）的数值解；了解常微分方程数值解。 任务：通过具体的问题，利用MATLAB软件来计算问题的结果，分析问题的结论。 2.实习要求 能够从案例的自然语言描述中，抽象出其中的数学模型；能够熟练应用所学的数值解计算方法；能够熟练使用MATLAB软件；对常微分方程数值解有所认识，包括对不同算法有所认识和对步长有所认识。 3.实习地点 学生宿舍 4.主要仪器设备 计算机、 Microsoft Windows XP Matlab 7.6 5.实习内容 5.1 数值求解Lorenz方程，模拟混沌现象 Lorenz方程：，其中a=10,b=8/3,c=28. 5.1.1 数值求解方程 程序代码： 建立函数文件 function xdot = lorenzeq(t,x) a=10;r=28;b=8/3; xdot=[a*(x(2)-x(1)); r*x(1)-x(2)-x(1)*x(2); x(1)*x(2)-b*x(3)]; 调用函数 global a b c a=10;b=8/3;c=28;%定义参数 lorenz=@(t,x)[a*(x(2)-x(1));c*x(1)-x(1)*x(3)-x( 2);x(1)*x(2)-b*x(3)]; % 定义函数 [T,X]=ode45(lorenz,[0,20],[0;1;2]);%数值法解微分方程 Data=[T,X] plot3(X(:,1),X(:,2),X(:,3),m)%绘图 view(-20,60);%设置视角 xlabel(x);ylabel(y);zlabel(z);%标记坐标轴 Data = 0 0 1.0000 2.0000 0.0001 0.0012 0.9999 1.9994 0.0002 0.0025 0.9998 1.9987 0.0004 0.0037 0.9996 1.9980 0.0005 0.0050 0.9995 1.9974 ................................... ................................... 19.9735 -10.0007 0.2594 38.1095 19.9801 -9.3215 0.8813 37.4048 19.9867 -8.6489 1.4130 36.6808 19.9934 -7.9881 1.8611 35.9479 20.0000 -7.3435 2.2331 35.2142 运行则得到参数a=10，b=8/3，c=28自变