2013-02-18 36 views
1

破坏这个公式当指数曲线拟合噪声

y=a*exp(-x/b) 

x被固定x=[13 26 39 52 65 78 91]做曲线拟合时,我遇到了一个问题。 y是输入。 a和b不知道。 b是输出。我使用LSQ估计来进行曲线拟合,并为输出b:b添加一个约束,范围为[0,1000]

目前,该系统的工作原理是这样的:当我有这样

y=[460 434 288 218 164 114 89] 

输入序列输出为b=51.46,这是很好的。

如果输入序列是

y=[599 640 592 609 550 588 573 626] 

的估计结果是b=1000。这也很好。没问题。

但是,当我输入一个纯粹的噪声序列:

y=[24 19 31 5 27 31 17] 

结果我从我的曲线拟合算法得到的是b=1000。在这种情况下,输出b是一个非常高的信号,这对于系统来说是不可接受的。我预计输出一个低值b,比如b = 0

我试图在Y添加一个门槛,说

if y<50 then b=0 

但是该系统不是很稳定。噪音水平会随时变化。还有其他方法可以解决这个问题吗?先谢谢你。

回答

2

首先,请注意,这类问题通常以logistic growth model(或参见here)的形式出现在文献中。我相信你的具体问题应该在包含固定效应和随机效应的统计模型Mixed Model中考虑。

更具体地说,你可以在其统计工具箱中使用Matlab的nlmefit。 在此ppt中可以找到nlme的鸟瞰图。