 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1524 | 回復(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è)小程序快把我折磨死了。現(xiàn)在出來NLINFIT did NOT converge. Returning results from last iteration. beta = 0.0271 -8.1892 按道理，-8.1892不合理。出來的參數(shù)應(yīng)該和我預(yù)估的差不多。大家看看，這是怎樣回事？ |
回復(fù)此樓» 猜你喜歡
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有178人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有79人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 用parstat2273電化學(xué)工作站測正極材料的循環(huán)伏安曲線、交流阻抗的參數(shù)怎么設(shè)置呀？
已經(jīng)有4人回復(fù)
- 高手幫忙擬合一組曲線，解決另加分！
已經(jīng)有19人回復(fù)
- 求助粉晶擬合，歡迎各位高手！
已經(jīng)有19人回復(fù)
- origin 曲線擬合求助，高手請進(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高手 非線性擬合問題
已經(jīng)有6人回復(fù)
- 用1stopt進(jìn)行數(shù)據(jù)擬合時(shí)，無法運(yùn)行sharedmodel ,求助高手
已經(jīng)有7人回復(fù)
- 【求助】Arcgis與GS+在地統(tǒng)計(jì)分析中的區(qū)別
已經(jīng)有39人回復(fù)
lidaxue
木蟲 (正式寫手)
之乎者也
- 應(yīng)助: 10 (幼兒園)
- 金幣: 1952.1
- 散金: 612
- 紅花: 2
- 帖子: 527
- 在線: 856.1小時(shí)
- 蟲號: 259962
- 注冊: 2006-06-17
- 性別: GG
- 專業(yè): 幾何學(xué)
4樓2011-05-31 18:47:03
lidaxue
木蟲 (正式寫手)
之乎者也
- 應(yīng)助: 10 (幼兒園)
- 金幣: 1952.1
- 散金: 612
- 紅花: 2
- 帖子: 527
- 在線: 856.1小時(shí)
- 蟲號: 259962
- 注冊: 2006-06-17
- 性別: GG
- 專業(yè): 幾何學(xué)
2樓2011-05-31 10:06:39
vs570588
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3112.1
- 散金: 456
- 紅花: 1
- 帖子: 923
- 在線: 420.6小時(shí)
- 蟲號: 822119
- 注冊: 2009-08-04
- 性別: GG
- 專業(yè): 環(huán)境工程
|
謝謝你回復(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ù)。 |
3樓2011-05-31 16:42:18
vs570588
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3112.1
- 散金: 456
- 紅花: 1
- 帖子: 923
- 在線: 420.6小時(shí)
- 蟲號: 822119
- 注冊: 2009-08-04
- 性別: GG
- 專業(yè): 環(huán)境工程
|
你好，我用你給介紹的，參考別人寫的程序，用數(shù)值解求參數(shù)，程序?qū)懙暮芊爆�，你能幫我改改嗎？另外，現(xiàn)在運(yùn)行不下去，提示說divided by zero.你能給看看，怎樣把數(shù)據(jù)處理就能好些？ S=dsolve(‘Dy=-k1*y*214.63/(y+k2)’,’y(0)= 255.55’) simplify(S) %微分方程積分，求出來式子相當(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))); |
5樓2011-06-02 15:08:43
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854AI CV方向招收調(diào)劑 +3 | 章小魚567 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6 | 小羊睡著了? 2026-03-23 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +5 | 林之夕 2026-03-24 | 5/250 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|