24小时热门版块排行榜

返回列表

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

最自私的女人

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 136.5
红花: 1
帖子: 24
在线: 7小时
虫号: 2922685
注册: 2014-01-09
专业: 化学反应工程

[求助] 动力学方程组参数求解问题已有3人参与

本人动力学方程如下：
dx1/dt=k21x2+k31x3+k41x4+k51x5+k61x6+k71x7+k81x8
dx2/dt=-k21x2+k32x3+k42x4+k52x5+k62x6+k72x7+k82x8
dx3/dt=-k31x3-k32x3+k43x4+ k53x5+k63x6+k73x7+k83x8
dx4/dt=-k41x4-k42x4-k43x4+ k54x5+k64x6+k74x7+k84x8
dx5/dt=-k51x5-k52x5-k53x5-k54x5-k56x52-k57x52+k65x6+k75x7+k85x8
dx6/dt=-k61x6-k62x6-k63x6-k64x6-k65x6-k67x62+k56x52+k76x7+k86x8
dx7/dt=-k71x7-k72x7-k73x7-k74x7-k75x7-k76x7-k78x72+k57x52+k67x62+k87x8
dx8/dt=-k81x8-k82x8-k83x8-k84x8-k85x8-k86x8-k87x8-k89x82+k78x72
dx9/dt=k89x82
方程中有多个缩合反应为表观二级（即平方），已上传图片供参考
k为动力学参数，k大于0，共33个
x为质量分数，x＞0，九个x之和为1
已知数据如下：
温度时间 x1 x2 x3 x4 x5 x6 x7 x8 x9
430 0.00 0.0192 0.0236 0.1006 0.3284 0.1647 0.1327 0.1981 0.0319 0.0008
430 20.00 0.0287 0.0503 0.1313 0.3156 0.1391 0.1206 0.1755 0.0366 0.0022
430 40.00 0.0349 0.0858 0.1827 0.3058 0.1203 0.0805 0.1476 0.0384 0.0039
430 60.00 0.0423 0.1250 0.2086 0.2950 0.1005 0.0624 0.1180 0.0428 0.0053
本人想编写一个MATLAB程序实现以下功能：
设k初值均为0，求解方程中的33个参数，将求得的参数代入方程中，根据t=0min的数据求得t=20、40、60min的x计算值，将x计算值与x实验值比较，如果误差＜5%，得到的参数就是最终参数，如果误差大于5%，就将得到的参数最为初值循环上述计算，直到误差＜5%。
哪位大侠能帮忙编写一下，小女子不胜感激呢~

动力学方程组.jpg

回复此楼