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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 有没有人会用1stopt或者matlab拟合这样数据，急求啊
51/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

铁杆木虫 (职业作家)
引用回帖:
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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