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

24小时热门版块排行榜    

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

馨香诺言

新虫 (小有名气)

[求助] 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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
走好每一步！
1楼 2016-03-26 15:50:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

皓小天

木虫之王 (职业作家)
引用回帖:
20楼: Originally posted by Mr__Right at 2016-05-05 09:31:57


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

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

专门的基础知识的私人辅导都是要收费的。...

这么简单的问题，你有什么好洋气的，你自己看下面图结果，我就用你给出的方法。指数是1.5时，曲线根本没有0点，所以此方程无解的；而你用指数为0.5算才会有两个解，但指数为0.5直接用matlab的bvp方法就可以得到，初值选1/2*sin(x)和1/2*cos(x)。
matlab求解二阶微分方程-2
untitled.jpg


matlab求解二阶微分方程-3
untitled1.jpg


matlab求解二阶微分方程-4
result.jpg
赞一下(2人) 回复此楼
请不要站内找我要书，如果需要请到书籍板块求助
21楼2016-05-05 10:54:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
4楼: Originally posted by 馨香诺言 at 2016-03-27 11:14:56
谢谢你的回答，我使用的这个公式是一篇文献上推导出的无量纲方程，我使用边值法一直出错。
你说的初值问题，我一开始想用四阶龙格库塔法求，但当时认为y'(0)不知道，就放弃了。你能详细说一下你指的这个迭代法怎么 ...

打靶法。
设定y'（0）的初始值，然后计算由此得到的y(1)的值和1的差，调整y'(0)初始值，直到得到的y(1)充分接近1
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
5楼2016-03-27 17:17:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓小天

木虫之王 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
方程无解，系数矩阵非常不好
赞一下 回复此楼
请不要站内找我要书，如果需要请到书籍板块求助
6楼2016-03-27 21:18:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)
引用回帖:
5楼: Originally posted by cooooldog at 2016-03-27 17:17:54
打靶法。
设定y'（0）的初始值，然后计算由此得到的y(1)的值和1的差，调整y'(0)初始值，直到得到的y(1)充分接近1...

谢谢你的回答，我按照你的方法再试一下。
赞一下 回复此楼
走好每一步！
7楼2016-03-28 09:09:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

馨香诺言

新虫 (小有名气)
引用回帖:
6楼: Originally posted by 皓小天 at 2016-03-27 21:18:04
方程无解，系数矩阵非常不好

这是我看的一篇文献上的，他说他按照这个方程解出来，不知道真假。
赞一下 回复此楼
走好每一步！
8楼2016-03-28 09:10:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓小天

木虫之王 (职业作家)
引用回帖:
8楼: Originally posted by 馨香诺言 at 2016-03-28 09:10:45
这是我看的一篇文献上的，他说他按照这个方程解出来，不知道真假。...

你可以拟合它根据微分方程画出的曲线找个合适的原函数，代入微分方程进行检验。
赞一下 回复此楼
请不要站内找我要书，如果需要请到书籍板块求助
9楼2016-03-28 09:13:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
8楼: Originally posted by 馨香诺言 at 2016-03-28 09:10:45
这是我看的一篇文献上的，他说他按照这个方程解出来，不知道真假。...

具体是哪个文献？中文的还是外文的？
作者用的什么方法有没有提到？
这个方程直接求解的确不容易。
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
10楼2016-03-28 09:30:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 馨香诺言 的主题更新
251/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 ，六级已过，有科研经历 +14 曦熙兮 2026-03-15 14/700 2026-03-21 19:12 by ColorlessPI
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
关闭
请填处理意见
关闭