24小时热门版块排行榜    

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

洛羽

木虫 (小有名气)

[求助] 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

非常感谢!
回复此楼

» 猜你喜欢

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

笑笑而已
已阅   回复此楼   关注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的回帖
查看全部 8 个回答

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的回帖

xxz903

金虫 (小有名气)

算了一下  还是不是很好
5楼2013-11-04 21:58:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见