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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt拟合问题求助
81/1返回列表
查看: 1484  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

洛羽

木虫 (小有名气)

[求助] 1stopt拟合问题求助

用1stopt 1.5拟合，得到全零解或者拟合度很差的解，想请有高版本的帮忙运行下，看看是软件的问题还是方程设计的问题，代码如下：
Title "iPP-SL Equation";
Parameters D0,P0,T0;
Variable T,D,P;
Function P=P0*(-(D/D0)^2-T/T0*(LN(1-D/D0)+(1-8314.5*T0*D0/52000/P)*D/D0));
Data;
443.8        784.6829881        0.1
453.6        763.4170547        0.1
473.5        752.2190462        0.1
493.5        741.4547342        0.1
513.4        730.4068366        0.1
533.6        719.5279896        0.1
553.7        708.2654579        0.1
573.6        696.961249        0.1
443.8        805.3475074        10
453.6        772.797527        10
473.5        762.7183281        10
493.5        752.8419785        10
513.4        742.8868583        10
533.6        732.976618        10
553.7        722.9612493        10
573.6        712.9616427        10
443.8        823.0452675        20
453.6        781.6774799        20
473.5        772.3200494        20
493.5        763.1257631        20
513.4        753.9772299        20
533.6        744.712541        20
553.7        735.6727728        20
573.6        726.4802034        20
443.8        864.7526807        50
453.6        803.5355564        50
473.5        795.1653944        50
493.5        787.3395796        50
513.4        779.4839816        50
533.6        771.8431615        50
553.7        764.0586797        50
573.6        756.3724378        50
443.8        900.4952724        100
453.6        830.7027745        100
473.5        823.6553826        100
493.5        816.7932696        100
513.4        810.1758082        100
533.6        803.4709947        100
553.7        797.1938776        100
573.6        790.6388362        100
443.8        919.0331771        150
453.6        852.151683        150
473.5        845.5944529        150
493.5        839.419122        150
513.4        833.2638947        150
533.6        827.2667108        150
553.7        821.3552361        150
573.6        815.7272208        150

非常感谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
笑笑而已
1楼 2013-10-25 16:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dingd: 金币+1, 谢谢相助。 2013-10-25 23:20:02
洛羽: 金币+30, ★★★★★最佳答案, Thx！ 2013-10-26 08:37:05
洛羽: 金币+10, ★★★★★最佳答案, 谢了哈 2013-11-05 08:30:00
5.0的运算结果

迭代数: 16
计算用时(时:分:秒:毫秒): 00:00:05:955
计算中止原因: 达到收敛判定标准
优化算法: 通用全局优化法(UGO1)
目标函数值(最小): 58383392.4877615
均方差(RMSE): 48.0442103004379
残差和(SSE): 110795.81488285
相关系数(R): 0.488662637237418
决定系数(DC): 0.201253232358661

参数        最佳估算
----------        -------------
d0        -2.77787802237206
p0        -0.000463828584826673
t0        -0.563656106249702

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

目标p        计算p
0.1        0.40410567441027
0.1        0.394358307221224
0.1        0.380577948472944
0.1        0.372127307528996
0.1        0.363569400130792
0.1        0.355903224888798
0.1        0.349110103455253
0.1        0.343253646743498
10        64.4567453353209
10        65.5251730522155
10        69.496642865531
10        73.4206671302621
10        77.1608254892335
10        80.8464186725441
10        84.4042235274274
10        87.7822080787686
20        64.9285322093494
20        65.7385790640006
20        69.8227060687401
20        73.884560424076
20        77.7448585842808
20        81.5302809995688
20        85.192293762852
20        88.7327391189727
50        66.0575791958963
50        66.5711350827895
50        70.8652952654182
50        75.1160106781415
50        79.2335046056866
50        83.3091526370455
50        87.2290531905943
50        91.079799122838
100        67.0838569831947
100        67.5625758741098
100        72.0877039163822
100        76.5808409555551
100        80.975773555208
100        85.3349383219211
100        89.5912248773237
100        93.7252973209497
150        67.5002348670369
150        68.2952674713412
150        73.0033399654263
150        77.670453732329
150        82.2182762696684
150        86.7648139404564
150        91.2244806288243
150        95.5579799554526

====== 计算结束 ======
赞一下 回复此楼
2楼2013-10-25 21:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)
应该是模型问题~~你再检查一下
赞一下 回复此楼
3楼2013-10-25 21:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

洛羽

