MathGroup Archive 2005

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: global fit

  • To: mathgroup at smc.vnet.net
  • Subject: [mg57482] Re: [mg57435] global fit
  • From: Chris Chiasson <chris.chiasson at gmail.com>
  • Date: Sun, 29 May 2005 01:03:42 -0400 (EDT)
  • References: <200505280938.FAA21631@smc.vnet.net>
  • Reply-to: Chris Chiasson <chris.chiasson at gmail.com>
  • Sender: owner-wri-mathgroup at wolfram.com

If you look in the advanced documentation for find minimum (I think it
is called unconstrained optimization), it shows how root finding
functions may not necessarily find the root nearest their starting
point. I think there are options to establish a weighted trust region
about the starting point that can help reduce the step size. Maybe
those could help you?

On 5/28/05, Pascal <Pascal.Plaza at ens.fr> wrote:
> Hi,
> 
> Here is a fascinating problem:
> 
> B is a large data matrix, typically 500x60. Each line represents the
> variations of a given physical quantity as a function of time. There
> are here 60 different times.
> 
> The problem is to fit all the lines of B with a linear combination of a
> reduced number (say 5) of base functions: f1[t], f2[t]...
> If A is the 5x60 matrix A[[i,j]] = fi[tj], I have to find an X matrix
> (500x5) such that X.A is as least-squares close to B as possible.
> At this point the solution is just X = B.PseudoInverse[A]
> 
> It gets more interesting when the base functions depend non-linearly on
> a few parameters: a, b, c...
> Now the problem is to find X and a, b, c... such that X.A[a,b,c...]is
> the best fit to the data matrix B.
> 
> I built a function of a, b, c... returning the 2-norm of X.A-B , after
> calculating A with parameters a, b, c... and estimating X with
> B.PseudoInverse[A]
> Then I tried to use FindMinimum on this function. I entered a very good
> starting point but the Kernel rapidly uses all the memory available and
> shuts down.
> 
> Does anyone has a suggestion, maybe with a completely diffrent method?
> Thanks.
> Pascal.
> 
> 


-- 
Chris Chiasson
http://chrischiasson.com/
1 (810) 265-3161


  • References:
  • Prev by Date: Re: Sorting nested list
  • Next by Date: Re: Complex Oddity
  • Previous by thread: global fit
  • Next by thread: Re: global fit