 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 598 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】修改程序
|
|||
|
目前下面的程序能得到r=kPA^ n PB^ m中的k、n、m值，但是我想把程序修改成能直接求r=k0e(-Ea/RT)PA^nPB^m，中的k0、Ea、n和m。應(yīng)該如何修改程序，求大神賜教，不甚感激。 數(shù)據(jù)如下： T=310+273.15K時 PA0 = [3.7449 3.7500 3.7548 3.7568]; PB0 = [0.0807 0.0835 0.0861 0.0871]; r = [0.0019 0.0022 0.0024 0.0025]; T=330+273.15K時 PA0 = [3.7138 3.7215 3.7448 3.7572]; PB0 = [0.04 0.0458 0.0618 0.071]; r = [0.0041 0.0097 0.0109 0.0111]; Ps：R是常數(shù)R=8.3145 r=kPA^nPB^m的程序 function KineticsEst3 % 動力學(xué)參數(shù)辨識: 用微分法進(jìn)行反應(yīng)速率分析得到速率常數(shù)k和反應(yīng)級數(shù)n % Reaction of the type -- rate = kCA^order % order - reaction order % rate -- reaction rate vector % CA -- concentration vector for reactant A % T -- vector of reaction time % N -- number of data points % k- reacion rate constant clear all clc global negr0m PA0 PB0 PA0 = [3.7449 3.7500 3.7548 3.7568]; PB0 = [0.0807 0.0835 0.0861 0.0871]; % negr0m = [0.5 0.63 0.83 1.0 1.28 0.33 0.8 1.5 2.21 3.33]; % Problems ... negr0m = [0.0019 0.0022 0.0024 0.0025]; % 用多變量線性回歸方法估計(jì)動力學(xué)參數(shù) y = log(negr0m); x1 = log(PA0); x2 = log(PB0); y = y'; X = [ones(size(y)) x1' x2']; [b,bint] = regress(y,X,0.1); k = exp(b(1)) n = b(2) m = b(3) % 用多變量非線性回歸方法（以線性回歸的結(jié)果作為非線性回歸的初值） beta0 = [k n m]; lb = [0 0 0]; ub = [1 3 3]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,beta0,lb,ub); ci = nlparci(beta,resid,jacobian) % 殘差關(guān)于擬合值的殘差圖 negrA0c = Rate(beta,PA0,PB0); plot(negrA0c,resid,'*') xlabel('反應(yīng)速率擬合值, torr s^-^1') ylabel('殘差R, torr s^-^1') refline(0,0) % 參數(shù)辨識結(jié)果 fprintf('Estimated Parameters:\n') fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tn = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2)) fprintf('\tm = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tThe sum of the squares is: %.1e\n\n',resnorm) % ------------------------------------------------------------------ function f = ObjFunc(beta) global negr0m PA0 PB0 f = negr0m - Rate(beta,PA0,PB0); % ------------------------------------------------------------------ function negrA = Rate(beta,PA,PB) negrA = beta(1)*PA.^beta(2).*PB.^beta(3); % k=beta(1),n=beta(2),m=beta(3); |
回復(fù)此樓» 猜你喜歡
290分調(diào)劑求助
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有12人回復(fù)
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有14人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
070300化學(xué)354求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有10人回復(fù)
-
本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分
已經(jīng)有5人回復(fù)
-
313求調(diào)劑
已經(jīng)有5人回復(fù)
-
362求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有27人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
齊魯工業(yè)大學(xué)招收輕工、材料、化工、環(huán)境類研究生
+1/100
-
資源利用與植物保護(hù)專碩
+1/86
-
福建師范大學(xué)招收2026年化學(xué)、材料碩士3-4名
+1/84
-
上海電力大學(xué)先進(jìn)儲能電池技術(shù)課題組招生
+1/78
-
資源與環(huán)境專碩
+1/72
-
上海交大化院鄧楠楠課題組招收26年秋季入學(xué)申請考核博士
+1/45
-
考研調(diào)劑!
+1/41
-
歡迎化學(xué)、生物與醫(yī)藥、藥學(xué)，生物、食品等相關(guān)專業(yè)的同學(xué)
+1/37
-
中原食品實(shí)驗(yàn)室 & 鄭州輕工業(yè)大學(xué)2026年碩士研究生聯(lián)合培養(yǎng)調(diào)劑公告
+1/37
-
天津科技大學(xué)-（院士團(tuán)隊(duì)）先進(jìn)纖維與紙基功能材料團(tuán)隊(duì) 招收碩士生
+1/37
-
武漢工程大學(xué)化環(huán)學(xué)院程老師課題組招收化學(xué)類學(xué)碩、化學(xué)工程與環(huán)境工程專碩調(diào)劑
+2/36
-
南京信息工程大學(xué) 化材院接受調(diào)劑研究生
+1/36
-
【輕松準(zhǔn)入】【生物醫(yī)藥】2026青島大學(xué)招博士生
+1/19
-
長安大學(xué)智能爆破課題組招收2026級博士和碩士研究生
+1/18
-
寧夏大學(xué)化學(xué)化工學(xué)院魏逸彬研究員課題組招聘啟事
+1/7
-
江西科技師范大學(xué) 材料與能源招收調(diào)劑研究生
+1/6
-
湖北師范大學(xué)招調(diào)劑啦
+1/6
-
北京雙一流高校--環(huán)境科學(xué)與工程專業(yè)招收調(diào)劑
+1/5
-
2026年云南大學(xué)青拔課題組招收地理學(xué)，大氣科學(xué)、資源與環(huán)境等專業(yè)研究生
+1/4
-
湖北師范大學(xué) 材料科學(xué)與工程學(xué)院接收2026考研調(diào)劑
+1/2
robert2020:建議使用“引用回復(fù)該帖”，不然對方收不到你的信息。 2010-12-16 09:15:35
|
clear all clc global negr0m PA0 PB0 RT PA0 = [3.7377 3.7449 3.7500 3.7548 3.7568 3.5690 3.5859 3.5970 3.6046 3.6119 3.6161]; PB0 = [0.0769 0.0807 0.0835 0.0861 0.0871 0.0389 0.0468 0.0520 0.0554 0.0581 0.0601]; RT1=8.3145*(310+273.15)*ones(1,5); RT2=8.3145*(330+273.15)*ones(1,6); RT=[RT1 RT2]; % negr0m = [0.5 0.63 0.83 1.0 1.28 0.33 0.8 1.5 2.21 3.33]; % Problems ... negr0m = [0.0023 0.0024 0.0024 0.0025 0.0025 0.0040 0.0042 0.0044 0.0046 0.0048 0.0049]; % 用多變量線性回歸方法估計(jì)動力學(xué)參數(shù) y = log(negr0m); x1 = log(PA0); x2 = log(PB0); x3 =log(RT); y = y'; X = [ones(size(y)) x1' x2' x3']; [b,bint] = regress(y,X,0.1); k0 = exp(b(1)) n = b(2) m = b(3) Ea = b(4) % 用多變量非線性回歸方法（以線性回歸的結(jié)果作為非線性回歸的初值） beta0 = [k0 n m Ea]; lb = [0 0 0 -10]; ub = [1 3 3 10]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,beta0,lb,ub); ci = nlparci(beta,resid,jacobian) % 殘差關(guān)于擬合值的殘差圖 negrA0c = Rate(beta,PA0,PB0,RT); plot(negrA0c,resid,'*') xlabel('反應(yīng)速率擬合值, torr s^-^1') ylabel('殘差R, torr s^-^1') refline(0,0) % 參數(shù)辨識結(jié)果 fprintf('Estimated Parameters:\n') fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tn = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2)) fprintf('\tm = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tEa = %.2f ± %.2f\n',beta(4),ci(4,2)-beta(4)) fprintf('\tThe sum of the squares is: %.1e\n\n',resnorm) %------------------------------------------------------------------ function f = ObjFunc(beta) global negr0m PA0 PB0 RT f = negr0m - Rate(beta,PA0,PB0,RT); % ------------------------------------------------------------------ function negrA = Rate(beta,PA,PB,RT) negrA = beta(1)*PA.^beta(2).*PB.^beta(3).*exp(-beta(4)./RT); % k=beta(1),n=beta(2),m=beta(3);Ea=beta(4); 為什么最終的結(jié)果和 lb = [0 0 0 -10]; ub = [1 3 3 10]; 有很大關(guān)系，該如何解決 [ Last edited by 小扶搖 on 2010-12-15 at 13:35 ] |
3樓2010-12-15 13:30:40
★ ★
小扶搖(金幣+10):太謝謝了，我先研究下，有問題的話還要麻煩下 2010-12-11 09:38:01
zzuwangshilei(金幣+2):積極參與,辛苦了 2010-12-11 10:36:53
小扶搖(金幣+10):太謝謝了，我先研究下，有問題的話還要麻煩下 2010-12-11 09:38:01
zzuwangshilei(金幣+2):積極參與,辛苦了 2010-12-11 10:36:53
|
我把你的程序修改了一下，能運(yùn)行了，但是需要增加一組實(shí)驗(yàn)數(shù)據(jù)，如果不懂，請給我短消息~ function KineticsEst3_1 % 動力學(xué)參數(shù)辨識: 用微分法進(jìn)行反應(yīng)速率分析得到速率常數(shù)k和反應(yīng)級數(shù)n % Reaction of the type -- rate = kCA^order % order - reaction order % rate -- reaction rate vector % CA -- concentration vector for reactant A % T -- vector of reaction time % N -- number of data points % k- reacion rate constant clear all clc global negr0m PA0 PB0 RT PA0 = [3.7449 3.7500 3.7548 3.7568 3.7588 ]; PB0 = [0.0807 0.0835 0.0861 0.0871 0.0891]; RT=8.3*(310+273)*ones(1,5); % negr0m = [0.5 0.63 0.83 1.0 1.28 0.33 0.8 1.5 2.21 3.33]; % Problems ... negr0m = [0.0019 0.0022 0.0024 0.0025 0.0028]; % 用多變量線性回歸方法估計(jì)動力學(xué)參數(shù) y = log(negr0m); x1 = log(PA0); x2 = log(PB0); x3 =log(8.3*(310+273))*ones(1,5); y = y'; X = [ones(size(y)) x1' x2' x3']; [b,bint] = regress(y,X,0.1); k0 = exp(b(1)) n = b(2) m = b(3) Ea = b(4) % 用多變量非線性回歸方法（以線性回歸的結(jié)果作為非線性回歸的初值） beta0 = [k0 n m Ea]; lb = [0 0 0 -10]; ub = [1 3 3 10]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,beta0,lb,ub); ci = nlparci(beta,resid,jacobian) % 殘差關(guān)于擬合值的殘差圖 negrA0c = Rate(beta,PA0,PB0,RT); plot(negrA0c,resid,'*') xlabel('反應(yīng)速率擬合值, torr s^-^1') ylabel('殘差R, torr s^-^1') refline(0,0) % 參數(shù)辨識結(jié)果 fprintf('Estimated Parameters:\n') fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tn = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2)) fprintf('\tm = %.2f ± %.2f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tEa = %.2f ± %.2f\n',beta(4),ci(4,2)-beta(4)) fprintf('\tThe sum of the squares is: %.1e\n\n',resnorm) %------------------------------------------------------------------ function f = ObjFunc(beta) global negr0m PA0 PB0 RT f = negr0m - Rate(beta,PA0,PB0,RT); % ------------------------------------------------------------------ function negrA = Rate(beta,PA,PB,RT) negrA = beta(1)*PA.^beta(2).*PB.^beta(3).*exp(beta(4)./RT); % k=beta(1),n=beta(2),m=beta(3);3.75883.7588 |
2樓2010-12-11 09:28:21
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +16 | land xuxu 2026-03-26 | 16/800 |
|
|---|---|---|---|---|
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|
|
[考研] 085602化工求調(diào)劑（331分） +8 | 111@127 2026-03-30 | 8/400 |
|
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 12/600 |
|
|
|
[考研] 289求調(diào)劑 +16 | 新時代材料 2026-03-27 | 16/800 |
|
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|