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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
81/1返回列表
查看: 1164  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者azraelgaolei将赠送您 20 个金币

azraelgaolei

铁虫 (小有名气)

[求助] 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。

小弟需要解一个4阶非线性常微分方程，用mathematica一直没法得到需要结果，而且一直说步长小于精度，出现警告，搞了好久好久好久都没搞出来。据说1stopt简单，上手快，但是网上下了个1.5，按教程复制的程序，点了F9也没反应，也许是语法问题吧。我的问题如下，请哪位大哥空闲的话帮我跑一下，因为是微观的问题，所以尺度非常小，纳米级，看能能不能算，感激涕零。
方程是   ODEFuntion y''''=-1.20237*10^19*(y')^2*y''-6.2819*10^21
变量范围   x=[0,1.74*10^-7]
边界条件有4个  y=[0]=0, y[1.74*10^-7]=0, y'[0]=0, y'[1.74*10^-7]=0
小弟想知道x范围中点处的y值。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
帅得惊动裆中央。
1楼2012-07-12 21:41:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
dbb627: 金币+1, 感谢应助 2012-07-21 22:34:58
数据相差量级也太大了，下面不知对否，供参考：
CODE:
ODEStep = 1.74*10^(-7)/10;
Variable x,y,y',
ODEFunction y''''=-1.20237*10^19*(y')^2*y''-6.2819*10^21;
Data;
0,0,0
1.74*10^(-7),0,0

参数                  最佳估算
--------------------        -------------
y''初值        11995971.7330153
y'''初值        67927.458477625

No        x        目标y        计算y        目标y'        计算y'
1        1E-8        NAN        5.47047804059444E-10        NAN        0.0913097707467425
2        2E-8        NAN        1.09261736443608E-9        NAN        -0.0014427838744353
3        3E-8        NAN        4.92582333734155E-10        NAN        -0.0981845660250561
4        4E-8        NAN        -6.11168342806935E-11        NAN        0.00300253947914656
5        5E-8        NAN        4.71438462863779E-10        NAN        0.083380296077392
6        6E-8        NAN        9.4995567942877E-10        NAN        -0.00601770541774702
7        7E-8        NAN        2.83556771692151E-10        NAN        -0.103649819779204
8        8E-8        NAN        -2.48263437258767E-10        NAN        0.00812475646552469
9        9E-8        NAN        2.73915892295962E-10        NAN        0.0738683414753312
10        1E-7        NAN        6.47410579732159E-10        NAN        -0.0163892027996882
11        1.1E-7        NAN        -1.2455710142929E-10        NAN        -0.104911525261711
12        1.2E-7        NAN        -5.69277358685573E-10        NAN        0.0176540829734727
13        1.3E-7        NAN        -5.20803008701728E-11        NAN        0.0607494957855953
14        1.4E-7        NAN        1.37257022458506E-10        NAN        -0.0394864648151568
15        1.5E-7        NAN        -7.84216562091161E-10        NAN        -0.0914689201111017
16        1.6E-7        NAN        -1.02279018874811E-9        NAN        0.0324114075575231
17        1.7E-7        NAN        -5.37939651149806E-10        NAN        0.0381175630542217
18        1.74E-7        0        -4.55651976657919E-10        0        3.98986399474666E-17
赞一下 回复此楼
2楼2012-07-13 10:12:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

azraelgaolei

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2012-07-13 10:12:09
数据相差量级也太大了，下面不知对否，供参考：

ODEStep = 1.74*10^(-7)/10;
Variable x,y,y',
ODEFunction y''''=-1.20237*10^19*(y')^2*y''-6.2819*10^21;
Data;
0,0,0
1.74*10^(-7),0,0


参数       ...

非常感谢你给我的帮助！这个结果和我用mathematica做的类似，都是震荡的，肯定是有问题的，因为结果中的y在我文章中是位移，而且全部是正的。这个尺度太小确实有问题，我也考虑过尺度放大再做，但是怕数量级搞错，我若确定能放大麻烦再帮帮我！
赞一下 回复此楼
帅得惊动裆中央。
3楼2012-07-13 10:52:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

azraelgaolei

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2012-07-13 10:12:09
数据相差量级也太大了，下面不知对否，供参考：

