dd_rice2

Sum of two 3D-Rice distributions

Syntax

info = dd_rice2()
P = dd_rice2(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 R(r,\nu_1,\sigma_1) + (1-a_1) R(r,\nu_2,\sigma_2)

R(r,\nu,\sigma) = \frac{\nu^{n/2-1}}{\sigma^2}r^{n/2}\exp\left(-\frac{(r^2+\nu^2)}{2\sigma^2}\right)I_{n/2-1}\left(\frac{r\nu}{\sigma^2} \right)

where n=3 and I_{n/2-1}(x) is the modified Bessel function of the first kind with order n/2-1. This is a three-dimensional non-central chi distribution, the 3D generalization of the 2D Rice distribution.

Variable Symbol Default Lower Upper Description
param(1) \nu_1 2.5 1.0 10 center, 1st component
param(2) \sigma_1 0.7 0.1 5 width, 1st component
param(3) a_1 0.5 0 1 amplitude, 1st component
param(4) \nu_2 4.0 1.0 10 center, 2nd component
param(5) \sigma_2 0.7 0.1 5 width, 2nd component

Example using default parameters:

../_images/model_dd_rice2.png

Description

info = dd_rice2()

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_rice2(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.