DescriptionDOG vs MHF.png	Comparison of Difference of Gaussians and Mexican Hat Function, intended to illustrate that DoG is a reasonable approximation to MHF.
Date	6 décembre 2005, 04:58 (UTC)
Source	Originally uploaded to en.wikipedia. Generated using Mathematica 5.0.   Ce diagramme a été créé avec Mathematica
Auteur	JustinWick

Cette œuvre a été placée dans le domaine public par son auteur, JustinWick sur Wikipédia anglais. Ceci s’applique dans le monde entier. Dans certains pays, ceci peut ne pas être possible ; dans ce cas : JustinWick accorde à toute personne le droit d’utiliser cette œuvre dans n’importe quel but, sans aucune condition, sauf celles requises par la loi.Public domainPublic domainfalsefalse

Notes on Generation

To reproduce this plot in mathematica, use the following code:

<< "Statistics`NormalDistribution`";
DOG[mu_, s1_, s2_][x_] := PDF[NormalDistribution[mu, s1], x] - PDF[NormalDistribution[mu, s2], x]; 
MHF[x_] := ((2/(Sqrt[3]*Pi^(1/4)))*(1 - x^2))/E^(x^2/2);
Show[Plot[c*DOG[0, a, b][x] /. {a -> 1.3383326406221043,  b -> 0.7477302449634321, c -> -3.812773094633841}, {x, -5, 5}, 
PlotRange -> {-1,  1}, Frame -> True, Axes -> False, PlotStyle -> {RGBColor[0, 0, 1]}, PlotLabel -> "Difference of Gaussians"], 
Plot[MHF[x], {x, -5, 5}, PlotRange -> {-1, 1}, Frame -> True, Axes -> False, 
PlotStyle -> RGBColor[1, 0, 0], PlotLabel -> "Mexican Hat Function"], PlotLabel -> "Difference of Gaussians and Mexican Hat Function"]

Note, to fit the DoG best to the MHF, I utilized the following code:

NMinimize[{Hold[NIntegrate[Abs[c*DOG[0, a, b][x] - MHF[x]], {x, -5, 5}]], {a > 0, b > 0}}, {a, b, c}]

This uses least squares fitting, which may be slightly suboptimal.

File history on en.wikipedia

Description page history

File versions

• 04:58, 6 December 2005 . . JustinWick . . 486x900 (10415 bytes) (Comparison of Difference of Gaussians Difference of Gaussians and Mexican Hat Function. This file was created at 19:02, 18 March 2006 (UTC) by JustinWick, using Mathematica 5.0.)