24小时热门版块排行榜    

查看: 1164  |  回复: 10

yanjia119

铁虫 (初入文坛)

[求助] 1stopt,我用最佳参数带回模型里面运算得不到运算结果的计算值啊?, 已有1人参与

程序如下,我把得到的最佳参数abcd带到方程 z=a-b*log(x+exp((a-(c-d*log(y+0.5)))/b))里面去计算z得不到最下面结果输出里面的计算值z,差距很大,我算出来很多负数。

Variable   x,y,z;
Parameters   a,b,c,d;
Function    z=a-b*log(x+exp((a-(c-d*log(y+0.5)))/b));
data ;
0        0        1
10        1        0.810000000000000
100        1        0.600000000000000
1000        1        0.345000000000000
10        10        0.418000000000000
100        10        0.360000000000000
1000        10        0.291000000000000
10        100        0.229000000000000
100        100        0.217000000000000
1000        100        0.200000000000000
10        1000        0.109000000000000
100        1000        0.103000000000000
1000        1000        0.100000000000000

迭代数: 44
计算用时(时:分:秒:微秒): 00:00:01:999
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0760567087005865
残差平方和(SSE): 0.0752000981987563
相关系数(R): 0.959037618081535
相关系数之平方(R^2): 0.919753152895504
决定系数(DC): 0.919753152895504
卡方系数(Chi-Square): 0.15024761250399
F统计(F-Statistic): 126.077036630187

参数        最佳估算
----------        -------------
a                  1.27367539826881
b                  0.319548736974564
c                  0.209323754540026
d                  0.584710951477248

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

No        实测值z        计算值z
1        1        0.8878759
2        0.81        0.7347398
3        0.6        0.5892871
4        0.345        0.3097744
5        0.418        0.5446592
6        0.36        0.4911163
7        0.291        0.2919232
8        0.229        0.3017430
9        0.217        0.2908335
10        0.2        0.2126957
11        0.109        0.0493632
12        0.103        0.0475338
13        0.1        0.0304548
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
1\可以得到
No        实测值z        计算值z
1        1        0.8878759
2        0.81        0.7347398
3        0.6        0.5892871
4        0.345        0.3097744
5        0.418        0.5446592
6        0.36        0.4911163
7        0.291        0.2919232
8        0.229        0.3017430
9        0.217        0.2908335
10        0.2        0.2126957
11        0.109        0.0493632
12        0.103        0.0475338
13        0.1        0.0304548
2\公式与数据匹配一般,拟合效果一般。
一团和气
2楼2017-04-24 09:57:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
3楼2017-08-18 04:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
4楼2017-08-18 19:03:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2017-08-19 03:15:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
6楼2017-08-19 08:55:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
7楼2017-08-20 07:08:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
8楼2017-08-20 18:34:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
9楼2017-08-21 01:49:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
10楼2017-08-21 02:29:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yanjia119 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:42 by 4IW0sJvtEqX8
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[教师之家] 课题组招2027级博士 上海工程技术大学 激光智能制造方向 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:37 by oxidpl
[有机交流] 有机合成以后会不会被AI改变?做科研的虫友怎么看 +5 zapen 2026-09-20 8/400 2026-09-23 17:21 by tianxiaoxian
信息提示
请填处理意见