24СʱÈÈÃŰæ¿éÅÅÐаñ    

±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1276  |  »Ø¸´: 3

liqiang0915

ľ³æ (³õÈëÎÄ̳)

[ÇóÖú] ÄÄλ´óÉñ°ïæ¿´Ò»ÏÂMATLABÇó½â΢·Ö·½³ÌµÄÎÊÌ⣿ ÒÑÓÐ1È˲ÎÓë

clear all
clc
global  L Ct Ac rhoB Cp H1 H2 U dt TJ U TJ  us E1 E2 R us
L = 4;                  % ·´Ó¦¹Ü³¤, m
rhoB = 770;             % ´ß»¯¼Á¶Ñ»ýÃܶÈ, kg/m3
rhog =691.82            % ÒºÌå»ìºÏÎïÃܶÈ, kg/m3
CA0=4.136              % 2M1BĦ¶ûÁ÷ÂÊ, moles/hr
CB0=37.22             % 2M2BĦ¶ûÁ÷ÂÊ, moles/hr
CC0 = 41.36;                  % MeOHĦ¶ûÁ÷ÂÊ, moles/hr
Cp = 131.33;                  % ±ÈÈÈ, kJ/kmol K
U = 0.105;                         % ´«ÈÈϵÊý, kJ/m2 hr K
dt = 0.025;                    % ¹Ü¾¶, m  
TJ = 333.15;                      % ÀäȴζÈ, K
T0 = 343.15;                 % ÎïÁϽø¿ÚζÈ(³õʼζÈ), K
H1 = 33600;                  % ·´Ó¦2M1B¡úTµÄ·´Ó¦ÈÈ, kJ/kmol
H2 = 26800;             % ·´Ó¦2M2B¡úTµÄ·´Ó¦ÈÈ, kJ/kmol
us = 3.11               % ¿ÕËÙ£¬h-1

% »î»¯ÄÜ, kJ/kmol
E1 =89830  ;           
E2 =120010 ;                                      

R = 8.314;                      % ÀíÏëÆøÌå³£Êý, kJ/kmol K

Ac = pi*(dt/2)^2;       % ·´Ó¦¹ÜµÄºá½ØÃæ»ý, m2


[z, y] = ode45(@Equations, [0 L], [CA0 CB0 CC0 T0])
CA = y(:, 1);
CB = y(:, 2);
CC = y(:, 3);
xA = (CA0-CA)./CA0;                        % AµÄת»¯ÂÊ
xB = (CB0-CB)./CB0;                    % BµÄת»¯ÂÊ
xB = [0; xB]
xC = (CC0-CC)./CC0;                    % CµÄת»¯ÂÊ
xC = [0; xC]

% ͼÐÎÊä³ö
plot(z, y(:, 4))        % ζȷֲ¼
xlabel('z')
ylabel('T (K)')
figure
plot(z, xA, 'r-')       % ת»¯ÂÊ·Ö²¼
xlabel('z')
ylabel('x_A')
figure
plot(z, CA, 'r-', z, CB, 'k--', z, CC, 'b:')    % Ũ¶È·Ö²¼
xlabel('z')
ylabel('C_A, C_B, C_C')
legend('C_A', 'C_B', 'C_C')

function [rA, rB, rC, k1, k2,K1,K2] = Rates(y(1), y(2), y(3), T)    % ·´Ó¦¶¯Á¦Ñ§
global E1 E2 R CC0

% Ëٶȳ£Êý, kmol/kg catalyst hr
k1 = exp(-E1/(R*T) + 35.35);
k2 = exp(-E2/(R*T) + 42.47);
% ·´Ó¦Æ½ºâ³£Êý
K1= exp(-8.3881+4041.2/T);
K2= exp(-8.2473+3225.3/T);

% ·´Ó¦ËÙ¶È, kmol/kg catalyst hr
     
