24小时热门版块排行榜    

查看: 2393  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:45 by o3EMxHUAQDtS
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:40 by o3EMxHUAQDtS
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +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
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 基金不中,共勉 +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
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 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
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见