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

登錄注冊

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

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

返回列表

pangteng1990

金蟲 (正式寫手)

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

[求助] 濃度一時(shí)間比法確定反應(yīng)級數(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ì)算測定數(shù)據(jù)中毎兩組數(shù)據(jù)的時(shí)間比
      tb=[tb,b];
end
end
st=sum(tb); %計(jì)算總時(shí)間比
tn=[];
for n=-3:0.001:3; %逐步逼近搜索，從-3級到3級，步長為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)級數(shù)n=1時(shí)的濃度比
         else
            w=((1-x(1,j)).^(1-n)-1)./((1-x(1,i)).^(1-n)-1); %反應(yīng)級數(shù)n≠1時(shí)的濃度比
         end
         xn=[xn,w];
      end
end
tn=[tn,n];
sn=[sn,sum(xn)]; %建立不同級數(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)級數(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)級數(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)級數(shù)n≠1時(shí)的速率常數(shù)
      kc=[kc,p];
end
end
k=sum(kc)./(N-1); %計(jì)算速率常數(shù)的平均值
n
k

回復(fù)此樓

» 本帖附件資源列表

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

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有180人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

一個(gè)人的武林

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

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

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

MATLAB應(yīng)用四階龍格庫塔法求解反應(yīng)動(dòng)力學(xué)問題? 已經(jīng)有21人回復(fù)
如何求這個(gè)matlab程序畫出的legend曲線的法矢量？已經(jīng)有6人回復(fù)
微分方程ode45求解，最小二乘法優(yōu)化微分方程參數(shù)，程序運(yùn)行求助已經(jīng)有5人回復(fù)
Matlab同時(shí)擬合活化能，反應(yīng)級數(shù)，速率常數(shù) 已經(jīng)有26人回復(fù)
求用牛頓迭代法解非線性方程組的matlab程序已經(jīng)有4人回復(fù)
MATLAB無法計(jì)算出這個(gè)含【級數(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編寫最大熵法計(jì)算概率密度的程序已經(jīng)有8人回復(fù)

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

科研從小木蟲開始，人人為我，我為人人

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