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

vs570588

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 3112.1
散金: 456
紅花: 1
帖子: 923
在線: 420.6小時(shí)
蟲(chóng)號(hào): 822119
注冊(cè): 2009-08-04
性別: GG
專業(yè): 環(huán)境工程

[求助] 求高手，擬合求參數(shù)

function M=Monod(c,Y)
M= -c(1).*Y./(Y+c(2))

Y=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31 0];
x=[-0.78 -2.2268 -5.2033 -6.1377 -8.6137 -8.6428 -8.4792 -8.1692 -7.7128 -7.11 -6.3608 -1.9];
x=x/214.63;
c0=[0.03 0.3];beta=nlinfit(Y,x ,’Monod’,c0);
為了參數(shù)c(1),c(2)，這個(gè)小程序快把我折磨死了�，F(xiàn)在出來(lái)NLINFIT did NOT converge. Returning results from last iteration.
beta =

0.0271
-8.1892
按道理，-8.1892不合理。出來(lái)的參數(shù)應(yīng)該和我預(yù)估的差不多。大家看看，這是怎樣回事？

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有150人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

用parstat2273電化學(xué)工作站測(cè)正極材料的循環(huán)伏安曲線、交流阻抗的參數(shù)怎么設(shè)置呀？已經(jīng)有4人回復(fù)
高手幫忙擬合一組曲線，解決另加分！已經(jīng)有19人回復(fù)
求助粉晶擬合，歡迎各位高手！已經(jīng)有19人回復(fù)
origin 曲線擬合求助，高手請(qǐng)進(jìn) 已經(jīng)有11人回復(fù)
利用origin 擬合曲線求參數(shù)，程序老出現(xiàn)錯(cuò)誤，急死了，求各位高手指點(diǎn)！已經(jīng)有7人回復(fù)
origin 擬合公式求助已經(jīng)有4人回復(fù)
幫幫忙看看我的難題已經(jīng)有15人回復(fù)
求助origin高手非線性擬合問(wèn)題已經(jīng)有6人回復(fù)
用1stopt進(jìn)行數(shù)據(jù)擬合時(shí)，無(wú)法運(yùn)行sharedmodel ,求助高手已經(jīng)有7人回復(fù)
【求助】Arcgis與GS+在地統(tǒng)計(jì)分析中的區(qū)別已經(jīng)有39人回復(fù)

1樓 2011-05-27 20:14:00

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

vs570588

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 3112.1
散金: 456
紅花: 1
帖子: 923
在線: 420.6小時(shí)
蟲(chóng)號(hào): 822119
注冊(cè): 2009-08-04
性別: GG
專業(yè): 環(huán)境工程

引用回帖:

Originally posted by dbb627 at 2011-06-02 21:50:44:
計(jì)算的圖見(jiàn)附件

你好，謝謝你了。你用maple做出的解析解。我matlab是個(gè)菜鳥(niǎo)，怎樣用matlab實(shí)現(xiàn)，我不會(huì)。你能幫我看看嗎？另外，這里還有幾篇英文文獻(xiàn)，處理類似的問(wèn)題。
第一篇：
Non-linear least-square error minimization was used to estimate best-fit values for the kinetic parameters (Sa´ez and Rittmann, 1992). In this technique, the modeling equations are solved numerically, and parameters are selected to minimize the sum of the relative least-square residuals. The equations were solved by finite differences in a Microsoft Excel spreadsheet.
第二篇：
2.7.5. Data fitting
AQUASIM version 2.1f (Reichert, 1995) was used to fit kinetic parameters. AQUASIM estimates kinetic parameters by minimizing
the sum of the squares of the weighted deviations between actual data and results of the calculated model. The calculation step size was 0.01 days. The secant method was used with a maximum iteration number of 100.
第三篇：
The fitting method adopted is detailed as follows. Eqs. (1) and（8) were solved numerically by the finite-difference method with
a finite-difference of dt= 0.00625 h and with an initial guess of qmax and Ks. The optimal qmax and Ks were then obtained by
changing their values in Microsoft Excel Solver to reach the minimum SSE between the model-calculated and observed data.
我希望我出來(lái)的擬合圖形不應(yīng)該是折線圖，而是如附件所示。我的數(shù)據(jù)和第三篇中的數(shù)據(jù)很像，也是一式三組試驗(yàn)，其中一組和另外兩組數(shù)據(jù)有區(qū)別。第三篇擬合求參方法如上所示。
最后，還是要謝謝你。

贊一下

回復(fù)此樓

9樓2011-06-04 18:32:27

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

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

lidaxue

木蟲(chóng) (正式寫(xiě)手)

之乎者也

