| 查看: 275 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】有人能帮忙求解这个微分方程吗组？
|
|||
|
dx/dt = -3*sin(y)-0.5*sqrt(1+0.01*x)*(cos(y))^2; dy/dt = -3*cos(y)+0.5*sqrt(1+0.01*x)*(sin(y))^2; 其中：t = 0:10:800 x(800)=191 y(8000)=0.029 龙哥库塔法处理的是t=0的初值，这里却是终点处的值，不知该如何处理才好…… |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
3楼2009-09-11 18:42:29
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 11859.7
- 散金: 200
- 帖子: 642
- 在线: 678.9小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
★ ★ ★
chenyouzhi(金币+2,VIP+0):先送两个，能在解释一下原理吗，很希望向你学习！不管解释与否，明天奉送另外两个…… 9-11 18:45
kuhailangyu(金币+1,VIP+0):欢迎参与 9-11 20:09
chenyouzhi(金币+2,VIP+0):先送两个，能在解释一下原理吗，很希望向你学习！不管解释与否，明天奉送另外两个…… 9-11 18:45
kuhailangyu(金币+1,VIP+0):欢迎参与 9-11 20:09
|
我有个问题，你求解的步长是不是太大了？ 相关程序如下： clc;close all;clear; df=inline(['[-3*sin(x(2))-0.5*sqrt(1+0.01*x(1))*cos(x(2))^2;',... '-3*cos(x(2))+0.5*sqrt(1+0.01*x(1))*sin(x(2))^2]'],'t','x'); tspan=800:-10:0; x0=[191,0.029]; [t,x]=ode45(df,tspan,x0); whos plot(tspan,x(:,1),'k');hold on plot(tspan,x(:,2),'r'); legend('x','y') 下面是输出图形： |
2楼2009-09-11 18:28:39
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 11859.7
- 散金: 200
- 帖子: 642
- 在线: 678.9小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
4楼2009-09-11 18:51:07
5楼2009-09-11 18:55:12