24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cuphunter

金虫 (知名作家)

应助: 2 (幼儿园)
金币: 1195.5
散金: 1440
红花: 13
沙发: 26
帖子: 6312
在线: 612.1小时
虫号: 794813
注册: 2009-06-15
性别: GG
专业: 化学反应工程

[求助] matlab拟合动力学参数

本人为matlab菜鸟，希望有高手相助，给出matlab代码，本人万分感谢。这里给出的数据只是一个温度下的，我想先学习一下，看自己能不能弄别的温度，再拟合活化能。这些等之后再交流。

数据：
初始值： L1(0) L2(0) L3(0) L4 (0) L5(0) L6(0)
X0 4.29 32.77 15.50 21.41 5.80 20.23

实验值：
编号 L1 L2 L3 L4 L5 L6
1 6.95 36.54 13.13 17.56 3.86 21.96
2 6.47 34.98 11.95 17.09 4.07 25.45
3 9.45 40.82 11.87 11.80 3.03 23.03
4 8.62 50.21 10.01 7.49 1.71 21.96
5 10.33 46.08 9.11 9.63 2.55 22.3
6 9.13 40.62 11.38 12.01 3.29 23.56
7 7.26 37.41 13.64 16.65 3.46 21.57
8 7.61 37.20 13.22 15.21 3.46 23.3
9 8.07 37.86 12.28 14.67 3.58 23.54
10 7.32 36.66 14.39 14.58 3.72 23.54

动力学模型为：

回复此楼

珍爱生命，远离木虫.............................

1楼2011-12-22 15:50:29

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

dingd

铁杆木虫 (职业作家)

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

1stOpt求解很方便，只是你的数据拟合结果不是太好！；

CODE:

Variable x,a1,a2,a3,a4,a5,a6;

ODEFunction a1'=k31*a3;

            a2'=k42*a4;

            a3'=-(k31+k34+k35)*a3;

            a4'=-((k42+k45)*a4-k34*a3);

            a5'=-(k56*a5-k35*a3-k45*a4);

            a6'=k56*a5;

Data;

x        a1        a2        a3        a4        a5        a6

0        4.29        32.77        15.50        21.41        5.80        20.23

1        6.95        36.54        13.13        17.56        3.86        21.96

2        6.47        34.98        11.95        17.09        4.07        25.45

3        9.45        40.82        11.87        11.80        3.03        23.03

4        8.62        50.21        10.01        7.49        1.71        21.96

5        10.33        46.08        9.11        9.63        2.55        22.3

6        9.13        40.62        11.38        12.01        3.29        23.56

7        7.26        37.41        13.64        16.65        3.46        21.57

8        7.61        37.20        13.22        15.21        3.46        23.3

9        8.07        37.86        12.28        14.67        3.58        23.54

10        7.32        36.66        14.39        14.58        3.72        23.54

参数最佳估算
-------------------- -------------
k31 0.0419690845909743
k42 0.0662882491041956
k34 0.413188420544208
k35 -0.42242795564006
k45 0.375861403372349
k56 0.0940392220979913