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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

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

hustwzh

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 576.7
散金: 236
紅花: 4
帖子: 204
在線: 129.9小時
蟲號: 1873788
注冊: 2012-06-30
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

[求助] 求大神幫忙解決一下，用MATLAB求解動力學(xué)數(shù)據(jù)總是出錯~ 已有1人參與

如題，本人是個Matlab方面的小白，臨時要用到，但自己糾結(jié)了好久還是沒弄會，只好來向諸位大師求助了，代碼及問題如下，多謝了

CODE:

function KineticsEst5 % 動力學(xué)ODE方程模型的參數(shù)估計

 %

 %

 % 

% The variables y here are y(1)=xB, y(2)=xoNB, y(3)=xmNB,y(4)=xpNB,y(5)=xDNB . 

clear all 

clc 

k0 = [5 5 5 5 5]; % 參數(shù)初值 

lb = [0 0 0 0 0]; % 參數(shù)下限 

ub = [+inf +inf +inf +inf +inf]; % 參數(shù)上限 

x0 = [0 0 0 0 0 0]; 

KineticsData; 

yexp = ExpData(:,2:6); % yexp: 實驗數(shù)據(jù)[x1 x4 x5 x6] 

% 使用函數(shù)fmincon()進(jìn)行參數(shù)估計 

[k,fval,flag] = fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp); 

fprintf('\n使用函數(shù)fmincon()估計得到的參數(shù)值為:\n')

fprintf('\tk1 = %.4f\n',k(1)) 

fprintf('\tk2 = %.4f\n',k(2))

fprintf('\tk3 = %.4f\n',k(3)) 

fprintf('\tk4 = %.4f\n',k(4)) 

fprintf('\tk5 = %.4f\n',k(5)) 

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,[],x0,yexp); 

 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,[],x0,yexp); 

ci = nlparci(k,residual,jacobian); 

fprintf('\n\n以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:\n') 

Output

 % ------------------------------------------------------------------ 

function f = ObjFunc4Fmincon(k,x0,yexp) 

tspan = [0.00 : 5.03 : 15.09]; 

[t x] = ode45(@KineticEqs,tspan,x0,[],k); 

y(:,1) = x(:,1); 

y(:,2:4) = x(:,4:6); 

f = sum((y(:,1)-yexp(:,1)).^2) + sum((y(:,2)-yexp(:,2)).^2) ... 

     + sum((y(:,3)-yexp(:,3)).^2) + sum((y(:,4)-yexp(:,4)).^2); 

% ------------------------------------------------------------------ 

function f = ObjFunc4LNL(k,x0,yexp) 

tspan = [0.00 : 5.03 : 15.09]; 

[t x] = ode45(@KineticEqs,tspan,x0,[],k); 

y(:,1) = x(:,1); 

y(:,2:4) = x(:,4:6); 

f1 = y(:,1) - yexp(:,1); 

f2 = y(:,2) - yexp(:,2); 

f3 = y(:,3) - yexp(:,3); 

f4 = y(:,4) - yexp(:,4); 

f = [f1; f2; f3; f4]; 

% ------------------------------------------------------------------ 

function dxdt = KineticEqs(t,x,k) 

M=1.003; 

dxdt = ... 

   [ ( (k(1)+k(2)+k(3))*(1-x(1))*M*(1-x(2)-x(3)-x(4)-2*x(5)) ) 

     ( M*(1-x(2)-x(3)-x(4)-2*x(5))*(k(1)*(1-x(1))-k(4)*x(2)) ) 

     ( k(2)*(1-x(1))*M*(1-x(2)-x(3)-x(4)-2*x(5)) ) 

     ( M*(1-x(2)-x(3)-x(4)-2*x(5))*(k(3)*(1-x(1))-k(5)*x(4)) ) 

     ( M*(1-x(2)-x(3)-x(4)-2*x(5))*(k(4)*x(2)+k(5)*x(4)) ) 

];

錯誤使用 odearguments (line 92)
KINETICEQS 返回的矢量的長度為 5，但初始條件矢量的長度為 6。KINETICEQS 返回的矢量和初始條件矢量的元素數(shù)目必須相同。

出錯 ode45 (line 113)
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

出錯 KineticsEst>ObjFunc4Fmincon (line 48)
[t x] = ode45(@KineticEqs,tspan,x0,[],k);

出錯 fmincon (line 564)
initVals.f = feval(funfcn{3},X,varargin{:});

出錯 KineticsEst (line 19)
[k,fval,flag] = fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp);

原因:
Failure in initial user-supplied objective function evaluation. FMINCON cannot continue.

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : KineticsData.m

2017-02-18 14:07:42, 498 bytes

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有120人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(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)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有64人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

