24小时热门版块排行榜    

查看: 1114  |  回复: 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的回帖

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的回帖

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的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lgd224: 金币+20, ★★★★★最佳答案 2014-04-01 11:02:23
CODE:
Constant R=8.314;
VarConstant E=[0.01,0.1,1];
ParameterDomain = [0,];
Variable T,Y;
Function y=ArcSinh(power(E/exp(-Q/(R*T))/A,1/n))/a1;
Data;
1123.15        270.246
1223.15        151.858
1323.15        83.005
1423.15        45.71
1523.15        26.393
Data;
1123.15        359.754
1223.15        204.59
1323.15        107.459
1423.15        78.443
1523.15        63.934
Data;
1123.15        379.235
1223.15        284.426
1323.15        192.35
1423.15        112.022
1523.15        87.158

均方差(RMSE):13.018674791844
残差平方和(SSE):2542.28840003693
相关系数(R): 0.994905726348786
相关系数之平方(R^2): 0.989837404321606
确定系数(DC): 0.985581325706425
F统计(F-Statistic): 24.7915890860313

参数                  最佳估算
--------------------        -------------
q        473044.570077009
a        6.02336806702476E17
n        5.06819838070165
a1        0.00652970729120344

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

文件: 数据文件-1
No        目标 y        计算 y
1        270.246        265.53439728803
2        151.858        156.714944848118
3        83.005        87.8252550235811
4        45.71        50.177292211177
5        26.393        30.2351106945208

文件: 数据文件-2
No        目标 y        计算 y
1        359.754        332.270564685145
2        204.59        214.682587105772
3        107.459        129.896721026133
4        78.443        77.1510400555287
5        63.934        47.183363204705

文件: 数据文件-3
No        目标 y        计算 y
1        379.235        400.657013280714
2        284.426        278.76023492921
3        192.35        183.146119215513
4        112.022        115.470459190782
5        87.158        72.7345737366438
13楼2014-04-01 10:27:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgd224 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 75批博士后面上项目公布了 +4 ishow1989 2024-06-28 5/250 2024-06-30 00:58 by ishow1989
[基金申请] 2024国自然会评日期 +3 王兄弟 2024-06-29 4/200 2024-06-29 23:39 by yatou
[访问学者] 国家公派访问学者出结果了 +10 326lhpqk 2024-06-29 18/900 2024-06-29 21:55 by shi01587
[公派出国] 今天公派访问学者项目会出结果吗? +8 326lhpqk 2024-06-28 10/500 2024-06-29 20:09 by shi01587
[基金申请] 感觉地区基金太难了,申请了多年不中 +24 123为wqw 2024-06-27 31/1550 2024-06-29 14:14 by Siuhong Yiu
[基金申请] 统计一波计分系数 +5 lyfbangong 2024-06-28 6/300 2024-06-29 13:30 by lddscut
[基金申请] 能看出是否上会了吗 +11 articlefan 2024-06-23 17/850 2024-06-29 13:28 by articlefan
[找工作] offer选择 +13 羡渔樵 2024-06-26 13/650 2024-06-28 23:59 by 高敖曹
[考研] 刚当完兵回来想考研 +8 五官在线 2024-06-23 26/1300 2024-06-28 19:20 by 五官在线
[基金申请] 博后面上今天下午会公布吗?大家有无消息? +11 地理学1995 2024-06-24 16/800 2024-06-28 14:54 by 地理学1995
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +9 babu2015 2024-06-25 9/450 2024-06-28 09:11 by newfuzzy1
[论文投稿] 为何全世界都用JCR分区,咱们要弄个中科院分区?意义何在? +6 意得辑_editage 2024-06-25 7/350 2024-06-27 20:59 by keyaner23
[基金申请] 国基在研影响申请结果吗 +14 WOWO159357 2024-06-26 23/1150 2024-06-27 20:27 by wspglt
[基金申请] 基金评审时,评审专家还回去看申请人代表作之外的文章吗? +12 孤单12站 2024-06-26 18/900 2024-06-27 18:30 by 孤单12站
[论文投稿] Nano energy 投稿一个月with editor 到今天状态没有变过 5+3 甜甜甜。 2024-06-26 3/150 2024-06-27 18:00 by nono2009
[基金申请] 今年什么时候会评啊 +9 lancet0903 2024-06-24 9/450 2024-06-27 13:22 by yrzhu1984
[基金申请] 青年基金E02口青基去年几个函评专家? +6 他山攻玉之石 2024-06-25 9/450 2024-06-26 15:09 by 他山攻玉之石
[基金申请] 今天能不能出来名单 +8 地理学1995 2024-06-25 10/500 2024-06-26 09:46 by msjy
[有机交流] 求助析晶问题 20+4 dengdawang 2024-06-24 5/250 2024-06-24 21:22 by cc116
[有机交流] 生成亚胺的反应怎么能进行完全 +3 1369836 2024-06-23 3/150 2024-06-23 18:44 by hwqMSE
信息提示
请填处理意见