24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

xiaofenger01

金蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 709.2
帖子: 44
在線: 92.7小時(shí)
蟲(chóng)號(hào): 2318750
注冊(cè): 2013-03-05
性別: MM
專業(yè): 化學(xué)反應(yīng)工程

[求助] Matlab ode45的運(yùn)用

下面這段程序中，T是常數(shù)，所以比較簡(jiǎn)單，如果我想讓T隨著時(shí)間變化，比如寫(xiě)成T=50+0.5*tspan，這樣T和tspan一樣是1*301的向量了，下面函數(shù)MassEquations里的k也是T的函數(shù)了，這樣的話應(yīng)該如何改寫(xiě)？
謝謝

function Cha5demo4
T = 220;
R = 8.315;
k0 = [5.78052E+10 3.92317E+12 1.64254E+4 6.264E+8];
Ea = [124670 150386 77954 111528];
C0 = [1 0 0 0 0];
tspan = [0:1:300];
[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--');

function dCdt = MassEquations(t,C,k0,Ea,R,T)
k = k0.*exp(-Ea/(R*(T+273.15)));
rA = -(k(1)+k(2))*C(1);
rB = k(1)*C(1)-k(3)*C(2);
rC = k(2)*C(1)-k(4)*C(3);
rD = k(3)*C(2)-k(4)*C(4);
rE = k(4)*C(3);
dCdt = [rA; rB; rC; rD; rE];

回復(fù)此樓

» 猜你喜歡

277跪求調(diào)劑已經(jīng)有11人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有18人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué) 已經(jīng)有14人回復(fù)
面上5B能上會(huì)嗎？已經(jīng)有7人回復(fù)
0703本科鄭州大學(xué)求調(diào)劑已經(jīng)有7人回復(fù)
286分調(diào)劑已經(jīng)有3人回復(fù)
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可已經(jīng)有15人回復(fù)
317分一志愿南理工材料工程本科湖工大求調(diào)劑已經(jīng)有11人回復(fù)
332求調(diào)劑已經(jīng)有12人回復(fù)

1樓 2015-04-10 10:54:03

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛(ài)小木蟲(chóng)

找到一些相關(guān)的精華帖子，希望有用哦~

微分方程ode45求解，最小二乘法優(yōu)化微分方程參數(shù)，程序運(yùn)行求助已經(jīng)有5人回復(fù)
求助！！對(duì)于調(diào)用龍格庫(kù)塔法如何對(duì)結(jié)果判斷反饋已經(jīng)有6人回復(fù)
matlab 運(yùn)行ode45出錯(cuò)，不知道什么原因已經(jīng)有4人回復(fù)
使用matlab最優(yōu)化方法擬合獲得多個(gè)動(dòng)力學(xué)參數(shù)中的問(wèn)題已經(jīng)有4人回復(fù)
求助：用matlab或其它軟件進(jìn)行如下計(jì)算模擬，多謝已經(jīng)有3人回復(fù)
matlab編程 ODE45 已經(jīng)有4人回復(fù)
求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率已經(jīng)有13人回復(fù)
Matlab中ode45出現(xiàn)大誤差，關(guān)于用ODE45解偏微分方程的疑問(wèn)，有待解決，大家共同討論！已經(jīng)有4人回復(fù)
請(qǐng)問(wèn)什么是分段龍格庫(kù)塔的方法？需要用matlab軟件采用這種方法求二階微分方程的解已經(jīng)有3人回復(fù)
三元二次的三角方程如何用MATLAB求解？？已經(jīng)有10人回復(fù)
matlab的lsqnonlin函數(shù)怎么用已經(jīng)有7人回復(fù)
如何求解描述振動(dòng)的二階微分方程已經(jīng)有21人回復(fù)
請(qǐng)幫忙用matlab求矩陣特征值已經(jīng)有8人回復(fù)
求助利用matlab的ODE45求解微分方程已經(jīng)有16人回復(fù)
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復(fù)
【求助】MATLAB導(dǎo)入數(shù)據(jù)精度問(wèn)題。。。已經(jīng)有4人回復(fù)
【求助】急求各位大俠，關(guān)于用ode45解微分方程已經(jīng)有8人回復(fù)
【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問(wèn)題已經(jīng)有13人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xiaofenger01 的主題更新