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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者h(yuǎn)zd250將贈送您 10 個金幣

hzd250

新蟲 (小有名氣)

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

[求助] matlab擬合反應(yīng)動力學(xué) 已有2人參與

剛?cè)腴Tmatlab的小白，最近想做反應(yīng)的動力學(xué)，按照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行的實驗數(shù)據(jù)，勞煩大佬幫我瞅瞅哪里需要改動，萬分感謝。

function Kinetics
%反應(yīng)一：A+B=C+M
%r=k*XA*XB-K*XC*XM
%反應(yīng)二：A+C=D+M
%r=K*XA*XC-K*XD*XM
%反應(yīng)三：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];%參數(shù)初值
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);%動力學(xué)計算得到的點進(jìn)行樣條插值
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 實驗點數(shù)和自由度分別為 Ne = %d和 Np = %d\n',Ne,Np)
fprintf('\t 決定性指標(biāo)ρ^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: 決定性指標(biāo)
F = (sy - s)*(Ne-Np)/(Np*s); %F：F比

回復(fù)此樓

» 猜你喜歡

德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生已經(jīng)有0人回復(fù)
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試已經(jīng)有11人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有206人回復(fù)
2026-CJ開始申報了已經(jīng)有1人回復(fù)
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士已經(jīng)有13人回復(fù)
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有1人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有2人回復(fù)

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

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

hzlhm

至尊木蟲 (著名寫手)

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

【答案】應(yīng)助回帖

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

贊一下(1人)

回復(fù)此樓

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

hzd250

QQ：2120156492

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

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

hzd250

新蟲 (小有名氣)

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

送紅花一朵

引用回帖:

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

有的，太謝謝你了

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

贊一下

回復(fù)此樓

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

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

hzd250

新蟲 (小有名氣)

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

引用回帖:

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

...

這個是簡化的數(shù)據(jù)

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

贊一下

回復(fù)此樓

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

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

hzlhm

至尊木蟲 (著名寫手)

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

【答案】應(yīng)助回帖

引用回帖:

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

...

每列數(shù)據(jù)對應(yīng)的變量是什么？

贊一下(1人)

回復(fù)此樓

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

hzd250

QQ：2120156492

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

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

hzd250

新蟲 (小有名氣)

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

每列數(shù)據(jù)對應(yīng)的變量就是五中物質(zhì)的濃度變化：x1 x2 x3 x4 x5

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

贊一下

回復(fù)此樓

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

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

hzd250

新蟲 (小有名氣)

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

送紅花一朵

引用回帖:

5樓: Originally posted by hzlhm at 2021-05-18 18:39:58
每列數(shù)據(jù)對應(yīng)的變量是什么？...

每列數(shù)據(jù)對應(yīng)的變量就是五中物質(zhì)的濃度變化：x1 x2 x3 x4 x5，就是后面的微分速率方程里面的x1 x2 x3 x4 x5

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

贊一下

回復(fù)此樓

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

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

dingd

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

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

【答案】應(yīng)助回帖

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

引用回帖:

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

...

參考下：

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

贊一下

回復(fù)此樓

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

hzd250

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

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

hzd250

新蟲 (小有名氣)

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

送紅花一朵

引用回帖:

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 ...

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

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

贊一下

回復(fù)此樓

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

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

hzd250

新蟲 (小有名氣)

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

送紅花一朵

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

返回列表

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

普通表情龍兔虎貓

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

[考研] 一志愿中南大學(xué)理學(xué)化學(xué) +4	15779376950 2026-03-01	6/300	2026-03-04 09:03 by 無際的草原
[考研] 0703化學(xué)306調(diào)劑 +4	26要上岸 2026-03-03	4/200	2026-03-03 20:08 by yoyohj0427
[考研] 290求調(diào)劑 +9	ErMiao1020 2026-03-02	9/450	2026-03-03 18:03 by linlonghao
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 材料270求調(diào)劑 6+6	Eiiiio 2026-03-01	11/550	2026-03-03 15:04 by Phoebelala
[考研] 材料學(xué)碩318求調(diào)劑 +11	February_Feb 2026-03-01	11/550	2026-03-03 14:24 by NUAAZXWS
[考研] 0857調(diào)劑 +6	一ll半 2026-02-28	8/400	2026-03-03 09:30 by 熱情沙漠
[考研] 324求調(diào)劑 +4	wxz2 2026-03-03	5/250	2026-03-03 09:25 by 楊楊楊紫
[考研] 338求調(diào)劑 +5	18162027187 2026-03-02	6/300	2026-03-03 07:52 by njzyff
[考研] 0856調(diào)劑 +10	劉夢微 2026-02-28	10/500	2026-03-02 23:42 by hcy618
[考研] 材料化工調(diào)劑 +13	今夏不夏 2026-03-01	16/800	2026-03-02 22:11 by sunny81
[考研] 求調(diào)劑 +7	repeatt?t 2026-02-28	7/350	2026-03-02 20:34 by hypershenger
[考研] 材料復(fù)試調(diào)劑 +5	學(xué)材料的點 2026-03-01	6/300	2026-03-02 20:01 by hypershenger
[考研] 材料085601調(diào)劑 +5	多多子. 2026-03-02	5/250	2026-03-02 19:15 by zhukairuo
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3	shs1083 2026-03-02	3/150	2026-03-02 17:27 by houyaoxu
[考研] 303求調(diào)劑 +5	今夏不夏 2026-03-01	5/250	2026-03-02 15:01 by 向上的胖東
[考博] 誠招農(nóng)業(yè)博士 +3	心欣向榮 2026-02-28	3/150	2026-03-02 13:33 by 時間不狗
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6	禮堂丁真258 2026-02-28	9/450	2026-03-02 12:04 by 52hz~~
[考研] 272求調(diào)劑 +6	田智友 2026-02-28	6/300	2026-03-01 21:40 by 公瑾逍遙
[考研] 尋找調(diào)劑 +4	LYidhsjabdj 2026-02-28	4/200	2026-03-01 10:56 by sunny81