版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

shidoudou

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲號(hào): 1410410
注冊(cè): 2011-09-21
專業(yè): 生物化工與食品化工

[求助] matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取

本人matlab不通，勉強(qiáng)做了一個(gè)程序能夠達(dá)到預(yù)期目的，可是老師想要中間運(yùn)算的導(dǎo)數(shù)值，也就是dxdt的值，不知道如何填補(bǔ)表達(dá)式實(shí)現(xiàn)，望各位高手出手相助，不勝感激。
主程序：
format short
global Umax Ks Kp Y1 Y2 Uo Ko D Ux Kx Ku Ki;
Ki=0.0001;Ku=0.0006;
k0=[0.2 0.5 0.1 10 10 0.01 0.01 0.1 0.01 0.01]; %參數(shù)的初始值
x0=[0.0514 0.250 0]; %菌體濃度、底物濃度和產(chǎn)物濃度的初始值
t1=[0 6 12 18 24 30 36 42 48 54 60 66]';%發(fā)酵時(shí)間的實(shí)驗(yàn)數(shù)據(jù)
tspan=[0 6 12 18 24 30 36 42 48 54 60 66 ]';%時(shí)間間隔
%yexp實(shí)驗(yàn)數(shù)據(jù)[菌體濃度x1、底物濃度x2和產(chǎn)物濃度x3]
yexp=[[0.0514 0.0711 0.0997 0.1887 0.3478 0.6634 0.8711 1.1023 1.2336 1.4239 1.5849 1.6739];[0.25 0.2447 0.2366 0.2258 0.2196 0.203 0.1872 0.1698 0.1453 0.1124 0.0855 0.0608];[0 0.0012 0.0021 0.0039 0.0088 0.0137 0.0183 0.0268 0.0361 0.0438 0.0557 0.0643]]';
lb=[0 0 0 0 0 0.0216 0.001 0 0.001 0.001];ub=[1 10 10 50 50 2 1 1 0.1 0.1];%參數(shù)的下、上限
[k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL103,k0,lb,ub,[],x0,yexp);%非線性擬合
y1=[yexp(:,1)];y2=[yexp(:,2)];y3=[yexp(:,3)]';
[t4plot,x4plot]=ode45(@kineticseqs103,[0 100],x0,[],k);
plot(t1,y1,'b*',t4plot,x4plot(:,1),'k-'),xlabel('T(h)'),ylabel('Cx/（g/L）');
figure
plot(t1,y2,'g*',t4plot,x4plot(:,2),'k-'),xlabel('T(h)'),ylabel('Cs/(mmol/L)');
figure
plot(t1,y3,'r*',t4plot,x4plot(:,3),'k-'),xlabel('T(h)'),ylabel('Cp/(mmol/L)');
disp(k)
子程序：
1）
function f=ObjFunc4LNL103(k,x0,yexp)
tspan=[0 6 12 18 24 30 36 42 48 54 60 66];
[t1,x]=ode45(@kineticseqs103,tspan,x0,[],k);
y(:,1)=x(:,1);y(:,2)=x(:,2);y(:,3)=x(:,3);
f1=y(: ,1)-yexp(: ,1);f2=y(: ,2)-yexp(: ,2);f3=y(: ,3)-yexp(: ,3);
f=[f1*1 f2*5 f3*10];
2）
function dxdt=kineticseqs103(t,x,k) %模型方程
global Umax Ks Kp Y1 Y2 Uo Ko D Ux Kx Ku Ki
Umax=k(1); Ks=k(2); Kp=k(3); Y1=k(4); Y2=k(5); Ux=k(6); Kx=k(7); D=k(8); Uo=k(9); Ko=k(10);
  dxdt1=Umax*x(2)*x(1)*(1-Kp*x(3))/(Ks+x(2));
if dxdt1<0
   dxdt1=0;
end
  dxdt2=-1/Y1*dxdt1-Ux*x(2)*x(1)/(Kx+x(2))*(Ki/(Ki+x(3)));
  dxdt3=D*(-dxdt2-1/Y2*dxdt1-Uo*x(2)*x(1)/(Ko+x(2))*(Ku/(Ku+x(3))));
if dxdt3<0
   dxdt3=0;
  end
  dxdt=[dxdt1;dxdt2;dxdt3];
只要能顯示對(duì)應(yīng)時(shí)間t1,這些點(diǎn)的dxdt值就可以了。謝謝

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

matlab

» 本帖已獲得的紅花（最新10朵）

shidoudou

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有281人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(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)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有52人回復(fù)

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

matlab擬合擬均相動(dòng)力學(xué)參數(shù) 已經(jīng)有6人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
吸附的動(dòng)力學(xué)方程擬合已經(jīng)有18人回復(fù)
matlab可否求離散三維函數(shù)（即一組（x，y，z）數(shù)據(jù)）的偏導(dǎo)數(shù)（如：Dz/Dy）？已經(jīng)有5人回復(fù)
【求助】擬合動(dòng)力學(xué)方程求助已經(jīng)有13人回復(fù)
origin，sigmaplot等軟件自定義曲線擬合（主要是一些動(dòng)力學(xué)方程的模擬）已經(jīng)有2人回復(fù)
【求助】光催化動(dòng)力學(xué)方程擬合已經(jīng)有7人回復(fù)
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問題已經(jīng)有7人回復(fù)
【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問題已經(jīng)有13人回復(fù)

1樓 2012-02-09 14:46:40

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

感謝參與，應(yīng)助指數(shù) +1
dbb627(金幣+1): 歡迎交流 2012-02-10 18:25:29

最好用非Matlab代碼的文本（圖片）方式把問題描述一下。

