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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

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)          % 目標函數(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年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有159人回復(fù)
哈爾濱理工大學物理系招收物理學考研調(diào)劑已經(jīng)有5人回復(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ù)
光學工程學碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有53人回復(fù)

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

matlab里擬合后，知道y坐標求x坐標已經(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)化動力學參數(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)去除一部分數(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的回帖

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

月只藍

主管區(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è): 宇宙學
管轄: 計算模擬區(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)           % 目標函數(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的回帖

普通回帖

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
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學
管轄: 計算模擬區(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ù)的準確度要求不是那么的高，不要差的太多就好了！

贊一下

回復(fù)此樓

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

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

dingd

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

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

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

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

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

贊一下

回復(fù)此樓

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

已閱回復(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的回帖

lf0612

木蟲 (小有名氣)

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

贊一下

回復(fù)此樓

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

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

lf0612

木蟲 (小有名氣)

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

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

贊一下

回復(fù)此樓

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

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

小牧的影子

鐵蟲 (小有名氣)

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

不錯不錯，學習了

回復(fù)此樓

學習，學習是立足之本。

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

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

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

返回列表

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

[考研] 0703化學調(diào)劑，求各位老師收留 +9	秋有木北 2026-03-14	9/450	2026-03-18 19:49 by macy2011
[考研] 0703化學 305求調(diào)劑 +3	FY_yy 2026-03-14	3/150	2026-03-18 19:40 by macy2011
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 化學工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 0703化學調(diào)劑，六級已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +6	薛云鵬 2026-03-13	6/300	2026-03-18 14:14 by 脫穎而出
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 302求調(diào)劑 +4	小賈同學123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學123
[考研] 285化工學碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] [導(dǎo)師推薦]西南科技大學國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 326求調(diào)劑 +4	諾貝爾化學獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學獎覬?/a>
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 081200-11408-276學碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

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

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

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