| 查看: 2390 | 回复: 15 | ||
lgd224木虫 (正式写手)
|
[求助]
求大神用matlab给我处理一组数据，有人没有啊？
|
|
温度 T ε=0.01 ε=0.1 ε =1 1123.15 277.88 374.69 1223.15 185.93 227.18 280.02 1323.15 103.56 143.35 196.46 1423.15 72.107 93.341 128.58 1523.15 38.036 59.045 86.149 ε=A[sinh（ασ）]^n*exp(-Q/RT) 其中σ为277.88 374.69 185.93 227.18 280.02 103.56 143.35 196.46 72.107 93.341 128.58 38.036 59.045 86.149 R=8.314，求大神求出A，n，α，Q的值。  |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有14人回复
-
为什么国自然不能直接公布
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有3人回复
-
系统查不到
已经有9人回复
-
科研孤儿太难了
已经有22人回复
-
怎么查啊
已经有4人回复
-
国合里面能看到了
已经有7人回复
-
国际合作可查了，中了面上
已经有18人回复
-
能否退出参与的面上项目解除限项
已经有25人回复
-
系统进不去
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何实现数据流聚类过程模拟的
已经有10人回复
- 关于matlab进行傅里叶变换和逆变换的程序，求大神指导！
已经有9人回复
- matlab解决重心法选址问题，运行不了啊，求助各位大神
已经有10人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
- Pd XPS数据处理，不会啊，求大神
已经有11人回复
- MATLAB作图求助：怎样在一个圆中用颜色渐变表示数据的分布情况
已经有4人回复
- Matlab或origin 数据处理 求零点处的值
已经有4人回复
- 【求助MATLAB大神】MATLAB画多行图时，每行的高度不太一样
已经有14人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- MATLAB关于lsqcurvefit程序几个错误，求大神指导一下
已经有6人回复
- matlab数据拟合问题求助
已经有7人回复
- 求大神教这个图怎么画
已经有16人回复
- matlab作业哪位大神麻烦给做一下,十万火急,周五之间要交
已经有3人回复
- MATLAB如何随机分组数据
已经有9人回复
- 如何用matlab处理txt文件里的数据
已经有7人回复
- 怎么处理matlab中的数据？
已经有7人回复
- 【求助】一组光谱数据的多元散射校正（MSC），求具体的MATLAB程序及用法
已经有6人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
chouchouNUDT
铜虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 918.6
- 散金: 117
- 红花: 3
- 沙发: 1
- 帖子: 593
- 在线: 115.5小时
- 虫号: 2275887
- 注册: 2013-02-06
- 专业: 信号理论与信号处理
2楼2013-11-20 22:51:40
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
★ ★
jjdg: 金币+2, 感谢参与 2013-11-21 00:12:51
jjdg: 金币+2, 感谢参与 2013-11-21 00:12:51
|
将σ视作因变量，T为自变量，ε用E表示，α用a1表示，原公式“ε=A[sinh（ασ）]^n*exp(-Q/RT)”可以变为： y=ArcSinh(power(E/exp(-Q/(R*T))/A,1/n))/a1 这样就变成一个多数据的拟合问题。用1stOpt求解，代码如下： 可得稳定唯一解： 均方差(RMSE):3.62215747674245 残差平方和(SSE):183.680347008497 相关系数(R): 0.999153030194409 相关系数之平方(R^2): 0.998306777746669 确定系数(DC): 0.998014412702569 F统计(F-Statistic): 625.625393341874 参数 最佳估算 -------------------- ------------- q 532019.492416723 a 8.7094212625832E18 n 5.59413166119775 a1 0.00817508910031064 ====== 结果输出 ====== 文件: 数据文件-1 No 目标 y 计算 y 1 277.88 277.612887666796 2 185.93 180.962712452548 3 103.56 110.089822992317 4 72.107 65.2588659639294 5 38.036 39.6542003250407 文件: 数据文件-2 No 目标 y 计算 y 1 227.18 227.782578218994 2 143.35 149.450943505574 3 93.341 93.772222805066 4 59.045 58.6306425200283 文件: 数据文件-3 No 目标 y 计算 y 1 374.69 377.253442161494 2 280.02 276.482215698101 3 196.46 193.930707772309 4 128.58 129.94332959974 5 86.149 84.9260713830243  E=0.01  E=0.1  E=1 |
3楼2013-11-20 23:08:33
lgd224
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 3210.1
- 散金: 6
- 红花: 7
- 帖子: 584
- 在线: 98.6小时
- 虫号: 2084021
- 注册: 2012-10-24
- 性别: GG
- 专业: 金属基复合材料
4楼2013-11-20 23:10:06
lgd224
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 3210.1
- 散金: 6
- 红花: 7
- 帖子: 584
- 在线: 98.6小时
- 虫号: 2084021
- 注册: 2012-10-24
- 性别: GG
- 专业: 金属基复合材料
5楼2013-11-20 23:14:35
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
jjdg: 金币+2, 感谢参与 2013-11-21 00:13:02
感谢参与，应助指数 +1
jjdg: 金币+2, 感谢参与 2013-11-21 00:13:02
|
你·公式中“exp(-Q/RT)”是否应该是“exp(-Q/(R*T))”?，若是，就是上面的结果，若不是，参考下面结果： 均方差(RMSE):6.05927259398182 残差平方和(SSE):514.006981154508 相关系数(R): 0.998695771969003 相关系数之平方(R^2): 0.997393244948764 确定系数(DC): 0.994741899984704 F统计(F-Statistic): 130.402251286684 参数 最佳估算 -------------------- ------------- q -0.33474848568755 a 1.76661899323681E-24 n 7.42615822684346 a1 0.00526155750103373 ====== 结果输出 ====== 文件: 数据文件-1 No 目标 y 计算 y 1 277.88 267.960601187546 2 185.93 183.131483206876 3 103.56 116.343445052685 4 72.107 70.339425391132 5 38.036 41.5101766354818 文件: 数据文件-2 No 目标 y 计算 y 1 227.18 229.879170433833 2 143.35 152.000645570673 3 93.341 94.2077240999639 4 59.045 56.2269281097231 文件: 数据文件-3 No 目标 y 计算 y 1 374.69 377.740955119525 2 280.02 281.213018627753 3 196.46 194.368631063173 4 128.58 124.648474340735 5 86.149 75.7665218221402 |
6楼2013-11-20 23:30:07
lgd224
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 3210.1
- 散金: 6
- 红花: 7
- 帖子: 584
- 在线: 98.6小时
- 虫号: 2084021
- 注册: 2012-10-24
- 性别: GG
- 专业: 金属基复合材料
7楼2013-11-20 23:42:13
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lgd224: 金币+20, ★★★★★最佳答案, 大神谢谢你，但是需要接近，别正好450000啊，能不能给我计算一下Q=472180时候其他参数的值啊！ 2013-11-21 13:53:49
lgd224: 金币+20, ★★★★★最佳答案, 大神谢谢你，但是需要接近，别正好450000啊，能不能给我计算一下Q=472180时候其他参数的值啊！ 2013-11-21 13:53:49
|
基于你的模型和数据，从数学计算来说，q=532019.492416723就是最优值；如果非要强调q=450000，那可以令q为常数450000，此时就只有另外三个参数了： 均方差(RMSE):6.78769969975805 残差平方和(SSE):645.020140997337 相关系数(R): 0.999115394404998 相关系数之平方(R^2): 0.998231571337055 确定系数(DC): 0.99395943457163 F统计(F-Statistic): 323.554807613202 参数 最佳估算 -------------------- ------------- a 5.65197707433542E15 n 4.72728071332255 a1 0.00847329748896421 ====== 结果输出 ====== 文件: 数据文件-1 No 目标 y 计算 y 1 277.88 265.736461617687 2 185.93 172.574542392554 3 103.56 104.559735056992 4 72.107 61.8048866739181 5 38.036 37.4988826530263 文件: 数据文件-2 No 目标 y 计算 y 1 227.18 226.131220942749 2 143.35 149.768333067832 3 93.341 94.74327628421 4 59.045 59.5098630637393 文件: 数据文件-3 No 目标 y 计算 y 1 374.69 379.581710788027 2 280.02 282.026824187704 3 196.46 201.485452678126 4 128.58 137.793949165241 5 86.149 91.5507811431674 |
8楼2013-11-21 09:38:51
lgd224
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 3210.1
- 散金: 6
- 红花: 7
- 帖子: 584
- 在线: 98.6小时
- 虫号: 2084021
- 注册: 2012-10-24
- 性别: GG
- 专业: 金属基复合材料
9楼2013-11-21 13:55:30
lgd224
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 3210.1
- 散金: 6
- 红花: 7
- 帖子: 584
- 在线: 98.6小时
- 虫号: 2084021
- 注册: 2012-10-24
- 性别: GG
- 专业: 金属基复合材料
10楼2013-11-21 14:05:15
