Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#2 |
高级会员
注册日期: 2008-11-07
住址: 湖南长沙
帖子: 233
声望力: 21 ![]() |
![]()
代码我写好了:
%这就是一个例子,初值我都自己定义了,你要计算可以修改。把上面的X1'=AX0+(B+1)X1+X1^2*X2;X2'=X0+BX1-X1^2*X2这个方程组写成一个函数文件,存为myfunfun.m t_final=100; x0=[1;1]; [t,x]=ode45('myfunfun',[0,t_final],x0); plot(t,x) 可惜了,因为有非线性项,我电脑没有在短时间内算出来。明天我去实验室用大型机给你算一下,再找个时间把结果发给你。 |
![]() |
![]() |