【答案】应助回帖

» 本帖附件资源列表

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

2016-03-26 15:54:55, 561 bytes

» 猜你喜欢

三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复

走好每一步！

1楼2016-03-26 15:50:02

馨香诺言

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 141.5
散金: 50
红花: 2
帖子: 190
在线: 130.8小时
虫号: 3742322
注册: 2015-03-16
专业: 电化学分析

引用回帖:

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

查看全部 25 个回答

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

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

方程自己无解。

应该看原始问题的具体背景，方程是如何得到和建立的？尤其两个边值条件？

ส็็็็็็็็็็็็็็็็็็็็

2楼2016-03-26 21:42:06

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

【答案】应助回帖

感谢参与，应助指数 +1

原常微分方程的数值解不唯一，这个问题可以改成初值问题来做。
即，用迭代的方法

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

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

不过为何边值方法不能求解？这个的确很奇怪。

ส็็็็็็็็็็็็็็็็็็็็

3楼2016-03-27 09:43:21

馨香诺言

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 141.5
散金: 50
红花: 2
帖子: 190
在线: 130.8小时
虫号: 3742322
注册: 2015-03-16
专业: 电化学分析

引用回帖:

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