MATLAB爱好者论坛-LabFans.com - 查看单个帖子

Marthay · 2009-12-24, 21:55

用龙格-库塔4 阶算法，ode45即可
function defg=group(t,x)
defg=[-x(1)*exp(1-t)+0.8*x(2)*(t-0.1);x(1)-x(2)^3];

clc
clear
x0=[0 2]
[t,x]=ode45('group',[0,40],x0);
x1=x(:,1);
x2=x(:,2);
plot(t,x1,'r')
hold on
plot(t,x2,'g')

2009-12-24, 21:55	#2
Marthay 初级会员注册日期: 2009-02-16 年龄: 42 帖子: 1 声望力: 0	回复: 一个很简短的微分方程组，进来看看吧！用龙格-库塔4 阶算法，ode45即可 function defg=group(t,x) defg=[-x(1)exp(1-t)+0.8x(2)*(t-0.1);x(1)-x(2)^3]; clc clear x0=[0 2] [t,x]=ode45('group',[0,40],x0); x1=x(:,1); x2=x(:,2); plot(t,x1,'r') hold on plot(t,x2,'g')