Simulates Fractional Brownian field on unit disk, with Hurst parameter 'H';
Note that the covariance function is isotropic, see reference below.
- 'H' is the Hurst parameter of the Gaussian process
- 'n' is the number of grid points, where 'n' is a power of 2;
if the 'n' supplied is not a power of two,
then we set n=2^ceil(log2(n)); default is n=2^8;
- two statistically independent fields 'field1' and 'field2'
over unit disk; if not output requested, then function
outputs a figure of one of the fields
- vectors 'tx' and 'ty' so that the field is plotted via
Kroese, D. P., & Botev, Z. I. (2015). Spatial Process Simulation.
In Stochastic Geometry, Spatial Statistics and Random Fields(pp. 369-404)
Springer International Publishing, DOI: 10.1007/978-3-319-10064-7_12
How can we measure these fractals using box counting dimension method.
Thanks for sharing this file. I am trying to run the file, but Matlab says that is necessary to have a function rho. Could you provide a rho function as an example? You have done it in your other code (stationary_Gaussian_process) and it was great to have a start point.
- updated reference
Create scripts with code, output, and formatted text in a single executable document.