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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解二阶微分方程
252/3返回列表上一页123下一页
查看: 5546  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 生物学071000 329分求调剂 +5 我爱生物生物爱� 2026-03-17 5/250 2026-03-22 16:42 by tcx007
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 生物学一志愿985，分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
信息提示
关闭
请填处理意见
关闭