zhl717

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 332
帖子: 60
在線: 28.7小時
蟲號: 3705129
注冊: 2015-03-02
性別: GG
專業(yè): 微生物遺傳育種學(xué)

[求助] matlab擬合微分方程組中的參數(shù) 已有2人參與

最近一直在看matlab求解微分方程組中參數(shù)問題，但是由于本人以前沒有接觸過matlab，所以還是不太懂，請懂的大神幫忙寫一下，謝謝��！

有高版本1st opt的大神也可以幫忙跑一下，上回發(fā)了一個帖子，沒人回應(yīng)，謝謝各位��！

dx/dt=a*x*（1-x/b)*(1+s/c)^(-1)
dp/dt=d*dx/dt+e*[s/(s+f)]*x
ds/dt=g*dx/dt+h*dp/dt+i*x

t為自變量，x、p、s為因變量，a-i  均為參數(shù)

數(shù)據(jù)如下：
t, X， P,    S
0 0.35  0    99.8
12  2.75  0.55 88.96
24  4.6 3.72 73.49
36  6.27  8.19 61.57
48  8.07  14.12  47.38
60  9.73 19.37  33.63
72  10.41 22.54  24.65
84  10.7  27.61  11.19
96  10.53 32.49  3.2
108 10.56 35.62  0
120 10.59 37.58  0
132 10.4  38.44  0
144 10.77 39.88  0

回復(fù)此樓

» 本帖@通知

@月只藍(lán)

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有283人回復(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ù)
北工大宋曉艷教授團隊招收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ù)

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

matlab擬合反應(yīng)速率常數(shù) 已經(jīng)有5人回復(fù)
Matlab擬合動力學(xué)參數(shù)問題，求指導(dǎo)！已經(jīng)有9人回復(fù)
求助大神，如何用MATLAB擬合微分方程的參數(shù)，關(guān)鍵里面有兩個自變量用三組數(shù)據(jù)擬合？已經(jīng)有1人回復(fù)
求助各位大神，如何用MATLAB擬合微分方程的參數(shù)？已經(jīng)有0人回復(fù)
MATLAB擬合動力學(xué)問題已經(jīng)有6人回復(fù)
matlab常微分方程多參數(shù)擬合的參數(shù)之間cross correlations 求助@月只藍(lán)以及matlab高手已經(jīng)有8人回復(fù)
求助修改該matlab擬合帶參數(shù)的時滯微分方程程序已經(jīng)有1人回復(fù)
呼叫版主，在線緊急求助，關(guān)于matlab中微分方程組參數(shù)擬合得問題！已經(jīng)有12人回復(fù)
求助1stopt或者matlab 微分方程擬合已經(jīng)有0人回復(fù)
求助，matlab自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
請問怎么擬合微分方程中的常數(shù) 已經(jīng)有3人回復(fù)
關(guān)于matlab微分、及非線性擬合問題已經(jīng)有5人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
matlab參數(shù)擬合已經(jīng)有6人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
matlab-常微分方程參數(shù)估計已經(jīng)有12人回復(fù)
50金幣求一MATLAB的擬合程序已經(jīng)有11人回復(fù)

加油�。�！

1樓 2015-08-31 10:57:30

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

月只藍(lán)

