 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 459 | 回復(fù): 0 | ||
zhang66163銀蟲 (正式寫手)
|
[求助]
Matlab 程序求助
|
|
各位大俠，用matlab求擬合動(dòng)力學(xué)數(shù)據(jù)時(shí)，寫了以下的程序。結(jié)果出現(xiàn)錯(cuò)誤是：Too many output arguments. 請(qǐng)教以下如何對(duì)程序進(jìn)行改進(jìn)。謝謝 clear all clc T = 150 + 273.15; R = 8.31434; a=0.1; b=0.2; c=60; d=0.257; e=-0.691; f=-0.465; g=-0.238; h=1.7; i=2.5; % Arrhenius constant, 1/s k0 = [5.39E+13 9.55 2.04 5.80 3.76E+07 110.39]; % Activation energy, kJ/kmol Ea = [141.3E+3, 49.9E+3, 42.6E+3, 46.3E+3, 79.9E+3, 52.7E+3]; % 初始濃度C0(i), mol/m^3 C0 = [4.96 0 0 0 0 0 0]; tspan = [0 1e4]; [t,C] = ode45(@MassEquations, tspan, C0,[],k0,Ea,R,T) % 繪圖 plot(t,C(:,1),'r-',t,C(:,2),'k:',t,C(:,3),'b-.',t,C(:,4),'k--',t,C(:,5),'g*',t,C(:,5),'r-',t,C(:6),'k:',t,C(:7),'b*'); xlabel('Time (s)'); ylabel('Concentration (mol/m^3)'); legend('A','B','C','D','E','F','G') CDmax = max(D(:,4)); % CDmax: the maximum concentration of B, kmol/m^3 yDmax = CDmax/3/C0(1) % yBmax: the maximum yield of B index = find(C(:,4)==CDmax); t_opt = t(index) % t_opt: the optimum batch time, s % ------------------------------------------------------------------ function dCdt = MassEquations(t,C,k0,Ea,R,T) % Reaction rate constants, 1/s k = k0.*exp(-Ea/(R*T)); % Reaction rates, kmoles/m3 s rA = -((a^d)*k(1)+a^e*b*k(2)*c)*C(1); rB = (a^d)*k(1)*C(1)+(a^d)*2*k(5)*(C(5)^f)-(a^g)*b*k(3)*C(3)*c-a^d*k(6)*(C(3)^i); rC = a^d*k(1)*C(1)-a^f*b*k(4)*C(5)*c-a^d*k(5)*(C(5)^f); rD = (a^g)*b*k(3)*C(3)*c; rE = (a^f)*b*k(4)*C(5)*c; rF = (a^e)*b*k(2)*C(1)*c; rG = (a^d)*k(6)*(C(3)^i); % Mass balances dCdt = [rA; rB; rC; rD; rE; rF; rG]; |
回復(fù)此樓» 猜你喜歡
焦慮
已經(jīng)有13人回復(fù)
-
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有9人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有5人回復(fù)
-
285化工學(xué)碩求調(diào)劑（081700）
已經(jīng)有11人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有9人回復(fù)
-
274求調(diào)劑
已經(jīng)有6人回復(fù)
-
求材料調(diào)劑
已經(jīng)有10人回復(fù)
-
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑
已經(jīng)有4人回復(fù)
-
國(guó)自科面上基金字體
已經(jīng)有8人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhang66163 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +9 | rare12345 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +7 | iwuli 2026-03-12 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|