版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(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
專(zhuān)業(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é)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專(zhuān)博名額均有】！已經(jīng)有1人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有160人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有12人回復(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
專(zhuān)業(yè): 環(huán)境工程

引用回帖:

Originally posted by lidaxue at 2011-05-31 18:47:03:
樓主你好，看了你的問(wèn)題，其實(shí)不是很難，請(qǐng)樓主參考我給你的ppt，里面有個(gè)問(wèn)題和你的問(wèn)題比較相似，時(shí)間挺緊，還請(qǐng)樓主多多努力！

你好，我用你給介紹的，參考別人寫(xiě)的程序，用數(shù)值解求參數(shù)，程序?qū)懙暮芊爆崳隳軒臀腋母膯幔苛硗�，現(xiàn)在運(yùn)行不下去，提示說(shuō)divided by zero.你能給看看，怎樣把數(shù)據(jù)處理就能好些？
S=dsolve(‘Dy=-k1*y*214.63/(y+k2)’,’y(0)= 255.55’)
simplify(S)                                     %微分方程積分，求出來(lái)式子相當(dāng)繁瑣

function monodfit2
clear all;
t= [0 2 7 9 19 22 24 26 28 30 32 40]’;
c=[255.55 246.44 237.28 228.36 136.08 114 99.16 82.33 69.4 56.94 42.31 0]’;
[y_row,y_col]=size(c);
beta0=[0.03,0.3];
c0=255.55;
lb=[0 0];ub=[inf inf];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@seqfun,beta0,lb,ub,[],t,c,y_col,c0);
ci = nlparci(beta,residual,jacobian);
function y = seqfun(beta,t,c,y_col,c0)    % Objective function
tspan = [0  max(t)];
[tt yy] = ode45(@modeleqs,tspan,c0,[],beta);
for col = 1:y_col
yc(:,col) = spline(tt,yy(:,col),t);
end
y=[c(:,1)-yc(:,1)];

function dydt = modeleqs(t,y,beta)    % Model equation
dydt=beta(2)*lambertw(1/beta(2)*exp(-1/100*(21463*t*beta(1)-25555-100*beta(2)*log(19)-100* beta(2)*log(269)+200* beta(2)*log(2)+100* beta(2)*log(5))/ beta(2)));

贊一下

回復(fù)此樓

5樓2011-06-02 15:08:43

已閱回復(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
專(zhuān)業(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
專(zhuān)業(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
專(zhuān)業(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的回帖

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

普通表情龍兔虎貓

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

[考研] 0710生物學(xué)求調(diào)劑！ +4	敘述文 2026-03-31	4/200	2026-03-31 14:32 by 記事本2026
[考研] 070300化學(xué)354求調(diào)劑 +13	101次希望 2026-03-28	13/650	2026-03-31 13:38 by anny19840123
[考研] 266分，求材料相關(guān)專(zhuān)業(yè)調(diào)劑 +10	哇呼哼呼哼 2026-03-30	12/600	2026-03-31 11:00 by 熊一刀
[論文投稿] chinese chemical letters英文版投稿求助 130+3	Yishengeryi 2026-03-30	3/150	2026-03-31 10:10 by 北京萊茵潤(rùn)色
[考研] 物理學(xué)調(diào)劑 +3	小羊36 2026-03-30	3/150	2026-03-31 09:20 by jp9609
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[考研] 材料科學(xué)與工程 317求調(diào)劑 +7	JKSOIID 2026-03-26	7/350	2026-03-30 20:10 by 滴滴上岸呀
[考研] 359求調(diào)劑 +5	王了個(gè)楠 2026-03-25	5/250	2026-03-30 19:36 by 源_2020
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +8	26要上岸 2026-03-28	8/400	2026-03-30 16:39 by herarysara
[碩博家園] 求調(diào)劑有機(jī)化學(xué)考研356分 +10	Nadiums 2026-03-25	11/550	2026-03-30 11:56 by yjolah
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 321求調(diào)劑 +7	璞玉~~ 2026-03-25	8/400	2026-03-29 06:41 by 544594351
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 化學(xué)調(diào)劑 +4	愛(ài)吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 305求調(diào)劑 +5	哇盧卡庫(kù) 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 調(diào)劑求收留 +7	果然有我 2026-03-26	7/350	2026-03-27 00:26 by wxiongid
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947