【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1055  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 钟声入翠微 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 想问一下有没有需要文章但是没时间带学生的研究生导师。 +9 lekinna 2024-04-19 18/900 2024-04-20 01:33 by 斯文赫
[公派出国] CSC在职申请求问 +7 Amber0919 2024-04-16 7/350 2024-04-19 21:37 by 小狐_小胡
[育儿交流] 小孩7岁,上一年级,不自信,学习较差 10+7 ljaxida 2024-04-16 34/1700 2024-04-19 18:49 by 水可南阳
[找工作] 事业单位还是大学好? +17 青萍之沫 2024-04-16 18/900 2024-04-19 17:45 by charles-c
[基金申请] 基金和生小孩 +33 Ausy 2024-04-15 35/1750 2024-04-19 16:21 by feng6531
[有机交流] 傅克酰基化,产率大于百分之一百,求解,很急 90+5 hsn991013 2024-04-15 11/550 2024-04-19 14:49 by scdxyouji
[基金申请] 特别资助审核状态 +7 pantray 2024-04-17 7/350 2024-04-19 13:51 by Breezy恋
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +4 lin_2018 2024-04-17 4/200 2024-04-19 08:46 by TsolmonZ
[无机/物化] H3PO2在130度附近分解为磷化氢(剧毒),怎样防护,求指点。 1(金币+5) +3 大力2010 2024-04-13 9/450 2024-04-17 21:55 by 大力2010
[考研] 工学322,英语六级求调剂 +3 domax 2024-04-17 10/500 2024-04-17 21:25 by 东方花园小万
[基金申请] 迟国泰通过向学生发放劳务费再回收的方式套取科学基金重点项目 +6 babu2015 2024-04-13 7/350 2024-04-16 20:32 by sundiv
[考研] 求调剂 +4 鹿萌月yy 2024-04-15 5/250 2024-04-16 16:16 by chenweiwade
[考研] 294求调剂 +3 694062003 2024-04-15 4/200 2024-04-16 15:01 by 邹邹哈哈
[考研] 334求调剂 +4 学药救人 2024-04-14 4/200 2024-04-15 15:05 by hunanzang
[考研] 322求调剂 +6 本己上岸 2024-04-15 7/350 2024-04-15 14:19 by mthwyj
[材料综合] 请问哪里可以测试低压吸氢的PCT曲线 +3 yunshengcd 2024-04-14 5/250 2024-04-15 13:59 by zqdsb
[考研] 287求调剂 +6 南沨 2024-04-14 6/300 2024-04-14 23:08 by lincunhui
[考研] 323求调剂 +15 啊Q精神~ 2024-04-13 17/850 2024-04-14 12:44 by qjhawk
[考研] 338求调剂 +3 18280338551 2024-04-14 5/250 2024-04-14 10:03 by tcni
[考研] 086000生医没考数学304分,一志愿西农求调剂 +3 王有田有学上 2024-04-13 5/250 2024-04-13 16:18 by coco1981
信息提示
请填处理意见