vs570588

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3112.1
散金: 456
紅花: 1
帖子: 923
在線: 420.6小時(shí)
蟲號: 822119
注冊: 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)在出來NLINFIT did NOT converge. Returning results from last iteration.
beta =

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

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)

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

關(guān)于報(bào)到證改派的問題，請過來人幫幫忙已經(jīng)有4人回復(fù)
幾個(gè)物理化學(xué)問題　請大家?guī)蛶兔?/a> 已經(jīng)有26人回復(fù)

請問VSM磁滯回線的超順磁性問題？高手幫幫忙解答···· 已經(jīng)有15人回復(fù)
關(guān)于配合物固體熒光問題，大家?guī)蛶兔Π?/a> 已經(jīng)有9人回復(fù)

磁控濺射出問題了大蝦們幫幫忙已經(jīng)有6人回復(fù)
【討論】無溶劑加熱遇到了難題，大家?guī)蛶兔?/a> 已經(jīng)有10人回復(fù)

【求助/交流】qPCR問題求助，請大家?guī)蛶兔�，謝謝已經(jīng)有5人回復(fù)
【求助】關(guān)于出國讀博的弱弱問題，請幫幫忙已經(jīng)有6人回復(fù)
【求助】保研聯(lián)系導(dǎo)師的問題，各位幫幫忙吧…… 已經(jīng)有24人回復(fù)
【求助】聚碳酸酯二醇PCDL合成的問題！�。�！大家?guī)蛶兔Π?/a> 已經(jīng)有5人回復(fù)

1樓 2011-05-28 07:58:00

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

文俊點(diǎn)點(diǎn)

木蟲 (著名寫手)

應(yīng)助: 30 (小學(xué)生)
金幣: 2529.5
散金: 1796
紅花: 1
帖子: 1151
在線: 500.1小時(shí)
蟲號: 1219128
注冊: 2011-03-02
性別: GG
專業(yè): 文藝美學(xué)

★
dubo(金幣+1): 謝謝參與 2011-05-30 19:02:30

表示看不懂，短程序應(yīng)該邏輯沒有問題，是細(xì)節(jié)出現(xiàn)錯(cuò)誤。

贊一下(1人)

回復(fù)此樓

這是你所擁有的時(shí)間·····這是你所能改變的生活·········

2樓2011-05-28 17:52:47

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

木蟲pursuer

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1988.3
散金: 3
紅花: 2
帖子: 457
在線: 56.4小時(shí)
蟲號: 1260049
注冊: 2011-04-09
性別: GG
專業(yè): 電化學(xué)

★
dubo(金幣+1): 謝謝參與 2011-05-30 19:03:40

只略懂些c++，頂起來，讓高手進(jìn)

贊一下(1人)

回復(fù)此樓

努力提高自己，努力再努力

3樓2011-05-28 18:19:14

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

vs570588

木蟲 (正式寫手)

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

高手在哪里，大家看是數(shù)據(jù)本身的問題嗎？

贊一下

回復(fù)此樓

4樓2011-05-28 19:32:02

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

wmc_1979

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1232.8
散金: 30
帖子: 128
在線: 110.9小時(shí)
蟲號: 1243434
注冊: 2011-03-24
專業(yè): 運(yùn)籌學(xué)

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

★
jjdg(金幣+1): 感謝參與 2011-05-30 04:22:31

我?guī)湍闼懔艘幌?br /> beta =

0.0249 -14.4864

贊一下(1人)

回復(fù)此樓

5樓2011-05-29 22:01:36

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

vs570588

木蟲 (正式寫手)

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

引用回帖:

Originally posted by wmc_1979 at 2011-05-29 22:01:36:
我?guī)湍闼懔艘幌?br /> beta =

0.0249 -14.4864

謝謝，第二個(gè)參數(shù)是負(fù)的就沒有意義，能在0到1之間取到最佳解嗎

贊一下

回復(fù)此樓

6樓2011-05-30 16:36:11

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

wmc_1979

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1232.8
散金: 30
帖子: 128
在線: 110.9小時(shí)
蟲號: 1243434
注冊: 2011-03-24
專業(yè): 運(yùn)籌學(xué)

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

