版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

lf0612

木蟲 (小有名氣)

應助: 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)          % 目標函數(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];

在此感謝！

回復此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復
哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有6人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有203人回復
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復
求調(diào)劑已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復
光學工程學碩調(diào)劑信息已經(jīng)有30人回復
歡迎加入課題組已經(jīng)有0人回復
散金幣，求好運，祝面上順利！已經(jīng)有64人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經(jīng)有1人回復

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

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

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

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

回帖支持 ( 顯示支持度最高的前 50 名 )

月只藍

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

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

【答案】應助回帖

今晚比較忙，今天才著手做，不知道來得及不。程序已經(jīng)調(diào)通

CODE:

function k1k2k3k4k5

format long

clear all

clc

tspan = [0 1 2 3 4 5 6 7 8 9 10];

x0=[0.012 0.005 0.037 0.008];

%k0 = [-0.5 -0.8 -0.8 -1.5 2];  

k0 = [-0.87 -0.9 48.11 -47.80 -0.142];

lb = -[1 1 1 1 1]*1e9;

ub = [1 1 1 1 1]*1e5;

data=[0.073,0.089,0.092,0.316,0.498,1.810,2.314,0.679,0.063,0.044;

    0.084,0.093,0.270,0.406,0.642,0.814,1.026,0.422,0.053,0.014;

    0.638,0.969,2.057,20.419,39.655,69.107,34.099,5.657,0.697,0.346;

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

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)           % 目標函數(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];

參數(shù)比較多，你再調(diào)調(diào)初值k0吧。
多參數(shù)的擬合，還是用1stopt軟件為好，在此推薦專家dingd。

贊一下(1人)

回復此樓

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

5樓2014-06-25 09:55:08

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

普通回帖

lf0612

木蟲 (小有名氣)

應助: 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

贊一下

回復此樓

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

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

月只藍

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

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

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

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

lf0612

木蟲 (小有名氣)

應助: 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ù)的準確度要求不是那么的高，不要差的太多就好了！

贊一下

回復此樓

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

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

dingd

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

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

【答案】應助回帖

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

公式數(shù)據(jù)對不對��？效果很差的。

贊一下

回復此樓

6樓2014-06-25 15:49:35

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

lf0612

木蟲 (小有名氣)

應助: 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ù)，找了個最好的帶進去了！再次感謝！

贊一下

回復此樓

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

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

lf0612

木蟲 (小有名氣)

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

引用回帖:

5樓: Originally posted by 月只藍 at 2014-06-25 09:55:08
今晚比較忙，今天才著手做，不知道來得及不。程序已經(jīng)調(diào)通

function k1k2k3k4k5
format long
clear all
clc
tspan = ;
x0=;
%k0 = ;
k0 = ;
lb = -*1e9;
ub = *1e5;

data=';
size(data);
yexp ...

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

贊一下

回復此樓

8樓2014-07-04 16:31:57

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

lf0612

木蟲 (小有名氣)

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

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

贊一下

回復此樓

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

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

小牧的影子

鐵蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 462.4
散金: 25
帖子: 181
在線: 38.1小時
蟲號: 2650509
注冊: 2013-09-12
性別: GG
專業(yè): 光學

不錯不錯，學習了

回復此樓

學習，學習是立足之本。

10樓2015-06-04 21:47:45

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

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

返回列表

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

[考研] 085700資源與環(huán)境308求調(diào)劑 +8	墨墨漠 2026-03-18	8/400	2026-03-20 08:59 by lbsjt
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +8	步川酷紫123 2026-03-13	8/400	2026-03-20 08:47 by xingguangj
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +3	lemonzzn 2026-03-17	3/150	2026-03-19 23:58 by 23Postgrad
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 274求調(diào)劑 +6	S.H1 2026-03-18	6/300	2026-03-19 09:34 by 花店相見
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 331求調(diào)劑（0703有機化學 +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 考研化學學碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 080500，材料學碩302分求調(diào)劑學校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike

24小時熱門版塊排行榜

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

» 猜你喜歡

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

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

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

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