木虫 (小有名气)
★
fegg7502: 金币+1, 鼓励交流 2013-11-05 08:54:44
引用回帖:
3楼: Originally posted by xxz903 at 2013-10-25 21:26:25
应该是模型问题~~你再检查一下

能再帮忙运行下程序吗，谢谢啦，这次是能拟合出结果，但解不稳定，想知道高版本能不能得到稳定解。
Title "iPP-SL Equation";
Parameters D0,P0,T0;
Variable T,D,P;
Function D=D0*(1-exp(-((D/D0)^2+P/P0)/T*T0-(1-8314.5*T0*D0/52000/P0)*D/D0));
Data;
443.8        784.6829881        0.1
453.6        763.4170547        0.1
473.5        752.2190462        0.1
493.5        741.4547342        0.1
513.4        730.4068366        0.1
533.6        719.5279896        0.1
553.7        708.2654579        0.1
573.6        696.961249        0.1
443.8        805.3475074        10
453.6        772.797527        10
473.5        762.7183281        10
493.5        752.8419785        10
513.4        742.8868583        10
533.6        732.976618        10
553.7        722.9612493        10
573.6        712.9616427        10
443.8        823.0452675        20
453.6        781.6774799        20
473.5        772.3200494        20
493.5        763.1257631        20
513.4        753.9772299        20
533.6        744.712541        20
553.7        735.6727728        20
573.6        726.4802034        20
443.8        864.7526807        50
453.6        803.5355564        50
473.5        795.1653944        50
493.5        787.3395796        50
513.4        779.4839816        50
533.6        771.8431615        50
553.7        764.0586797        50
573.6        756.3724378        50
443.8        900.4952724        100
453.6        830.7027745        100
473.5        823.6553826        100
493.5        816.7932696        100
513.4        810.1758082        100
533.6        803.4709947        100
553.7        797.1938776        100
573.6        790.6388362        100
443.8        919.0331771        150
453.6        852.151683        150
473.5        845.5944529        150
493.5        839.419122        150
513.4        833.2638947        150
533.6        827.2667108        150
553.7        821.3552361        150
573.6        815.7272208        150
443.8        932.2270905        200
453.6        872.2958828        200
473.5        863.8562543        200
493.5        858.0744809        200
513.4        852.442247        200
533.6        846.9551961        200
553.7        841.4675194        200
573.6        836.1903169        200
赞一下 回复此楼
笑笑而已
4楼2013-11-04 17:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)
算了一下  还是不是很好
赞一下 回复此楼
5楼2013-11-04 21:58:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zee2004

