应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » Matlab数值求解二阶常微分方程
101/1返回列表
查看: 3319  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

815292578

木虫 (著名写手)

[求助] Matlab数值求解二阶常微分方程

现在遇到一个二阶微分方程，需要用Matlab 求解
方程如下：
d/dy(dy/dr)=(2y/a-1/r)dy/dr+y/r^2
其中a为常数。
首先这个方程是什么类型的方程啊？查阅过一些书籍，都没有与之匹配的类型啊，方程含有因变量y，我想差分后迭代求解。
请各位大侠指教啊.......
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-08 11:15:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher
d/dy(dy/dr)这里打错没有
赞一下 回复此楼
2楼2013-09-08 12:53:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

815292578

木虫 (著名写手)
引用回帖:
2楼: Originally posted by 锐利的碎片 at 2013-09-08 12:53:20
d/dy(dy/dr)这里打错没有

哈，太粗心了啊。不是d/dy(dy/dr)而是d/dr(dy/dr),是y对r的二阶导数。
方程应为：
d/dr(dy/dr)=(2y/a-1/r)dy/dr+y/r^2
还有边界条件：
当r=0时，y=0
当r=r0时，y=b
其中，a、r0、b都是常数。
赞一下 回复此楼
3楼2013-09-08 14:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher
引用回帖:
3楼: Originally posted by 815292578 at 2013-09-08 14:41:15
哈，太粗心了啊。不是d/dy(dy/dr)而是d/dr(dy/dr),是y对r的二阶导数。
方程应为：
d/dr(dy/dr)=(2y/a-1/r)dy/dr+y/r^2
还有边界条件：
当r=0时，y=0
当r=r0时，y=b
其中，a、r0、b都是常数。...

当r=0时，y=0
r不能为0啊
赞一下 回复此楼
4楼2013-09-08 15:21:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

815292578

木虫 (著名写手)
引用回帖:
4楼: Originally posted by 锐利的碎片 at 2013-09-08 15:21:34
当r=0时，y=0
r不能为0啊...

这个方程表示一段圆弧曲线。其中y代表夹角西塔（曲线切线与竖直方向的夹角）；r代表半径
边界条件表示：r=0时，夹角西塔（即y）=0
                    r=r0时，夹角西塔（即y）=b（临界夹角）
方程曲线图如下所示
Matlab数值求解二阶常微分方程
360截图20130908202833640.jpg
赞一下 回复此楼
5楼2013-09-08 20:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

710506937

铜虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
ben_ladeng: 金币+1, 欢迎常来 2013-09-09 22:08:58
ts=[0,10];%求解区间
y=[0,b];%初始条件
ode45('D2y=2*y/a-1/r.*Dy+y/r^2',ts,y）
赞一下(1人) 回复此楼
天道酬勤！永不言弃！
6楼2013-09-08 20:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

815292578

木虫 (著名写手)
引用回帖:
6楼: Originally posted by 710506937 at 2013-09-08 20:38:26
ts=;%求解区间
y=;%初始条件
ode45('D2y=2*y/a-1/r.*Dy+y/r^2',ts,y）

十分感谢....
能不能把具体的Matlab程序写出来？我感觉没有那么简单，我曾套用很多例子去解，都没有得出结果。
赞一下 回复此楼
7楼2013-09-09 09:21:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

710506937

铜虫 (小有名气)

【答案】应助回帖

★
ben_ladeng: 金币+1, 欢迎常来 2013-09-09 22:09:12
这样看行不行！
syms  a,b，r;
ts=[0,10];%求解区间
y=[0,b];%初始条件
ode45('D2y=（2*y/a - 1/r）.*Dy+y/r^2',ts,y）
赞一下 回复此楼
天道酬勤！永不言弃！
8楼2013-09-09 09:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

★
感谢参与，应助指数 +1
ben_ladeng: 金币+1, 欢迎常来 2013-09-09 22:09:22
你给的边界条件是奇点，看起来比较麻烦。matlab的bvp4c处理不了，mathematica也解不出来。另外ode45是解初值问题的，对于二阶方程需要一点的值和导数作初始值，和你要的边值问题不一样。
赞一下(1人) 回复此楼
9楼2013-09-09 10:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

815292578

木虫 (著名写手)
引用回帖:
9楼: Originally posted by 锐利的碎片 at 2013-09-09 10:50:07
你给的边界条件是奇点，看起来比较麻烦。matlab的bvp4c处理不了，mathematica也解不出来。另外ode45是解初值问题的，对于二阶方程需要一点的值和导数作初始值，和你要的边值问题不一样。

是的，感觉是比较麻烦，很长时间一直没有解决。
现在想先做一下差分，看能不能化成三对角矩阵形式，这样就可以用迭代、追赶法求解。
最后再次感谢你的建议....
赞一下 回复此楼
10楼2013-09-09 19:06:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 815292578 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天基金会出结果吗？20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] 只有每年这种时候来逛逛小木虫 +4 yaoyewhu2008 2026-08-20 4/200 2026-08-20 08:01 by 201120518022
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