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

zhl717

銅蟲 (小有名氣)

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

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

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

有高版本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論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有207人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？已經(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ù)問題，求指導(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é)問題已經(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ù)來(lái)求方程中的參數(shù) 已經(jīng)有4人回復(fù)
請(qǐng)問怎么擬合微分方程中的常數(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的回帖

zhl717

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 332
帖子: 60
在線: 28.7小時(shí)
蟲號(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ù)，這種問題還是推薦用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的回帖

查看全部 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í)
蟲號(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ù)，這種問題還是推薦用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小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

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

皓小天

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

應(yīng)助: 19 (小學(xué)生)
貴賓: 0.1
金幣: 76016.6
散金: 502
紅花: 137
帖子: 3801
在線: 869小時(shí)
蟲號(hào): 1864745
注冊(cè): 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ù)此樓

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

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小時(shí)
蟲號(hào): 1786201
注冊(cè): 2012-04-29
性別: MM
專業(yè): 管理學(xué)

引用回帖:

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

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

回復(fù)此樓

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

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

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

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

[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4	zxts12 2026-03-21	7/350	2026-03-22 09:57 by zxts12
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 323求調(diào)劑 +3	洼小桶 2026-03-18	3/150	2026-03-20 22:54 by JourneyLucky
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518