| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
 |
| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
|
2009-12-24, 21:55
|
#2 |
|
初级会员
注册日期: 2009-02-16
年龄: 42
帖子: 1
声望力: 0  |
回复: 一个很简短的微分方程组，进来看看吧！
用龙格-库塔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') |
|
|
树形模式