| 查看: 2542 | 回复: 8 | ||||
| 本帖产生 1 个 仿真EPI ，点击这里进行查看 | ||||
bdwhq新虫 (著名写手)
|
[求助]
请各位前辈指点一下，如何用MATLAB求解如下的微分方程，谢谢！
|
|||
|
请各位前辈指点一下，如何用MATLAB求解如下的微分方程，谢谢！ 我想用matlab求解下面的微分方程组： x'_1=x_2+x_1*sin(x_1); x'_2=w+x_1*x_2; w'=|v'|*(0.5*v-w)+0.3*v'; 其中V是可以给出的函数：v=-5*(x_2-u)-x_1*x_2,v函数中的u=-10*x_1-x_1*sin(x_1). 另外，x'_1,x'_2,w',v'表示各个变量对时间t的导数，X_1(0),x_2(0),w(0),v(0),也就是各个变量在t=0时的初值是给定的，比如X_1(0)=0.1,x_2(0)=0.1,w(0)=0,v(0)=0,也可以是其他的值。想画出t从0到20秒x_1，x_2的图像，敬请大家给予指点，谢谢！ 有的前辈可能说把v及u的导数具体算出来，然后利用ode45来求解，可是这样来算是特别麻烦，并且结果还不理想。烦请各位看看能否好的办法，小弟在此谢谢各位了！！！ 若悬赏金币少，可追加。 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】用pdepe求解偏微分方程
已经有10人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】Python 中 求解微分方程 odeint 的用法
已经有11人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
★ ★
感谢参与，应助指数 +1
bdwhq: 金币+1, ★有帮助 2012-03-28 19:01:54
臭水沟: 金币+1, 谢谢交流~~ 2012-03-31 12:30:33
感谢参与，应助指数 +1
bdwhq: 金币+1, ★有帮助 2012-03-28 19:01:54
臭水沟: 金币+1, 谢谢交流~~ 2012-03-31 12:30:33
|
本帖内容被屏蔽 |
2楼2012-03-28 12:24:19
★ ★ ★
bdwhq: 金币+2, ★有帮助 2012-03-28 19:01:40
臭水沟: 金币+1, 谢谢交流~~ 2012-03-31 12:30:41
bdwhq: 金币+2, ★有帮助 2012-03-28 19:01:40
臭水沟: 金币+1, 谢谢交流~~ 2012-03-31 12:30:41
|
本帖内容被屏蔽 |
3楼2012-03-28 12:25:14
bdwhq
新虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 1175.1
- 散金: 5117
- 红花: 1
- 帖子: 1123
- 在线: 639.1小时
- 虫号: 415789
- 注册: 2007-06-29
- 性别: GG
- 专业: 控制论中的数学方法
4楼2012-03-28 18:58:35
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
caemechanics: 仿真EPI+1, 感谢应助 2012-03-28 21:48:58
bdwhq: 金币+7, ★★★★★最佳答案, 非常感谢dbb627 的热心帮助，谢谢谢谢！！！！！！！ 2012-03-29 08:32:26
感谢参与，应助指数 +1
caemechanics: 仿真EPI+1, 感谢应助 2012-03-28 21:48:58
bdwhq: 金币+7, ★★★★★最佳答案, 非常感谢dbb627 的热心帮助，谢谢谢谢！！！！！！！ 2012-03-29 08:32:26
|
>> x1=sym('x1(t)');x2=sym('x2(t)'); >> vv=-5*(x2+10*x1+x1*sin(x1))-x1*x2 vv = - 50*x1(t) - 5*x2(t) - x1(t)*x2(t) - 5*sin(x1(t))*x1(t) >> syms t >> dv=diff(vv,t) dv = - x1(t)*diff(x2(t), t) - x2(t)*diff(x1(t), t) - 5*sin(x1(t))*diff(x1(t), t) - 50*diff(x1(t), t) - 5*diff(x2(t), t) - 5*cos(x1(t))*x1(t)*diff(x1(t), t) ================================ v‘=-x1*x2'-x2*x1'-5*sin(x1)*x1'-50*x1'-5*x2'-5*cos(x1)*x1*x1' 微分方程组变为 x1=x(1) x2=x(2) v=x(3) w=x(4) =========== x1’=x2+x1*sin(x1); x2‘=w+x1*x2; v‘=-x1*x2'-x2*x1'-5*sin(x1)*x1'-50*x1'-5*x2'-5*cos(x1)*x1*x1' w'=|v'|*(0.5*v-w)+0.3*v'; matlab 代码 |
5楼2012-03-28 21:15:34
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
6楼2012-03-28 21:16:11
|
本帖内容被屏蔽 |
7楼2012-03-30 14:02:08
 |
8楼2012-03-31 09:30:57
|
本帖内容被屏蔽 |
9楼2012-04-04 12:14:57
25