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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

lf0612

木蟲 (小有名氣)

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

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

參考版主“月只藍(lán)”以前帖子中的代碼，寫了如下關(guān)于微分方程組中參數(shù)擬合得代碼，進(jìn)行了修改，結(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ù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有271人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊招收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ù)

» 本主題相關(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ìn)行擬合，自己擬合總是出錯！已經(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ù)進(jìn)行擬合求解參數(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ù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(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的回帖

月只藍(lán)

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

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.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的回帖

查看全部 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的回帖

lf0612

木蟲 (小有名氣)

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

引用回帖:

3樓: Originally posted by 月只藍(lán) 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的回帖

月只藍(lán)

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

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

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

今晚比較忙，今天才著手做，不知道來得及不。程序已經(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)           % 目標(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];

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

贊一下(1人)

回復(fù)此樓

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

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

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

查看全部 13 個回答

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

[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	13/650	2026-03-27 22:42 by Micky11223
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 一志愿211院校 344分東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5	丶風(fēng)雪夜歸人丶 2026-03-26	8/400	2026-03-27 19:22 by 丶風(fēng)雪夜歸人丶
[考研] 291求調(diào)劑 +12	hhhhxn.. 2026-03-23	18/900	2026-03-27 17:41 by ditto77778
[考研] 308求調(diào)劑 +6	墨墨漠 2026-03-27	6/300	2026-03-27 15:20 by caszguilin
[考研] 調(diào)劑推薦 +5	清酒714 2026-03-26	6/300	2026-03-27 11:12 by 不吃魚的貓
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +4	舉個栗子oi 2026-03-24	5/250	2026-03-26 23:15 by 不吃魚的貓
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4	柳溪邊 2026-03-26	4/200	2026-03-26 21:14 by XPU李慶
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6	jiajunX 2026-03-22	6/300	2026-03-25 23:05 by licg0208
[考研] 一志愿南航 335分 \| 0856材料化工 \| GPA 4.07 \| 有科研經(jīng)歷 +6	cccchenso 2026-03-23	6/300	2026-03-25 22:25 by 544594351
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 化學(xué)調(diào)劑 +6	yzysaa 2026-03-21	6/300	2026-03-25 09:27 by aa331100
[考研] B區(qū)考研調(diào)劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[有機交流] 有機合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch
[考研] 求助 +5	夢里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO

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

24小時熱門版塊排行榜

lf0612

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

» 猜你喜歡

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

月只藍(lán)

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

lf0612

lf0612

月只藍(lán)

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

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

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