MathGroup Archive 2007

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

Search the Archive

Re: Fit

  • To: mathgroup at smc.vnet.net
  • Subject: [mg76173] Re: Fit
  • From: dh <dh at metrohm.ch>
  • Date: Thu, 17 May 2007 06:00:14 -0400 (EDT)
  • References: <f2ejm8$d2$1@smc.vnet.net>


Hi Mathieu,

it depend what you want to do. If it is important that the approximation 

is good near a point, e.g. 0, one can expand both function in a series 

and equate the coefficients:

Series[f[x,y],{x,0,2},{y,0,2}]//Normal//N

Series[Gaussian2D[x,y,rad,amp],{x,0,2},{y,0,2}]//Normal//N

in your case we would euqte the zero and quadratic (x^2,y^2 have the 

same factor) terms.

On the otherhand, if the approximation is to be used over a region, one 

can do it numerical. Choose enough points in the region, calculate the 

function values at these points and use e.g. FindFit to fit the Gaussion.

hope this helps, Daniel



Mathieu G wrote:

> Hello,

> I have an expression of the form:

> 

> 8000000 Sqrt[2/\[Pi]] (-Erf[(-300 + x)/(4000 Sqrt[2])] +

>     Erf[(300 + x)/(4000 Sqrt[2])]) (-1/2 Sqrt[\[Pi]/2]

>       Erf[(-300 + y)/(4000 Sqrt[2])] +

>     1/2 Sqrt[\[Pi]/2] Erf[(300 + y)/(4000 Sqrt[2])])

> 

> how can I make a fit using a 2D Gaussian function of the form

> 

> Gaussian2D[x_, y_, Radius_, Amplitude_] :=

>    Amplitude Exp[-1/2 ((x/Radius)^2 + (y/Radius)^2)];

> 

> Thank you for your help!

> MG

> 




  • Prev by Date: Re: Why Sort[{"AX", "!D", "EX"}] -> {"AX", "!D", "EX"} ?
  • Next by Date: Re: Re: Compatibility woes
  • Previous by thread: Fit
  • Next by thread: Re: Fit