[Date Index]
[Thread Index]
[Author Index]
Re: maximum entropy method for deconvolution
*To*: mathgroup at smc.vnet.net
*Subject*: [mg75420] Re: maximum entropy method for deconvolution
*From*: Roman <rschmied at gmail.com>
*Date*: Sun, 29 Apr 2007 03:12:51 -0400 (EDT)
*References*: <f0v6f3$929$1@smc.vnet.net>
Dan,
Here's a list of books:
http://www.amazon.de/s?ie=UTF8&index=books-de-intl-us&field-keywords=Maximum%20entropy%20method&page=1
I had used the one named "Maximum Entropy in Action", but took quite
some time to implement something even vaguely correct. I had found all
explanations quite opaque. About a year ago I searched for quite a
while for a free implementation of MEM, but could not find any. I'll
think about how to do your problem, since I guess it is of very
general interest. One of the central parts is the question on how to
define the entropy of the resulting F. I imagine you could use
something of the sort of
S = Integrate[-F(w)*Log(F(w)), {w,-Inf,Inf}]
More later, now it's off to the mountains for some hiking.
Roman.
On Apr 28, 12:08 pm, dantimatter <dantimat... at gmail.com> wrote:
> hello all,
>
> first off, many thanks to 'Roman' et al for all the previous help with
> my inversion problem.
>
> i have a convolution function G which is the convolution of F and p (G
> = F**p). i know G and i know p, and i'd like to extract F. i can do
> this by taking the Fourier transform of G, dividing by the Fourier
> transform of p, and inverting the result to get F. the problem is
> that p is a step function (p = UnitStep[t]*UnitStep[33.6-t]) which has
> a lot of zeros in frequency space, and thus it is difficult to get at
> F via inversion. Mathematica is happily doing the inversion but the
> results are very noticeably wrong.
>
> i understand from my conversations with some of you and much time
> spent in the library that this is in general a difficult problem, but
> there are some methods that are known to make this type of problem
> tractable, such as the maximum entropy method (MEM) for inversion. is
> anyone aware of an implementation of a MEM algorithm in Mathematica?
> i have read Numerical Recipes a couple of times and i am unable to get
> my head around the relevant chapter. If there isn't a Mathematica
> implementation, perhaps someone could offer some advice on where else
> to look? if it exists, a "for dummies" type book with step-by-step
> instructions would be the best resource for me...
>
> cheers,
> dan
Prev by Date:
**Re: RE: Re: Sometimes <space> means multiple**
Next by Date:
**Re: Sometimes <space> means multiple**
Previous by thread:
**Re: maximum entropy method for deconvolution**
Next by thread:
** Loading packages from the internet?**
| |