| 5 | 1/1 | 返回列表 |
| 查看: 1528 | 回復(fù): 13 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
vs570588木蟲 (正式寫手)
|
[求助]
求高手,擬合求參數(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)在出來NLINFIT did NOT converge. Returning results from last iteration. beta = 0.0271 -8.1892 按道理,-8.1892不合理。出來的參數(shù)應(yīng)該和我預(yù)估的差不多。大家看看,這是怎樣回事? |
榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +4 |

木蟲 (正式寫手)
之乎者也

木蟲 (正式寫手)
|
謝謝你回復(fù),這是我的問題,你看看,有啥辦法能解決?ds/dt = -q*S*X/(k+S)這里未知參數(shù)是q和K, q是比最大降解速率,K是半飽和常數(shù),X是污泥濃度214.63,這個(gè)值是定值。S是污染物的濃度, t肯定就是時(shí)間了。我具體試驗(yàn)是隔一段時(shí)間,取一個(gè)樣品測出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ù)來擬合出上面微分方程里中的未知參數(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í)在想不出好辦法。也有人說用有限差分法,求出數(shù)值解,再代入,求最優(yōu)化參數(shù)。 |
木蟲 (正式寫手)
之乎者也

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華北電力大學(xué)能動(dòng)專碩,293,求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +7 | Li李魚 2026-03-26 | 7/350 |
|
|
[考研] 求化學(xué)調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 11/550 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 一志愿國科過程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|