24小时热门版块排行榜    

Znn3bq.jpeg
查看: 5591  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

馨香诺言

新虫 (小有名气)

[求助] matlab求解二阶微分方程 已有3人参与

我想求解这样一个二阶微分方程(下图1),我使用了matlab边值函数,按着书上编写的程序如下(下图2),可为什么求不出来。提示错误:错误使用 bvp4c (line 251)
无法求解排列方程式 - 遇到了奇异性雅可比行列式。
有做过的大神帮忙给看看哪里出错了,不胜感激。

matlab求解二阶微分方程
1.jpg


matlab求解二阶微分方程-1
2.jpg
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 程序.m
  • 2016-03-26 15:54:55, 561 bytes

» 猜你喜欢

走好每一步!
已阅   回复此楼   关注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

这两个初值分别是:
-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求解二阶微分方程-3
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
13楼2016-03-28 14:00:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

cooooldog

铁杆木虫 (著名写手)

ส็็็

问题可能出在方程的边值条件

方程自己无解。

应该看原始问题的具体背景,方程是如何得到和建立的?尤其两个边值条件?
ส็็็็็็็็็็็็็็็็็็็็
2楼2016-03-26 21:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

感谢参与,应助指数 +1
原常微分方程的数值解不唯一,这个问题可以改成初值问题来做。
即,用迭代的方法

y[0]=1;
y'[0]=u ; (设置不同的u,直到得到数值解满足 y[1]=1为止)

得到的满足条件的数值解不唯一。

不过为何边值方法不能求解?这个的确很奇怪。
ส็็็็็็็็็็็็็็็็็็็็
3楼2016-03-27 09:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)

引用回帖:
3楼: Originally posted by cooooldog at 2016-03-27 09:43:21
原常微分方程的数值解不唯一,这个问题可以改成初值问题来做。
即,用迭代的方法

y=1;
y'=u ; (设置不同的u,直到得到数值解满足 y=1为止)

得到的满足条件的数值解不唯一。

不过为何边值方法不能求解?这 ...

谢谢你的回答,我使用的这个公式是一篇文献上推导出的无量纲方程,我使用边值法一直出错。
你说的初值问题,我一开始想用四阶龙格库塔法求,但当时认为y'(0)不知道,就放弃了。你能详细说一下你指的这个迭代法怎么求吗?
走好每一步!
4楼2016-03-27 11:14:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 300求调剂 +9 橙a777 2026-04-15 9/450 2026-04-16 06:47 by fenglj492
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 复试调剂 +21 积极向上; 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
信息提示
请填处理意见