24小时热门版块排行榜    

查看: 2397  |  回复: 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的值。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
查看全部 16 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见