vs570588

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3112.1
散金: 456
紅花: 1
帖子: 923
在線: 420.6小時
蟲號: 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)，這個小程序快把我折磨死了�，F(xiàn)在出來NLINFIT did NOT converge. Returning results from last iteration.
beta =

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

回復(fù)此樓

» 猜你喜歡

能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)
296求調(diào)劑已經(jīng)有5人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有8人回復(fù)
268求調(diào)劑已經(jīng)有8人回復(fù)
材料，紡織，生物（0856、0710），化學(xué)招生啦已經(jīng)有7人回復(fù)
301求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有7人回復(fù)
材料與化工求調(diào)劑已經(jīng)有6人回復(fù)
299求調(diào)劑已經(jīng)有4人回復(fù)
277調(diào)劑已經(jīng)有6人回復(fù)

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

關(guān)于報到證改派的問題，請過來人幫幫忙已經(jīng)有4人回復(fù)
幾個物理化學(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ī)蛶兔Γx謝已經(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的回帖

文俊點點

木蟲 (著名寫手)

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

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

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

贊一下(1人)

回復(fù)此樓

這是你所擁有的時間·····這是你所能改變的生活·········

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小時
蟲號: 1260049
注冊: 2011-04-09
性別: GG
專業(yè): 電化學(xué)

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

只略懂些c++，頂起來，讓高手進

贊一下(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小時
蟲號: 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小時
蟲號: 1243434
注冊: 2011-03-24
專業(yè): 運籌學(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小時
蟲號: 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

謝謝，第二個參數(shù)是負的就沒有意義，能在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小時
蟲號: 1243434
注冊: 2011-03-24
專業(yè): 運籌學(xué)

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

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

引用回帖:

Originally posted by vs570588 at 2011-05-30 16:36:11:
謝謝，第二個參數(shù)是負的就沒有意義，能在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

金蟲 (著名寫手)

程序強帖: 4
應(yīng)助: 23 (小學(xué)生)
貴賓: 0.779
金幣: 569.2
散金: 3220
紅花: 31
沙發(fā): 1
帖子: 1821
在線: 349.5小時
蟲號: 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小時
蟲號: 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];%時間采樣點
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ù)不好。

贊一下(1人)

回復(fù)此樓

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

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

dubo

金蟲 (著名寫手)

程序強帖: 4
應(yīng)助: 23 (小學(xué)生)
貴賓: 0.779
金幣: 569.2
散金: 3220
紅花: 31
沙發(fā): 1
帖子: 1821
在線: 349.5小時
蟲號: 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];%時間采樣點
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;
plot(x,y,'o');
數(shù)據(jù)點沒規(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];%時間采樣點
% 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;
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ā)表

[考研] 301求調(diào)劑 +5	yy要上岸呀 2026-03-17	5/250	2026-03-17 20:20 by peike
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6	SUICHILD 2026-03-12	6/300	2026-03-17 09:24 by 霧散后相遇lc
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 材料專碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4	一氧二氮 2026-03-11	4/200	2026-03-13 22:35 by JourneyLucky
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 考研調(diào)劑 +4	芬達46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu