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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解二阶微分方程
251/3返回列表123下一页
查看: 5585  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854求调剂 +12 门路摸摸 2026-04-15 13/650 2026-04-16 00:14 by Art1977
[考研] 273求调剂 +6 白居不易. 2026-04-09 8/400 2026-04-15 22:02 by wooluyong
[考研] 291求调剂 +7 关忆北. 2026-04-14 7/350 2026-04-15 14:39 by 黑科技花岗岩
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 求调剂288 +7 ioodiiij 2026-04-10 9/450 2026-04-13 08:33 by Hayaay
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[找工作] 山东高校教师考核超级无底线，员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 277 数一104，学硕，求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
信息提示
关闭
请填处理意见
关闭