rA = -k1*CA*CC + k1/K1*(CC0-CC)         % 2M1BµÄÏûºÄËÙÂÊ
rB = -k2*CB*CC + k2/K2*(CC0-CC);        % 2M2BµÄÏûºÄËÙÂÊ
rC = -k1*CA*CC + k1/K1*(CC0-CC)-k2*CB*CC + k2/K2*(CC0-CC);   % ¼×´¼µÄ×Ü·´Ó¦ËÙ¶È

function dydz = Equations(z, y)                         % Ä£ÐÍ·½³Ì×é
global   Ct Ac rhoB Cp H1 H2 U dt TJ U
CA = y(1);                                    % A---2M1B
CB = y(2);                                    % B---2M2B
CC = y(3);                                    % C---MeOH
T = y(4);

% ·´Ó¦ËÙ¶È
[rA, rB, rC, k1, k2,K1,K2] = Rates(y(1), y(2), y(3), T)

% ÎïÁÏÆ½ºâ
dCAdz = rhoB*rA/us;
dCBdz = rhoB*rB/us;
dCCdz = rhoB*rC/us;

% ÈÈÁ¿ºâËã
dTdz = ( rhoB*((-H1*k1*y(1)*y(3))-(-H2*k2*y(2)*y(3)))-4*U*(T-TJ)/dt )/(us*rhog*Cp);
dydz = [dCAdz; dCBdz; dCCdz; dTdz];
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

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

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
ͼÐÎÈ總ͼ1¡£
CODE:
function test33
clear all;clc
L = 4;                  % ·´Ó¦¹Ü³¤, m
CA0=4.136;              % 2M1BĦ¶ûÁ÷ÂÊ, moles/hr
CB0=37.22  ;           % 2M2BĦ¶ûÁ÷ÂÊ, moles/hr
CC0 = 41.36;                  % MeOHĦ¶ûÁ÷ÂÊ, moles/hr
Cp = 131.33;                  % ±ÈÈÈ, kJ/kmol K
U = 0.105;                         % ´«ÈÈϵÊý, kJ/m2 hr K
dt = 0.025;                    % ¹Ü¾¶, m  
TJ = 333.15;                      % ÀäȴζÈ, K
T0 = 343.15;  
[z, y] = ode45(@Equations,[0 L],[CA0 CB0 CC0 T0]);
[z y(:,1) y(:,2) y(:,3) ]
CA = y(:, 1);
CB = y(:, 2);
CC = y(:, 3);
xA = (CA0-CA)./CA0;                        % AµÄת»¯ÂÊ
xB = (CB0-CB)./CB0;                    % BµÄת»¯ÂÊ
xB = [0; xB];
xC = (CC0-CC)./CC0;                    % CµÄת»¯ÂÊ
xC = [0; xC];

% ͼÐÎÊä³ö
figure(1)
plot(z, y(:, 4))        % ζȷֲ¼
xlabel('z')
ylabel('T (K)')
figure(2)
plot(z, xA, 'r-')       % ת»¯ÂÊ·Ö²¼
xlabel('z')
ylabel('x_A')
figure(3)
plot(z, CA, 'ro-', z, CB, 'm>-', z, CC, 'b*-')    % Ũ¶È·Ö²¼
xlabel('z')
ylabel('C_A, C_B, C_C')
legend('C_A', 'C_B', 'C_C')



function dydz = Equations(z, y)

rhoB = 770;             % ´ß»¯¼Á¶Ñ»ýÃܶÈ, kg/m3
rhog =691.82   ;         % ÒºÌå»ìºÏÎïÃܶÈ, kg/m3
CA0=4.136;              % 2M1BĦ¶ûÁ÷ÂÊ, moles/hr
CB0=37.22  ;           % 2M2BĦ¶ûÁ÷ÂÊ, moles/hr
CC0 = 41.36;                  % MeOHĦ¶ûÁ÷ÂÊ, moles/hr
Cp = 131.33;                  % ±ÈÈÈ, kJ/kmol K
U = 0.105;                         % ´«ÈÈϵÊý, kJ/m2 hr K
dt = 0.025;                    % ¹Ü¾¶, m  
TJ = 333.15;                      % ÀäȴζÈ, K
T0 = 343.15;                 % ÎïÁϽø¿ÚζÈ(³õʼζÈ), K
H1 = 33600;                  % ·´Ó¦2M1B¡úTµÄ·´Ó¦ÈÈ, kJ/kmol
H2 = 26800;             % ·´Ó¦2M2B¡úTµÄ·´Ó¦ÈÈ, kJ/kmol
us = 3.11   ;            % ¿ÕËÙ£¬h-1

