很多小伙伴问MATLAB 解偏微分方程组的数值解法编程问题，这里提供2种常用的MATLAB数值解法及源代码（ODE45和ODE15s）。

洛伦兹混沌

clc;clear;close all

tspan=[0,100];

cond=[1.3 4 2];

opts=odeset('RelTol',1e-5,'AbsTol',1e-5);

[t,x]=ode45(@chaos,tspan,cond,opts)

plot(x(:,1),x(:,2))

plot(x(:,1),x(:,3))

plot(x(:,2),x(:,3))

plot3(x(:,1),x(:,2),x(:,3))

opts=odeset('RelTol',1e-5,'AbsTol',1e-5);

[t,x]=ode15s(@chaos,tspan,cond,opts)

plot(x(:,1),x(:,2))

plot(x(:,1),x(:,3))

plot(x(:,2),x(:,3))

plot3(x(:,1),x(:,2),x(:,3))

function dx=chaos(t,x)

ru=10;

b=8/3;

r=28;

dx=[ru*(x(2)-x(1));r*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-b*x(3)];

end

得到的图形为：