How do you fit non-negative exponential decay that is biased with non-uniform noise over time?
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Eric Diaz
el 18 de Nov. de 2015
Respondida: Eric Diaz
el 30 de Dic. de 2015
We often make assumptions when we model exponential decays.
The most common assumptions that I know of are, 1) the signal is contaminated with zero-mean, gaussian distributed noise and 2) the noise distribution is uniform across time with the decay.
What do we do when we know these assumptions are both false. In my particular case, the signal is biased by a positive, non-zero mean, noise distribution that is dependent on the signal intensity, and thus the noise distribution changes with time as the signal decays.
0 comentarios
Respuesta aceptada
John D'Errico
el 18 de Nov. de 2015
Editada: John D'Errico
el 18 de Nov. de 2015
Ok, so postulate some sort of behavior for that distribution as a function of time. Use a couple of parameters to model the behavior, choosing some logical candidate distribution. It need not be perfect, just as good as you can do.
Then use maximum likelihood estimation to find the parameters for the noise model, as well as the exponential decay model. This is just an optimization, although you will probably need constraints on the parameters. And most of the time, you will need to optimize the log of the likelihood function to make things numerically tractable.
So, for example, you might postulate some sort of gamma distribution for the noise model, where the gamma parameters are a function of time. Or, pick some other distribution, anything from lognormal, to beta, even uniform. Whatever makes sense in context.
Sorry, I won't/can't write it for you, since you need to do the work up front to choose these models. But MLE is simple in concept. It ends up as a product of probabilities (thus sum of logs when you log it), which you need to optimize over.
5 comentarios
John D'Errico
el 10 de Dic. de 2015
Sorry, but no. If all that you know is the standard deviation, there are potentially many possible different distributions that have the same standard deviations. Each of those potential distributions will have very different characteristics in the modeling.
So many people assume that since they know the standard deviation, then they know everything important about the distribution under study. This probably comes from the prevalence of the Normal (Gaussian) distribution, since then knowing the mean and variance of a Normal distribution tells you all that you need to know. In this case, your distribution is VERY non-normal.
At best, you might ASSUME the noise distribution follows some specific distribution. Pick a gamma distribution, or a lognormal, for example. Once the variance of that distribution is known, then you can make inferences about the actual distribution parameters of that distribution. Depending on the distribution you choose, this may not be sufficient information to completely compute all the parameters of the distribution.
Más respuestas (1)
Ver también
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!