24小时热门版块排行榜

返回列表

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
独孤神宇: 金币+1 2018-07-05 14:30:10
l1003785517: 金币+20, ★★★很有帮助 2018-07-05 15:28:26

CODE:

Title "kinetic values";

Constant R=8.314;

VarConstant T=[275,282,287,292];

ConstStr a1=k1*exp(-Ea1/R/T),a2=k2*exp(-Ea2/R/T);

Variables x(1:3),y(1:2);

Parameters a(3:5);

ParameterDomain = [0,];

ShareModel;

Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);

         y2=(a2*x2-a1*x1)/(1+a3*x1+a4*x2+a5*x2);

Data; //T=275

6.31317        0.029302 0.007644 0.0574 0.0345

6.29412        0.044590 0.011466 0.0407 0.0298

6.27507        0.059878 0.015288 0.0241 0.0251

6.26237        0.069433 0.018473 0.0278 0.0203

6.24713        0.080899 0.022295 0.0315 0.0208

6.23316        0.091091 0.026117 0.0268 0.0194

6.22109 0.099372 0.029939 0.0220 0.0160

6.21093        0.107016 0.032487 0.0172 0.0127

Data; //T=282

6.14807        0.133133        0.069433        0.2894        0.1089

6.062345        0.186004        0.102557        0.1519        0.0714

6.021705        0.202566        0.126763        0.0144        0.0338

5.991225        0.212758        0.147147        0.0425        -0.0038

5.968365        0.220402        0.162435        0.0705        0.0217

5.937885        0.234416        0.178997        0.0575        0.0297

5.91566        0.244608        0.1911        0.0445        0.0201

5.892165        0.254163        0.205114        0.0316        0.0105

Data; //T=287

6.0452        0.187278        0.118482        0.4615        0.1918

5.89026        0.277095        0.184093        0.273        0.1397

5.80136        0.326781        0.223587        0.0845        0.0876

5.74675        0.356083        0.249067        0.0829        0.0356

5.69976        0.377741        0.274547        0.0813        0.0243

5.67182        0.386022        0.294294        0.0677        0.0176

5.648325        0.390481        0.313404        0.0541        0.0152

5.614035        0.400673        0.33761        0.0404        0.0129

Data; //T=292

5.90677        0.220402        0.224224        0.6236        0.1761

5.741035        0.303212        0.307671        0.4216        0.1584

5.559425        0.383474        0.409591        0.2197        0.1408

5.434965        0.437619        0.480298        0.2059        0.1231

5.311775        0.497497        0.543998        0.1921        0.0677

5.26669        0.514696        0.572026        0.1219        0.0332

5.229225        0.527436        0.596869        0.0817        0.0195

5.201285        0.53508        0.617253        0.0585        0.0058

结果不好，自己好好检查下公式和数据：

Root of Mean Square Error (RMSE): 0.0829320200375007
Sum of Squared Residual: 0.440174076640027
Correlation Coef. (R): 0.0189691245411682
R-Square: 0.000359827685858348
Adjusted R-Square: 0.51886698386634
Determination Coef. (DC): -2.88389955816035
F-Statistic: 0.226197729277559

Parameter Best Estimate
-------------------- -------------
k1 4.60347494177944E35
ea1 166179.343552407
k2 9.73280643547253E23
ea2 93748.2902053327
a3 65737.9703181681
a4 4.29079677286443E-8
a5 49386056.7270752

赞一下(1人)

回复此楼

11楼2018-07-04 23:05:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

专家经验: +65
应助: 122 (高中生)
金币: 4841.3
红花: 43
帖子: 722
在线: 435.7小时
虫号: 3929201
注册: 2015-06-17
性别: GG
专业: 流体力学
管辖: 计算模拟

引用回帖:

6楼: Originally posted by l1003785517 at 2018-07-04 12:55:03
是的，作为一个求助者真的不应该出现这种错误，很抱歉。
我在代码中给出的公式，Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
y2=(a2*x2-a1*x1)/(1+a3*x1+a4 ...

你是四组数据对应四组a,k和ea吧。

赞一下

回复此楼

一团和气

12楼2018-07-05 11:21:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

l1003785517

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 11.4
散金: 30
帖子: 59
在线: 22.1小时
虫号: 4326948
注册: 2016-01-01
专业: 能源化工

引用回帖:

11楼: Originally posted by dingd at 2018-07-04 23:05:52
Title "kinetic values";
Constant R=8.314;
VarConstant T=;
ConstStr a1=k1*exp(-Ea1/R/T),a2=k2*exp(-Ea2/R/T);
Variables x(1:3),y(1:2);
Parameters a(3:5);
ParameterDomain = ;
...

非常感谢您的帮助，我想问一下：我把您给的代码在1.5版上运行的结果和您给的结果为什么不一样呢？而且每一次运行结果都不同，是因为软件版本太低的原因吗？
您说的结果不好应该是公式的问题，我自己推了几个可能的公式。您能分别帮我算一下这四个公式吗？谢谢！

（1）：Function y1=a1*x1/(1+a3*x1+a4*x2+a5*x3);
                        y2=a2*x1/(1+a3*x1+a4*x2+a5*x3);
（2）：Function y1=a1*a3*x1/(1+a3*x1+a4*x2+a5*x3)^4;
                        y2=a2*a3*x1/(1+a3*x1+a4*x2+a5*x3)^6;
（3）：Function y1=a1*a3*x1/(1+a3*x1+a4*x2+a5*x3)^7;
                        y2=a2*a3*x1/(1+a3*x1+a4*x2+a5*x3)^11;
（4）：Function y1=a1*a4*x2/(1+a3*x1+a4*x2+a5*x3);
                        y2=a2*a4*x3/(1+a3*x1+a4*x2+a5*x3);

赞一下

回复此楼

13楼2018-07-05 15:41:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