| 查看: 431 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 wangyanjun7392 的 5 个金币 | |||
wangyanjun7392银虫 (小有名气)
|
[交流]
【求助】如何用MATLAB求解这个微分方程
|
||
|
如何用MATLAB解微分方程0.0469y''+2y'+2.45cosy=1 非常感谢高手指点 |
回复此楼» 猜你喜欢
上海大学实验技术岗位非升即走
已经有8人回复
-
今年审到国自然15份，谈谈感受
已经有14人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
-
半夜喝咖啡
已经有5人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
★ ★ ★ ★ ★ ★ ★
nono2009(金币+2,VIP+0):谢谢应助！ 12-7 10:58
nono2009(金币+0,VIP+0):欢迎去投票，寻找更多同行，同时领取金币。http://emuch.net/bbs/viewthread.php?tid=1649328 12-7 10:58
wangyanjun7392(金币+5,VIP+0): 12-7 16:10
nono2009(金币+2,VIP+0):谢谢应助！ 12-7 10:58
nono2009(金币+0,VIP+0):欢迎去投票，寻找更多同行，同时领取金币。http://emuch.net/bbs/viewthread.php?tid=1649328 12-7 10:58
wangyanjun7392(金币+5,VIP+0): 12-7 16:10
|
这个简单， 编写求导数的函数并以wyj存盘： function dy = wyj(x, y) dy = [y(2); (1-2.45*cos(y(1))-2*y(2))/0.0469]; 然后在Matlab工作窗口输入： >>tspan=[0 10]; >>y0=[0 1]; >>options=odeset('RelTol', 1e-4, 'AbsTol', [1e-4 1e-4]); >>[x, y]=ode45(@wyj, tspan, y0, options); 就差不多了。 |
2楼2009-12-07 10:53:06
wangyanjun7392
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 298.6
- 散金: 2
- 帖子: 182
- 在线: 10.8小时
- 虫号: 833617
- 注册: 2009-08-24
- 性别: MM
- 专业: 控制理论与方法
3楼2009-12-07 16:11:28
30