版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

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

返回列表

pangteng1990

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 695.8
散金: 19
紅花: 2
帖子: 312
在線(xiàn): 72.3小時(shí)
蟲(chóng)號(hào): 2328422
注冊(cè): 2013-03-08
性別: GG
專(zhuān)業(yè): 無(wú)機(jī)化工

[求助] 濃度一時(shí)間比法確定反應(yīng)級(jí)數(shù)和計(jì)算速率常數(shù) 附帶MATLAB程序出錯(cuò)，求大神教我

edit t.m %輸入時(shí)間數(shù)據(jù)
edit c.m %輸入濃度數(shù)據(jù)、
%數(shù)據(jù)處理程序 fyjs.m
load t.m %調(diào)入時(shí)間數(shù)據(jù)
load c.m %調(diào)入濃度數(shù)據(jù)
NO=max(size(t));
CO=c(1,1);
x=[];
tb=[];
for i=2:1:NO
a=(c(1,1)-c(1,i))./c(1,1); %將濃度數(shù)據(jù)轉(zhuǎn)化為轉(zhuǎn)化率
x=[x,a];
end
t(1:1)=[];
N=max(size(t));
sn=[];
for i=1:1:N-1
for j=i+1:1:N
      b=t(1,j)./t(1,i); %計(jì)算測(cè)定數(shù)據(jù)中毎兩組數(shù)據(jù)的時(shí)間比
      tb=[tb,b];
end
end
st=sum(tb); %計(jì)算總時(shí)間比
tn=[];
for n=-3:0.001:3; %逐步逼近搜索，從-3級(jí)到3級(jí)，步長(zhǎng)為0.001
xn=[];
for i=1:1:N-1
      for j=i+1:1:N
         if n==1
            w=log(1-x(1,j))./log(1-x(1,i)); %反應(yīng)級(jí)數(shù)n=1時(shí)的濃度比
         else
            w=((1-x(1,j)).^(1-n)-1)./((1-x(1,i)).^(1-n)-1); %反應(yīng)級(jí)數(shù)n≠1時(shí)的濃度比
         end
         xn=[xn,w];
      end
end
tn=[tn,n];
sn=[sn,sum(xn)]; %建立不同級(jí)數(shù)下總濃度比集合
end
d=[];
for k=1:1:max(size(sn))
d=abs([d,sn(1,k)-st]); %建立總濃度比和總時(shí)間比的差的集合
end
for k=1:1:max(size(sn))
if abs(sn(1,k)-st)==min(d); %搜索最小差值的位置
      n=tn(1,k); %記錄最小差值的反應(yīng)級(jí)數(shù)
end
end
kc=[];
for k=1:1:NO-1
if n==1
      p=-log(1-x(1,k))./t(1,k); %計(jì)算反應(yīng)級(jí)數(shù)n=1時(shí)的速率常數(shù)
else
      p=CO.^(1-n).*((1-x(1,k)).^(1-n)-1)./(n-1)./t(1,k); %計(jì)算反應(yīng)級(jí)數(shù)n≠1時(shí)的速率常數(shù)
      kc=[kc,p];
end
end
k=sum(kc)./(N-1); %計(jì)算速率常數(shù)的平均值
n
k

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 濃度_時(shí)間比法確定反應(yīng)級(jí)數(shù)和計(jì)算速率常數(shù)_粟智.caj

2015-06-16 07:50:38, 143.42 K

附件 2 : fanyingjishu.m

2015-06-16 07:52:29, 1.33 K

» 猜你喜歡

求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有84人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)

一個(gè)人的武林

1樓 2015-06-16 07:51:47

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

智能機(jī)器人

Robot (super robot)

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

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

MATLAB應(yīng)用四階龍格庫(kù)塔法求解反應(yīng)動(dòng)力學(xué)問(wèn)題? 已經(jīng)有21人回復(fù)
如何求這個(gè)matlab程序畫(huà)出的legend曲線(xiàn)的法矢量？已經(jīng)有6人回復(fù)
微分方程ode45求解，最小二乘法優(yōu)化微分方程參數(shù)，程序運(yùn)行求助已經(jīng)有5人回復(fù)
Matlab同時(shí)擬合活化能，反應(yīng)級(jí)數(shù)，速率常數(shù) 已經(jīng)有26人回復(fù)
求用牛頓迭代法解非線(xiàn)性方程組的matlab程序已經(jīng)有4人回復(fù)
MATLAB無(wú)法計(jì)算出這個(gè)含【級(jí)數(shù)和超幾何函數(shù)】的式子已經(jīng)有6人回復(fù)
求助用matlab三次樣條法算一個(gè)反應(yīng)的反應(yīng)速率已經(jīng)有13人回復(fù)
求助一個(gè)最小二乘法matlab擬合程序。已經(jīng)有7人回復(fù)
求助��！求一程序，用matlab程序做，用蒙特卡洛方法模擬已經(jīng)有10人回復(fù)
【求助】催化反應(yīng)動(dòng)力學(xué)matlab計(jì)算各基元反應(yīng)的速率常數(shù)時(shí)，該如何避免較小量被忽略？已經(jīng)有3人回復(fù)
【求助】MATLAB編寫(xiě)最大熵法計(jì)算概率密度的程序已經(jīng)有8人回復(fù)

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

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

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