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

lf0612

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1912.8
散金: 100
帖子: 289
在線: 120.1小時(shí)
蟲號(hào): 964589
注冊(cè): 2010-03-08
專業(yè): 計(jì)算機(jī)軟件

[求助] 呼叫版主，在線緊急求助，關(guān)于matlab中微分方程組參數(shù)擬合得問題！已有2人參與

參考版主“月只藍(lán)”以前帖子中的代碼，寫了如下關(guān)于微分方程組中參數(shù)擬合得代碼，進(jìn)行了修改，結(jié)果有問題，因?yàn)槭切率郑a的含義不全懂，加之時(shí)間急著要，所以特來求助，往各位高手多多幫助，在此先謝謝了！

function k1k2k3
format long
clear all
clc
tspan = [0 1 2 3 4 5 6 7 8 9 10 ];
x0=[1 1 1 1 ];
k0 = [1 1 1 1 1 1 1 1 1 1 1 ]*1e-4;
lb = [0 0 0 0 0 0 0 0 0 0 0];
ub = [1 1 1 1 1 1 1 1 1 1 1 ]*1e5;

data=[0.012,0.073,0.089,0.092,0.316,0.498,1.810,2.314,0.679,0.063,0.044;
0.005,0.084,0.093,0.270,0.406,0.642,0.814,1.026,0.422,0.053,0.014;
0.037,0.638,0.969,2.057,20.419,39.655,69.107,34.099,5.657,0.697,0.346;
0.008,0.224,0.234,0.347,1.860,2.854,3.786,1.883,0.601,0.157,0.079;
]';
size(data);
yexp = data;

[k,resnorm,residual,exitflag,output,lambda,jacobian] =lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n')
fprintf('\tk1 = %.16f \n',k(1))
fprintf('\tk2 = %.16f \n',k(2))
fprintf('\tk3 = %.16f \n',k(3))
fprintf('\tk4 = %.16f \n',k(4))
fprintf('\tk5 = %.16f \n',k(5))
fprintf('\tk5 = %.16f \n',k(6))
fprintf('\tk5 = %.16f \n',k(7))
fprintf('\tk5 = %.16f \n',k(8))
fprintf('\tk5 = %.16f \n',k(9))

ts=0:1:max(tspan);

[ts ys]=ode45(@KineticsEqs,ts,x0,[],k);
[ttt XXsim] = ode45(@KineticsEqs,tspan,x0,[],k);

figure(1),
plot(ts,ys(:,1),'b',tspan(2:end),yexp(:,1),'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best');

figure(2),
plot(ts,ys(:,2),'b',tspan(2:end),yexp(:,2),'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best');

figure(3),
plot(ts,ys(:,3),'b',tspan(2:end),yexp(:,3),'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best');

figure(4),
plot(ts,ys(:,4),'b',tspan(2:end),yexp(:,4),'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best');

function f = ObjFunc(k,tspan,x0,yexp)          % 目標(biāo)函數(shù)
[t Xsim] = ode45(@KineticsEqs,tspan,x0,[],k);
Xsim1=Xsim(:,1);
Xsim2=Xsim(:,2);
Xsim3=Xsim(:,3);
Xsim4=Xsim(:,4);
ysim(:,1) = Xsim1(2:end);
ysim(:,2) = Xsim2(2:end);
ysim(:,3) = Xsim3(2:end);
ysim(:,4) = Xsim4(2:end);
size(ysim(:,1));
size(ysim(:,2));
size(ysim(:,3));
size(ysim(:,4));
size(yexp(:,1));
size(yexp(:,2));
size(yexp(:,3));
size(yexp(:,4));
f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3))  (ysim(:,4)-yexp(:,4)) ];

function dEdt = KineticsEqs(t,C,k)
% ODE模型方程 %g(1)=k(4);g(2)=k(5);g(3)=k(6);m(1)=k(7);m(2)=k(8);m(3)=k(9);
CA=C(1);CB=C(2);CC=C(3);CD=C(4);
%dCAdt = k(1)*CA*(CB/2-1)-k(2)*CB+k(3)*CC;
%dCBdt=-g(1)*CA+g(2)*CB.*(1-CB/5)+g(3)*CC.*(1-CC/10);
%dCCdt=m(1)*CA.*(CA/2-1)-m(2)*CB-m(3)*CC;
dCAdt= -k(1)*CA-k(2)*CA;
dCBdt=k(1)*CA-k(3)*CB-k(4)*CB;
dCCdt=k(3)*CB-k(5)*CC;
dCDdt=k(2)*CA+k(4)*CB+k(5)*CC;
dEdt=[dCAdt;dCBdt;dCCdt;dCDdt];

在此感謝！

回復(fù)此樓

» 猜你喜歡

請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有131人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有78人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)

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

