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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

zhaoshazhu

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

[求助] Matlab指導(dǎo)

這是我編的程序代碼，但不知哪地方有問題，請高手幫忙指點一下。
function KineticsEstHydrogenation
%

clear all
clc
format long
%    KH KDMM KMe KHPM KPDO k1    k2    k3
k0 = [2.780  4.049 3.225  1.540  2.456 27.967 42.249 10];
lb = [0 0 0 0 0 0 0 0];                % 參數(shù)下限
ub = [+inf  +inf  +inf  +inf +inf +inf +inf  +inf]; % 參數(shù)上限
% 動力學(xué)數(shù)據(jù):
% T =195
%  t    H2in DMMin Mein HPMin PDOin NPAin HPMout PDOout NPAout
data=...
[2.05 0.00258    0.92932 0.06810    0 0 0 3.1360E-05 1.1642E-03 0.0011365;
1.37 0.00155    0.93029 0.06817    0 0 0 7.8438E-05 1.3477E-03 0.0005733;
1.62 0.00173    0.91758 0.08068    0 0 0 1.0444E-04 1.7271E-03 0.0004754;
1.98 0.00187    0.89928 0.09884    0 0 0 2.8942E-04 2.3631E-03 0.0004436;
1.21 0.00158    0.91772 0.08070    0 0 0 4.1373E-04 1.8968E-03 0.0002797;
1.49 0.00224    0.89895 0.09881    0 0 0 6.4058E-04 2.0193E-03 0.0002921;
2.70 0.00531  0.81544 0.17925    0 0 0 1.3040E-03 2.3203E-03 0.0003679;
1.02 0.00254    0.92936 0.06810    0 0 0 2.0611E-04 1.7761E-03 0.0003612;
0.82 0.00273    0.92918 0.06809    0 0 0 4.6464E-04 1.6084E-03 0.0001887;
0.97 0.00213    0.91722 0.08065    0 0 0 7.9002E-04 1.3764E-03 0.0001489;
1.19 0.00379    0.89756 0.09865    0 0 0 9.7162E-04 1.7297E-03 0.0002101;
2.16 0.00663    0.81436 0.17902    0 0 0 1.7438E-03 2.1190E-03 0.0003488
];

% 使用函數(shù)fmincon()進(jìn)行參數(shù)估計
[k,fval,flag] = fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[]);
fprintf('\n使用函數(shù)fmincon()估計得到的參數(shù)值為:\n')
fprintf('\tk1(KH) = %.4f\n',k(1))
fprintf('\tk2(KDMM) = %.4f\n',k(2))
fprintf('\tk3(KMe) = %.4f\n',k(3))
fprintf('\tk4(KHPM) = %.4f\n',k(4))
fprintf('\tk5(KPDO) = %.4f\n',k(5))
fprintf('\tk6(k1) = %.4f\n',k(6))
fprintf('\tk7(k2) = %.4f\n',k(7))
fprintf('\tk8(k3) = %.4f\n',k(8))
fprintf('  The sum of the squares is: %.1e\n\n',fval)
k_fmincon = k;

% 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[]);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')
Output

