當(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)懂的大神幫忙寫(xiě)一下，謝謝��！

有高版本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)

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有274人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(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的回帖

皓小天

木蟲(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)找我要書(shū)，如果需要請(qǐng)到書(shū)籍板塊求助

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

已閱回復(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
金幣: 68127.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的回帖

narian

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

總司令

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

引用回帖:

高手，有問(wèn)題請(qǐng)教

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

回復(fù)此樓

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

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

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

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

[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +5	幸運(yùn)哩哩 2026-03-22	9/450	2026-03-22 18:34 by 幸運(yùn)哩哩
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	17/850	2026-03-22 15:18 by 脫穎而出
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢(mèng)人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*