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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 有没有人会用1stopt或者matlab拟合这样数据，急求啊
52/1返回列表
查看: 1381  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