24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 857  |  回复: 8

茶禅子

新虫 (正式写手)

[求助] 帮忙用1stOpt跑一下数据~~~MATLAB出错跑不出来已有1人参与

数学模型为y=a*[sinh(bx)]^b   a、b、c为系数,模型是基于双曲正弦函数的,数据如下。用matlab拟合不出来,总是出错。求帮下忙用1stopt分别跑一下下面三组数据得出a、b、c值。。。万分感谢~~~!!!
第一组:
  x                   y
200                65.23447
459                60.36114
997                59.05684
2000        58.05068
5015        55.4592
9078        56.00643
第二组:
x                                     y
25.00002                        68.15676
57.375                        65.99296
124.62498                61.14024
250.0002                        56.03819
626.87524                51.89158
1134.75061                50.46947
1875.00598                49.12768
2500.00195                48.95017
3124.99805                49.30959
3749.9939                49.52093
第三组:
x                                              y
3.125                                34.06138
7.17188                                58.20715
15.57812                                71.6804
31.25002                                81.03972
78.35941                                91.04919
141.84383                        60.65113
234.37575                        56.37441
312.50024                        54.54849
390.62476                        53.16828
468.74924                        51.08537
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

怎么 又是你 ,你matlab 提示 什么错误?
数值计算
2楼2018-12-12 20:11:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
茶禅子: 金币+33, ★★★★★最佳答案 2018-12-12 20:25:26
第一组 :均方差(RMSE): 1.03387792236279
残差平方和(SSR): 6.41342135009515
相关系数(R): 0.947853347298768
相关系数之平方(R^2): 0.898425967985478
修正R平方(Adj. R^2): 0.830709946642464
确定系数(DC): 0.898392774282132
卡方系数(Chi-Square): 0.0532666001387005
F统计(F-Statistic): 13.4247261315101

参数        最佳估算
----------        -------------
a        14.79528351424
b        5.66088918227076E-19
c        -0.0397738716506555

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

No.        实测值 y        计算值 y
1        65.23447        63.7330900479526
2        60.36114        61.6608774656741
3        59.05684        59.7876456011801
4        58.05068        58.155043030452
5        55.4592        56.0670221504062
6        56.00643        54.7592147520807

第二组:

均方差(RMSE): 1.36599478654719
残差平方和(SSR): 18.6594175687411
相关系数(R): 0.98104001414782
相关系数之平方(R^2): 0.962439509359155
修正R平方(Adj. R^2): 0.951707940604628
确定系数(DC): 0.962407772851148
卡方系数(Chi-Square): 0.174860579780483
F统计(F-Statistic): 90.6312777054339

参数        最佳估算
----------        -------------
a        4.68584011883617
b        4.04862941341875E-18
c        -0.0727187067456189

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

No.        实测值 y        计算值 y
1        68.15676        68.2179633487173
2        65.99296        64.2204563151755
3        61.14024        60.6976600481208
4        56.03819        57.7013036048222
5        51.89158        53.97016375581
6        50.46947        51.6907496988387
7        49.12768        49.8370843729887
8        48.95017        48.8053494675879
9        49.30959        48.0197863405008
10        49.52093        47.3873360915942

第三组:
均方差(RMSE): 14.6251920369
残差平方和(SSR): 2138.96242116202
相关系数(R): 0.315050166243831
相关系数之平方(R^2): 0.0992566072502655
修正R平方(Adj. R^2): -0.158098647821087
确定系数(DC): 0.0992066961226141
卡方系数(Chi-Square): 17.7026969550528
F统计(F-Statistic): 0.368607471400822

参数        最佳估算
----------        -------------
a        66.1751724444395
b        1.51568429187171
c        -0.000319653899523198

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

No.        实测值 y        计算值 y
1        34.06138        66.0896993778328
2        58.20715        65.960243250757
3        71.6804        65.6921482607219
4        81.03972        65.1952401009426
5        91.04919        63.7240617227339
6        60.65113        61.793883105336
7        56.37441        59.084769513297
8        54.54849        56.8901551916159
9        53.16828        54.7770557005547
10        51.08537        52.7424451430111

第三组 效果不好,检查一下数据。
数值计算
3楼2018-12-12 20:17:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶禅子

新虫 (正式写手)

