理力-龙格-张新意.docVIP

龙格-库塔法解洛伦兹方程 结构01 张新意 一．洛伦兹方程 美国科学家Lorenz在研究大气对流模型时得到了以下方程组 X=σY-XY=ρX-Y-XZZ=-βZ+XY 称为Lorenz方程。 现对采用常用具体参数σ=10，β=8/3，ρ=28的Lorenz方程进行求解，描绘出X、Y、Z随时间变化的关系图形和三维图形，并验证该方程的初值敏感性质。 二．求解过程 （1）建立函数 用三维向量dx代表X、Y、Z，程序代码为： function dx=Lorenz(t,x) dx=zeros(3,1); dx(1)=10*(-x(1)+x(2)); dx(2)=28*x(1)-x(2)-x(1)*x(3); dx(3)=x(1)*x(2)-8*x(3)/3; 保存为Lorenz.m （2）用四五阶龙格-库塔法建立X、Y、Z关于时间的函数关系并绘制图形 程序代码： clc clear [t,x]=ode45(lorenz,[0 20],[10,10,10]); figure(1) hold on title(X(t)图像) plot(t,x(:,1)) xlabel(t) ylabel(X) figure(2) hold on title(Y(t)图像) plot(t,x(:,2)) xlabel(t) ylabel(Y) figure(3) hold on tit