ODEStep = 1.74*10^(-7)/10;
Variable x,y,y',
ODEFunction y''''=-1.20237*10^19*(y')^2*y''-6.2819*10^21;
Data;
0,0,0
1.74*10^(-7),0,0


参数       ...

而且您看，这个算出的数据在0和1.74×10^-7边界处，都没法满足y=0的，所以尺度影响很大。
赞一下 回复此楼
帅得惊动裆中央。
4楼2012-07-13 10:53:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★
dbb627: 金币+1, 感谢应助 2012-07-14 18:53:51
在x=0处y=0，y'=0是作为初值的，肯定是满足起始边界的。

x=1.74×10^-7时，y=-4.55651976657919E-10，y'= 3.98986399474666E-17

也基本满足为0的边界条件了。
赞一下 回复此楼
5楼2012-07-13 11:11:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

azraelgaolei

铁虫 (小有名气)
引用回帖:
5楼: Originally posted by dingd at 2012-07-13 11:11:55
在x=0处y=0，y'=0是作为初值的，肯定是满足起始边界的。

x=1.74×10^-7时，y=-4.55651976657919E-10，y'= 3.98986399474666E-17

也基本满足为0的边界条件了。

x=1.74×10^-7时，y=-4.55651976657919E-10
结果中的y的量级就是E-10，所以-4.55太大了
赞一下 回复此楼
帅得惊动裆中央。
6楼2012-07-13 11:15:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

太微观了，先想办法调整下吧。
赞一下 回复此楼
7楼2012-07-13 11:42:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

azraelgaolei

铁虫 (小有名气)
引用回帖:
7楼: Originally posted by dingd at 2012-07-13 11:42:33
太微观了，先想办法调整下吧。

您好！我已经将尺度改小，变为[0,1000]，方程变为：
2.671972 (w'[x])^2 w''[x] + 7.34` w''''[x] + 5.568`*^-9 == 0,
w[0] == 0, w[1000] == 0, w'[0] == 0, w'[1000] == 0

我能否麻烦您再帮我跑一下，谢谢！感激不尽
赞一下 回复此楼
帅得惊动裆中央。
8楼2012-07-18 14:03:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 azraelgaolei 的主题更新
81/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[基金申请] 刚刚收到科研之友邮件 +16 olivermiaoer 2024-06-19 22/1100 2024-06-20 11:07 by 小木虫吴彦祖分�
[基金申请] 2024国社科通讯评审 +21 qsd10086 2024-06-13 44/2200 2024-06-20 09:22 by 全球漂
[催化] 镍负载氧化铝的保存问题 8+3 lwn0130 2024-06-15 6/300 2024-06-20 09:00 by lwn0130
[论文投稿] 投稿求助 6+3 小短腿安琪拉 2024-06-19 3/150 2024-06-20 08:59 by nono2009
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中！ +10 恐龙爸爸 2024-06-14 10/500 2024-06-20 08:42 by jsguo4127
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
[论文投稿] Scientific Reports二审90天了，催稿也没反应怎么办？ +3 智晨君 2024-06-18 6/300 2024-06-19 20:08 by licheng78
[基金申请] F口信息学部拿面上，大概需要什么样的成果 +7 _奋黎_ 2024-06-16 18/900 2024-06-19 20:02 by _奋黎_
[基金申请] 面上博士后 +3 jsqy 2024-06-19 3/150 2024-06-19 17:58 by bxwisdom_886
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
[访问学者] 国家公派访问学者申请结果出了吗？ +4 65syn 2024-06-13 4/200 2024-06-19 16:40 by 海洋之心168
[基金申请] Nature 11日发文，中国著名学者们称造假迫不得已 +8 babu2015 2024-06-14 8/400 2024-06-19 15:25 by 风今25
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀？ +7 65syn 2024-06-14 11/550 2024-06-19 10:19 by 65syn
[论文投稿] 审稿人含糊拒稿，还需要回复吗？如何回复？ 20+4 BruceChum 2024-06-15 22/1100 2024-06-19 08:00 by kanyechris
[有机交流] 跑板能跑开，过柱过不纯怎么办 +4 小胡在努力 2024-06-18 6/300 2024-06-18 15:20 by long8811
[论文投稿] 投稿求助 10+3 2022_灵魂工程师 2024-06-16 8/400 2024-06-18 12:12 by 投必得科研顾问
[基金申请] 博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[基金申请] 博士后基金需要结题吗？ +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
信息提示
关闭
请填处理意见
关闭