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

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

» 猜你喜歡

【新加坡】納米電子器件項(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)有161人回復(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é)智能制造與未來能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(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的回帖

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的回帖

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

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的回帖

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的回帖

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

[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +11	Micky11223 2026-03-25	11/550	2026-03-31 04:05 by fmesaito
[考研] 一志愿鄭大材料工程290求調(diào)劑 +12	Youth_ 2026-03-30	12/600	2026-03-31 03:34 by 蒙奇奇521
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +7	真名有冰 2026-03-29	8/400	2026-03-30 21:08 by maaj02
[考研] 286分調(diào)劑 +10	Faune 2026-03-30	11/550	2026-03-30 20:35 by 啊李999
[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +7	蘿卜燉青菜 2026-03-28	8/400	2026-03-30 19:59 by 無際的草原
[考研] 085600 材料與化工 329分求調(diào)劑 +18	Mr. Z 2026-03-25	19/950	2026-03-30 19:01 by Wang200018
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7	Shern—- 2026-03-28	7/350	2026-03-30 16:31 by nothing投稿中
[考研] 考研調(diào)劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 070300化學(xué)354求調(diào)劑 +7	101次希望 2026-03-28	7/350	2026-03-30 13:22 by Iveryant
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7	nhj_ 2026-03-25	7/350	2026-03-30 12:44 by fangnagu
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6	王一一依依 2026-03-30	8/400	2026-03-30 11:52 by yjolah
[考研] 求調(diào)劑 +10	張zz111 2026-03-27	11/550	2026-03-30 09:17 by 無際的草原
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂呵呵的追夢(mèng)人
[考研] 一志愿吉大071010，316分求調(diào)劑 +3	xgbiknn 2026-03-27	3/150	2026-03-27 10:36 by guoweigw
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky