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

zhl717

銅蟲 (小有名氣)

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

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

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

有高版本1st opt的大神也可以幫忙跑一下，上回發(fā)了一個(gè)帖子，沒人回應(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論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有156人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

matlab擬合反應(yīng)速率常數(shù) 已經(jīng)有5人回復(fù)
Matlab擬合動力學(xué)參數(shù)問題，求指導(dǎo)！已經(jīng)有9人回復(fù)
求助大神，如何用MATLAB擬合微分方程的參數(shù)，關(guān)鍵里面有兩個(gè)自變量用三組數(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ù)的時(shí)滯微分方程程序已經(jīng)有1人回復(fù)
呼叫版主，在線緊急求助，關(guān)于matlab中微分方程組參數(shù)擬合得問題！已經(jīng)有12人回復(fù)
求助1stopt或者matlab 微分方程擬合已經(jīng)有0人回復(fù)
求助，matlab自定義微分方程擬合實(shí)驗(yàn)數(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ì) 已經(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

木蟲 (著名寫手)

總司令

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

引用回帖:

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

function k1k2k3
format long
cle ...

高手，有問題請教

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

回復(fù)此樓

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

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

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

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

月只藍(lán)

主管區(qū)長 (職業(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í)
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

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

要擬合九個(gè)參數(shù)，這種問題還是推薦用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小問題、普通問題請發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

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

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

zhl717

銅蟲 (小有名氣)

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

引用回帖:

好的，非常感謝

回復(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小時(shí)
蟲號: 1864745
注冊: 2012-06-19
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

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

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

我記得已經(jīng)回復(fù)了你這個(gè)問題，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的回帖

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

[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 269專碩求調(diào)劑 +5	金恩貝 2026-03-21	5/250	2026-03-21 22:37 by zhyzzh
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 311求調(diào)劑 +13	冬十三 2026-03-15	14/700	2026-03-21 22:10 by peike
[考研] 化學(xué)工程321分求調(diào)劑 +18	大米飯！ 2026-03-15	22/1100	2026-03-21 20:20 by HH領(lǐng)袖
[基金申請] 山東省面上項(xiàng)目限額評審 +3	石瑞0426 2026-03-19	3/150	2026-03-21 18:59 by RRRKKK
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