24小时热门版块排行榜    

查看: 1793  |  回复: 12

lgd224

木虫 (正式写手)

[求助] 有没有人会用1stopt或者matlab拟合这样数据,急求啊 已有1人参与

温度        

     T          ε=0.01      ε=0.1     ε =1
1123.15        270.246        359.754        379.235
1223.15        151.858        204.59        284.426
1323.15        83.005        107.459        192.35
1423.15        45.71        78.443        112.022
1523.15        26.393        63.934        87.158


ε=A[sinh(ασ)]^n*exp(-Q/RT)
(σ为上表中3组数据)
R=8.314,求大神求出A,n,α,Q的值。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

lgd224

木虫 (正式写手)

2楼2014-03-31 21:58:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)

这是大神你上次用的代码,就是数据变了一下。@dingd
将σ视作因变量,T为自变量,ε用E表示,α用a1表示,原公式“ε=A[sinh(ασ)]^n*exp(-Q/RT)”可以变为:

y=ArcSinh(power(E/exp(-Q/(R*T))/A,1/n))/a1

这样就变成一个多数据的拟合问题。用1stOpt求解
3楼2014-03-31 22:00:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
3楼: Originally posted by lgd224 at 2014-03-31 22:00:12
这是大神你上次用的代码,就是数据变了一下。dingd
将σ视作因变量,T为自变量,ε用E表示,α用a1表示,原公式“ε=A^n*exp(-Q/RT)”可以变为:

y=ArcSinh(power(E/exp(-Q/(R*T))/A,1/n))/a1

这样就变成一个 ...

自己先用1stOpt试试,有问题再问。
4楼2014-03-31 22:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)

内容已删除
5楼2014-03-31 22:35:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)

引用回帖:
4楼: Originally posted by dingd at 2014-03-31 22:09:48
自己先用1stOpt试试,有问题再问。...

我把上次你弄的那个代码完整复制到代码本,然后执行,出来结果不对,Q才5000多,不知道是不是哪儿不对,是不是操作不对?
6楼2014-03-31 22:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
第一组:
均方差(RMSE):0.55317791225107
残差平方和(SSE):1.53002901301226
相关系数(R): 0.999980853458417
相关系数之平方(R^2): 0.999961707283424
确定系数(DC): 0.999960903658602
卡方系数(Chi-Square): 0.0142617105566451
F统计(F-Statistic): 8525.94696622839

参数        最佳估算
----------        -------------
q        -12384.7168982259
a        0.0029745312581722
n        -0.12908019461298
a1        0.00597845424309054

老版本似乎也能算的,你自己试试。
7楼2014-03-31 23:03:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)

引用回帖:
7楼: Originally posted by dingd at 2014-03-31 23:03:16
第一组:
均方差(RMSE):0.55317791225107
残差平方和(SSE):1.53002901301226
相关系数(R): 0.999980853458417
相关系数之平方(R^2): 0.999961707283424
确定系数(DC): 0.999960903658602
卡方系数(Chi-Square ...

我不会设置算法里面的那些,这个结果应该不对吧,q不可能是负值啊,q应该算的是443176.053 左右吧 麻烦您再帮我看看
8楼2014-03-31 23:06:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

参数有范围,你应该事先就说出来啊,他人哪里知道,另外你的公式也有问题,有过拟合现象,会导致参数值结果不唯一。
9楼2014-03-31 23:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)

引用回帖:
9楼: Originally posted by dingd at 2014-03-31 23:23:47
参数有范围,你应该事先就说出来啊,他人哪里知道,另外你的公式也有问题,有过拟合现象,会导致参数值结果不唯一。

对不起,这个公式还是上次你给我计算用的公式啊!我自己是用笨办法算的,但是老师要求要有误差系数,那就设定α为0.098,Q=443176.053.麻烦你了
10楼2014-03-31 23:29:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgd224 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
请填处理意见