主管區(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
zhl717: 金幣+5, ★★★很有幫助 2015-08-31 17:03:01

要擬合九個參數(shù)，這種問題還是推薦用1stOpt，參數(shù)個數(shù)太多，用MATLAB比較困難。
下面給出MATLAB參考代碼，代碼中初值k0可以試著調(diào)一下，我初步調(diào)了一下初值，擬合效果不好。

CODE:

function k1k2k3

format long

clear all

clc

tspan = 0:12:144;

x0 = [0.35  0      99.8];

k0 = [0.5 8 1.5 2 -0.5 0.1 0.5 5 -0.01];   %k的初值，最需要調(diào)節(jié)

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

ub = [1 1 1 1 1 1 1 1 1]*inf;

data=...

    [

12  2.75  0.55   88.96

24  4.6   3.72   73.49

36  6.27  8.19   61.57

48  8.07  14.12  47.38

60  9.73 19.37  33.63

72  10.41 22.54  24.65

84  10.7  27.61  11.19

96  10.53 32.49  3.2

108 10.56 35.62  0

120 10.59 37.58  0

132 10.4  38.44  0

144 10.77 39.88  0

];

yexp = data(:,2:4);

options=optimset('MaxFunEvals',1500);

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

    lsqnonlin(@ObjFunc,k0,lb,ub,options,tspan,x0,yexp);       

ci = nlparci(k,residual,jacobian);

[t uu]=ode45(@KineticsEqs,tspan,x0,[],k);

uexp=[0   0.35  0      99.8

12  2.75  0.55   88.96

24  4.6   3.72   73.49

36  6.27  8.19   61.57

48  8.07  14.12  47.38

60  9.73 19.37  33.63

72  10.41 22.54  24.65

84  10.7  27.61  11.19

96  10.53 32.49  3.2

108 10.56 35.62  0

120 10.59 37.58  0

132 10.4  38.44  0

144 10.77 39.88  0

];

figure(1)

size(uexp(:,2));

size(uexp(:,3));

size(uexp(:,4));

size(uu(:,1));

size(uu(:,2));

size(uu(:,3));

R1_Square=1-sum((uexp(2:end,2)-uu(2:end,1)).^2)./sum((uexp(2:end,2)-mean( uu(2:end,1)  )).^2);

R2_Square=1-sum((uexp(2:end,3)-uu(2:end,2)).^2)./sum((uexp(2:end,3)-mean( uu(2:end,2)  )).^2);

R3_Square=1-sum((uexp(2:end,4)-uu(2:end,3)).^2)./sum((uexp(2:end,4)-mean( uu(2:end,3)  )).^2);

plot(tspan,uexp(:,2),'or',tspan,uu(:,1),'-r',tspan,uexp(:,3),'>b',tspan,uu(:,2),'b-',tspan,uexp(:,4),'*k',tspan,uu(:,3),'k-')

fprintf('\n\tk1=%.9f',k(1))

fprintf('\n\tk2=%.9f',k(2))

fprintf('\n\tk3=%.9f',k(3))

fprintf('\n\tk4=%.9f',k(4))

fprintf('\n\tk5=%.9f',k(5))

fprintf('\n\tk6=%.9f',k(6))

fprintf('\n\tk7=%.9f',k(7))

fprintf('\n\tk8=%.9f',k(8))

fprintf('\n\tk9=%.9f',k(9))

fprintf('\n\t方程一的決定系數(shù)R1_Square=%.6f',R1_Square);

fprintf('\n\t方程二的決定系數(shù)R2_Square=%.6f',R2_Square);

fprintf('\n\t方程三的決定系數(shù)R3_Square=%.6f',R3_Square);

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(:,2);

ysim(:,1) = Xsim1(2:end);

ysim(:,2) = Xsim2(2:end);

ysim(:,3) = Xsim3(2:end);

f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2)) (ysim(:,3)-yexp(:,3))];

function dCdt = KineticsEqs(t,u,k)  

x=u(1);

p=u(2);

s=u(3);

a=k(1);

b=k(2);

c=k(3);

d=k(4);

e=k(5);

f=k(6);

g=k(7);

h=k(8);

i=k(9);

dxdt = a*x*(1-x/b).*(1+s/c).^(-1);

dpdt=d*dxdt+e*(s./(s+f)).*x;

dsdt=g*dxdt+h*dpdt+i*x;

dCdt = [dxdt;dpdt;dsdt];

計算結(jié)果：
k1=-24.867032689
k2=5.968347685
k3=0.525742903
k4=11.267708389
k5=6.039602224
k6=-1.614606972
k7=-15.309885669
k8=6.596877645
k9=0.664727164
方程一的決定系數(shù)R1_Square=-0.002145
方程二的決定系數(shù)R2_Square=0.359295
方程三的決定系數(shù)R3_Square=-0.035619

附圖1.png

贊一下

回復(fù)此樓

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

2樓2015-08-31 16:48:55

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

zhl717

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 332
帖子: 60
在線: 28.7小時
蟲號: 3705129
注冊: 2015-03-02
性別: GG
專業(yè): 微生物遺傳育種學(xué)

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2015-08-31 16:48:55
要擬合九個參數(shù)，這種問題還是推薦用1stOpt，參數(shù)個數(shù)太多，用MATLAB比較困難。
下面給出MATLAB參考代碼，代碼中初值k0可以試著調(diào)一下，我初步調(diào)了一下初值，擬合效果不好。

function k1k2k3
format long
cle ...

好的，非常感謝

回復(fù)此樓

加油�。�！

3樓2015-08-31 17:03:35

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

皓小天

木蟲之王 (職業(yè)作家)

應(yīng)助: 19 (小學(xué)生)
貴賓: 0.1
金幣: 76016.6
散金: 502
紅花: 137
帖子: 3801
在線: 869小時
蟲號: 1864745
注冊: 2012-06-19
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

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

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

我記得已經(jīng)回復(fù)了你這個問題，http://www.gaoyang168.com/bbs/viewthread.php?tid=9255377&authorid=1864745

贊一下

回復(fù)此樓

請不要站內(nèi)找我要書，如果需要請到書籍板塊求助

4樓2015-08-31 22:03:03

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

narian

木蟲 (著名寫手)

總司令

應(yīng)助: 6 (幼兒園)
金幣: 12.3
散金: 1943
紅花: 62
帖子: 1158
在線: 234.4小時
蟲號: 1786201
注冊: 2012-04-29
性別: MM
專業(yè): 管理學(xué)

引用回帖:

高手，有問題請教

[ 發(fā)自小木蟲客戶端 ]

回復(fù)此樓

回不去的青春，那是我一生一次的認(rèn)真！

5樓2015-09-01 06:30:47

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

narian

木蟲 (著名寫手)

總司令

應(yīng)助: 6 (幼兒園)
金幣: 12.3
散金: 1943
紅花: 62
帖子: 1158
在線: 234.4小時
蟲號: 1786201
注冊: 2012-04-29
性別: MM
專業(yè): 管理學(xué)

引用回帖:

4樓: Originally posted by 皓小天 at 2015-08-31 22:03:03
我記得已經(jīng)回復(fù)了你這個問題，http://www.gaoyang168.com/bbs/viewthread.php?tid=9255377&authorid=1864745

高手，私信請教你一些問題

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

回不去的青春，那是我一生一次的認(rèn)真！

6樓2015-09-01 06:32:16

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

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

[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +8	舉個栗子oi 2026-03-24	9/450	2026-03-29 13:08 by peike
[考研] 0703化學(xué) +11	妮妮ninicgb 2026-03-27	11/550	2026-03-29 06:45 by 544594351
[考研] 332求調(diào)劑 +4	@MZB382400 2026-03-28	4/200	2026-03-28 21:02 by 唐沐兒
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +7	小小甜面團 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 求調(diào)劑 +3	QiMing7 2026-03-25	4/200	2026-03-28 14:30 by QiMing7
[考研] 求化學(xué)調(diào)劑 +4	wulanna 2026-03-28	4/200	2026-03-28 13:37 by 唐沐兒
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4	來財aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 085701環(huán)境工程求調(diào)劑 +9	多久上課 2026-03-27	9/450	2026-03-28 03:58 by fmesaito
[考研] 085600材料與化工調(diào)劑 +10	A-哆啦Z夢 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4	@taotao 2026-03-26	5/250	2026-03-27 08:10 by hypershenger
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 一志愿北化315 求調(diào)劑 +3	akrrain 2026-03-24	3/150	2026-03-24 19:35 by 了了了了。。
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait