龙格—库塔法分析Lorenz方程课程设计.docVIP

《数值分析》课程设计实验报告 龙格—库塔法分析Lorenz方程 问题叙述 考虑著名的Lorenz方程 其中s，r，b为变化区域内有一定限制的实参数，该方程形式简单，表面上看并无惊人之处，但由该方程揭示出的许多现象，促使“混沌”成为数学研究的崭新领域，在实际应用中也产生了巨大的影响。 问题分析 Lorenz方程实际上是一个四元一阶常微分方程，用解析法精确求解是不可能的，只能用数值计算，最主要的有欧拉法、亚当法和龙格- 库塔法等。为了得到较高精度的，我们采用经典四阶龙格—库塔方法求解该问题。 实验程序及注释 （1）算法程序 function [T]=Runge_Kutta(f,x0,y0,h,n) ％定义算法，其中f为待解方程组，x0是初始自变量，y0是初始函数值，h是步长，n为步数 if nargin5 n=100; ％如果输入参数个数小于5，则步数n=100 end r=size(y0);r=r(1); ％返回初始输出矩阵的行列数，并将值赋给r(1) s=size(x0);s=s(1); ％返回初始输入矩阵的行列数，并将值赋给s(1) r=r+s; T