Znn3bq.jpeg
²é¿´: 835  |  »Ø¸´: 1

zyj8119

ľ³æ (ÖøÃûдÊÖ)

[ÇóÖú] Ò»¸öMATLAB³ÌÐò£¬ÔËÐв»³öÀ´

CODE:
function IsothermTR
clear all;clc
global keq1 keq2 k1 k2 P
stoptime=1;dt=0.01;
T=573.9;keq1=0.312;keq2=0.480;
k1=1.496e7*exp(-8.44e3/T);k2=8.67e6*exp(-8.44e3/T);
P=1;F=128.2;
CA0=1/22.4;
CA0=CA0/453.6/(1e-3*3.2808^3);
v0=FA0/CA0;
x0=[0 0];t=[0,0.005,0.01,0.02,0.04,0.06,0.08,0.10,0.12,0.14,0.18,0.22,0.26,0.30,0.40];
[t,x]=ode45(@Equations,t,x0);
tau=CA0*t;SV=1./tau;
xBt=x(:,1)+x(:,2);
t42p=spline(xBt,t,0.42);
VR=t42p*FA0;
PB=1-x(:,1)-x(:,2);PD=x(:,1)/2-x(:,2);PT=x(:,2);
tmax=t(find(PD==max(PD)));
tau_max=tau(find(PD==max(PD)));
SV_max=SV(find(PD==max(PD)));
%Êä³ö½á¹û
disp('result:')
fprintf('\n ËùÐèµÄ·´Ó¦Æ÷Ìå')
fprintf('t x1 x2,xBt PB PD PT')
disp([t x(:,1:2) xBt PB PD PT]),plot(t,xBt),xlabel('V/F_A_0,ft^3.hr/lbmol'),ylabel('x_B_t')
figure,plot(t,PB,'b--',t,PD,'m-',t,PT,'r-'),xlabel('V/F_A_0,ft^3.hr/lbmol'),ylabel('P,atm')
legend('P_B','P_D','P_T')
%----------------------------------------------
function dxdt=Equations(t,x)
global keq1 keq2 k1 k2 P
PB=P*(1-x(1)-x(2));PD=P*(x(1)/2-x(2));PH=P*(x(1)/2+x(2));PT=P*x(2);
r1=k1*(PB*PB-PD*PH/keq1);r2=k2*(PB*PD-PH*PT/keq2);
%---------------------------------------------
dx1dt=r1;dx2dt=r2;
dxdt=[dx1dt;dx2dt];

»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ºÃºÃѧϰ£¬ÌìÌìÏòÉÏ¡£
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

change0618

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

·½ÕÉ´óʦ

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xiegangmai: ½ð±Ò+2, ллӦÖú£¡ 2012-05-27 10:07:46
P=1;F=128.2;
¸Ä³É
P=1;FA0=128.2;
2Â¥2012-05-27 08:30:03
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zyj8119 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁϹ¤³Ì085601£¬270Çóµ÷¼Á +35 @ASDF1234 2026-04-08 37/1850 2026-04-13 01:18 by ÐÒÃâ ..
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ´óѧ 0855 »úе 286 Çóµ÷¼Á +10 ²»»á³ÔÈâ 2026-04-12 10/500 2026-04-12 22:51 by ÁãÁã¶þ
[¿¼ÑÐ] ²ÄÁÏ085601µ÷¼Á +27 ºÎÈó²É123 2026-04-10 29/1450 2026-04-12 22:31 by fqwang
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +5 ×ÓľÄÅ 2026-04-12 5/250 2026-04-12 22:15 by paopaotu326
[¿¼ÑÐ] 086000µ÷¼Á +6 Ê®Æßsa 2026-04-07 6/300 2026-04-12 11:05 by ´óÁ¦Ë®ÊÖÁ¦´óÎÞÇ
[¿¼ÑÐ] 308Çóµ÷¼Á +5 VvvvL 2026-04-10 5/250 2026-04-12 10:17 by babysonlkd
[»ù½ðÉêÇë] ɽ¶«Ê¡»ù½ð2026 +5 jerry681 2026-04-08 6/300 2026-04-12 08:33 by kudofaye
[¿¼ÑÐ] 359Çóµ÷¼Á +5 θ¾·ÂÎÀÛÁË 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[¿¼ÑÐ] 085400 328·Ö Çóµ÷¼Á +10 ιÄãÒ»¸ö´ó³È×Ó 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[¿¼ÑÐ] 352 Çóµ÷¼Á +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by Ã÷Ô´ËʱÓÐ
[¿¼ÑÐ] 085404 293Çóµ÷¼Á +9 ÓÂÔ¶¿â°®314 2026-04-06 10/500 2026-04-11 10:36 by ×ÏêØ×ÏÆå
[¿¼ÑÐ] 087100³õÊÔ311Çóµ÷¼Á +4 ÈÎÑÅÇÙ 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[¿¼ÑÐ] 368Çóµ÷¼Á +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by µç×Ӱµ±ë
[¿¼ÑÐ] ±¾¿ÆÎ÷¹¤´ó 324Çóµ÷¼Á +4 wysyjs25 2026-04-10 4/200 2026-04-10 20:00 by À´¿´Á÷ÐÇÓê10
[¿¼ÑÐ] ²ÄÁÏרҵ344Çóµ÷¼Á +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[¿¼ÑÐ] 368»¯Ñ§Çóµ÷¼Á +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[¿¼ÑÐ] 353Çóµ÷¼Á +8 Çç¿ÕÍòÀïair 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[¿¼ÑÐ] 304Çóµ÷¼Á +10 ËØÄê¼ÀÓï 2026-04-06 17/850 2026-04-08 09:05 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +15 ÐܶþÏëÉϰ¶ 2026-04-06 15/750 2026-04-08 04:53 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 333Çóµ÷¼Á +6 ºÏ³ËÑîϰϦ 2026-04-06 6/300 2026-04-07 09:44 by Öí»á·É
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û