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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解二阶微分方程
252/3返回列表上一页123下一页
查看: 5591  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

馨香诺言

新虫 (小有名气)
引用回帖:
10楼: Originally posted by cooooldog at 2016-03-28 09:30:41
具体是哪个文献？中文的还是外文的？
作者用的什么方法有没有提到？
这个方程直接求解的确不容易。...

是一篇中文文献，我上传了，红线标出了求解方程。他不是用边值法，他用的是四阶龙格库塔法。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 液桥自由表面变形的二维数值模拟_彭岚.pdf
    • 2016-03-28 10:32:45, 286.19 K
走好每一步！
11楼2016-03-28 10:32:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)
引用回帖:
9楼: Originally posted by 皓小天 at 2016-03-28 09:13:19
你可以拟合它根据微分方程画出的曲线找个合适的原函数，代入微分方程进行检验。...

这个方程求解不出来该怎么拟合？
赞一下 回复此楼
走好每一步！
12楼2016-03-28 10:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
11楼: Originally posted by 馨香诺言 at 2016-03-28 10:32:51
是一篇中文文献，我上传了，红线标出了求解方程。他不是用边值法，他用的是四阶龙格库塔法。...

用Mathematica求解发现，当前的边值条件下可以得到两个完全符合条件的数值解。

非线性的微分方程其解的唯一性本来就难以确定，这看上去也很正常。这个问题因为简单，所以，两个解容易发现，复杂的问题就难说。

首先要把问题变成初值形式，比如y'[0]=某个给定的值，然后求微分方程，打靶、迭代直至给定的值求出微分方程数值解y[x]中对于的y[1]=1

给定了y'[0]=v 的情况下，求y[1]的函数写成这样：
CODE:
obj[v_] :=
Module[{yy1, y, x},
  yy1[x_] =
   y[x] /. NDSolve[{y''[x] ==
        204/100 (1 - x) (1 + y'[x]^2)^(1/2) + (1 + y'[x]^2)/y[x],
       y[0] == 1, y'[0] == v}, y[x], {x, 0, 2}][[1]];
  yy1[1]
  ]

绘图可以发现，这个函数等于1 (也就是y[1]=1）有两种初值v都能满足条件
matlab求解二阶微分方程

这两个初值分别是：
-2.58238958868605
和
-3.59253658799395

这用牛顿法对上面的函数等于1的方程，从-2 , -4初值出发迭代就可以分别得到
CODE:
FindRoot[obj[x] == 1, {x, -2}, Evaluated -> False]
FindRoot[obj[x] == 1, {x, -4}, Evaluated -> False]

然后可以求两个不同的非线性常微分方程的初值问题得到两个都满足条件的数值解。
matlab求解二阶微分方程-1
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
13楼2016-03-28 14:00:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)
引用回帖:
13楼: Originally posted by Mr__Right at 2016-03-28 14:00:25
用Mathematica求解发现，当前的边值条件下可以得到两个完全符合条件的数值解。

非线性的微分方程其解的唯一性本来就难以确定，这看上去也很正常。这个问题因为简单，所以，两个解容易发现，复杂的问题就难说。
...

谢谢你的解答，你是用Mathematica这个软件求得整个过程吗？这个软件从没接触过。
赞一下 回复此楼
走好每一步！
14楼2016-03-28 15:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)
引用回帖:
14楼: Originally posted by 馨香诺言 at 2016-03-28 15:14:07
谢谢你的解答，你是用Mathematica这个软件求得整个过程吗？这个软件从没接触过。...

求解思路相似，用matlab也完全可以做到。

mathematica解决这个问题的好处是，

语法独特，代码简练。

matlab如果四阶荣格库塔，应该ode45
不过还要把微分方程改写成常微分方程组的形式，要单独写成一个函数，……
这些在mathematica中不是必须的；

而且对一个求解微分方程的过程作迭代的话，matlab中需要嵌套函数，
写法更加繁琐。没有mathematica的notebook演算方式自由。

听说R2016a的matlab中这方面功能有所增强，不妨试试去。
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
15楼2016-03-28 18:28:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)
引用回帖:
15楼: Originally posted by Mr__Right at 2016-03-28 18:28:53
求解思路相似，用matlab也完全可以做到。

mathematica解决这个问题的好处是，

语法独特，代码简练。

matlab如果四阶荣格库塔，应该ode45
不过还要把微分方程改写成常微分方程组的形式，要单独写成一个函 ...

你好，如果用mathematica的话，就是用你写得上面那两段小程序就可以了吗？
赞一下 回复此楼
走好每一步！
16楼2016-03-28 18:48:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

郭琪琪1028

木虫 (正式写手)
这个方程的解，貌似求不出来的……

发自小木虫Android客户端
赞一下 回复此楼
17楼2016-03-28 19:06:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)
引用回帖:
17楼: Originally posted by 郭琪琪1028 at 2016-03-28 19:06:57
这个方程的解，貌似求不出来的……

好像是。。。。。我再试试
回复此楼
走好每一步！
18楼2016-03-29 08:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓小天

木虫之王 (职业作家)
引用回帖:
13楼: Originally posted by Mr__Right at 2016-03-28 14:00:25
用Mathematica求解发现，当前的边值条件下可以得到两个完全符合条件的数值解。

非线性的微分方程其解的唯一性本来就难以确定，这看上去也很正常。这个问题因为简单，所以，两个解容易发现，复杂的问题就难说。
...

楼主提的问题的指数不是0.5而是1.5，你换成1.5试试？mathematica也报错吧
赞一下 回复此楼
请不要站内找我要书，如果需要请到书籍板块求助
19楼2016-05-05 09:03:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)
引用回帖:
19楼: Originally posted by 皓小天 at 2016-05-05 09:03:09
楼主提的问题的指数不是0.5而是1.5，你换成1.5试试？mathematica也报错吧...



你说试试我就去试试，那也太没面子了吧。

你是自己不会用软件、或对求解ODE/PDE问题不熟悉而已，不关我事。

专门的基础知识的私人辅导都是要收费的。
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
20楼2016-05-05 09:31:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 馨香诺言 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料299专硕求调剂 +17 +21 2026-04-09 17/850 2026-04-15 17:20 by zs92450
[考研] 085500求调剂材料 +11 易11122 2026-04-09 11/550 2026-04-14 17:59 by lhj2009
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 335求调剂 +19 想上岸呀！！ 2026-04-12 21/1050 2026-04-14 16:23 by Art1977
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 人工智能320调剂08工类还有机会吗 +18 振—TZ 2026-04-10 19/950 2026-04-14 10:34 by screening
[考研] 求调剂 +3 我爱高数高数爱� 2026-04-12 3/150 2026-04-14 01:00 by 王珺璞
[考研] 332求调剂 +15 蕉蕉123 2026-04-10 15/750 2026-04-13 23:12 by pies112
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 346，工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 一志愿北理工298英一数二已上岸，感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
信息提示
关闭
请填处理意见
关闭