dd_gauss2

Sum of two Gaussian distributions

Syntax

info = dd_gauss2()
P = dd_gauss2(r,param)
Parameters
  • r - Distance axis (N-array)
  • param - Model parameters
Returns
  • P - Distance distribution (N-array)
  • info - Model information (struct)

Model

P(r) = a_1\sqrt{\frac{2}{\pi}}\frac{1}{\sigma_1}\exp\left(-\frac{(r-\left<r_1\right>)^2}{\sigma_1^2}\right) + (1 - a_1)\sqrt{\frac{2}{\pi}}\frac{1}{\sigma_2}\exp\left(-\frac{(r-\left<r_2\right>)^2}{\sigma_2^2}\right)

with \sigma_i = \mathrm{FWHM}_i/\sqrt{2ln(2)}

Variable Symbol Default Lower Upper Description
param(1) \left<r_1\right> 2.5 1.0 20 1st Gaussian center distance
param(2) \mathrm{FWHM}_1 0.5 0.2 5 1st Gaussian FWHM
param(3) a_1 0.5 0 1 1st Gaussian relative amplitude
param(4) \left<r_2\right> 3.5 1.0 20 2nd Gaussian center distance
param(5) \mathrm{FWHM}_2 0.5 0.2 5 2nd Gaussian FWHM

Example using default parameters:

../_images/model_dd_gauss2.png

Description

info = dd_gauss2()

Returns an info structure containing the specifics of the model:

  • info.model - Full name of the parametric model.
  • info.nparam - Total number of adjustable parameters.
  • info.parameters - Structure array with information on individual parameters.
P = dd_gauss2(r,param)

Computes the distance distribution model P from the axis r according to the parameters array param. The required parameters can also be found in the info structure.