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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求大神用matlab给我处理一组数据，有人没有啊？
161/2返回列表12下一页
查看: 1630  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者lgd224将赠送您 20 个金币

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-11-20 22:23:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chouchouNUDT

铜虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-11-21 00:12:37
这个恐怕不能直接计算，只能通过大数据来逼近一个近似值

[ 发自小木虫客户端 ]
赞一下 回复此楼
It will be better tomorrow.
2楼2013-11-20 22:51:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★ ★
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求解，代码如下：
CODE:
Constant R=8.314;
VarConstant E=[0.01,0.1,1];
Variable T,y;
Function y=ArcSinh(power(E/exp(-Q/(R*T))/A,1/n))/a1;
Data;
1123.15        277.88
1223.15        185.93
1323.15        103.56
1423.15        72.107
1523.15        38.036
Data;
1223.15        227.18
1323.15        143.35
1423.15        93.341
1523.15        59.045
Data;
1123.15        374.69
1223.15        280.02
1323.15        196.46
1423.15        128.58
1523.15        86.149

可得稳定唯一解：

均方差(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
求大神用matlab给我处理一组数据，有人没有啊？
E=0.01


求大神用matlab给我处理一组数据，有人没有啊？-1
E=0.1


求大神用matlab给我处理一组数据，有人没有啊？-2
E=1
赞一下 回复此楼
3楼2013-11-20 23:08:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)
那也行啊！帮帮弄弄啊！
赞一下 回复此楼
4楼2013-11-20 23:10:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)
引用回帖:
3楼: Originally posted by dingd at 2013-11-20 23:08:33
将σ视作因变量，T为自变量，ε用E表示，α用a1表示，原公式“ε=A^n*exp(-Q/RT)”可以变为：

y=ArcSinh(power(E/exp(-Q/(R*T))/A,1/n))/a1

这样就变成一个多数据的拟合问题。用1stOpt求解，代码如下：

Co ...

谢谢大神啊，不过这个Q求出来之后有点大啊，正常的事450000左右啊，怎么调整一下，可以求出q比较接近这个数啊？
赞一下 回复此楼
5楼2013-11-20 23:14:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)
引用回帖:
6楼: Originally posted by dingd at 2013-11-20 23:30:07
你·公式中“exp(-Q/RT)”是否应该是“exp(-Q/(R*T))”?，若是，就是上面的结果，若不是，参考下面结果：

均方差(RMSE):6.05927259398182
残差平方和(SSE):514.006981154508
相关系数(R): 0.998695771969003
...

是的，就应该是exp(-Q/(R*T))，大神能不能给我修改一下 然后让Q 接近450000啊！
赞一下 回复此楼
7楼2013-11-20 23:42:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)
引用回帖:
8楼: Originally posted by dingd at 2013-11-21 09:38:51
基于你的模型和数据，从数学计算来说，q=532019.492416723就是最优值；如果非要强调q=450000，那可以令q为常数450000，此时就只有另外三个参数了：

均方差(RMSE):6.78769969975805
残差平方和(SSE):645.0201409 ...

最好能把代码也给我发一下，我又这个1stopt的软件，但是只会最简单拟合，求代码以及说明啊
赞一下 回复此楼
9楼2013-11-21 13:55:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgd224

木虫 (正式写手)
引用回帖:
8楼: Originally posted by dingd at 2013-11-21 09:38:51
基于你的模型和数据，从数学计算来说，q=532019.492416723就是最优值；如果非要强调q=450000，那可以令q为常数450000，此时就只有另外三个参数了：

均方差(RMSE):6.78769969975805
残差平方和(SSE):645.0201409 ...

大神有没有别的联系方式，教教我怎么用1stopt，或者给我教程啊，小弟不胜感激啊！
赞一下 回复此楼
10楼2013-11-21 14:05:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgd224 的主题更新
161/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 每到中夜，情难自抑 +32 sioc-sunj 2024-05-28 56/2800 2024-06-01 16:55 by 佩佩姐618
[基金申请] 数理的人才答辩通知发了吗 +7 yzy3327 2024-05-30 7/350 2024-06-01 09:15 by 撒哈拉西部
[教师之家] 双非本科毕业论文，气人 +6 河西夜郎 2024-05-27 7/350 2024-06-01 08:50 by Quakerbird
[考博] 求生物学的博导收留 +4 申博成功功成博� 2024-06-01 4/200 2024-06-01 08:01 by 起名字哎
[教师之家] 选择 +6 459582015 2024-05-31 6/300 2024-06-01 07:32 by otani
[硕博家园] 哈工大硕博招生!博士每月入学! +4 nailooo 2024-05-30 5/250 2024-06-01 06:47 by anevay
[考研] 研0二导师分到新来的博士后靠谱吗 +6 sone9 2024-05-31 6/300 2024-06-01 06:35 by sjianglei
[基金申请] 讨论贴：好本子一定会被打招呼的淘汰吗 +36 ioe3937 2024-05-27 46/2300 2024-05-31 22:06 by wanghuawei
[考博] 广东以理材料系碳点与功能材料课题组 — 2博士名额 / 科研助理 +4 小城夜很美 2024-05-27 11/550 2024-05-31 21:26 by 小城夜很美
[找工作] 找工作如此之难 +7 探123 2024-05-25 7/350 2024-05-31 15:58 by 平南小刀
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠 +10 瞬息宇宙 2024-05-27 11/550 2024-05-31 15:38 by 烟 火
[高分子] MMA预聚体光固化发雾问题求助 +3 惠亚金总 2024-05-29 10/500 2024-05-31 14:59 by 惠亚金总
[电化学] 有锂金属软包电池代工公司代做软包嘛 +5 yyywj995 2024-05-26 6/300 2024-05-30 08:33 by 那年的旧城
[论文投稿] 审稿专家比较坚定的让补充实验，但实在没法补充实验，修回还有希望吗？ (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
[有机交流] 液相纯度高，但产品析不太出来 10+4 cui19236 2024-05-27 9/450 2024-05-30 07:45 by yuanjijoy
[论文投稿] 《控制理论与应用》期刊收版面费吗？ +4 ygj2015 2024-05-29 4/200 2024-05-30 06:59 by nono2009
[论文投稿] 有没有老师需要发表论文 +4 金老师论文助理- 2024-05-29 4/200 2024-05-29 16:51 by liuyupu132
[基金申请] 信息学部函评结束了吗？ +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[论文投稿] 真急着毕业，CPB主编终审17天了，邮件催稿了两次，就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[基金申请] 面上基金会评专家，有回避机制吗？ +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
信息提示
关闭
请填处理意见
关闭