[MATLAB基础] 一个很简短的微分方程组，进来看看吧！ - MATLAB爱好者论坛-LabFans.com

wangwei058043 · 2009-12-23, 08:02

dx/dt=-x*e^(1-t)+0.8y(t-0.1)。 x(0)=0。
dy/dt=x-y^3。 y(0)=2。
上面给定的是微分方程组及其初始条件，要求作出两个图形来。

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')

wangwei058043 · 2009-12-25, 08:38

引用:

作者: wangwei058043

dx/dt=-x*e^(1-t)+0.8y(t-0.1)。 x(0)=0。
dy/dt=x-y^3。 y(0)=2。
上面给定的是微分方程组及其初始条件，要求作出两个图形来。

在这里我做点补充：0.8y(t-0.1)中的不是y*(t-0.1)而是y(t)向右移了0.1个单位即y(t-0.1)

laosam280 · 2009-12-25, 13:29

上面的结果是不对的，时滞微分方程怎么可以直接调用龙格库塔的程序？

2009-12-23, 08:02	#1
wangwei058043 初级会员注册日期: 2009-12-01 帖子: 5 声望力: 0	一个很简短的微分方程组，进来看看吧！ dx/dt=-x*e^(1-t)+0.8y(t-0.1)。 x(0)=0。 dy/dt=x-y^3。 y(0)=2。上面给定的是微分方程组及其初始条件，要求作出两个图形来。

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')

2009-12-25, 13:29	#4
laosam280 高级会员注册日期: 2008-11-07 住址: 湖南长沙帖子: 233 声望力: 21	回复: 一个很简短的微分方程组，进来看看吧！上面的结果是不对的，时滞微分方程怎么可以直接调用龙格库塔的程序？ __________________ 坚持就是胜利，努力就有奇迹。