24小时热门版块排行榜    

查看: 1282  |  回复: 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 的主题更新
信息提示
请填处理意见