24小时热门版块排行榜    

查看: 1461  |  回复: 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

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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的回帖
查看全部 13 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] (出分啦!!!)26 考研调剂 20w+调剂信息火速更新中!【捡漏】qq 群:1087823172 +3 db58 2026-02-28 79/3950 2026-02-28 17:50 by nbf68
[考研] 085600材料工程一志愿中科大总分312求调剂 +4 吃宵夜1 2026-02-28 6/300 2026-02-28 17:38 by 吃宵夜1
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[考研] 高分子化学与物理调剂 +3 好好好1233 2026-02-28 4/200 2026-02-28 17:13 by sandychj
[考研] 化工专硕348,一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856材料求调剂 +8 hyf hyf hyf 2026-02-28 9/450 2026-02-28 16:27 by etapple
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 290求调剂 +4 材料专硕调剂; 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见