| 查看: 521 | 回复: 0 | ||
vs570588木虫 (正式写手)
|
[求助]
大家帮帮忙给看看吧？
|
|
我把我要拟合求参的问题给你说说。 ds/dt = -q*S*X/(k+S)这里未知参数是q和K, q是比最大降解速率，K是半饱和常数，X是污泥浓度213.63，S是污染物的浓度, t肯定就是时间了。我试验是各一段时间，取一个样品测出S,所以我最原始数据是 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]; 就那这一组数据来拟合出上面微分方程里中的未知参数。你看能用啥好办法？另外，我也看宋新山《matlab在环境科学中的应用》，上面也用个例子，但是有个例子直接给出了一系列ds/dt的值，并且这些值呈递增。但你也知道，实际试验不会出现这种理想情况。所以我求ds/dt值是用多项式拟合，求各个点的导数，肯定这样误差大。但我实在想不出好办法。也有人说用有限差分法，求出数值解，再代入，求最优化参数。 这个问题在几个论坛都发过贴了，望大家给点处理方法。急死我了，包括改改数据。 我也重新参考别人写的程序，用数值解求参数，程序写的很繁琐，你能帮我改改吗？另外，现在运行不下去，提示说divided by zero.你能给看看，怎样把数据处理就能好些？ 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))); |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有18人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有22人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有22人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
大家帮忙看看吧，求高人指点
已经有4人回复
- 3'RACE做了几次，为什么老是这样子，大家看一下帮帮忙吧，苦闷死了！~
已经有10人回复
- ICASS 2012，征集艺术类的稿件，大家帮忙看看吧，想参加。
已经有6人回复
- 小弟研一新生，导师对我做的PPT不满意，大家帮忙看看，教教小弟吧
已经有20人回复
- 快硕士毕业了，大家帮帮忙，看我该如何做出选择
已经有8人回复
- 大家来帮帮忙，急，看看这是什么仪器
已经有3人回复
- 大家帮忙看看俺的评审意见，给些建议吧
已经有8人回复
- 关于波生坦中间体的合成，请大家帮忙看看吧
已经有14人回复
- 请大家帮帮忙，看看这CASTEP结构优化怎么失败了
已经有11人回复
- 大家帮忙看看我的文章适合投什么杂志呢？给点建议吧！十分感谢！
已经有18人回复
- 【求助】外导要推荐信，学校申请时间快到了，大家帮忙看看吧
已经有7人回复
- 【求助】外导要面谈 大家帮帮忙 看看是几点
已经有3人回复
- 【求助】审稿人意见，不是很懂，大家帮帮忙！
已经有21人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20