版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

liqiang0915

木蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1718
帖子: 18
在線: 352.6小時
蟲號: 1455044
注冊: 2011-10-22
專業(yè): 化學反應工程

[求助] 哪位大神幫忙看一下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];

回復此樓

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

1樓 2014-02-17 16:13:25

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經驗: +1059
計算強帖: 8
應助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

感謝參與，應助指數 +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];

附圖1.jpg

贊一下(1人)

回復此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2014-02-17 21:02:13

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

liqiang0915

木蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1718
帖子: 18
在線: 352.6小時
蟲號: 1455044
注冊: 2011-10-22
專業(yè): 化學反應工程

引用回帖:

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 ...

算的時候還是出現錯誤，煩勞您再看看。
??? 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發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經驗: +1059
計算強帖: 8
應助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(qū)

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liqiang0915: 金幣+50, ★★★★★最佳答案, 3Q 2014-02-18 15:18:45

引用回帖:

3樓: Originally posted by liqiang0915 at 2014-02-18 10:30:35
算的時候還是出現錯誤，煩勞您再看看。
??? 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小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2014-02-18 11:07:41

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 liqiang0915 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求調劑 +6	十三加油 2026-03-21	6/300	2026-03-22 17:00 by i_cooler
[考研] 尋找調劑 +4	倔強芒? 2026-03-21	4/200	2026-03-22 16:14 by 木托莫露露
[考研] 319求調劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 298求調劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 291 求調劑 +3	化工2026屆畢業(yè)?/a> 2026-03-21	3/150	2026-03-22 14:26 by ColorlessPI
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 一志愿華中科技大學071000，求調劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 廣西大學材料導師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 297求調劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學員8dgXkO
[考研] 一志愿深大，0703化學，總分302，求調劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學員8dgXkO
[考研] 0805材料320求調劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 301求調劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 求調劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 材料與化工專碩調劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調劑） +3	123阿標 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 288求調劑，一志愿華南理工大學071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 085601材料工程專碩求調劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 286求調劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 一志愿南京大學，080500材料科學與工程，調劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 283求調劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

liqiang0915

[求助] 哪位大神幫忙看一下MATLAB求解微分方程的問題？ 已有1人參與

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

月只藍

【答案】應助回帖

liqiang0915

月只藍

【答案】應助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 哪位大神幫忙看一下MATLAB求解微分方程的問題？已有1人參與

» 本主題相關價值貼推薦，對您同樣有幫助: