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

客戶端APP下載

論文輔導(dǎo)

調(diào)劑小程序

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(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)有208人回復(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)有26人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有34人回復(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的回帖

回帖置頂 ( 共有1個 )

yanyl09

禁蟲 (初入文壇)

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

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

已閱關(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的回帖

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ā)表

[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 268求調(diào)劑 +3	簡單點0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 299求調(diào)劑 +4	△小透明* 2026-03-17	4/200	2026-03-17 20:09 by peike
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 274求調(diào)劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎覬?/a>
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky

24小時熱門版塊排行榜

» 猜你喜歡

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

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

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

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

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

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

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

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