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

zhl717

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

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

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

有高版本1st opt的大神也可以幫忙跑一下，上回發(fā)了一個(gè)帖子，沒(méi)人回應(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)

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有107人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收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ù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

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

narian

木蟲(chóng) (著名寫手)

總司令

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

引用回帖:

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

高手，私信請(qǐng)教你一些問(wèn)題

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

贊一下

回復(fù)此樓

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

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

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

查看全部 6 個(gè)回答

月只藍(lán)

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

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

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
zhl717: 金幣+5, ★★★很有幫助 2015-08-31 17:03:01

要擬合九個(gè)參數(shù)，這種問(wèn)題還是推薦用1stOpt，參數(shù)個(gè)數(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];

計(jì)算結(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小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

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

zhl717

銅蟲(chóng) (小有名氣)

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

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2015-08-31 16:48:55
要擬合九個(gè)參數(shù)，這種問(wèn)題還是推薦用1stOpt，參數(shù)個(gè)數(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的回帖

皓小天

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

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

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

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

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

贊一下

回復(fù)此樓

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

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

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

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

[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +11	Micky11223 2026-03-25	11/550	2026-03-31 04:05 by fmesaito
[考研] 307求調(diào)劑 +18	超級(jí)伊昂大王 2026-03-24	20/1000	2026-03-31 01:16 by BruceLiu320
[考研] 281求調(diào)劑 +5	亞克西good 2026-03-26	7/350	2026-03-30 20:42 by dophin1985
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 279求調(diào)劑 +12	j的立方 2026-03-29	12/600	2026-03-30 20:30 by dick_runner
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +12	辰沐5211314 2026-03-26	12/600	2026-03-30 19:28 by Wang200018
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +9	XX小鄧 2026-03-29	9/450	2026-03-30 17:18 by limeifeng
[考研] 071010 323 分求調(diào)劑 +3	Baekzhy 2026-03-27	3/150	2026-03-30 14:24 by andresqi
[考研] 考研調(diào)劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25	RichLi_ 2026-03-25	25/1250	2026-03-30 11:02 by klasasda
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無(wú)際的草原
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 0856求調(diào)劑 +13	zhn03 2026-03-25	14/700	2026-03-29 08:13 by fmesaito
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 求調(diào)劑 +7	爭(zhēng)取九點(diǎn)睡 2026-03-28	8/400	2026-03-28 21:07 by 爭(zhēng)取九點(diǎn)睡
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 279 分求調(diào)劑 +4	睡個(gè)好覺(jué)_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風(fēng)里
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-