;+ ; NAME: ; pinknoise ; ; ; PURPOSE: ; return a pink (1/f) noise spectra ; ; ; CATEGORY: ; stats, noise ; ; ; OPTIONAL INPUTS: ; n: number of samples to generate, default=1000 ; ; ; KEYWORD PARAMETERS: ; RANGE: range over which to generate samples, default=[5,100] ; SEED: random number seed ; ; ; OUTPUTS: ; array of size N distributed by a red-noise spectra (f^(-1) power law) ; ; ; EXAMPLE: ; histoplot, /xlog, /ylog, /fill, yrange=[1,100000], pinknoise(10000, range=[5,100]) ; ; ; MODIFICATION HISTORY: ; ; Tue Mar 24 12:21:57 2009, Brian Larsen ; written and tested ; ;- FUNCTION pinknoise, n, $ RANGE = range, $ SEED = seed SetDefaultValue, n, 1000 SetDefaultValue, range, [5, 100] randomp, ans, -1, n, range_x = range, SEED = seed RETURN, ans END