| 5 | 1/1 | 返回列表 |
| 查看: 1698 | 回復(fù): 15 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
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),這個小程序快把我折磨死了,F(xiàn)在出來NLINFIT did NOT converge. Returning results from last iteration. beta = 0.0271 -8.1892 按道理,-8.1892不合理。出來的參數(shù)應(yīng)該和我預(yù)估的差不多。大家看看,這是怎樣回事? |
木蟲 (正式寫手)
|
%非線性擬合濃度變化過程 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];%每時刻濃度觀測值 x=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31];%找出濃度的初始值,方程右邊的值就用它 i=1 length(t)-1);t1=t(i+1)-t(i);%求t間隔 s1=s(i+1)-s(i);%求s間隔 y=s1./t1; w=213.63;%污泥濃度 y=y./213.63; %-------------------------------------------------------------------------- Monod=inline('-beta(1).*x./(x+beta(2))','beta','x'); c0=[0.03 0.3]; beta=nlinfit(x,y,Monod,c0); m=beta(1); n=beta(2); xx=min(x):max(x); yy=-m.*xx./(xx+n); plot(x,y,'o',xx,yy,'r') 這個是別人幫我寫的,對著吧。但還是出現(xiàn)NLINFIT did NOT converge. Returning results from last iteration。有些人說是我擬合用的數(shù)學(xué)表達式,不應(yīng)該用雙曲線,但是我的模型確實是這樣表達式。難道是我原始數(shù)據(jù)不好。 |
木蟲 (著名寫手)

木蟲 (正式寫手)

木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +9 | 負心者當誅 2026-03-11 | 9/450 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 一志愿西南交大,材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|