lf0612

木蟲 (小有名氣)

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

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

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

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()估計得到的參數(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('計算值','實驗值','Location','best');

figure(2),
plot(ts,ys(:,2),'b',tspan(2:end),yexp(:,2),'or'),legend('計算值','實驗值','Location','best');

figure(3),
plot(ts,ys(:,3),'b',tspan(2:end),yexp(:,3),'or'),legend('計算值','實驗值','Location','best');

figure(4),
plot(ts,ys(:,4),'b',tspan(2:end),yexp(:,4),'or'),legend('計算值','實驗值','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ù)此樓

» 猜你喜歡

中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有135人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有1人回復(fù)

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

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對一方程進行擬合，自己擬合總是出錯！已經(jīng)有12人回復(fù)
matlab 擬合方法已經(jīng)有4人回復(fù)
matlab中參數(shù)改變后怎樣求得曲線？？？已經(jīng)有6人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有9人回復(fù)
使用lsqnonlin函數(shù)優(yōu)化動力學(xué)參數(shù)，總是得不到合理的結(jié)果已經(jīng)有6人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
求助，matlab自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
三維圖譜中，用matlab如何實現(xiàn)去除一部分?jǐn)?shù)據(jù)并實現(xiàn)曲面擬合已經(jīng)有10人回復(fù)
緊急求助，利用Matlab對實驗數(shù)據(jù)進行擬合求解參數(shù)。已經(jīng)有27人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）已經(jīng)有11人回復(fù)
微分方程組參數(shù)擬合的問題已經(jīng)有11人回復(fù)
關(guān)于matlab微分、及非線性擬合問題已經(jīng)有5人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
求高版本1stopt，解決一個問題，微分方程組參數(shù)擬合已經(jīng)有3人回復(fù)
求助擬合參數(shù)的置信區(qū)間問題已經(jīng)有6人回復(fù)
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較已經(jīng)有4人回復(fù)
matlab-常微分方程參數(shù)估計已經(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小時
蟲號: 964589
注冊: 2010-03-08
專業(yè): 計算機軟件

引用回帖:

6樓: Originally posted by dingd at 2014-06-25 15:49:35
公式數(shù)據(jù)對不對�。啃Ч懿畹�。

謝謝了，數(shù)據(jù)有點問題，不過已經(jīng)來不及了，所以最后手工試了幾個參數(shù)，找了個最好的帶進去了！再次感謝！

贊一下

回復(fù)此樓

7樓2014-07-04 16:31:45

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

查看全部 13 個回答

lf0612

木蟲 (小有名氣)

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

原始的微分方程組如下：
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的回帖

月只藍

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

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

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

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

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

贊一下

回復(fù)此樓

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

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小時
蟲號: 964589
注冊: 2010-03-08
專業(yè): 計算機軟件

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 085601材料工程找調(diào)劑 +14	oatmealR 2026-03-29	15/750	2026-03-29 21:48 by plmuchong
[考研] 349求調(diào)劑 +6	李木子啊哈哈 2026-03-25	6/300	2026-03-29 12:47 by 無際的草原
[考研] 332求調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 0703化學(xué) +11	妮妮ninicgb 2026-03-27	11/550	2026-03-29 06:45 by 544594351
[考研] 322求調(diào)劑 +5	舊吢 2026-03-24	5/250	2026-03-28 13:26 by Iveryant
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4	來財aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 求調(diào)劑推薦材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月
[考研] 材料考研調(diào)劑生 +3	黃粱一夢千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪！？！