24小时热门版块排行榜    

查看: 1669  |  回复: 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不合理。出来的参数应该和我预估的差不多。大家看看,这是怎样回事?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)

【答案】应助回帖

★ ★
vs570588(金币+1): 谢谢 2011-06-10 16:08:17
jjdg(金币+2): 辛苦了 2011-06-11 00:31:56
vs570588(金币+1): 谢谢你了,你看能拟合出来是个直线吗 2011-06-13 14:14:21
引用回帖:
Originally posted by vs570588 at 2011-05-31 16:46:03:
真的很感谢你了,另外我就把我要拟合求参的问题给你说说。
ds/dt  =  -q*S*X/(k+S)这里未知参数是q和K, q是比最大降解速率,K是半饱和常数,X是污泥浓度214.63,这个值是定值。S是污染物的浓度, t肯定就是时间 ...

楼主在最开始给的那些数据不应该用那个函数进行拟合,所以拟合出来有问题。

对于楼主给的这个微分方程,倒是可以直接拟合,但是很奇怪,怎么和楼主的数据不太符合呢,看下图,

楼主给的数据像是凸函数,而直接求解微分方程得到的典型曲线确实凹函数,所以我拟合的结果就不敢放上来了。
14楼2011-06-09 22:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

文俊点点

木虫 (著名写手)


dubo(金币+1): 谢谢参与 2011-05-30 19:02:30
表示看不懂,短程序应该逻辑没有问题,是细节出现错误。
这是你所拥有的时间·····这是你所能改变的生活·········
2楼2011-05-28 17:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木虫pursuer

木虫 (正式写手)


dubo(金币+1): 谢谢参与 2011-05-30 19:03:40
只略懂些c++,顶起来,让高手进
努力提高自己,努力再努力
3楼2011-05-28 18:19:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vs570588

木虫 (正式写手)

高手在哪里,大家看是数据本身的问题吗?
4楼2011-05-28 19:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见