24小时热门版块排行榜    

查看: 1080  |  回复: 6

钟声入翠微

铜虫 (小有名气)

[求助] 1stOpt非线性拟合自定义函数(gamma函数)已有1人参与

显示函数表达式错误,多次修改未成功,函数如附件所示
Title "Type your title here";
Parameters a,b;
Variable t,y;
Function y= t^(a-1) * exp(-t/b)*b^(-a) /(int(t^(a-1)*exp(-t),0,inf)  ;
Data;
0
12        -6.08
79        -7.32
120        -8.27
144        -9.38
172        -7.26
206        -9.03
230        -8.19
255        -7.00
351        -5.10
395        -5.70
425        -6.15
445        -6.66
486        -7.28
516        -7.77
545        -7.31
573        -7.02
606        -8.52
632        -8.07
665        -4.24
699        -3.24
723        -4.84
752        -5.97
789        -6.06
815        -7.00
851        -7.22
881        -7.65
912        -6.96
940        -6.86
966        -7.80
988        -6.78

1stOpt非线性拟合自定义函数(gamma函数)
3.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
钟声入翠微: 金币+30, ★★★★★最佳答案 2021-05-25 10:58:13
Parameters a,b;
Variable t,y;
Function y= t^(a-1)*exp(-t/b)/b^a/Gamma(a);
Data;
0
12        -6.08
......


***********************
均方差(RMSE): 6.51708182055846
残差平方和(SSR): 1316.64301913146
相关系数(R): 0.00732870945702214
相关系数之平方(R^2): 5.37099823054458E-5
修正R平方(Adj. R^2): -0.0713710250189585
确定系数(DC): -12.2545266162522
卡方系数(Chi-Square): -171.789745491722
F统计(F-Statistic): 26.8504257160124

参数        最佳估算
----------        -------------
a        -0.01320043171917
b        1.2315680815176E308

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

No.        实测值 y        计算值 y
1        0        0
2        -6.08        -12.3230025991806
3        -7.32        -1.82585740146398
4        -8.27        -1.19540786357906
5        -9.38        -0.993778589151446
6        -7.26        -0.830051534835844
7        -9.03        -0.691404461863923
8        -8.19        -0.618357711796463
9        -7        -0.556975246943513
10        -5.1        -0.402937024657334
11        -5.7        -0.357495143350713
12        -6.15        -0.331939278614345
13        -6.66        -0.316828278804579
14        -7.28        -0.289762657873544
15        -7.77        -0.272700287427244
16        -7.31        -0.258003339305755
17        -7.02        -0.245233609702048
18        -8.52        -0.231707974303466
19        -8.07        -0.22205251255645
20        -4.24        -0.210891625758833
21        -3.24        -0.20050164587487
22        -4.84        -0.19375962591301
23        -5.97        -0.186190829461552
24        -6.06        -0.177346970499099
25        -7        -0.171615811595814
26        -7.22        -0.164262167534759
27        -7.65        -0.158596128921691
28        -6.96        -0.153135329168843
29        -6.86        -0.148514555472775
30        -7.8        -0.144465229327475
数值计算
2楼2021-05-24 20:47:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2021-05-24 20:47:15
Parameters a,b;
Variable t,y;
Function y= t^(a-1)*exp(-t/b)/b^a/Gamma(a);
Data;
0
12        -6.08
......


***********************
均方差(RMSE): 6.51708182055846
残差平方和(SSR): 1316.6430 ...

谢谢大神!还有一个小问题,我该怎么定义参数a和b的范围呢?他们都是大于0的,
3楼2021-05-25 10:01:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 钟声入翠微 at 2021-05-25 10:01:07
谢谢大神!还有一个小问题,我该怎么定义参数a和b的范围呢?他们都是大于0的,...

模型和数据对不上,可能模型不适合。


Parameters a,b;
ParameterDomain = [0, ];
Variable t,y;
Function y= t^(a-1)*exp(-t/b)/b^a/Gamma(a);
Data;
0         0
.......

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

数值计算
4楼2021-05-25 10:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by 独孤神宇 at 2021-05-25 10:21:14
模型和数据对不上,可能模型不适合。


Parameters a,b;
ParameterDomain = ;
Variable t,y;
Function y= t^(a-1)*exp(-t/b)/b^a/Gamma(a);
Data;
0         0
..........

好的,谢谢大神!那么在拟合过程中,拟合出现了参数值、相关性和均方根误差,但是拟合的y值都为0,是不是就是说明模型不适合,没有拟合成功呢?
5楼2021-05-25 10:24:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 钟声入翠微 at 2021-05-25 10:24:23
好的,谢谢大神!那么在拟合过程中,拟合出现了参数值、相关性和均方根误差,但是拟合的y值都为0,是不是就是说明模型不适合,没有拟合成功呢?...

是的,拟合效果就是看统计指标

发自小木虫Android客户端
数值计算
6楼2021-05-25 10:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)

好的,谢谢大神!
7楼2021-05-25 10:57:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 钟声入翠微 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 想被211以上高校课题组接收 +11 风起沧澜 2024-05-16 13/650 2024-05-22 14:37 by keyaner23
[论文投稿] 为什么有的影响因子高的期刊分区不高呢? +4 安处一室 2024-05-21 4/200 2024-05-22 14:19 by cqupenghao
[硕博家园] 超过35岁的海外博士还可以进高校吗? +4 905452934 2024-05-20 5/250 2024-05-22 14:13 by xmc_ues_2015
[教师之家] 经常觉得挺累的 +3 zylfront 2024-05-22 3/150 2024-05-22 14:00 by 学酥一枚吖
[基金申请] 又有两个中科院三、四区计算机领域SCI/SSCI期刊爆出问题了! +4 瞬息宇宙 2024-05-19 6/300 2024-05-22 12:37 by yp_nupt
[论文投稿] 投稿journal of energy chemistry +3 小可爱嘿呀 2024-05-21 3/150 2024-05-22 11:11 by xs74101122
[考博] 化工三篇sci只想有书读,急寻博导! +3 8139173 2024-05-21 5/250 2024-05-22 10:56 by 安塔瓦拉多
[考博] 大扩招 +9 version211 2024-05-21 11/550 2024-05-22 09:59 by lu112175
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +3 第二种态度 2024-05-18 6/300 2024-05-21 15:37 by 安小樱
[论文投稿] 论文一审意见回来后发现实验程序编错了论证分析部分可能要大改 5+4 hshhenb 2024-05-20 5/250 2024-05-21 11:03 by bnullh
[论文投稿] Angew拒稿重投经历大修返回又送审了 48+4 YOONAZHA 2024-05-18 6/300 2024-05-21 10:35 by jurkat.1640
[基金申请] 基金委也是打工的 +25 wsjing 2024-05-16 34/1700 2024-05-21 08:31 by mrKiller
[基金申请] 申请基金代表性成果 +14 lancet0903 2024-05-17 20/1000 2024-05-21 00:23 by dxcharlary
[找工作] 绍兴文理学院怎么样?有没有坑啊 +9 zhaojiang427 2024-05-16 23/1150 2024-05-20 21:44 by 高敖曹
[有机交流] 除DMSO +5 Spiralup 2024-05-15 6/300 2024-05-20 17:49 by 刘洪振
[教师之家] 普通高校有编制讲师与211高校无编制专职科研岗,应该怎么选 +19 1137025804 2024-05-15 25/1250 2024-05-20 13:32 by Zhangsanxi
[论文投稿] 求推荐期刊 20+3 好困好困a 2024-05-18 4/200 2024-05-19 11:30 by nono2009
[基金申请] 大佬们,哪位大神知道E09有没有送啊? +4 牛棕色左右 2024-05-15 7/350 2024-05-18 22:48 by 牛棕色左右
[电化学] 锂离子电池石墨负极用 1M LiPF6 in DEC:EC=1:1 Vol% 可以吗? 50+3 fffhhhhh 2024-05-15 8/400 2024-05-17 14:57 by 多点关心多点i
[基金申请] 系统导出申报书中参与人那一页的格式,跨入两页,要通过调整行间距调整到一页吗? +3 danisha 2024-05-15 5/250 2024-05-16 18:21 by danisha
信息提示
请填处理意见