Generate spatial data

Generate 1/f noise spatial data.

Description

The script generates spatial data with a scale-invariant power spectrum (1/f noise) and a normal error distribution.

The spectral density of the data is proportional to f^BETA, where f is the frequency and BETA is the spectral exponent (BETA=0 is white noise, BETA=-2 is Brown noise).

To generate Brown noise on a 10x10 grid you type

spatialPattern([10,10],-2)

which gives a 10x10 matrix containihng the data.

1/f noise is not stationary, and so doesn't have a well defined variance. A variogram of the data show variance linearly increasing with increasing lag.

Time series data can also be generated. For example, a pink noise time series of lenth 1000 is produced by

spatialPattern([1000,1],-1)

• MATLAB release: MATLAB 6.5 (R13)

[Gomez Lopez]

 

[Madhukar U S]

Its of great use.

[Stanley Bileschi]

fantastic, works just like you want it to. It's a quick edit to make the default parameters a bit more user friendly, like so that Dim can be one number, and that automatically BETA sets to -2, for 1/f noise

[Md Nurujjaman]

 

[Terence Picton]

fine program - easy to use for pink noise for auditory purposes, haven't tried the spatial noises.

[Yuval Cohen]

I use it for audio noise generation. It works great. It would be beneficial to eliminate non-audible components (below 20Hz) in order to improve the power density of the signal.

[Gabriel Akira Schreiber]

Great, thank you very much.