% 以函數(shù)fmincon()估計得到的結(jié)果為初值，使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計
k0 = k_fmincon;
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[]);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n')
Output
function f = ObjFunc4Fmincon(k)
%  t    H2in DMMin Mein HPMin PDOin NPAin HPMout PDOout NPAout
data=...
[2.05 0.00258    0.92932 0.06810    0 0 0 3.1360E-05 1.1642E-03 0.0011365;
1.37 0.00155    0.93029 0.06817    0 0 0 7.8438E-05 1.3477E-03 0.0005733;
1.62 0.00173    0.91758 0.08068    0 0 0 1.0444E-04 1.7271E-03 0.0004754;
1.98 0.00187    0.89928 0.09884    0 0 0 2.8942E-04 2.3631E-03 0.0004436;
1.21 0.00158    0.91772 0.08070    0 0 0 4.1373E-04 1.8968E-03 0.0002797;
1.49 0.00224    0.89895 0.09881    0 0 0 6.4058E-04 2.0193E-03 0.0002921;
2.70 0.00531  0.81544 0.17925    0 0 0 1.3040E-03 2.3203E-03 0.0003679;
1.02 0.00254    0.92936 0.06810    0 0 0 2.0611E-04 1.7761E-03 0.0003612;
0.82 0.00273    0.92918 0.06809    0 0 0 4.6464E-04 1.6084E-03 0.0001887;
0.97 0.00213    0.91722 0.08065    0 0 0 7.9002E-04 1.3764E-03 0.0001489;
1.19 0.00379    0.89756 0.09865    0 0 0 9.7162E-04 1.7297E-03 0.0002101;
2.16 0.00663    0.81436 0.17902    0 0 0 1.7438E-03 2.1190E-03 0.0003488
];
yexp = [data(:,8),data(:,9),data(:,10)];
for i=1:12
[t x] = ode45(@KineticEqs,[0,data(:,1)],[data(:,2),data(:,3),data(:,4),data(:,5),data(:,6),data(:,7)],[],k);
y(:,1:3) = x(:,4:6)
f = sum((y(:,1)-yexp(:,1)).^2+(y(:,2)-yexp(:,2)).^2+(y(:,3)-yexp(:,3)).^2);
end
% ------------------------------------------------------------------
function f = ObjFunc4LNL(k)
%  t    H2in DMMin Mein HPMin PDOin NPAin HPMout PDOout NPAout
data=...
[2.05 0.00258    0.92932 0.06810    0 0 0 3.1360E-05 1.1642E-03 0.0011365;
1.37 0.00155    0.93029 0.06817    0 0 0 7.8438E-05 1.3477E-03 0.0005733;
1.62 0.00173    0.91758 0.08068    0 0 0 1.0444E-04 1.7271E-03 0.0004754;
1.98 0.00187    0.89928 0.09884    0 0 0 2.8942E-04 2.3631E-03 0.0004436;
1.21 0.00158    0.91772 0.08070    0 0 0 4.1373E-04 1.8968E-03 0.0002797;
1.49 0.00224    0.89895 0.09881    0 0 0 6.4058E-04 2.0193E-03 0.0002921;
2.70 0.00531  0.81544 0.17925    0 0 0 1.3040E-03 2.3203E-03 0.0003679;
1.02 0.00254    0.92936 0.06810    0 0 0 2.0611E-04 1.7761E-03 0.0003612;
0.82 0.00273    0.92918 0.06809    0 0 0 4.6464E-04 1.6084E-03 0.0001887;
0.97 0.00213    0.91722 0.08065    0 0 0 7.9002E-04 1.3764E-03 0.0001489;
1.19 0.00379    0.89756 0.09865    0 0 0 9.7162E-04 1.7297E-03 0.0002101;
2.16 0.00663    0.81436 0.17902    0 0 0 1.7438E-03 2.1190E-03 0.0003488
];
for i=1:12
yexp = [data(:,8),data(:,9),data(:,10)];
[t x] = ode45(@KineticEqs,[0,data(:,1)],[data(:,2),data(:,3),data(:,4),data(:,5),data(:,6),data(:,7)],[],k);
y(:,1:3) = x(:,4:6);
f1 = y(:,1) - yexp(:,1);
f2 = y(:,2) - yexp(:,2);
f3 = y(:,3) - yexp(:,3);
f = [f1; f2; f3];
end

% ------------------------------------------------------------------
function dxdt = KineticEqs(t,x,k)
% 動力學(xué)方程，氫氣與DMO均解離吸附，表面反應(yīng)為控制步驟。
r1 = k(6).*((k(1).*x(1).*k(2).*x(2)).^0.5)./...
(1+2.*(k(2).*x(2)).^0.5+2.*(k(4).*x(4)).^0.5+(k(1).*x(1)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;

r2 = k(7).*((k(1).*x(1).*k(4).*x(4)).^0.5)./...
(1+2.*(k(2).*x(2)).^0.5+2.*(k(4).*x(4)).^0.5+(k(1).*x(1)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;

r3 = k(8).*k(5).*k(2).*x(5).*x(2)./...
(1+2.*(k(1).*x(1)).^0.5+2.*(k(4).*x(4)).^0.5+(k(2).*x(2)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;

dxdt=...
[-2.*r1-2.*r2-r3    % H2
-r1                % DMM
r1+r2                % ME
r1-r2                % HPM
r2-r3                % PDO
r3                   % NPA
   ];

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有162人回復(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ù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

1樓 2014-12-26 16:13:44

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

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

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷
[考研] 北科281學(xué)碩材料求調(diào)劑 +6	tcxiaoxx 2026-03-20	6/300	2026-03-22 20:23 by edmund7
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6	皮皮皮123456 2026-03-22	6/300	2026-03-22 16:52 by i_cooler
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

zhaoshazhu

[求助] Matlab指導(dǎo)

» 猜你喜歡

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频