24小时热门版块排行榜    

查看: 1808  |  回复: 7

sic_cyrus

铁虫 (小有名气)

[求助] 急:求救Origin8.5自定义非线性模拟的问题

各位:本人有一套数据,想利用文献中的公式模拟,在Origin8.5中自定义了如下公式:y=10573.04211/((1+L*(8.85*10^(-12)*10573.04211)^3*x^2)^(1/3))+4*P*P*V/(1.38*10^(-23)*298.15)/( 8.85*10^(-12)*(exp(2*P*V*x/(1.38*10^(-23)*298.15))+exp(-2*P*V*x /(1.38*10^(-23)*298.15))+2)),该公式中,L、P、V均是模拟后需要得到的参数,初始值设为3.0e+10、0.001、6.4e-26,但是始终无法模拟成功,要么是不收敛,要么是得不到正确的模拟曲线(模拟得到的曲线明显不对)。原以为三个参数的初始值设置的不对,但当我用这三个值去代替公式中的三个参数,然后根据原x值画出曲线后,与原始数据y却很接近,不知为何?请高手指教(对高手来讲应该比较容易),谢谢!文献中已经模拟出来了。
已经模拟了两周了,还没搞定,很急,十分感谢!
附上原始数据:
x的值
0
100000
200000
300000
400000
500000
600000
700000
800000
对应的y
10573.04211
9934.99353
8639.69837
7612.95107
6606.19154
5898.41003
5162.38239
4730.58808
4275.64697
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sic_cyrus

铁虫 (小有名气)

送鲜花一朵
仍旧谢谢,虽不能解决问题。
2楼2012-09-25 10:28:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sic_cyrus

铁虫 (小有名气)

怎么回复我的那条帖子不见了?
3楼2012-09-25 10:31:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
建议用1stopt 来试试  它的效果不错
4楼2012-09-25 10:35:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sic_cyrus

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by xxz903 at 2012-09-25 10:35:49
建议用1stopt 来试试  它的效果不错

不好意思,我不知道1stopt ,可否详细指明,谢谢!
5楼2012-09-25 10:49:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)

【答案】应助回帖

现在在宿舍里~~软件在实验室  明天去给你代码试一下
6楼2012-09-25 20:18:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)

【答案】应助回帖


dbb627: 金币+1, 感谢应助 2012-09-28 15:13:15
你和效果非常好
迭代数: 111
计算用时(时:分:秒:微秒): 00:00:16:375
优化算法: 通用全局优化法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 74.2281592974862
残差平方和(SSE): 49588.3766942369
相关系数(R): 0.999499585535417
相关系数之平方(R^2): 0.998999421485471
决定系数(DC): 0.998812026002177
卡方系数(Chi-Square): 4.43440469554311
F统计(F-Statistic): 2525.30779755917

参数        最佳估算
----------        -------------
l        23967668573.1521
p        13727.7872635405
v        -6.25812615721438E-39

====== 结果输出 =====

No        实测值y        计算值y
1        10573.04211        10540.6538047424
2        9934.99353        9927.30271919463
3        8639.69837        8682.95529425934
4        7612.95107        7498.57436590191
5        6606.19154        6551.32678871753
6        5898.41003        5815.88889375
7        5162.38239        5239.02078578504
8        4730.58808        4777.50406459364
9        4275.64697        4400.67737305877

» 本帖已获得的红花(最新10朵)

7楼2012-09-26 09:01:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sic_cyrus

铁虫 (小有名气)

送鲜花一朵
引用回帖:
7楼: Originally posted by xxz903 at 2012-09-26 09:01:29
你和效果非常好
迭代数: 111
计算用时(时:分:秒:微秒): 00:00:16:375
优化算法: 通用全局优化法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 74.2281592974862
残差平方和(SSE): 49588.3766942369
...

先谢谢你这么热心!
不过,得到的参数值是不合适的。三个值都须大于0,都在我给的初始值附近,差不多1~3个数量级范围内吧。
8楼2012-09-28 08:26:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sic_cyrus 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 17:45 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 14/700 2026-08-30 17:01 by hZiFeudZyoGR
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 13:54 by GkCVyODlx4mR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 07:05 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
信息提示
请填处理意见