版塊導(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í)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

shidoudou

新蟲(chóng) (初入文壇)

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

[求助] matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過(guò)程中導(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ù)此樓

» 收錄本帖的淘帖專(zhuān)輯推薦

matlab

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

shidoudou

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有85人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專(zhuān)博名額均有】！已經(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ù)

» 本主題相關(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é)方程問(wèn)題已經(jīng)有7人回復(fù)
【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問(wèn)題已經(jīng)有13人回復(fù)

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

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

shidoudou

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲(chóng)號(hào): 1410410
注冊(cè): 2011-09-21
專(zhuān)業(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的回帖

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

dingd

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

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

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

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

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

贊一下

回復(fù)此樓

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

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

shidoudou

新蟲(chóng) (初入文壇)

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

送鮮花一朵

非常感謝您哈。
我的本來(lái)的目的就是，建立了微生物發(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ù)值，我就不知道怎么樣加語(yǔ)句了。不知道我表達(dá)清楚了沒(méi)有～

贊一下

回復(fù)此樓

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

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

hhucy

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

DOE鍋爐工

應(yīng)助: 53 (初中生)
貴賓: 0.264
金幣: 7775.6
散金: 1861
紅花: 13
沙發(fā): 7
帖子: 2588
在線: 1275.5小時(shí)
蟲(chóng)號(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的回帖

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

普通表情龍兔虎貓

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

[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3	AZMK 2026-03-28	6/300	2026-03-28 17:29 by AZMK
[考研] 283求調(diào)劑 +3	A child 2026-03-28	3/150	2026-03-28 15:41 by ms629
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 081200-314 +3	LILIQQ 2026-03-27	4/200	2026-03-28 09:41 by 保護(hù)地球你我做?/a>
[考研] 材料與化工考研調(diào)劑 +17	孅華 2026-03-22	17/850	2026-03-28 08:35 by WYUMater
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 085600，材料與化工321分調(diào)劑 +4	大饞小子 2026-03-27	6/300	2026-03-27 14:11 by 松花缸1201
[考研] 324求調(diào)劑 +8	hanamiko 2026-03-26	10/500	2026-03-27 08:06 by hypershenger
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17	RichLi_ 2026-03-25	17/850	2026-03-26 19:44 by plmuchong
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 0854電子信息求調(diào)劑 324 +4	Promise-jyl 2026-03-23	4/200	2026-03-25 11:36 by Sugarlight
[考研] 318求調(diào)劑 +5	plum李子 2026-03-21	8/400	2026-03-25 09:26 by aa331100
[考研] 340求調(diào)劑 +5	話(huà)梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte