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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

yanyl09

禁蟲 (初入文壇)

本帖內(nèi)容被屏蔽

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有63人回復(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ù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有27人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有44人回復(fù)

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

origin 8.0 對數(shù)函數(shù)擬合已經(jīng)有13人回復(fù)
請問一個有理函數(shù)的曲線擬合過程如何做，給出數(shù)據(jù)，請求擬合，已經(jīng)有14人回復(fù)
origin8 雙曲函數(shù)擬合一問已經(jīng)有3人回復(fù)
【求助】催化反應(yīng)動力學(xué) 已經(jīng)有5人回復(fù)
【求助】matlab用lsqnonlin非線性擬合的問題【已解決】已經(jīng)有8人回復(fù)
【求助】如何同時擬合兩條曲線或者兩個函數(shù)！【已解決】已經(jīng)有7人回復(fù)
【求助】二元函數(shù)擬合已經(jīng)有12人回復(fù)
【求助】緊急求助非線性擬合問題已經(jīng)有3人回復(fù)
【討論】怎么對數(shù)據(jù)進行特定函數(shù)的擬合已經(jīng)有4人回復(fù)
origin7.5建立函數(shù)擬合為什么會這樣？已經(jīng)有11人回復(fù)
請教關(guān)于origin中自定義函數(shù)非線性回歸擬合問題！已經(jīng)有10人回復(fù)

1樓 2011-12-05 14:54:16

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

回帖支持 ( 顯示支持度最高的前 50 名 )

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

計算結(jié)果

贊一下(2人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

6樓2011-12-06 17:56:32

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

普通回帖

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

微分方程擬合推薦使用1stOpt，更強大方便些。

贊一下

回復(fù)此樓

2樓2011-12-05 17:01:18

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

yanyl09

禁蟲 (初入文壇)

本帖內(nèi)容被屏蔽

3樓2011-12-06 09:02:08

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

yanyl09

禁蟲 (初入文壇)

本帖內(nèi)容被屏蔽

4樓2011-12-06 15:58:49

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

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

yanyl09(金幣+10): 非常感謝，公式的修改對我啟發(fā)很大，我再試試，謝謝�。�！ 2011-12-06 22:37:09

我覺得你的模型不太對
我該了下，你試試，其實第一個生物量還是有問題，你沒有考慮死亡率，估計改后會好些

CODE:

function AA

k0=[0.1 2 2];

x0=[0.01 54.026 0.2435];

t1=0:6:48;

tspan=0:6:48;

yexp=load('sample_file.txt')';

lb=[-10 -10 -10];ub=[10 10 10];

[k,resnorm,resid,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL8,k0,[0 0 0],[],[],x0,yexp);

ci=nlparci(k,resid,jacobian);

y1=[yexp(:,1)]';

y2=[yexp(:,2)]';

y3=[yexp(:,3)]';

[t4plot x4plot]=ode45(@kineticseqs8,[tspan(1) tspan(end)],x0,[],k);

plot(t1,y1,'bo',t1,y2,'g+',t1,y3,'r*',t4plot,x4plot,'k-'),

legend('Exp-biomass','Exp-glucose','Exp-butanol','Model'),

xlabel('time(h)'),ylabel('biomass dry weight(g/l),glucose(g/l),butanol(g/l)')

fprintf('Estimated Parameter\n'),

fprintf('\tk=%.4f %.4f\n',k(1),ci(1,2)-k(1))

fprintf('\tk=%.4f %.4f\n',k(2),ci(2,2)-k(2))

fprintf('\tk=%.4f %.4f\n',k(3),ci(3,2)-k(3))

[t,x]=ode45(@kineticseqs8,tspan,x0,[],k);

figure,plot(x,resid,'*')

xlabel('biomass dry wight(g/l)、glucose(g/l)、butanol(u/ml)'),

ylabel('Residuals(g/l,u/ml)')

function dxdt=kineticseqs8(t,x,k)

% dxdt=[k(1)*(1-x(2)./12)*x(1)

%    -k(2)*x(1).^2

%     k(3)*k(1)*(1-x(2)./12)*x(1)];

dxdt=[k(1)*x(1)*(1-x(1)/3)

   -k(2)*x(1)*x(2)

    k(3)*x(2)*x(1)];

function f=ObjFunc4LNL8(k,x0,yexp)

x0=[0.01 54.026 0.2435];

tspan=[0:6:48]';

[t1 x]=ode45(@kineticseqs8,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 f2 f3];

贊一下

回復(fù)此樓

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

yanyl09

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5樓2011-12-06 17:55:00

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

dingd

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

計算強帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

1stOpt求解代碼：

CODE:

Variable t,x(3);

ODEFunction x1'=k1*x1*(1-x1/3);

            x2'=-k2*x1*x2;

            x3'=k3*x2*x1;

Data;

0        0.01        54.026        0.2435

6        0.37        51.6922        0.2578

12        1.41        45.8798        0.4849

18        2.24        37.4896        1.0704

24        2.59        28.0922        4.3128

30        1.94        15.5371        6.6818

36        1.98        3.6878        9.0772

42        0.89        0.3698        9.4994

48        0.001        0.4067        9.7617

參數(shù) 最佳估算
-------------------- -------------
k1 0.105907029853778
k2 0.622112009922133
k3 0.10538762549862

贊一下

回復(fù)此樓

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

yanyl09

7樓2011-12-06 20:44:02

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

yanyl09

禁蟲 (初入文壇)

送鮮花一朵

本帖內(nèi)容被屏蔽

8樓2011-12-06 22:38:44

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

yanyl09

禁蟲 (初入文壇)

送鮮花一朵

本帖內(nèi)容被屏蔽

9樓2012-03-07 11:24:48

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

邁阿密HEAT

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.5
帖子: 20
在線: 9.6小時
蟲號: 2013394
注冊: 2012-09-19
專業(yè): 運籌學(xué)

引用回帖:

9樓: Originally posted by yanyl09 at 2012-03-07 11:24:48
你好，我改了一下我的模型，考慮了死亡率，模型為
function dxdt=kineticseqs8(t,x,k)
dxdt=;
但是生物量的模擬曲線還是平的，模型改了好幾次模擬曲線都是平的，不知道問題出在哪里，希望版主賜教，謝謝！...

請問輸出來的結(jié)果是什么？是參數(shù)嗎？怎么區(qū)分？

Estimated Parameter
k=0.1059 0.0136
k=0.6221 0.1550
k=0.1054 0.0299
怎么對應(yīng)的？

贊一下

回復(fù)此樓

10樓2014-11-07 22:29:57

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

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

返回列表

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

[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	10/500	2026-03-18 09:21 by 咪嗚喵嗚
[考研] 268求調(diào)劑 +6	簡單點0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao