应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求高手，拟合求参数
51/1返回列表
查看: 1385  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

vs570588

木虫 (正式写手)

[求助] 求高手，拟合求参数

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);
为了参数c(1),c(2)，这个小程序快把我折磨死了。现在出来NLINFIT did NOT converge. Returning results from last iteration.
beta =

    0.0271
   -8.1892
按道理，-8.1892不合理。出来的参数应该和我预估的差不多。大家看看，这是怎样回事？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-27 20:14:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vs570588

木虫 (正式写手)
引用回帖:
Originally posted by lidaxue at 2011-05-31 10:06:39:
楼主的函数文件里面，好像没有涉及到x的啊，还有你的x=x/214.63;啥意思？

谢谢你回复，这是我的问题，你看看，有啥办法能解决？ds/dt  =  -q*S*X/(k+S)这里未知参数是q和K, q是比最大降解速率，K是半饱和常数，X是污泥浓度214.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值是用多项式拟合，求各个点的导数，肯定这样误差大。但我实在想不出好办法。也有人说用有限差分法，求出数值解，再代入，求最优化参数。
赞一下 回复此楼
3楼2011-05-31 16:42:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

lidaxue

木虫 (正式写手)

之乎者也

【答案】应助回帖

★
sunyang1988(金币+1): 谢谢交流 2011-06-01 18:34:14
楼主的函数文件里面，好像没有涉及到x的啊，还有你的x=x/214.63;啥意思？
赞一下(1人) 回复此楼
Comeon！
2楼2011-05-31 10:06:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lidaxue

木虫 (正式写手)

之乎者也

【答案】应助回帖

★
sunyang1988(金币+1): 谢谢交流 2011-06-01 18:34:23
vs570588(金币+1): 谢谢你了 2011-06-05 20:54:19
楼主你好，看了你的问题，其实不是很难，请楼主参考我给你的ppt，里面有个问题和你的问题比较相似，时间挺紧，还请楼主多多努力！
赞一下(1人) 回复此楼
Comeon！
4楼2011-05-31 18:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vs570588

木虫 (正式写手)
引用回帖:
Originally posted by lidaxue at 2011-05-31 18:47:03:
楼主你好，看了你的问题，其实不是很难，请楼主参考我给你的ppt，里面有个问题和你的问题比较相似，时间挺紧，还请楼主多多努力！

你好，我用你给介绍的，参考别人写的程序，用数值解求参数，程序写的很繁琐，你能帮我改改吗？另外，现在运行不下去，提示说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)));
赞一下 回复此楼
5楼2011-06-02 15:08:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