24小时热门版块排行榜    

查看: 2074  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 钟声入翠微 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +41 医学老男孩 2026-08-13 94/4700 2026-08-18 10:20 by 医学老男孩
[论文投稿] 投稿咨询 +3 wwm09 2026-08-17 4/200 2026-08-18 09:52 by xs74101122
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 5/250 2026-08-18 08:03 by 基诺咪客
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见