★
dubo(金幣+1): 謝謝參與 2011-05-30 19:03:56

引用回帖:

Originally posted by vs570588 at 2011-05-30 16:36:11:
謝謝，第二個(gè)參數(shù)是負(fù)的就沒有意義，能在0到1之間取到最佳解嗎

你是想建立y是x的函數(shù)還是x是y的函數(shù)呢?

贊一下(1人)

回復(fù)此樓

7樓2011-05-30 18:27:48

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

dubo

金蟲 (著名寫手)

程序強(qiáng)帖: 4
應(yīng)助: 23 (小學(xué)生)
貴賓: 0.779
金幣: 569.2
散金: 3220
紅花: 31
沙發(fā): 1
帖子: 1821
在線: 349.5小時(shí)
蟲號: 559371
注冊: 2008-05-17
性別: GG
專業(yè): 高分子組裝與超分子結(jié)構(gòu)

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

★ ★
余澤成(金幣+2): 謝謝參與應(yīng)助！ 2011-06-04 19:16:23

用matlab算的吧，
function M=Monod(c,Y)
M= -c(1).*Y./(Y+c(2))
函數(shù)定義不對
nlinfit函數(shù)在Matlab中是如何調(diào)用的，你可以參見
http://zhidao.baidu.com/question/25633638

贊一下(1人)

回復(fù)此樓

8樓2011-05-30 19:01:50

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

vs570588

木蟲 (正式寫手)

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

★
dubo(金幣+1): 謝謝參與 2011-06-13 19:30:40

引用回帖:

Originally posted by dubo at 2011-05-30 19:01:50:
用matlab算的吧，
function M=Monod(c,Y)
M= -c(1).*Y./(Y+c(2))
函數(shù)定義不對
nlinfit函數(shù)在Matlab中是如何調(diào)用的，你可以參見
http://zhidao.baidu.com/question/25633638

%非線性擬合濃度變化過程
t=[0 2 7 9 19 22 24 26 28 30 32 40];%時(shí)間采樣點(diǎn)
s=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31 0];%每時(shí)刻濃度觀測值
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')

這個(gè)是別人幫我寫的，對著吧。但還是出現(xiàn)NLINFIT did NOT converge. Returning results from last iteration。有些人說是我擬合用的數(shù)學(xué)表達(dá)式，不應(yīng)該用雙曲線，但是我的模型確實(shí)是這樣表達(dá)式。難道是我原始數(shù)據(jù)不好。

贊一下(1人)

回復(fù)此樓

9樓2011-05-30 19:13:36

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

dubo

金蟲 (著名寫手)

程序強(qiáng)帖: 4
應(yīng)助: 23 (小學(xué)生)
貴賓: 0.779
金幣: 569.2
散金: 3220
紅花: 31
沙發(fā): 1
帖子: 1821
在線: 349.5小時(shí)
蟲號: 559371
注冊: 2008-05-17
性別: GG
專業(yè): 高分子組裝與超分子結(jié)構(gòu)

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

★ ★ ★
余澤成(金幣+3): 謝謝參與應(yīng)助！ 2011-06-04 19:16:57
vs570588(金幣+2): 謝謝你了 2011-06-05 20:53:54
vs570588(金幣+1): 有勞了 2011-06-09 20:32:03

數(shù)據(jù)2：
t=[0 2 7 9 19 22 24 26 28 30 32 40];%時(shí)間采樣點(diǎn)
s=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31 0];%每時(shí)刻濃度觀測值
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;
plot(x,y,'o');
數(shù)據(jù)點(diǎn)沒規(guī)律

數(shù)據(jù)1：
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;
% t=[0 2 7 9 19 22 24 26 28 30 32 40];%時(shí)間采樣點(diǎn)
% s=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31 0];%每時(shí)刻濃度觀測值
% 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;
plot(x,y,'o');

數(shù)據(jù)像b*(y-a)^2+c=x

贊一下(1人)

回復(fù)此樓

10樓2011-05-30 19:43:12

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

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

[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 化學(xué)工程321分求調(diào)劑 +11	大米飯！ 2026-03-15	14/700	2026-03-17 17:11 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4	小丁想進(jìn)步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant