24小时热门版块排行榜    

查看: 1074  |  回复: 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普通老师 +4 timergoi 2024-05-14 4/200 2024-05-15 00:17 by crazy peng
[基金申请] 面上项目代表作必须有通讯作者文章吗? +25 zlxnnx111 2024-05-09 41/2050 2024-05-14 20:45 by qweaqq
[基金申请] 河南工业大学过年开学后,就不让用国基金报账,这合法合规吗? 100+8 知己而已 2024-05-09 36/1800 2024-05-14 20:37 by Math露珠
[基金申请] 这个教授的科研做得真好 +3 zju2000 2024-05-14 4/200 2024-05-14 18:10 by lwliwei1989
[论文投稿] 中文投稿 +5 liujinyu80 2024-05-13 5/250 2024-05-14 17:13 by icm639
[有机交流] 无水硫酸钠除水 80+3 桃桃PXS 2024-05-13 4/200 2024-05-14 16:49 by binzy
[硕博家园] 博导推荐,申博,博士招生 +3 大林挺好 2024-05-13 4/200 2024-05-14 16:10 by 安塔瓦拉多
[教师之家] 谁要做个人网站,可以交流啊! +6 hjc404 2024-05-14 12/600 2024-05-14 15:27 by wjykycg
[论文投稿] 投稿竟然倒在了Scientific Reports上。 +11 lizhengke06 2024-05-10 13/650 2024-05-14 13:17 by lyfbangong
[硕博家园] 哈工大硕博招生! +4 nailooo 2024-05-12 5/250 2024-05-13 22:04 by yuanjijoy
[硕博家园] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +9 Kelaizhang 2024-05-08 17/850 2024-05-13 19:42 by LittleBush
[基金申请] 再谈一作和二作 +19 jklily 2024-05-10 33/1650 2024-05-13 17:12 by mierbushiwo
[基金申请] 祈福第三次青基 +17 lylylyheart 2024-05-09 26/1300 2024-05-13 15:57 by lylylyheart
[论文投稿] Journal of Electrical Engineering&Technology Reviews Completed 快一周了 10+4 qweasd12345 2024-05-12 5/250 2024-05-13 10:00 by bear2007
[考博] 准研三生态毒理学,爱化学生物新能源碳中和,电池等,25博士岗位求捞啦! +4 wangchensi 2024-05-10 14/700 2024-05-13 00:26 by 大林挺好
[基金申请] 国社科系统提交后,下载申请书,里面的思路图不显示,是什么问题 +3 十三画福将 2024-05-09 3/150 2024-05-12 23:01 by 老虎当猫养
[基金申请] NSFC函评几月几日截止提交评审意见? +9 瞬息宇宙 2024-05-11 9/450 2024-05-12 19:53 by 20081002
[基金申请] 科研之友老是给我发消息 +6 问君611 2024-05-11 8/400 2024-05-12 17:24 by 淀粉搬运工
[硕博家园] 科研不是打打杀杀,科研就是人情世故 +10 杞天大圣 2024-05-08 10/500 2024-05-11 09:05 by 贪吃fish
[有机交流] 求助酰胺的合成 +4 洛悠灬 2024-05-08 5/250 2024-05-09 11:29 by 小肉干
信息提示
请填处理意见