木虫 (正式写手)
引用回帖:
4楼: Originally posted by 洛羽 at 2013-11-04 17:22:09
能再帮忙运行下程序吗，谢谢啦，这次是能拟合出结果，但解不稳定，想知道高版本能不能得到稳定解。
Title "iPP-SL Equation";
Parameters D0,P0,T0;
Variable T,D,P;
Function D=D0*(1-exp(-((D/D0) ...

"iPP-SL Equation"

迭代数: 2850
计算用时(时:分:秒:微秒): 00:00:00:938
优化算法: 遗传算法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 8.75660320527954E-6
残差平方和(SSE): 4.29397358290387E-9
相关系数(R): 1
相关系数之平方(R^2): 1
决定系数(DC): 0.999999999999973
卡方系数(Chi-Square): 2.67019528050438E-12
F统计(F-Statistic): 1.21597189939003E17

参数        最佳估算
----------        -------------
D0                 2512710778.98991
P0                 -594462605166871
T0                 0.237047461528293

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

No        实测值d        计算值d
1        784.6829881        784.6829914
2        763.4170547        763.4170612
3        752.2190462        752.2190541
4        741.4547342        741.4547436
5        730.4068366        730.4068476
6        719.5279896        719.5280019
7        708.2654579        708.2654717
8        696.961249        696.9612642
9        805.3475074        805.3475076
10        772.797527        772.7975320
11        762.7183281        762.7183347
12        752.8419785        752.8419865
13        742.8868583        742.8868675
14        732.976618        732.9766285
15        722.9612493        722.9612613
16        712.9616427        712.9616558
17        823.0452675        823.0452646
18        781.6774799        781.6774836
19        772.3200494        772.3200544
20        763.1257631        763.1257695
21        753.9772299        753.9772377
22        744.712541        744.7125499
23        735.6727728        735.6727831
24        726.4802034        726.4802147
25        864.7526807        864.7526705
26        803.5355564        803.5355565
27        795.1653944        795.1653960
28        787.3395796        787.3395823
29        779.4839816        779.4839856
30        771.8431615        771.8431665
31        764.0586797        764.0586860
32        756.3724378        756.3724450
33        900.4952724        900.4952554
34        830.7027745        830.7027703
35        823.6553826        823.6553795
36        816.7932696        816.7932677
37        810.1758082        810.1758074
38        803.4709947        803.4709949
39        797.1938776        797.1938787
40        790.6388362        790.6388384
41        919.0331771        919.0331562
42        852.151683        852.1516749
43        845.5944529        845.5944460
44        839.419122        839.4191162
45        833.2638947        833.2638899
46        827.2667108        827.2667069
47        821.3552361        821.3552333
48        815.7272208        815.7272190
49        932.2270905        932.2270668
50        872.2958828        872.2958709
51        863.8562543        863.8562441
52        858.0744809        858.0744716
53        852.442247        852.4422388
54        846.9551961        846.9551887
55        841.4675194        841.4675130
56        836.1903169        836.1903116
赞一下 回复此楼
MC
6楼2013-11-08 13:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zee2004

木虫 (正式写手)
"iPP-SL Equation"

迭代数: 2850
计算用时(时:分:秒:微秒): 00:00:00:938
优化算法: 遗传算法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 8.75660320527954E-6
残差平方和(SSE): 4.29397358290387E-9
相关系数(R): 1
相关系数之平方(R^2): 1
决定系数(DC): 0.999999999999973
卡方系数(Chi-Square): 2.67019528050438E-12
F统计(F-Statistic): 1.21597189939003E17

参数        最佳估算
----------        -------------
D0                 2512710778.98991
P0                 -594462605166871
T0                 0.237047461528293

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

No        实测值d        计算值d
1        784.6829881        784.6829914
2        763.4170547        763.4170612
3        752.2190462        752.2190541
4        741.4547342        741.4547436
5        730.4068366        730.4068476
6        719.5279896        719.5280019
7        708.2654579        708.2654717
8        696.961249        696.9612642
9        805.3475074        805.3475076
10        772.797527        772.7975320
11        762.7183281        762.7183347
12        752.8419785        752.8419865
13        742.8868583        742.8868675
14        732.976618        732.9766285
15        722.9612493        722.9612613
16        712.9616427        712.9616558
17        823.0452675        823.0452646
18        781.6774799        781.6774836
19        772.3200494        772.3200544
20        763.1257631        763.1257695
21        753.9772299        753.9772377
22        744.712541        744.7125499
23        735.6727728        735.6727831
24        726.4802034        726.4802147
25        864.7526807        864.7526705
26        803.5355564        803.5355565
27        795.1653944        795.1653960
28        787.3395796        787.3395823
29        779.4839816        779.4839856
30        771.8431615        771.8431665
31        764.0586797        764.0586860
32        756.3724378        756.3724450
33        900.4952724        900.4952554
34        830.7027745        830.7027703
35        823.6553826        823.6553795
36        816.7932696        816.7932677
37        810.1758082        810.1758074
38        803.4709947        803.4709949
39        797.1938776        797.1938787
40        790.6388362        790.6388384
41        919.0331771        919.0331562
42        852.151683        852.1516749
43        845.5944529        845.5944460
44        839.419122        839.4191162
45        833.2638947        833.2638899
46        827.2667108        827.2667069
47        821.3552361        821.3552333
48        815.7272208        815.7272190
49        932.2270905        932.2270668
50        872.2958828        872.2958709
51        863.8562543        863.8562441
52        858.0744809        858.0744716
53        852.442247        852.4422388
54        846.9551961        846.9551887
55        841.4675194        841.4675130
56        836.1903169        836.1903116
赞一下 回复此楼
MC
7楼2013-11-08 13:29:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

洛羽

木虫 (小有名气)
引用回帖:
6楼: Originally posted by zee2004 at 2013-11-08 13:29:10
"iPP-SL Equation"

迭代数: 2850
计算用时(时:分:秒:微秒): 00:00:00:938
优化算法: 遗传算法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 8.75660320527954E-6
残差平方和(SSE): 4.29397358 ...

1stopt给出的拟合值太大且有负值，似乎还不稳定，重新运行一遍得到的值完全不同，文献上给出的参数为885.6        300.7
690.6，不知1stopt能否限定拟合值的范围？
赞一下 回复此楼
笑笑而已
8楼2013-11-08 14:38:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 洛羽 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