| 查看: 803 | 回复: 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 ] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有11人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wangen994(金币+3,VIP+0):谢谢交流，欢迎常来 6-12 10:13
小木虫(金币+0.5):给个红包，谢谢回帖交流
wangen994(金币+3,VIP+0):谢谢交流，欢迎常来 6-12 10:13
|
P=yexp(:,1); W=yexp(:,2); F0=yexp(:,3); y(1)=yexp(:,4); y(2)=yexp(:,5); y(3)=yexp(:,6); y(4)=yexp(:,7); y(5)=yexp(:,8); y(6)=yexp(:,9); P1=P.*y(1); P2=P.*y(2);P3=P.*y(3);P4=P.*y(4);P5=P.*y(5);P6=P.*y(6); 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]; |
7楼2009-06-12 09:15:14
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.7小时
- 虫号: 242238
- 注册: 2006-04-11
- 性别: GG
- 专业: 分子生物学
4楼2009-06-06 13:24:14
