版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者hzd250將贈送您 10 個金幣

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

[求助] matlab擬合反應動力學已有2人參與

剛入門matlab的小白，最近想做反應的動力學，按照B站up主的視頻自己寫了一段代碼，但是運行總是出問題：

錯誤使用 odearguments (第 93 行)；FUNC 必須返回列向量，出錯 ode45 (第 115 行)，odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
出錯 Kinetics>fun (第 67 行)，[t,x]=ode45(@func,tspan,x0,[],k)；
出錯 lsqnonlin (第 218 行)，initVals.F = feval(funfcn{3},xCurrent,varargin{:});出錯 Kinetics (第 18 行)，lsqnonlin(@fun,k0,lb,ub,[],yexp);%非線性最小二乘法。原因:Failure in initial objective function evaluation. LSQNONLIN cannot continue.

下面是我寫的代碼，讀取的Excel表格里有5列*7行的實驗數據，勞煩大佬幫我瞅瞅哪里需要改動，萬分感謝。

function Kinetics
%反應一：A+B=C+M
%r=k*XA*XB-K*XC*XM
%反應二：A+C=D+M
%r=K*XA*XC-K*XD*XM
%反應三：B+C=E+M
%r=K*XB*XC-K*XE*XM
%XM=0.175
clc
clear all;
global a b
tspan=[0.5 1 4 6 8 12 16];
yexp=xlsread('reaction.xls');
k0=[0.1 0.01 0.01 0.001 0.001 0.001];%參數初值
lb=[0 0 0 0 0 0];%下邊界
ub=[+inf +inf +inf +inf +inf +inf];%上邊界
[k,resnorm,residual,exitflag,output,lambda,jacobian]=...
lsqnonlin(@fun,k0,lb,ub,[],yexp);%非線性最小二乘法
tspan=[0.5 1 4 6 8 12 16];
a=1;
b=a+6;
x0=yexp(a,

;%積分初值
[t,x]=ode45(@func,tspan,x0,[],k);
t1=linspace(0.5,16,200);
ya1=spline(t,x(:,1),t1);%動力學計算得到的點進行樣條插值
ya2=spline(t,x(:,2),t1);
ya3=spline(t,x(:,3),t1);
ya4=spline(t,x(:,4),t1);
ya5=spline(t,x(:,5),t1);
for m=1:7
for n=1:5
yy(a+m-1,n)=x(m,n);%每一次的值存入yy矩陣
end
end
figure(1)
plot(tspan,yexp(a:b,1),'k^',t1,ya1,'k-',tspan,yexp(a:b,2),'ro',t1,ya2,'r-',tspan,yexp(a:b,3),'bd',t1,ya3,'b-',...
tspan,yexp(a:b,4),'g*',t1,ya4,'g-',tspan,yexp(a:b,5),'yp',t1,ya5,'y-');
legend('','A濃度','','B濃度','','C濃度','','D濃度','','E濃度');
xlabel('t(h)');ylabel('濃度(mol/L)');title('170℃ 0.1wt%催化劑');
t1=linspace(0.5,16,200);
z1=spline(t,yy(1:7,1),t1);
h1=spline(t,yy(1:7,2),t1);
s1=spline(t,yy(1:7,3),t1);
b1=spline(t,yy(1:7,4),t1);
u1=spline(t,yy(1:7,5),t1);
xlswrite('result.xls',[t1' z1' h1' s1' b1' u1'],'sheet1');
xlswrite('result.xls',residual,'sheet2');
Ne = length(yexp(:,2)); %模型適定性判別
Np = length(k);
[rho2,F] = rho2_F(k,yexp,resnorm,Ne,Np);
ci=nlparci(k,residual,jacobian)
fprintf('\t k1,0=%.1f ± %.4f\n',k(1),ci(1,2)-k(1));
fprintf('\t k2,0=%.1f ± %.4f\n',k(2),ci(2,2)-k(2));
fprintf('\t k3,0=%.1f ± %.4f\n',k(3),ci(3,2)-k(3));
fprintf('\t k4,0=%.1f ± %.4f\n',k(4),ci(4,2)-k(4));
fprintf('\t k5,0=%.1f ± %.4f\n',k(5),ci(5,2)-k(5));
fprintf('\t 殘差平方和：%.3f\n',resnorm)
fprintf('\t 實驗點數和自由度分別為 Ne = %d和 Np = %d\n',Ne,Np)
fprintf('\t 決定性指標ρ^2: %.4f\n',rho2)
fprintf('\t F比: %.3f\n\n',F)
%=================================================================================
function f=fun(k,yexp)
f=[];
tspan=[0.5 1 4 6 8 12 16];
a=1;
x0=yexp(a,

;
[t,x]=ode45(@func,tspan,x0,[],k)
d=a+6;
yc1=x(:,1);
yc2=x(:,2);
yc3=x(:,3);
yc4=x(:,4);
yc5=x(:,5);
f11=yexp(a:d,1)-yc1;
f12=yexp(a:d,2)-yc2;
f13=yexp(a:d,3)-yc3;
f14=yexp(a:d,4)-yc4;
f15=yexp(a:d,5)-yc5;
ff=[f11 f12 f13 f14 f15];
f=[f;ff];
%=================================================================================
function dxdt=func(t,x,k)
r1=-k(1)*x(1)*x(2)-k(2)*x(1)*x(3)+k(4)*x(3)*0.175+k(5)*x(4)*0.175;
r2=-k(1)*x(1)*x(2)-k(3)*x(2)*x(3)+k(4)*x(3)*0.175+k(6)*x(5)*0.175;
r3=k(1)*x(1)*x(2)+k(5)*x(4)*0.175+k(6)*x(5)*0.175-k(2)*x(1)*x(3)-k(3)*x(2)*x(3)-k(4)*x(3)*0.175;
r4=k(2)*x(1)*x(3)-k(5)*x(4)*0.175;
r5=k(3)*x(2)*x(3)-k(6)*x(5)*0.175;
dxdt=[r1 r2 r3 r4 r5]
%=================================================================================
function [rho2,F] = rho2_F(k,yexp,s,Ne,Np)
y=yexp.^2;
sy = sum(y(

);
rho2 = 1 - s/sy; %rho2: 決定性指標
F = (sy - s)*(Ne-Np)/(Np*s); %F：F比

回復此樓

» 猜你喜歡

幫我的英語口語老師找學生已經有4人回復
計算三階力常數用thirdorder擴胞時參數如何選擇？臨近原子數怎么取值已經有7人回復
物理學I論文潤色/翻譯怎么收費? 已經有69人回復
2026年循環(huán)經濟功能材料國際會議（ICFMCE 2026）已經有0人回復
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經有0人回復
哈爾濱理工大學物理系招收物理學考研調劑已經有1人回復
0702一志愿吉大B區(qū)求調劑已經有2人回復
求調劑已經有0人回復

1樓 2021-05-12 19:59:54

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

hzlhm

至尊木蟲 (著名寫手)

應助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時
蟲號: 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動力系統(tǒng)

【答案】應助回帖

有具體的數據嗎？可以發(fā)給我嗎？，可以試一試幫你找問題？

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

hzd250

QQ：2120156492

2樓2021-05-17 22:00:00

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

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

送紅花一朵

引用回帖:

2樓: Originally posted by hzlhm at 2021-05-17 22:00:00
有具體的數據嗎？可以發(fā)給我嗎？，可以試一試幫你找問題？

有的，太謝謝你了

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

3樓2021-05-18 10:22:08

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

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

引用回帖:

3樓: Originally posted by hzd250 at 2021-05-18 10:22:08
有的，太謝謝你了

...

這個是簡化的數據

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

4樓2021-05-18 10:52:35

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

hzlhm

至尊木蟲 (著名寫手)

應助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時
蟲號: 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動力系統(tǒng)

【答案】應助回帖

引用回帖:

4樓: Originally posted by hzd250 at 2021-05-18 10:52:35
這個是簡化的數據

...

每列數據對應的變量是什么？

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

hzd250

QQ：2120156492

5樓2021-05-18 18:39:58

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

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

每列數據對應的變量就是五中物質的濃度變化：x1 x2 x3 x4 x5

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

6樓2021-05-18 21:31:33

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

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

送紅花一朵

引用回帖:

5樓: Originally posted by hzlhm at 2021-05-18 18:39:58
每列數據對應的變量是什么？...

每列數據對應的變量就是五中物質的濃度變化：x1 x2 x3 x4 x5，就是后面的微分速率方程里面的x1 x2 x3 x4 x5

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

7樓2021-05-18 21:34:00

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★ ★
獨孤神宇: 金幣+2, 鼓勵交流 2021-05-20 21:35:31

引用回帖:

4樓: Originally posted by hzd250 at 2021-05-18 10:52:35
這個是簡化的數據

...

參考下：

Root of Mean Square Error (RMSE): 0.0602102835008095
Sum of Squared Residual: 0.108758347177436
Correlation Coef. (R): 0.963380309380527
R-Square: 0.928101620502121

Parameter Best Estimate
-------------------- -------------
k1 0.0431976858691231
k2 0.0204925345429154
k4 8.27595775884273E-20
k5 3.49573492337918E-16
k3 0.0313036416446144
k6 2.17399315624263E-18

贊一下

回復此樓

» 本帖已獲得的紅花（最新10朵）

hzd250

8樓2021-05-19 20:54:58

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

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

送紅花一朵

引用回帖:

8樓: Originally posted by dingd at 2021-05-19 20:54:58
參考下：

Root of Mean Square Error (RMSE): 0.0602102835008095
Sum of Squared Residual: 0.108758347177436
Correlation Coef. (R): 0.963380309380527
R-Square: 0.928101620502121

Parameter ...

太謝謝啦！
大佬，輸出結果怎么才能用科學記數法表示呢，我這邊輸出的數據都是小數表示的

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

9樓2021-05-20 12:52:15

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

hzd250

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 227.1
紅花: 1
帖子: 124
在線: 8.5小時
蟲號: 14587529
注冊: 2019-03-13
專業(yè): 高分子合成化學

送紅花一朵

引用回帖:

大佬，我運行的k1 k2 k3值都和你給的一樣，但是后面三個k值都是2.22×10-4這是因為什么原因呢

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

10樓2021-05-20 13:09:46

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

相關版塊跳轉我要訂閱樓主 hzd250 的主題更新

返回列表

不應助 確定回帖應助 (注意：應助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

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

[考研] 0856材料與化工353分求調劑 +8	NIFFFfff 2026-03-09	8/400	2026-03-10 15:17 by houyaoxu
[考研] 材料與化工求調劑 +11	與冬清寧 2026-03-07	14/700	2026-03-10 14:23 by Linda Hu
[基金申請] PDF在線壓縮軟件 10+3	haxia 2026-03-08	4/200	2026-03-10 13:44 by lfy8008
[考研] 調劑 +4	調劑的考研學生 2026-03-09	4/200	2026-03-09 20:22 by wsxw
[考研] 0832食品科學與工程293調劑 +3	東東不東 2026-03-07	3/150	2026-03-09 15:58 by ls_13456
[基金申請] 第二個面上項目的獲批幾率如何？ +3	wawm9296 2026-03-04	5/250	2026-03-09 13:00 by xujie4132
[考研] 一志愿蘇大材料工程專碩293求調劑 +5	鋼鐵大炮 2026-03-04	5/250	2026-03-09 09:49 by EBSD
[考研] 材料與化工，291分，求調劑 +13	咕嚕咕嚕123123 2026-03-04	15/750	2026-03-09 09:32 by 30660438
[考研] 085701環(huán)境工程專碩求調劑 +5	幼稚園大班打手 2026-03-08	5/250	2026-03-09 07:07 by houyaoxu
[考研] 085701環(huán)境工程317分求調劑 +9	6汆尼9 2026-03-07	9/450	2026-03-08 06:41 by 劉兵
[考研] 求調劑 +4	呼呼？~+123456 2026-03-06	4/200	2026-03-06 23:11 by L135790
[考研] 306求調劑 +7	Bahati 2026-03-05	7/350	2026-03-06 22:11 by 星空星月
[考研] 材料專碩290求調劑 +8	杰尼龜aaa 2026-03-04	8/400	2026-03-06 16:06 by @颯颯颯颯
[考博] 2026申博自薦六級440電催化方向 +4	櫻落成影花成雙 2026-03-05	4/200	2026-03-06 07:57 by L135790
[考研] 化工282求調劑一志愿211 +5	NA0912 2026-03-05	6/300	2026-03-05 20:10 by 2735147993
[考研] 材料調劑 +4	L9370 2026-03-05	4/200	2026-03-05 19:36 by caszguilin
[考研] 一志愿西交化工專碩288專業(yè)課93求調劑 +6	好運好運接接 2026-03-04	6/300	2026-03-05 09:36 by houyaoxu
[考研] 0856材料專碩274能調劑去哪里？ +3	22735 2026-03-04	4/200	2026-03-05 09:06 by 斬魂滴兔子！
[論文投稿] 主編終審卡者不動 100+4	Stray2021 2026-03-03	4/200	2026-03-05 08:04 by bdhaoxiang
[考研] 085601 材料305分求助 +4	泡泡郵件 2026-03-03	6/300	2026-03-04 13:00 by 無際的草原

24小時熱門版塊排行榜

[求助] matlab擬合反應動力學 已有2人參與

» 猜你喜歡

【答案】應助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應助回帖

» 本帖已獲得的紅花（最新10朵）

[求助] matlab擬合反應動力學已有2人參與