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

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

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

返回列表

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

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

hzd250

新蟲(chóng) (小有名氣)

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

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

剛?cè)腴T(mén)matlab的小白，最近想做反應(yīng)的動(dòng)力學(xué)，按照B站up主的視頻自己寫(xiě)了一段代碼，但是運(yùn)行總是出問(wèn)題：

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

下面是我寫(xiě)的代碼，讀取的Excel表格里有5列*7行的實(shí)驗(yàn)數(shù)據(jù)，勞煩大佬幫我瞅瞅哪里需要改動(dòng)，萬(wàn)分感謝。

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);%動(dòng)力學(xué)計(jì)算得到的點(diǎn)進(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í)驗(yàn)點(diǎn)數(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ù)此樓

» 猜你喜歡

德國(guó)亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生已經(jīng)有0人回復(fù)
推薦一款可以AI輔助寫(xiě)作的Latex編輯器SmartLatexEditor，超級(jí)好用，推薦試試已經(jīng)有11人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有50人回復(fù)
2026-CJ開(kāi)始申報(bào)了已經(jīng)有1人回復(fù)
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士已經(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)有3人回復(fù)
考研交流已經(jīng)有0人回復(fù)

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

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

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

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

★ ★
獨(dú)孤神宇: 金幣+2, 鼓勵(lì)交流 2021-05-20 21:35:31

引用回帖:

4樓: Originally posted by hzd250 at 2021-05-18 10:52:35
這個(gè)是簡(jiǎn)化的數(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的回帖

查看全部 10 個(gè)回答

hzlhm

至尊木蟲(chóng) (著名寫(xiě)手)

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

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

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

贊一下(1人)

回復(fù)此樓

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

hzd250

QQ：2120156492

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

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

hzd250

新蟲(chóng) (小有名氣)

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

送紅花一朵

引用回帖:

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

有的，太謝謝你了

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

贊一下

回復(fù)此樓

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

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

hzd250

新蟲(chóng) (小有名氣)

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

引用回帖:

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

...

這個(gè)是簡(jiǎn)化的數(shù)據(jù)

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

贊一下

回復(fù)此樓

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

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

查看全部 10 個(gè)回答

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

普通表情龍兔虎貓

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

[考研] 080500材料科學(xué)與工程 +5	202114020319 2026-03-03	5/250	2026-03-04 09:33 by gaoxiaoniuma
[考研] 085600求調(diào)劑 +4	LRZZZZZZ 2026-03-02	6/300	2026-03-04 09:11 by 無(wú)際的草原
[考研] 289求調(diào)劑 +8	yang婷 2026-03-02	10/500	2026-03-04 09:10 by 無(wú)際的草原
[考研] 298求調(diào)劑一志愿中海洋 +3	lour. 2026-03-03	3/150	2026-03-03 20:41 by gxg2025
[考研] 理學(xué)，工學(xué)，農(nóng)學(xué)調(diào)劑，少走彎路，這里歡迎您！ +8	likeihood 2026-03-02	11/550	2026-03-03 19:39 by 一葉凡塵?
[考研] 085602化學(xué)工程350，調(diào)劑，有沒(méi)有211的 +5	利好利好. 2026-03-02	9/450	2026-03-03 17:06 by 利好利好.
[考研] 289求調(diào)劑 +7	BrightLL 2026-03-02	9/450	2026-03-03 16:02 by tgxtgxtgx9
[考研] 299求調(diào)劑 +5	kkcoco25 2026-03-02	9/450	2026-03-03 12:55 by 公瑾逍遙
[考研] 材料工程求調(diào)劑 +3	1431251 2026-03-03	3/150	2026-03-03 11:58 by EBSD
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5	awaystay 2026-03-02	6/300	2026-03-03 09:03 by houyaoxu
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5	蟲(chóng)蟲(chóng)蟲(chóng)蟲(chóng)蟲(chóng)7 2026-03-01	9/450	2026-03-02 21:33 by sunny81
[考研] 0856材料求調(diào)劑 +12	hyf hyf hyf 2026-02-28	13/650	2026-03-02 20:19 by hypershenger
[考研] 302材料工程求調(diào)劑 +5	Doleres 2026-03-01	6/300	2026-03-02 19:53 by 張曉芳0105
[考研] 材料與化工328求調(diào)劑 +3	。，。，。，。i 2026-03-02	3/150	2026-03-02 13:09 by houyaoxu
[考研] 274求調(diào)劑 +3	cgyzqwn 2026-03-01	7/350	2026-03-02 10:38 by lature00
[考研] 調(diào)劑 +3	13853210211 2026-03-02	4/200	2026-03-02 10:16 by 13853210211
[基金申請(qǐng)] 本子寫(xiě)完了，給DS兄弟看了，得了92分 +3	Doma 2026-03-01	7/350	2026-03-02 00:00 by jnzsy
[基金申請(qǐng)] 成果系統(tǒng)訪問(wèn)量大，請(qǐng)一小時(shí)后再嘗試。---NSFC啥時(shí)候好哦，已經(jīng)兩天這樣了 +4	NSFC2026我來(lái)了 2026-02-28	4/200	2026-03-01 22:37 by 鐵門(mén)栓
[考研] 化工299分求調(diào)劑一志愿985落榜 +5	嘻嘻(^ω^) 2026-03-01	5/250	2026-03-01 19:47 by 無(wú)際的草原
[考研] 328求調(diào)劑 +3	aaadim 2026-03-01	5/250	2026-03-01 17:29 by njzyff