应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 有没有人会用1stopt或者matlab拟合这样数据，急求啊
52/1返回列表
查看: 1804  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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的值。
回复此楼

» 猜你喜欢
1楼 2014-03-31 21:51:54
已阅   回复此楼   关注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

木虫 (正式写手)
@dingd
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的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:59 by zNcEhTcH7ihq
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-15 15:35 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考研] 售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,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 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.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