% »î»¯ÄÜ, kJ/kmol
E1 =89830  ;           
E2 =120010 ;                                       

R = 8.314;                      % ÀíÏëÆøÌå³£Êý, kJ/kmol K

Ac = pi*(dt/2)^2;       % ·´Ó¦¹ÜµÄºá½ØÃæ»ý, m2



CA = y(1);                                    % A---2M1B
CB = y(2);                                    % B---2M2B
CC = y(3);                                    % C---MeOH
T = y(4);
% Ëٶȳ£Êý, kmol/kg catalyst hr
k1 = exp(-E1/(R*T) + 35.35);
k2 = exp(-E2/(R*T) + 42.47);
% ·´Ó¦Æ½ºâ³£Êý
K1= exp(-8.3881+4041.2/T);
K2= exp(-8.2473+3225.3/T);
% ·´Ó¦ËÙ¶È, kmol/kg catalyst hr   
rA = -k1*CA*CC + k1/K1*(CC0-CC)  ;       % 2M1BµÄÏûºÄËÙÂÊ
rB = -k2*CB*CC + k2/K2*(CC0-CC);        % 2M2BµÄÏûºÄËÙÂÊ
rC = -k1*CA*CC + k1/K1*(CC0-CC)-k2*CB*CC + k2/K2*(CC0-CC);   % ¼×´¼µÄ×Ü·´Ó¦ËÙ¶È

                         % Ä£ÐÍ·½³Ì×é


% ·´Ó¦ËÙ¶È


% ÎïÁÏÆ½ºâ
dCAdz = rhoB*rA/us;
dCBdz = rhoB*rB/us;
dCCdz = rhoB*rC/us;

% ÈÈÁ¿ºâËã
dTdz = ( rhoB*((-H1*k1*y(1)*y(3))-(-H2*k2*y(2)*y(3)))-4*U*(T-TJ)/dt )/(us*rhog*Cp);
dydz = [dCAdz; dCBdz; dCCdz; dTdz];

ÄÄλ´óÉñ°ïæ¿´Ò»ÏÂMATLABÇó½â΢·Ö·½³ÌµÄÎÊÌ⣿
¸½Í¼1.jpg

MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
2Â¥2014-02-17 21:02:13
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

liqiang0915

ľ³æ (³õÈëÎÄ̳)

ÒýÓûØÌû:
2Â¥: Originally posted by ÔÂÖ»À¶ at 2014-02-17 21:02:13
ͼÐÎÈ總ͼ1¡£


function test33
clear all;clc
L = 4;                  % ·´Ó¦¹Ü³¤, m
CA0=4.136;              % 2M1BĦ¶ûÁ÷ÂÊ, moles/hr
CB0=37.22  ;           % 2M2BĦ¶ûÁ÷ÂÊ, moles/hr
CC0 = 41.3 ...

ËãµÄʱºò»¹ÊdzöÏÖ´íÎ󣬷³ÀÍÄúÔÙ¿´¿´¡£
??? Error using ==> feval
Undefined function or method 'Equations' for input arguments of type 'double'.

Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...
3Â¥2014-02-18 10:30:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÔÂÖ»À¶

Ö÷¹ÜÇø³¤ (Ö°Òµ×÷¼Ò)

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
liqiang0915: ½ð±Ò+50, ¡ï¡ï¡ï¡ï¡ï×î¼Ñ´ð°¸, 3Q 2014-02-18 15:18:45
ÒýÓûØÌû:
3Â¥: Originally posted by liqiang0915 at 2014-02-18 10:30:35
ËãµÄʱºò»¹ÊdzöÏÖ´íÎ󣬷³ÀÍÄúÔÙ¿´¿´¡£
??? Error using ==> feval
Undefined function or method 'Equations' for input arguments of type 'double'.