matlab里擬合后，知道y坐標(biāo)求x坐標(biāo) 已經(jīng)有4人回復(fù)
Matlab中ode45出現(xiàn)大誤差，關(guān)于用ODE45解偏微分方程的疑問，有待解決，大家共同討論！已經(jīng)有4人回復(fù)
]緊急求助復(fù)介電常數(shù)的模型擬合-Matlab or 1stOpt 已經(jīng)有9人回復(fù)
求助使用matlab對(duì)一方程進(jìn)行擬合，自己擬合總是出錯(cuò)！已經(jīng)有12人回復(fù)
matlab 擬合方法已經(jīng)有4人回復(fù)
matlab中參數(shù)改變后怎樣求得曲線？？？已經(jīng)有6人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有9人回復(fù)
使用lsqnonlin函數(shù)優(yōu)化動(dòng)力學(xué)參數(shù)，總是得不到合理的結(jié)果已經(jīng)有6人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
求助，matlab自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
三維圖譜中，用matlab如何實(shí)現(xiàn)去除一部分?jǐn)?shù)據(jù)并實(shí)現(xiàn)曲面擬合已經(jīng)有10人回復(fù)
緊急求助，利用Matlab對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合求解參數(shù)。已經(jīng)有27人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動(dòng)力系統(tǒng)參數(shù)估計(jì)）已經(jīng)有11人回復(fù)
微分方程組參數(shù)擬合的問題已經(jīng)有11人回復(fù)
關(guān)于matlab微分、及非線性擬合問題已經(jīng)有5人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
求高版本1stopt，解決一個(gè)問題，微分方程組參數(shù)擬合已經(jīng)有3人回復(fù)
求助擬合參數(shù)的置信區(qū)間問題已經(jīng)有6人回復(fù)
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
matlab-常微分方程參數(shù)估計(jì) 已經(jīng)有12人回復(fù)
lsqnonlin函數(shù)擬合微分方程組參數(shù)擬合問題已經(jīng)有10人回復(fù)
【求助】微分方程擬合求參數(shù) 已經(jīng)有7人回復(fù)

1樓 2014-06-24 21:08:21

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

lf0612

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1912.8
散金: 100
帖子: 289
在線: 120.1小時(shí)
蟲號(hào): 964589
注冊(cè): 2010-03-08
專業(yè): 計(jì)算機(jī)軟件

上網(wǎng)不方便，所以現(xiàn)在才回復(fù)二位，不好意思，謝謝了！

贊一下

回復(fù)此樓

9樓2014-07-04 16:32:28

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

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

lf0612

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1912.8
散金: 100
帖子: 289
在線: 120.1小時(shí)
蟲號(hào): 964589
注冊(cè): 2010-03-08
專業(yè): 計(jì)算機(jī)軟件

原始的微分方程組如下：
dS/dt=-k1*S-k2*S
dE/dt=k1*S-k3*E-k4*E
dI/dt=k3*E-k5*I
dR/dt=k2*S+-k4*E+k5*I

贊一下

回復(fù)此樓

2樓2014-06-24 21:17:46

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

月只藍(lán)

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
lf0612: 金幣+10, ★有幫助, 5 2014-07-04 16:29:05

如果你一定要用MATLAB的話，我可以幫你把程序調(diào)通。
不過對(duì)于常微分方程擬合，建議你找一下dingd專家，他有高版本的1stopt軟件，計(jì)算高效強(qiáng)大。也許MATLAB編了半天程序，擬合結(jié)果沒有1stopt花十分鐘計(jì)算的結(jié)果好。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

3樓2014-06-24 21:38:38

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

lf0612

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1912.8
散金: 100
帖子: 289
在線: 120.1小時(shí)
蟲號(hào): 964589
注冊(cè): 2010-03-08
專業(yè): 計(jì)算機(jī)軟件

引用回帖:

3樓: Originally posted by 月只藍(lán) at 2014-06-24 21:38:38
如果你一定要用MATLAB的話，我可以幫你把程序調(diào)通。
不過對(duì)于常微分方程擬合，建議你找一下dingd專家，他有高版本的1stopt軟件，計(jì)算高效強(qiáng)大。也許MATLAB編了半天程序，擬合結(jié)果沒有1stopt花十分鐘計(jì)算的結(jié)果好。

謝謝你了，幫忙把我把程序調(diào)通吧，我對(duì)參數(shù)的準(zhǔn)確度要求不是那么的高，不要差的太多就好了！

贊一下

回復(fù)此樓

4樓2014-06-24 21:47:56

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

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

[考研] 生物學(xué)學(xué)碩求調(diào)劑 +4	小羊睡著了? 2026-03-23	4/200	2026-03-23 21:42 by hyzs6688
[考研] 335分 \| 材料與化工專碩 \| GPA 4.07 \| 有科研經(jīng)歷 +3	cccchenso 2026-03-23	3/150	2026-03-23 19:55 by pswait
[考研] 求調(diào)劑材料學(xué)碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e