急�。�！求大神幫忙解決一下泛函方程的建立和壓力函數(shù)的求解已經(jīng)有0人回復(fù)
matlab----解微分方程組，求大神幫忙！！已經(jīng)有6人回復(fù)
求幫忙解一下，萬分感謝�。�！已經(jīng)有6人回復(fù)
求大神幫忙解下結(jié)構(gòu)，萬分感謝�。。。� 已經(jīng)有4人回復(fù)
跪求大神抽空幫忙解一下結(jié)構(gòu)，急急急！已經(jīng)有11人回復(fù)
matlab 中trim函數(shù)總出現(xiàn)Index exceeds matrix dimensions，求大神幫忙解決已經(jīng)有0人回復(fù)
VASP運行出現(xiàn)很多行不該有的文字，求大神幫忙解決！已經(jīng)有1人回復(fù)
模擬扣式電池的問題，求大神幫忙解決！已經(jīng)有6人回復(fù)
eps圖片插入word中遇到了一個很神奇的現(xiàn)象，求大神幫忙解決！已經(jīng)有2人回復(fù)
晶體解析有A類錯誤跪求大神幫忙修一下吧已經(jīng)有4人回復(fù)
求大神幫忙解決幾個A類錯誤已經(jīng)有1人回復(fù)
小弟不才，求大神們解決以下CIF問題已經(jīng)有2人回復(fù)

1樓 2017-02-18 13:59:20

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

hustwzh

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 576.7
散金: 236
紅花: 4
帖子: 204
在線: 129.9小時
蟲號: 1873788
注冊: 2012-06-30
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by chendequan at 2017-02-19 10:02:37
>> KineticsEst5

Local minimum possible. Constraints satisfied.

fmincon stopped because the size of the current step is less than
the default value of the step size tolerance and constrai ...

你好，你可以把修改后的程序發(fā)給我看一下嗎？

贊一下

回復(fù)此樓

5樓2017-02-19 15:45:40

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

查看全部 6 個回答

chendequan

鐵蟲 (小有名氣)

應(yīng)助: 152 (高中生)
金幣: 918.6
紅花: 17
帖子: 262
在線: 208.6小時
蟲號: 994474
注冊: 2010-04-11
專業(yè): 結(jié)構(gòu)工程

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

感謝參與，應(yīng)助指數(shù) +1

>> KineticsEst5

Local minimum possible. Constraints satisfied.

fmincon stopped because the size of the current step is less than
the default value of the step size tolerance and constraints are
satisfied to within the default value of the constraint tolerance.

<stopping criteria details>

使用函數(shù)fmincon()估計得到的參數(shù)值為:
k1 = 0.0358
k2 = 0.0021
k3 = 5.6484
k4 = 22.7533
k5 = 14.7665
The sum of the squares is: 1.5e+00

Local minimum possible.

lsqnonlin stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.

<stopping criteria details>

使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:
k1 = 0.8360
k2 = 0.0128
k3 = 3.0166
k4 = 38.9765
k5 = 7.6220

Local minimum possible.

lsqnonlin stopped because the size of the current step is less than
the default value of the step size tolerance.

<stopping criteria details>

以fmincon()的結(jié)果為初值，使用函數(shù)lsqnonlin()估計得到的參數(shù)值為:
k1 = 0.0358
k2 = 0.0021
k3 = 5.6484
k4 = 22.7533
k5 = 14.7665

贊一下

回復(fù)此樓

QQ：516477448，真心幫助解決MATLAB相關(guān)問題，提供詳細(xì)資料，Word文檔明確具體問題及要求，盡力而為！

2樓2017-02-19 10:02:37

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

hustwzh

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 576.7
散金: 236
紅花: 4
帖子: 204
在線: 129.9小時
蟲號: 1873788
注冊: 2012-06-30
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

你好，你可以把修改后的程序發(fā)給我看一下嗎？

贊一下

回復(fù)此樓

3樓2017-02-19 15:44:27

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

hustwzh

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 576.7
散金: 236
紅花: 4
帖子: 204
在線: 129.9小時
蟲號: 1873788
注冊: 2012-06-30
性別: GG
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

你好，你可以把修改后的程序發(fā)給我看一下嗎？

贊一下

回復(fù)此樓

4樓2017-02-19 15:45:13

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

查看全部 6 個回答

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

[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +5	adaie 2026-03-19	5/250	2026-03-20 10:44 by EBSD
[考研] 085600材料與化工求調(diào)劑 +17	enenenhui 2026-03-13	18/900	2026-03-20 08:55 by xingguangj
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 求調(diào)劑 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、

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

24小時熱門版塊排行榜

hustwzh

[求助] 求大神幫忙解決一下，用MATLAB求解動力學(xué)數(shù)據(jù)總是出錯~ 已有1人參與

» 本帖附件資源列表

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

hustwzh

chendequan

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

hustwzh

hustwzh

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

[求助] 求大神幫忙解決一下，用MATLAB求解動力學(xué)數(shù)據(jù)總是出錯~ 已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: