| 查看: 563 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】matlab程序
|
|||
|
% 动力学数据: % yCO yCO2 yH2 yH2O yCH3OH yDME % P W F0 y(1) y(2) y(3) y(4) y(5) y(6) % 1 2 3 4 5 6 7 8 9 % T=523.15K yexp= ... [ 4 2 5.5 0.3182 0.1663 0.3858 2.96E-04 9.83E-05 0.1293 4 2 8.3 0.3746 0.1084 0.4226 1.39E-03 3.51E-04 0.0927 4 2 11 0.4137 0.0766 0.4482 6.72E-04 6.70E-05 0.0607 5 2 6.9 0.2212 0.2619 0.3087 4.85E-04 1.27E-04 0.2076 5 2 10.3 0.3128 0.1648 0.3834 8.15E-04 5.97E-05 0.1381 5 2 13.8 0.3876 0.1065 0.4244 6.85E-04 8.56E-05 0.0807 ] 动力学数据如上贴 又知 P=[4;4;4;5;5;5]; F0 =[5.5; 8.3; 11; 6.9; 10.3; 13.8]; r2=1.503e-3.*(P5.^2./P4-P6/1.0e-013*0.4997); r1=8.101e-5.*P3.^(-1/2).*(P1.*P3.^2-P5./1.0e+061*7.6535)/(1+1.0e-006*0.476.*P2+(676.7713.*P3).^(1/2)).^2; r3=5.015e-3.*(85.3437.*P4-P2.*P3./P1)./(1+1.0e-006*0.476.*P2)+676.7713.^(1/2).*P1.^(1/2)); A=1+4*y(6)+2*y(5); dy1dW = A./F0*(r3+2*y(2)*r1); dy2dW = r3; dy3dW = r3-2*r1; dy4dW = r2; dy5dW = A./F0*(2*r1*y(6)+r2); dy6dW = A./F0*(r1+2*r1*y(5)-2*r2); dydW = [dy1dW; dy2dW; dy3dW; dy4dW; dy5dW; dy6dW]; P1=P.*y1; P2=P.*y2;P3=P.*y3;P4=P.*y4;P5=P.*y5;P6=P.*y6; y1 y2 y3 y4 y5 y6 分别是动力学数据的对应列 [ Last edited by myp5188 on 2009-6-3 at 21:45 ] |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有21人回复
-
导师想让我从独立一作变成了共一第一
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
申请2026年博士
已经有6人回复
8楼2009-06-14 11:33:31
2楼2009-06-04 15:36:13
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2009-06-05 21:45:50
shinevip
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2487.3
- 散金: 1292
- 红花: 4
- 帖子: 489
- 在线: 465.5小时
- 虫号: 242238
- 注册: 2006-04-11
- 性别: GG
- 专业: 分子生物学
4楼2009-06-06 13:24:14