應(yīng)助: 10 (幼兒園)
金幣: 1952.1
散金: 612
紅花: 2
帖子: 527
在線: 856.1小時(shí)
蟲(chóng)號(hào): 259962
注冊(cè): 2006-06-17
性別: GG
專業(yè): 幾何學(xué)

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

★
sunyang1988(金幣+1): 謝謝交流 2011-06-01 18:34:14

樓主的函數(shù)文件里面，好像沒(méi)有涉及到x的啊，還有你的x=x/214.63;啥意思？

贊一下(1人)

回復(fù)此樓

Comeon！

2樓2011-05-31 10:06:39

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

vs570588

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 3112.1
散金: 456
紅花: 1
帖子: 923
在線: 420.6小時(shí)
蟲(chóng)號(hào): 822119
注冊(cè): 2009-08-04
性別: GG
專業(yè): 環(huán)境工程

引用回帖:

Originally posted by lidaxue at 2011-05-31 10:06:39:
樓主的函數(shù)文件里面，好像沒(méi)有涉及到x的啊，還有你的x=x/214.63;啥意思？

謝謝你回復(fù)，這是我的問(wèn)題，你看看，有啥辦法能解決？ds/dt = -q*S*X/(k+S)這里未知參數(shù)是q和K, q是比最大降解速率，K是半飽和常數(shù)，X是污泥濃度214.63，這個(gè)值是定值。S是污染物的濃度, t肯定就是時(shí)間了。我具體試驗(yàn)是隔一段時(shí)間，取一個(gè)樣品測(cè)出S,所以我最原始數(shù)據(jù)是
t=[0 2 7 9 19 22 24 26 28 30 32 40];
S=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31 0];
就那這一組數(shù)據(jù)來(lái)擬合出上面微分方程里中的未知參數(shù)。你看能用啥好辦法？另外，我也看宋新山《matlab在環(huán)境科學(xué)中的應(yīng)用》，上面也用個(gè)例子，但是有個(gè)例子直接給出了一系列ds/dt的值，并且這些值呈遞增。但你也知道，實(shí)際試驗(yàn)不會(huì)出現(xiàn)這種理想情況。所以我求ds/dt值是用多項(xiàng)式擬合，求各個(gè)點(diǎn)的導(dǎo)數(shù)，肯定這樣誤差大。但我實(shí)在想不出好辦法。也有人說(shuō)用有限差分法，求出數(shù)值解，再代入，求最優(yōu)化參數(shù)。

贊一下

回復(fù)此樓

3樓2011-05-31 16:42:18

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

lidaxue

木蟲(chóng) (正式寫(xiě)手)

之乎者也

應(yīng)助: 10 (幼兒園)
金幣: 1952.1
散金: 612
紅花: 2
帖子: 527
在線: 856.1小時(shí)
蟲(chóng)號(hào): 259962
注冊(cè): 2006-06-17
性別: GG
專業(yè): 幾何學(xué)

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

★
sunyang1988(金幣+1): 謝謝交流 2011-06-01 18:34:23
vs570588(金幣+1): 謝謝你了 2011-06-05 20:54:19

樓主你好，看了你的問(wèn)題，其實(shí)不是很難，請(qǐng)樓主參考我給你的ppt，里面有個(gè)問(wèn)題和你的問(wèn)題比較相似，時(shí)間挺緊，還請(qǐng)樓主多多努力！

贊一下(1人)

回復(fù)此樓

Comeon！

4樓2011-05-31 18:47:03

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

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

[考研] 求調(diào)劑院校信息 +4	CX 330 2026-03-21	4/200	2026-03-21 23:48 by ms629
[考研] 354求調(diào)劑 +6	Tyoumou 2026-03-18	9/450	2026-03-21 20:47 by lbsjt
[考研] 化學(xué)工程321分求調(diào)劑 +18	大米飯！ 2026-03-15	22/1100	2026-03-21 20:20 by HH領(lǐng)袖
[考研] 326求調(diào)劑 +5	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	8/400	2026-03-21 19:33 by ColorlessPI
[考研] 工科0856求調(diào)劑 +3	沐析汀汀 2026-03-21	3/150	2026-03-21 18:30 by 學(xué)員8dgXkO
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +11	吃吃吃才有意義 2026-03-19	11/550	2026-03-21 18:23 by 學(xué)員8dgXkO
[考研] 317求調(diào)劑 +9	申子申申 2026-03-19	15/750	2026-03-21 17:31 by 學(xué)員8dgXkO
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中國(guó)石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 324求調(diào)劑 +5	lucky呀呀呀鴨 2026-03-20	5/250	2026-03-20 22:30 by 促天成
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz