24小时热门版块排行榜    

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

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值。
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

帅得惊动裆中央。
已阅   回复此楼   关注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的回帖
查看全部 8 个回答

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


参数       ...

而且您看,这个算出的数据在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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +45 医学老男孩 2026-08-13 101/5050 2026-08-18 13:40 by 医学老男孩
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KZyU1THERRVo 2026-08-18 4/200 2026-08-18 13:28 by vHWVV9D5Deb5
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见