| 查看: 1812 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
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不合理。出来的参数应该和我预估的差不多。大家看看，这是怎样回事？ |
回复此楼» 猜你喜欢
申博/考博
已经有8人回复
-
291求调剂
已经有12人回复
-
一志愿A区211，22408 321求调剂
已经有8人回复
-
295分求调剂
已经有6人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有13人回复
-
085600材料与化工调剂
已经有6人回复
-
期刊推荐
已经有5人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有4人回复
-
337求调剂
已经有3人回复
vs570588
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3112.1
- 散金: 456
- 红花: 1
- 帖子: 923
- 在线: 420.6小时
- 虫号: 822119
- 注册: 2009-08-04
- 性别: GG
- 专业: 环境工程
★
dubo(金币+1): 谢谢参与 2011-06-13 19:31:05
dubo(金币+1): 谢谢参与 2011-06-13 19:31:05
|
你好，我用重新参考别人写的程序，用数值解求参数，程序写的很繁琐，你能帮我改改吗？另外，现在运行不下去，提示说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))); |
12楼2011-06-02 15:07:50
文俊点点
木虫 (著名写手)
- 应助: 30 (小学生)
- 金币: 2529.5
- 散金: 1796
- 红花: 1
- 帖子: 1151
- 在线: 500.1小时
- 虫号: 1219128
- 注册: 2011-03-02
- 性别: GG
- 专业: 文艺美学
2楼2011-05-28 17:52:47
木虫pursuer
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1988.3
- 散金: 3
- 红花: 2
- 帖子: 457
- 在线: 56.4小时
- 虫号: 1260049
- 注册: 2011-04-09
- 性别: GG
- 专业: 电化学
3楼2011-05-28 18:19:14
vs570588
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3112.1
- 散金: 456
- 红花: 1
- 帖子: 923
- 在线: 420.6小时
- 虫号: 822119
- 注册: 2009-08-04
- 性别: GG
- 专业: 环境工程
4楼2011-05-28 19:32:02
5