嗯···又是我,matlab错误是  Inf computed by model function,fitting cnnot continue.Try using or tightening upper and lower bounds on coefficients。。。有时候能出来,但是参数值都是乱的 @beefly
4楼2018-12-12 20:27:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶禅子

新虫 (正式写手)

引用回帖:
3楼: Originally posted by 独孤神宇 at 2018-12-12 20:17:37
第一组 :均方差(RMSE): 1.03387792236279
残差平方和(SSR): 6.41342135009515
相关系数(R): 0.947853347298768
相关系数之平方(R^2): 0.898425967985478
修正R平方(Adj. R^2): 0.830709946642464
确定系数(DC ...

第三组可能第一个数据误差值太大了...
5楼2018-12-12 20:30:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

你用的cftool?多试几组初始值

发自小木虫Android客户端
数值计算
6楼2018-12-12 20:32:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶禅子

新虫 (正式写手)

引用回帖:
6楼: Originally posted by 独孤神宇 at 2018-12-12 20:32:55
你用的cftool?多试几组初始值

对用的CFTOOL ,试过不同初始值还是一直出不来,不知道是不是算法的问题, cftool 里面只有Levenberg-Marquardt和  Trust-Region可选,1stOpt是不是算法比MATLAB强大许多。。。
7楼2018-12-12 20:43:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶禅子

新虫 (正式写手)

引用回帖:
6楼: Originally posted by 独孤神宇 at 2018-12-12 20:32:55
你用的cftool?多试几组初始值

这出的结果完全就是错的。。。
帮忙用1stOpt跑一下数据~~~MATLAB出错跑不出来
1.jpg

8楼2018-12-12 20:44:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

你可以将1stOpt的结果作为初始值试一下。

发自小木虫Android客户端
数值计算
9楼2018-12-12 21:49:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 茶禅子 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金开始函评了吗? +14 wych1103 2024-04-25 14/700 2024-04-26 15:17 by bnullh
[论文投稿] 研二光催化6月底四篇二区什么水平 5+5 wjtab 2024-04-22 14/700 2024-04-26 15:02 by 莱茵润色
[论文投稿] with editor 两个月了,什么原因? +6 yiersan9 2024-04-24 14/700 2024-04-26 14:03 by topedit
[教师之家] 博士论文被抄袭 +25 和尚敲小木鱼 2024-04-22 42/2100 2024-04-26 13:55 by ZHONGWU_U
[硕博家园] 考研,求职还是考编? +15 xizj 2024-04-21 24/1200 2024-04-26 11:49 by Kan客
[考研] 没学上 +4 季向阳 2024-04-26 8/400 2024-04-26 08:52 by 18975788067
[考博] 取博导收留 5+4 zzb777888 2024-04-20 10/500 2024-04-26 08:52 by polymerfriend
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +9 pom戴墨镜 2024-04-24 16/800 2024-04-26 08:37 by lfq_198989
[论文投稿] Nature一直在编辑手里,考虑好几天了,是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[教师之家] 期末给学生划重点都是什么话术啊 +16 luokereng 2024-04-20 18/900 2024-04-25 15:46 by BusyGer
[硕博家园] 聊天 +11 暮色恋伊人 2024-04-22 12/600 2024-04-25 13:53 by UCTS
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才!!! +4 黑魔变身啾 2024-04-20 12/600 2024-04-25 11:18 by dodonaomi
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然,现在基金委发来调查函 +33 bjdxyxy 2024-04-19 53/2650 2024-04-25 08:05 by 大雁高飞
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿,黑贝
[有机交流] 紧急求助,有谁用过三甲基硅醇钾TMSOK吗? 1000+3 genius2008 2024-04-19 10/500 2024-04-24 12:37 by genius2008
[教师之家] 大家访学都是怎么找的啊? +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[考博] 青岛科技大学顾玉超教授课题组招收博士生 +3 yanqingwang 2024-04-19 3/150 2024-04-22 15:21 by 理科生yes
[考研] 问题已经解决。 +13 lekinna 2024-04-19 30/1500 2024-04-22 13:09 by 矛dei到
[硕博家园] 材料工程,二本硕士,值得读吗 +3 简风23 2024-04-20 6/300 2024-04-21 12:19 by 简风23
信息提示
请填处理意见