This package simulate trajectories of the comb model via Langevin equations driven by fractional Gaussian noises (long-range correlated).
If you have used comb-model
in a scientific publication, we would appreciate citations to the following reference [1]:
- Haroldo V. Ribeiro, Angel A. Tateishi, Luiz G. A. Alves, Rafael S. Zola, Ervin K Lenzi (2014). Investigating the interplay between mechanisms of anomalous diffusion via fractional Brownian walks on a comb-like structure. New Journal of Physics 16, 093050
For more detailed information about the methods functions in comb-model
, please
consult its documentation.
The comb-model
package is available on pypi and can be installed using pip
pip install comb-model
numpy
stochastic
tdqm
This package offers functions to simulate Fractional Brownian walks on a plane or Comb-like structures via via Langevin equations.
# Fractional Brownian walks on a plane
import matplotlib.pyplot as plt
from combmodel import langevin_fbn
x,y=langevin_fbn(x0=0,y0=0,betax=1,betay=1,hurstx=0.5,hursty=0.5,tmax=1000,disable_tqdm=True)
f, ax = plt.subplots(nrows=1, ncols=1, figsize=(9.1,7))
ax.plot(x,y,color='#D62728',markersize=0,linewidth=2)
ax.set_xlabel(r'Position, $x(t)$')
ax.set_ylabel(r'Position, $y(t)$')
# Fractional Brownian walks on a Comb-like structure
import matplotlib.pyplot as plt
from combmodel import comb_model
x, y = comb_model(x0=0,y0=0,betax=1,betay=1,hurstx=0.5,hursty=0.5,tmax=5000,eps=1,disable_tqdm=True)
f, ax = plt.subplots(nrows=1, ncols=1, figsize=(9.1,7))
ax.plot(x,y,color='#D62728',markersize=0,linewidth=2)
ax.set_xlabel(r'Position, $x(t)$')
ax.set_ylabel(r'Position, $y(t)$')
[1] | Haroldo V. Ribeiro, Angel A. Tateishi, Luiz G. A. Alves, Rafael S. Zola, Ervin K Lenzi (2014). Investigating the interplay between mechanisms of anomalous diffusion via fractional Brownian walks on a comb-like structure. New Journal of Physics 16, 093050 |