发布网友 发布时间:2023-07-13 12:26
共1个回答
热心网友 时间:2024-01-18 21:29
用matlab求解lorenz方程这段程序的是数值微分方程(组)的变步长四五阶龙格——库塔方法,即ode45()函数。 Lorenz.m function s = Lorenz(t,x) s = [-8/3*x(1)+x(2)*x(3); %Lorenz方程组 -10*x(2)+10*x(3); -x(1)*x(2)+28*x(2)-x(3)];主函数: x0 = [1e-17;1e-17;1e-17]; [t,x] = ode45('Lorenz',[0,150],x0); axis equal plot3(x(1:(end/2),1),x(1:(end/2),2),x(1:(end/2),3),'r') hold on plot3(x((end/2):end,1),x((end/2):end,2),x((end/2):end,3),'b') hold off追问这个程序运行出来是0、1序列么