 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 720 | 回復(fù): 2 | ||
happytan新蟲 (小有名氣)
|
[求助]
請那位幫我看看這個MATLAB程序那里有問題！先謝謝
|
|
請那位幫我看看這個MATLAB程序那里有問題！先謝謝 function RH_Cat_Kinetics %用數(shù)值積分法進(jìn)行反應(yīng)速率分析得到速率常數(shù)，并將結(jié)果與實驗值對比分析 clear all;clc global C kk0 t=[0:300:2400]; % t:時間，s C=[0.00 0.11 0.19 0.25 0.31 0.44 0.54 0.62 0.65; 0.00 0.10 0.18 0.24 0.28 0.39 0.48 0.52 0.53; 0.10 3.20 1.70 1.30 2.50 3.90 5.20 8.20 9.7]; C=C';%溫度150℃的動力學(xué)實驗數(shù)據(jù)C(1) &C(2): kmol/m^3, C(3):% %非線性擬合 kk0=[5.1095e-004; 2.8939;9.5453e+003]; %利用suresh模型算出同樣溫度條件下的速率常數(shù)作為其初始值 tspan=[0:300:2400]; %時間階梯 C0=[0;0;0.1]; % 0.1:氧氣在反應(yīng)體系氣相中的初始分率，考慮了環(huán)己烷的飽和蒸 %汽壓，在實驗條件下為50%左右，除去該蒸汽壓以及惰性氣體所 %占的比例后得出的上述數(shù)值 kk=lsqnonlin(@myfunc,kk0,0,inf,[ ],tspan,C0,C); %非線性最小二乘法擬合。調(diào)用函數(shù)myfunc ci=nlparci(beta,resid,jacobian); %回歸系數(shù)，殘差，雅閣比矩陣 %擬合效果圖(實驗值與擬合值的比較) M=[ 1 0 0 0; %解微分代數(shù)方程組時的系數(shù)矩陣 0 1 0 0; 0 0 0 0; 0 0 0 0]; options=odeset; %產(chǎn)生/改變參數(shù)結(jié)構(gòu) options.Mass=M; options.TolRol=le-7;%定義精度 tspan=[5:40]; C0=[0 0 0.1 ];%初值C P=1.2;%MPa G=1.5;%L/min kla=0.14;%1/s; apxl=input('Pls.input gas holdup you got according reaction conditions using hydro_mass_ c:','s'); [t_sim,C_sim]=ode45(@lj_ssm,tspan,CO,options,kk,P,G,kla,apxl); plot(t,C_exp(:,1),'-ks',t_sim,C_sim(:,1),':ko',t,C_exp(:,2),'-kd',t_sim,C_sim(:,2),':k+'); legend('Exp.CP','Model.CP','Exp.CI','Mode1.CI'),xlabel('時間t/s'),ylabel('濃度C/(kmol/m^3)'); %定義非線性優(yōu)化的目標(biāo)函數(shù) function f=myfunc(kk,tspan,C0,C) P=12; G=1.5; kla=0.14; apxl=0.2; [t,CC]=ode45(@lj_ssm,tspan,C0,[],kk,P,G,kla,apxl); f=CC-C; %定義待求的動力學(xué)方程 function dCdt=lj_ssm(t,C,kk,P,G,kla,apxl) G=G/(60*22.4);%改變氣體流量的單位，便于后面的等式成立，,mol/s yin=0.21; GI=G*(1-yin);%按照進(jìn)氣中氧氣含量21%計算惰性氣體的體積流量 Vm=22.4;%氣體的摩爾體積，,L Vl=0.3; % VL:液體體積，按300mL計算 VG1=Vl/(1-apxl); %反應(yīng)體系的總體積，L P=P*10; % P原來的單位是MPa，為了能夠利用后面的公式，將單位轉(zhuǎn)變?yōu)閎ar H=1.09e-2; % H:氧氣的Henry系數(shù) kk=kk0.*exp(-Ea/(R*T)); % T:K k0:指前因子Ea:活化能 dCdt=zeros(1,4);%先預(yù)分配空間給濃度C，這樣可以加快調(diào)用速度 dCPdt=kk(1)*kk(3)*C(1)*kk(4)/(kk(1)+kk(2)*C(1)+kk(3)*kk(4)); dC1dt=kk(3)*C(1)*kk(4)*(kk(1)-kk(2)*C(3))/(kk(1)+kk(2)*C(1)+kk(3)*kk(4)); gas_mass=Vm*(1-apxl)/(Vl*apxl)*(GI*(yin/(1-yin)-C(4)/(1-C(4)))-kla*(H*P*C(4)-kk(4))*VG1); liquid_mass=kla*(H*P*C(4)-C(2))-(kk(3)*C(1)*C(2)*(kk(1)+kk(2)*C(1))/(kk(1)+kk(2)*C(1)+kk(3)*C(2))); % gas mass & liquidse mass分別達(dá)到穩(wěn)態(tài)時氧氣在氣液相中的物料平衡方程左端 dCdt=[dCPdt;dC1dt;gas_mass;liquid_mass]; |
回復(fù)此樓» 猜你喜歡
290分調(diào)劑求助
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有12人回復(fù)
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有14人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有10人回復(fù)
-
本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分
已經(jīng)有5人回復(fù)
-
313求調(diào)劑
已經(jīng)有5人回復(fù)
-
362求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有27人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 幫忙看看這個平面圖用matlab怎么畫
已經(jīng)有19人回復(fù)
- 師哥師姐幫我看看我這個邀請函符合要求嗎~~先謝謝了
已經(jīng)有12人回復(fù)
- 請大家?guī)臀铱纯催@個核磁譜圖吧，實在是沒有辦法了
已經(jīng)有30人回復(fù)
- 【求助】請教一個matlab程序，謝謝幫忙��！
已經(jīng)有15人回復(fù)
- 【請教】大家?guī)臀铱纯催@種薄膜樣品的截面SEM照片是怎么獲得的，謝謝了
已經(jīng)有18人回復(fù)
- 【求助】急急！��！哪位懂xps，幫我看一下我的數(shù)據(jù)吧。謝謝
已經(jīng)有5人回復(fù)
- 【求助】大家?guī)臀铱纯催@個MATLAB程序？
已經(jīng)有6人回復(fù)
2樓2011-09-28 16:59:12
信彼南山
木蟲 (著名寫手)
- 應(yīng)助: 33 (小學(xué)生)
- 金幣: 4142.9
- 散金: 1221
- 紅花: 16
- 帖子: 1178
- 在線: 233.5小時
- 蟲號: 1133529
- 注冊: 2010-10-27
- 專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)
3樓2011-09-28 19:44:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 happytan 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +6 | 小黃鴨寶 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +7 | 小黃鴨寶 2026-03-30 | 7/350 |
|
|
|
[考研] 266分，求材料相關(guān)專業(yè)調(diào)劑 +9 | 哇呼哼呼哼 2026-03-30 | 11/550 |
|
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +11 | lllllos 2026-03-30 | 11/550 |
|
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25 | 稻妻小編 2026-03-26 | 25/1250 |
|
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|