matlab怎样求微分方程组特解(解偏微分方程组的数值解法)
很多小伙伴问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
得到的图形为:
,
免责声明:本文仅代表文章作者的个人观点,与本站无关。其原创性、真实性以及文中陈述文字和内容未经本站证实,对本文以及其中全部或者部分内容文字的真实性、完整性和原创性本站不作任何保证或承诺,请读者仅作参考,并自行核实相关内容。文章投诉邮箱:anhduc.ph@yahoo.com