贊一下

回復(fù)此樓

2樓2012-02-09 17:16:45

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

shidoudou

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲號(hào): 1410410
注冊(cè): 2011-09-21
專業(yè): 生物化工與食品化工

送鮮花一朵

非常感謝您哈。
我的本來的目的就是，建立了微生物發(fā)酵的動(dòng)力學(xué)方程，是三個(gè)微分方程，然后利用實(shí)驗(yàn)所得的數(shù)據(jù)就是yexp那組數(shù)據(jù)，非線性擬合微分方程里面的參數(shù)即可，以上的程序可以實(shí)現(xiàn)這個(gè)了�？墒牵蠋熥屛夷軌蝻@示這三個(gè)微分方程在對(duì)應(yīng)的時(shí)間點(diǎn)t1那幾個(gè)點(diǎn)處的導(dǎo)數(shù)值，我就不知道怎么樣加語句了。不知道我表達(dá)清楚了沒有～

贊一下

回復(fù)此樓

3樓2012-02-09 20:30:07

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

hhucy

實(shí)習(xí)版主 (著名寫手)

DOE鍋爐工

應(yīng)助: 53 (初中生)
貴賓: 0.264
金幣: 7775.6
散金: 1861
紅花: 13
沙發(fā): 7
帖子: 2588
在線: 1275.5小時(shí)
蟲號(hào): 1079650
注冊(cè): 2010-08-22
性別: GG
管轄: 能源

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

感謝參與，應(yīng)助指數(shù) +1
dbb627(金幣+1): 歡迎交流 2012-02-10 18:25:37

添加作圖的命令，就能直觀顯示

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

shidoudou

人生那么多不確定，你怕什么

4樓2012-02-09 20:40:21

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

shidoudou

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲號(hào): 1410410
注冊(cè): 2011-09-21
專業(yè): 生物化工與食品化工

送鮮花一朵

引用回帖:

樓: Originally posted by hhucy at 2012-02-09 20:40:21:
添加作圖的命令，就能直觀顯示

感謝哈。
能不能詳細(xì)點(diǎn)致命咧，呵呵，我對(duì)此還真不通呢，可是老師說想要數(shù)值看一看，還得交差哈

贊一下

回復(fù)此樓

5樓2012-02-09 20:53:56

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

change0618

鐵桿木蟲 (著名寫手)

方丈大師

計(jì)算強(qiáng)帖: 1
應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517
注冊(cè): 2008-01-19
專業(yè): 化學(xué)反應(yīng)工程

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

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

dxdt=kineticseqs103(t,x,k) %模型方程

這不就是濃度梯度啊？

贊一下

回復(fù)此樓

6樓2012-02-10 15:58:39

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

change0618

鐵桿木蟲 (著名寫手)

方丈大師

計(jì)算強(qiáng)帖: 1
應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時(shí)
蟲號(hào): 496517
注冊(cè): 2008-01-19
專業(yè): 化學(xué)反應(yīng)工程

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

dbb627:編輯內(nèi)容 2012-02-10 18:24
dbb627(金幣+2): 感謝應(yīng)助 2012-02-10 18:25:08

CODE:

[t,x]=ode45(@kineticseqs103,t1,x0,[],k);

n = length(x);

dxdt = zeros(size(x));

for i = 1:n

    dxdt(i,:) = kineticseqs103(t,x(i,:),k);

end

dxdt

[ Last edited by dbb627 on 2012-2-10 at 18:24 ]

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

shidoudou

7樓2012-02-10 16:30:07

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

715211229

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 3652.8
散金: 20
紅花: 6
帖子: 339
在線: 249.2小時(shí)
蟲號(hào): 1357002
注冊(cè): 2011-07-27
性別: GG
專業(yè): 傳動(dòng)機(jī)械學(xué)

matlab命令窗口輸入cftool，這個(gè)是擬合工具箱，擬合好后，最后一項(xiàng)分析里面有求一二階導(dǎo)倒數(shù)的按鈕

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

shidoudou

我是蝸牛

8樓2012-02-11 00:49:28

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

shidoudou

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲號(hào): 1410410
注冊(cè): 2011-09-21
專業(yè): 生物化工與食品化工

送鮮花一朵

引用回帖:

樓: Originally posted by change0618 at 2012-02-10 16:30:07:

CODE:

[t,x]=ode45(@kineticseqs103,t1,x0,[],k);

n = length(x);

dxdt = zeros(size(x));

for i = 1:n

    dxdt(i,:) = kineticseqs103(t,x(i,:),k);

end

dxdt

[ Last edited by dbb627 on ...

謝謝哈

回復(fù)此樓

9樓2012-02-12 09:08:14

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

shidoudou

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲號(hào): 1410410
注冊(cè): 2011-09-21
專業(yè): 生物化工與食品化工

送鮮花一朵

引用回帖:

樓: Originally posted by 715211229 at 2012-02-11 00:49:28:
matlab命令窗口輸入cftool，這個(gè)是擬合工具箱，擬合好后，最后一項(xiàng)分析里面有求一二階導(dǎo)倒數(shù)的按鈕

謝謝哈

回復(fù)此樓

10樓2012-02-12 09:08:30

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 shidoudou 的主題更新

返回列表

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

[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao

24小時(shí)熱門版塊排行榜

[求助] matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取

» 收錄本帖的淘帖專輯推薦

» 本帖已獲得的紅花（最新10朵）

» 猜你喜歡

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

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

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

» 本帖已獲得的紅花（最新10朵）

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

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

» 本帖已獲得的紅花（最新10朵）

» 本帖已獲得的紅花（最新10朵）

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