anisotropic_numerical_msd {AIUQ}R Documentation

Compute anisotropic numerical MSD

Description

Compute numerical mean squared displacement(MSD) based on particle trajectory for anisotropic processes in x,y-directions separately.

Usage

anisotropic_numerical_msd(pos, M, len_t)

Arguments

pos

position matrix for particle trajectory. See 'Details'.

M

number of particles

len_t

number of time steps

Details

Input pos should be the position matrix with dimension M\timeslen_t. See bm_particle_intensity, ou_particle_intensity, fbm_particle_intensity, fbm_ou_particle_intensity.

Value

A matrix of numerical MSD for given lag times in x,y-directions, dimension 2 by len_t.

Author(s)

Yue He [aut], Xubo Liu [aut], Mengyang Gu [aut, cre]

Examples

library(AIUQ)
# Simulate particle trajectory for BM
M = 10
len_t = 50
sigma = c(0.5,0.1)
pos0 = matrix(100/8+0.75*100*runif(M*2),nrow=M,ncol=2)
pos = anisotropic_bm_particle_intensity(pos0=pos0,M=M,len_t=len_t,sigma=sigma)

# Compute numerical MSD
(num_msd = anisotropic_numerical_msd(pos=pos, M=M, len_t=len_t))

[Package AIUQ version 0.5.2 Index]