Error in ==> odearguments at 110
f0 = feval( ...

ÎÒÔÚ¶þÂ¥Ìù³öµÄ´úÂ룬ÕûÌ叴֯£¬Õ³Ìù½øÒ»¸ömÎļþÖУ¬ÔËÐС£
ÎÒÔÚÉÏÎç11µã×óÓÒ²âÊÔ£¬ÎÞ´íÎó¡£
MATLAB¡¢MSСÎÊÌâ¡¢ÆÕͨÎÊÌâÇë·¢ÌûÇóÖú£¡Ê±¼ä¾«Á¦ÓÐÏÞ£¬Ë¡²»½ÓÊÜÎÞ³¥Ë½ÐÅÇóÖú¡£
4Â¥2014-02-18 11:07:41
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ liqiang0915 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ÇóÖú071001µ÷¼Á£¡£¡£¡ +3 »ÆÊØËÉ 2026-04-05 4/200 2026-04-05 21:01 by barlinike
[¿¼ÑÐ] Ò»Ö¾Ô¸ÉúÎïÓëÒ½Ò©£¬296·Ö£¬Çóµ÷¼Á +9 66¹ 2026-04-03 10/500 2026-04-05 20:11 by lys0704
[¿¼ÑÐ] 296Çóµ÷¼Á +3 Íô£¡£¿£¡ 2026-04-05 5/250 2026-04-05 17:38 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 324Çóµ÷¼Á +9 ÏëÉÏѧÇóµ÷ 2026-04-03 9/450 2026-04-04 23:57 by ¹û¶³´óÍõ
[¿¼ÑÐ] 323Çóµ÷¼Á +8 Àî¼ÑÀÖ1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[¿¼ÑÐ] 333Çóµ÷¼Á +9 °¢¿ÆÒÝ 2026-03-31 9/450 2026-04-04 18:25 by macy2011
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +6 ·¶¸ùÅà 2026-04-04 6/300 2026-04-04 14:27 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] Çóµ÷¼Á +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[¿¼ÑÐ] 266·Ö£¬Ò»Ö¾Ô¸µçÆø¹¤³Ì£¬±¾¿Æ²ÄÁÏ£¬Çó²ÄÁÏרҵµ÷¼Á +9 ÍÛºôºßºôºß 2026-04-02 9/450 2026-04-03 12:05 by 1753564080
[¿¼ÑÐ] Ò»Ö¾Ô¸Éî´ó085601²ÄÁϹ¤³Ìרҵ£¨×¨Ë¶£©300·Ö¿ÉÒÔµ÷¼ÁÈ¥ÄÄ +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[¿¼ÑÐ] 312Çóµ÷¼Á +4 ÉÞÔÂÉ« 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£© +6 ÏàÐűػá¹ââÍòÕ 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©¿Æ¼¼²ÄÁÏ¿ÆÑ§Ó빤³Ì288·Ö£¬Çóµ÷¼Á +14 Êdz½°¡ 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÄÏʦ·¶´óѧ-22408¼ÆËã»ú-292·Ö-Çó»ªÄÏʦ·¶´óѧµ÷¼Á +4 °®¶ÁÊéµÄСöùÓã 2026-04-02 4/200 2026-04-02 18:35 by Çóµ÷¼Ázz
[¿¼ÑÐ] 271Çóµ÷¼Á +15 ÀÕ²¼ÀÊ@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[¿¼ÑÐ] 270µ÷¼Á +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[¿¼ÑÐ] 310·ÖÇóµ÷¼Á +4 ³É¹¦Éϰ¶wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[¿¼ÑÐ] ²ÄÁÏרҵµ÷¼Á +5 À²À²À²¿Þ 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[¿¼ÑÐ] 326Çóµ÷¼Á +4 áÌáÌ×Ð 2026-03-31 4/200 2026-04-01 09:58 by ÎҵĴ¬Îҵĺ£
[¿¼ÑÐ] ÎïÀíѧµ÷¼Á +4 СÑò36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û