 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 大家?guī)蛶兔o看看吧？
| 1 | 1/1 | 返回列表 |
| 查看: 425 | 回復: 0 | ||
vs570588木蟲 (正式寫手)
|
[求助]
大家?guī)蛶兔o看看吧？
|
|
我把我要擬合求參的問題給你說說。 ds/dt = -q*S*X/(k+S)這里未知參數(shù)是q和K, q是比最大降解速率，K是半飽和常數(shù)，X是污泥濃度213.63，S是污染物的濃度, t肯定就是時間了。我試驗是各一段時間，取一個樣品測出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)境科學中的應用》，上面也用個例子，但是有個例子直接給出了一系列ds/dt的值，并且這些值呈遞增。但你也知道，實際試驗不會出現(xiàn)這種理想情況。所以我求ds/dt值是用多項式擬合，求各個點的導數(shù)，肯定這樣誤差大。但我實在想不出好辦法。也有人說用有限差分法，求出數(shù)值解，再代入，求最優(yōu)化參數(shù)。 這個問題在幾個論壇都發(fā)過貼了，望大家給點處理方法。急死我了，包括改改數(shù)據(jù)。 我也重新參考別人寫的程序，用數(shù)值解求參數(shù)，程序寫的很繁瑣，你能幫我改改嗎？另外，現(xiàn)在運行不下去，提示說divided by zero.你能給看看，怎樣把數(shù)據(jù)處理就能好些？ S=dsolve(‘Dy=-k1*y*214.63/(y+k2)’,’y(0)= 255.55’) simplify(S) %微分方程積分，求出來式子相當繁瑣 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))); |
回復此樓» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有3人回復
-
346求調(diào)劑[0856]
已經(jīng)有6人回復
-
281求調(diào)劑（0805）
已經(jīng)有18人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有9人回復
-
288求調(diào)劑
已經(jīng)有4人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有12人回復
-
一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有5人回復
-
085601材料工程專碩求調(diào)劑
已經(jīng)有7人回復
-
321求調(diào)劑
已經(jīng)有6人回復
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有7人回復
找到一些相關的精華帖子，希望有用哦~
大家?guī)兔纯窗�，求高人指點
已經(jīng)有4人回復
- 3'RACE做了幾次，為什么老是這樣子，大家看一下幫幫忙吧，苦悶死了！~
已經(jīng)有10人回復
- ICASS 2012，征集藝術類的稿件，大家?guī)兔纯窗桑雲(yún)⒓印?/a>
已經(jīng)有6人回復
- 小弟研一新生，導師對我做的PPT不滿意，大家?guī)兔纯�，教教小弟�?/a>
已經(jīng)有20人回復
- 快碩士畢業(yè)了，大家?guī)蛶兔�，看我該如何做出選擇
已經(jīng)有8人回復
- 大家來幫幫忙，急，看看這是什么儀器
已經(jīng)有3人回復
- 大家?guī)兔纯窗车脑u審意見，給些建議吧
已經(jīng)有8人回復
- 關于波生坦中間體的合成，請大家?guī)兔纯窗?/a>
已經(jīng)有14人回復
- 請大家?guī)蛶兔�，看看這CASTEP結構優(yōu)化怎么失敗了
已經(jīng)有11人回復
- 大家?guī)兔纯次业奈恼逻m合投什么雜志呢？給點建議吧！十分感謝！
已經(jīng)有18人回復
- 【求助】外導要推薦信，學校申請時間快到了，大家?guī)兔纯窗?/a>
已經(jīng)有7人回復
- 【求助】外導要面談 大家?guī)蛶兔?看看是幾點
已經(jīng)有3人回復
- 【求助】審稿人意見，不是很懂，大家?guī)蛶兔Γ?/a>
已經(jīng)有21人回復
點擊這里搜索更多相關資源科研從小木蟲開始，人人為我，我為人人
相關版塊跳轉
我要訂閱樓主 學員BmWXvC 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +9 | 陌の森林 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 085601材料工程專碩求調(diào)劑 +7 | 慕寒mio 2026-03-16 | 7/350 |
|
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
|
[考研] 0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +9 | rare12345 2026-03-18 | 9/450 |
|
|
|
[考研] 0703化學 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 0703 物理化學調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|