用matlab求解lorenz方程，求问这段程序用的是什么方法

用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-171e-171e-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

解决边界问题，可以用山纤凳matlab内置的ode函数。

解竖轮决边逗旅界-初值问题（偏微分方程），可以用MATLAB内置的pdepe函数。

觉得有帮助就采纳吧。

---