 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 721 | 回復(fù): 2 | ||
happytan新蟲 (小有名氣)
|
[求助]
請(qǐng)那位幫我看看這個(gè)MATLAB程序那里有問題！先謝謝
|
|
請(qǐng)那位幫我看看這個(gè)MATLAB程序那里有問題！先謝謝 function RH_Cat_Kinetics %用數(shù)值積分法進(jìn)行反應(yīng)速率分析得到速率常數(shù)，并將結(jié)果與實(shí)驗(yàn)值對(duì)比分析 clear all;clc global C kk0 t=[0:300:2400]; % t:時(shí)間，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℃的動(dòng)力學(xué)實(shí)驗(yàn)數(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]; %時(shí)間階梯 C0=[0;0;0.1]; % 0.1:氧氣在反應(yīng)體系氣相中的初始分率，考慮了環(huán)己烷的飽和蒸 %汽壓，在實(shí)驗(yàn)條件下為50%左右，除去該蒸汽壓以及惰性氣體所 %占的比例后得出的上述數(shù)值 kk=lsqnonlin(@myfunc,kk0,0,inf,[ ],tspan,C0,C); %非線性最小二乘法擬合。調(diào)用函數(shù)myfunc ci=nlparci(beta,resid,jacobian); %回歸系數(shù)，殘差，雅閣比矩陣 %擬合效果圖(實(shí)驗(yàn)值與擬合值的比較) M=[ 1 0 0 0; %解微分代數(shù)方程組時(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('時(shí)間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; %定義待求的動(dòng)力學(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%計(jì)算惰性氣體的體積流量 Vm=22.4;%氣體的摩爾體積，,L Vl=0.3; % VL:液體體積，按300mL計(jì)算 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)時(shí)氧氣在氣液相中的物料平衡方程左端 dCdt=[dCPdt;dC1dt;gas_mass;liquid_mass]; |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研生物與醫(yī)藥調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)296求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有7人回復(fù)
-
安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑
已經(jīng)有7人回復(fù)
-
英一數(shù)一總分334求調(diào)劑
已經(jīng)有4人回復(fù)
-
070300一志愿211，312分求調(diào)劑院校
已經(jīng)有10人回復(fù)
-
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
362求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 幫忙看看這個(gè)平面圖用matlab怎么畫
已經(jīng)有19人回復(fù)
- 師哥師姐幫我看看我這個(gè)邀請(qǐng)函符合要求嗎~~先謝謝了
已經(jīng)有12人回復(fù)
- 請(qǐng)大家?guī)臀铱纯催@個(gè)核磁譜圖吧，實(shí)在是沒有辦法了
已經(jīng)有30人回復(fù)
- 【求助】請(qǐng)教一個(gè)matlab程序，謝謝幫忙��！
已經(jīng)有15人回復(fù)
- 【請(qǐng)教】大家?guī)臀铱纯催@種薄膜樣品的截面SEM照片是怎么獲得的，謝謝了
已經(jīng)有18人回復(fù)
- 【求助】急急！��！哪位懂xps，幫我看一下我的數(shù)據(jù)吧。謝謝
已經(jīng)有5人回復(fù)
- 【求助】大家?guī)臀铱纯催@個(gè)MATLAB程序？
已經(jīng)有6人回復(fù)
2樓2011-09-28 16:59:12
信彼南山
木蟲 (著名寫手)
- 應(yīng)助: 33 (小學(xué)生)
- 金幣: 4142.9
- 散金: 1221
- 紅花: 16
- 帖子: 1178
- 在線: 233.5小時(shí)
- 蟲號(hào): 1133529
- 注冊(cè): 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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿a區(qū)211，085601-307分求調(diào)劑 +7 | 黨嘉豪 2026-03-31 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 347求調(diào)劑 +11 | 山頂見α 2026-03-25 | 11/550 |
|
|
|
[考研] 調(diào)劑310 +13 | 溫柔的晚安 2026-03-25 | 14/700 |
|
|
|
[考研] 286求調(diào)劑 +5 | 丟掉懶惰 2026-03-27 | 8/400 |
|
|
|
[考博] 材料專業(yè)申博 +5 | 杜雨婷dyt 2026-03-29 | 5/250 |
|
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +11 | lllllos 2026-03-30 | 11/550 |
|
|
|
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|