Misplaced Pages

Pink noise

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Fractal noise) Signal with equal energy per octave
Colors of noise
For the album by Laura Mvula, see Pink Noise (album).This article is about 1/f noise. For quantum 1/f noise, see quantum 1/f noise."Fractal noise" redirects here. For the novel, see Fractal Noise.
A two-dimensional pink noise grayscale image, generated with a computer program; some fields observed in nature are characterized by a similar power spectrum
A 3D pink noise image, generated with a computer program, viewed as an animation in which each frame is a 2D slice

Pink noise, 1⁄f noise, fractional noise or fractal noise is a signal or process with a frequency spectrum such that the power spectral density (power per frequency interval) is inversely proportional to the frequency of the signal. In pink noise, each octave interval (halving or doubling in frequency) carries an equal amount of noise energy.

Pink noise sounds like a waterfall. It is often used to tune loudspeaker systems in professional audio. Pink noise is one of the most commonly observed signals in biological systems.

The name arises from the pink appearance of visible light with this power spectrum. This is in contrast with white noise which has equal intensity per frequency interval.

Definition

Pink noise 10 seconds of pink noise, normalized to −1 dBFS peak amplitude
Problems playing this file? See media help.

Within the scientific literature, the term 1/f noise is sometimes used loosely to refer to any noise with a power spectral density of the form

S ( f ) 1 f α , {\displaystyle S(f)\propto {\frac {1}{f^{\alpha }}},}

where f is frequency, and 0 < α < 2, with exponent α usually close to 1. One-dimensional signals with α = 1 are usually called pink noise.

The following function describes a length N {\displaystyle N} one-dimensional pink noise signal (i.e. a Gaussian white noise signal with zero mean and standard deviation σ {\displaystyle \sigma } , which has been suitably filtered), as a sum of sine waves with different frequencies, whose amplitudes fall off inversely with the square root of frequency u {\displaystyle u} (so that power, which is the square of amplitude, falls off inversely with frequency), and phases are random:

h ( x ) = σ N 2 u χ u u sin ( 2 π u x N + ϕ u ) , χ u χ ( 2 ) , ϕ u U ( 0 , 2 π ) . {\displaystyle h(x)=\sigma {\sqrt {\frac {N}{2}}}\sum _{u}{\frac {\chi _{u}}{\sqrt {u}}}\sin({\frac {2\pi ux}{N}}+\phi _{u}),\quad \chi _{u}\sim \chi (2),\quad \phi _{u}\sim U(0,2\pi ).}

χ u {\displaystyle \chi _{u}} are iid chi-distributed variables, and ϕ u {\displaystyle \phi _{u}} are uniform random.

In a two-dimensional pink noise signal, the amplitude at any orientation falls off inversely with frequency. A pink noise square of length N {\displaystyle N} can be written as: h ( x , y ) = σ N 2 u , v χ u v u 2 + v 2 sin ( 2 π N ( u x + v y ) + ϕ u v ) , χ u v χ ( 2 ) , ϕ u v U ( 0 , 2 π ) . {\displaystyle h(x,y)={\frac {\sigma N}{\sqrt {2}}}\sum _{u,v}{\frac {\chi _{uv}}{\sqrt {u^{2}+v^{2}}}}\sin \left({\frac {2\pi }{N}}(ux+vy)+\phi _{uv}\right),\quad \chi _{uv}\sim \chi (2),\quad \phi _{uv}\sim U(0,2\pi ).}

General 1/f-like noises occur widely in nature and are a source of considerable interest in many fields. Noises with α near 1 generally come from condensed-matter systems in quasi-equilibrium, as discussed below. Noises with a broad range of α generally correspond to a wide range of non-equilibrium driven dynamical systems.

Pink noise sources include flicker noise in electronic devices. In their study of fractional Brownian motion, Mandelbrot and Van Ness proposed the name fractional noise (sometimes since called fractal noise) to describe 1/f noises for which the exponent α is not an even integer, or that are fractional derivatives of Brownian (1/f) noise.

Description

Spectrum of a pink noise approximation on a log-log plot; power density falls off at 10 dB/decade of frequency
Relative intensity of pink noise (left) and white noise (right) on an FFT spectrogram with the vertical axis being linear frequency

In pink noise, there is equal energy per octave of frequency. The energy of pink noise at each frequency level, however, falls off at roughly 3 dB per octave. This is in contrast to white noise which has equal energy at all frequency levels.

The human auditory system, which processes frequencies in a roughly logarithmic fashion approximated by the Bark scale, does not perceive different frequencies with equal sensitivity; signals around 1–4 kHz sound loudest for a given intensity. However, humans still differentiate between white noise and pink noise with ease.

Graphic equalizers also divide signals into bands logarithmically and report power by octaves; audio engineers put pink noise through a system to test whether it has a flat frequency response in the spectrum of interest. Systems that do not have a flat response can be equalized by creating an inverse filter using a graphic equalizer. Because pink noise tends to occur in natural physical systems, it is often useful in audio production. Pink noise can be processed, filtered, and/or effects can be added to produce desired sounds. Pink-noise generators are commercially available.

One parameter of noise, the peak versus average energy contents, or crest factor, is important for testing purposes, such as for audio power amplifier and loudspeaker capabilities because the signal power is a direct function of the crest factor. Various crest factors of pink noise can be used in simulations of various levels of dynamic range compression in music signals. On some digital pink-noise generators the crest factor can be specified.

Generation

The spatial filter which is convolved with a one-dimensional white noise signal to create a pink noise signal

Pink noise can be computer-generated by first generating a white noise signal, Fourier-transforming it, then dividing the amplitudes of the different frequency components by the square root of the frequency (in one dimension), or by the frequency (in two dimensions) etc. This is equivalent to spatially filtering (convolving) the white noise signal with a white-to-pink-filter. For a length N {\displaystyle N} signal in one dimension, the filter has the following form:

a ( x ) = 1 N [ 1 + 1 N / 2 cos π ( x 1 ) + 2 k = 1 N / 2 1 1 k cos 2 π k N ( x 1 ) ] . {\displaystyle a(x)={\frac {1}{N}}\left.}

Matlab programs are available to generate pink and other power-law coloured noise in one or any number of dimensions.

Properties

The autocorrelation (Pearson's correlation coefficient) of one-dimensional (top) and two-dimensional (bottom) pink noise signals, across distance d (in units of the longest wavelength comprising the signal); grey curves are the autocorrelations of a sample of pink noise signals (comprising discrete frequencies), and black is their average, red is the theoretically calculated autocorrelation when the signal comprises these same discrete frequencies, and blue assumes a continuum of frequencies

Power-law spectra

The power spectrum of pink noise is 1 f {\displaystyle {\frac {1}{f}}} only for one-dimensional signals. For two-dimensional signals (e.g., images) the average power spectrum at any orientation falls as 1 f 2 {\displaystyle {\frac {1}{f^{2}}}} , and in d {\displaystyle d} dimensions, it falls as 1 f d {\displaystyle {\frac {1}{f^{d}}}} . In every case, each octave carries an equal amount of noise power.

The average amplitude a θ {\displaystyle a_{\theta }} and power p θ {\displaystyle p_{\theta }} of a pink noise signal at any orientation θ {\displaystyle \theta } , and the total power across all orientations, fall off as some power of the frequency. The following table lists these power-law frequency-dependencies for pink noise signal in different dimensions, and also for general power-law colored noise with power α {\displaystyle \alpha } (e.g.: Brown noise has α = 2 {\displaystyle \alpha =2} ):

Power-law spectra of pink noise
dimensions avg. amp. a θ ( f ) {\displaystyle a_{\theta }(f)} avg. power p θ ( f ) {\displaystyle p_{\theta }(f)} tot. power p ( f ) {\displaystyle p(f)}
1 1 / f {\displaystyle 1/{\sqrt {f}}} 1 / f {\displaystyle 1/f} 1 / f {\displaystyle 1/f}
2 1 / f {\displaystyle 1/f} 1 / f 2 {\displaystyle 1/f^{2}} 1 / f {\displaystyle 1/f}
3 1 / f 3 / 2 {\displaystyle 1/f^{3/2}} 1 / f 3 {\displaystyle 1/f^{3}} 1 / f {\displaystyle 1/f}
d {\displaystyle d} 1 / f d / 2 {\displaystyle 1/f^{d/2}} 1 / f d {\displaystyle 1/f^{d}} 1 / f {\displaystyle 1/f}
d {\displaystyle d} , power α {\displaystyle \alpha } 1 / f α d / 2 {\displaystyle 1/f^{\alpha d/2}} 1 / f α d {\displaystyle 1/f^{\alpha d}} 1 / f 1 + ( α 1 ) d {\displaystyle 1/f^{1+(\alpha -1)d}}

Distribution of point values

Consider pink noise of any dimension that is produced by generating a Gaussian white noise signal with mean μ {\displaystyle \mu } and sd σ {\displaystyle \sigma } , then multiplying its spectrum with a filter (equivalent to spatially filtering it with a filter a {\displaystyle {\boldsymbol {a}}} ). Then the point values of the pink noise signal will also be normally distributed, with mean μ {\displaystyle \mu } and sd a σ {\displaystyle \lVert {\boldsymbol {a}}\rVert \sigma } .

Autocorrelation

Unlike white noise, which has no correlations across the signal, a pink noise signal is correlated with itself, as follows.

1D signal

The Pearson's correlation coefficient of a one-dimensional pink noise signal (comprising discrete frequencies k {\displaystyle k} ) with itself across a distance d {\displaystyle d} in the configuration (space or time) domain is: r ( d ) = k cos 2 π k d N k k 1 k . {\displaystyle r(d)={\frac {\sum _{k}{\frac {\cos {\frac {2\pi kd}{N}}}{k}}}{\sum _{k}{\frac {1}{k}}}}.} If instead of discrete frequencies, the pink noise comprises a superposition of continuous frequencies from k min {\displaystyle k_{\textrm {min}}} to k max {\displaystyle k_{\textrm {max}}} , the autocorrelation coefficient is: r ( d ) = Ci ( 2 π k max d N ) Ci ( 2 π k min d N ) log k max k min , {\displaystyle r(d)={\frac {{\textrm {Ci}}({\frac {2\pi k_{\textrm {max}}d}{N}})-{\textrm {Ci}}({\frac {2\pi k_{\textrm {min}}d}{N}})}{\log {\frac {k_{\textrm {max}}}{k_{\textrm {min}}}}}},} where Ci ( x ) {\displaystyle {\textrm {Ci}}(x)} is the cosine integral function.

2D signal

The Pearson's autocorrelation coefficient of a two-dimensional pink noise signal comprising discrete frequencies is theoretically approximated as: r ( d ) = k J 0 ( 2 π k d N ) k k 1 k , {\displaystyle r(d)={\frac {\sum _{k}{\frac {J_{0}({\frac {2\pi kd}{N}})}{k}}}{\sum _{k}{\frac {1}{k}}}},} where J 0 {\displaystyle J_{0}} is the Bessel function of the first kind.

Occurrence

Pink noise has been discovered in the statistical fluctuations of an extraordinarily diverse number of physical and biological systems (Press, 1978; see articles in Handel & Chung, 1993, and references therein). Examples of its occurrence include fluctuations in tide and river heights, quasar light emissions, heart beat, firings of single neurons, resistivity in solid-state electronics and single-molecule conductance signals resulting in flicker noise. Pink noise describes the statistical structure of many natural images.

General 1/f noises occur in many physical, biological and economic systems, and some researchers describe them as being ubiquitous. In physical systems, they are present in some meteorological data series, the electromagnetic radiation output of some astronomical bodies. In biological systems, they are present in, for example, heart beat rhythms, neural activity, and the statistics of DNA sequences, as a generalized pattern.

An accessible introduction to the significance of pink noise is one given by Martin Gardner (1978) in his Scientific American column "Mathematical Games". In this column, Gardner asked for the sense in which music imitates nature. Sounds in nature are not musical in that they tend to be either too repetitive (bird song, insect noises) or too chaotic (ocean surf, wind in trees, and so forth). The answer to this question was given in a statistical sense by Voss and Clarke (1975, 1978), who showed that pitch and loudness fluctuations in speech and music are pink noises. So music is like tides not in terms of how tides sound, but in how tide heights vary.

Precision timekeeping

Main article: Allan variance

The ubiquitous 1/f noise poses a "noise floor" to precision timekeeping. The derivation is based on.

A clock is most easily tested by comparing it with a far more accurate reference clock. During an interval of time τ, as measured by the reference clock, the clock under test advances by τy, where y is the average (relative) clock frequency over that interval.

Suppose that we have a timekeeping device (it could be anything from quartz oscillators, atomic clocks, and hourglasses). Let its readout be a real number x ( t ) {\displaystyle x(t)} that changes with the actual time t {\displaystyle t} . For concreteness, let us consider a quartz oscillator. In a quartz oscillator, x ( t ) {\displaystyle x(t)} is the number of oscillations, and x ˙ ( t ) {\displaystyle {\dot {x}}(t)} is the rate of oscillation. The rate of oscillation has a constant component x ˙ 0 {\displaystyle {\dot {x}}_{0}} and a fluctuating component x ˙ f {\displaystyle {\dot {x}}_{f}} , so x ˙ ( t ) = x ˙ 0 + x ˙ f ( t ) {\textstyle {\dot {x}}(t)={\dot {x}}_{0}+{\dot {x}}_{f}(t)} . By selecting the right units for x {\displaystyle x} , we can have x ˙ 0 = 1 {\displaystyle {\dot {x}}_{0}=1} , meaning that on average, one second of clock-time passes for every second of real-time.

The stability of the clock is measured by how many "ticks" it makes over a fixed interval. The more stable the number of ticks, the better the stability of the clock. So, define the average clock frequency over the interval [ k τ , ( k + 1 ) τ ] {\displaystyle } as y k = 1 τ k τ ( k + 1 ) τ x ˙ ( t ) d t = x ( ( k + 1 ) τ ) x ( k τ ) τ {\displaystyle y_{k}={\frac {1}{\tau }}\int _{k\tau }^{(k+1)\tau }{\dot {x}}(t)dt={\frac {x((k+1)\tau )-x(k\tau )}{\tau }}} Note that y k {\displaystyle y_{k}} is unitless: it is the numerical ratio between ticks of the physical clock and ticks of an ideal clock.

The Allan variance of the clock frequency is half the mean square of change in average clock frequency: σ 2 ( τ ) = 1 2 ( y k y k 1 ) 2 ¯ = 1 K k = 1 K 1 2 ( y k y k 1 ) 2 {\displaystyle \sigma ^{2}(\tau )={\frac {1}{2}}{\overline {(y_{k}-y_{k-1})^{2}}}={\frac {1}{K}}\sum _{k=1}^{K}{\frac {1}{2}}(y_{k}-y_{k-1})^{2}} where K {\displaystyle K} is an integer large enough for the averaging to converge to a definite value. For example, a 2013 atomic clock achieved σ ( 25000  seconds ) = 1.6 × 10 18 {\displaystyle \sigma (25000{\text{ seconds}})=1.6\times 10^{-18}} , meaning that if the clock is used to repeatedly measure intervals of 7 hours, the standard deviation of the actually measured time would be around 40 femtoseconds.

Now we have y k y k 1 = R g ( k τ t ) x ˙ f ( t ) d t = ( g x ˙ f ) ( k τ ) {\displaystyle y_{k}-y_{k-1}=\int _{\mathbb {R} }g(k\tau -t){\dot {x}}_{f}(t)dt=(g\ast {\dot {x}}_{f})(k\tau )} where g ( t ) = 1 [ 0 , τ ] ( t ) + 1 [ τ , 0 ] ( t ) τ {\displaystyle g(t)={\frac {-1_{}(t)+1_{}(t)}{\tau }}} is one packet of a square wave with height 1 / τ {\displaystyle 1/\tau } and wavelength 2 τ {\displaystyle 2\tau } . Let h ( t ) {\displaystyle h(t)} be a packet of a square wave with height 1 and wavelength 2, then g ( t ) = h ( t / τ ) / τ {\displaystyle g(t)=h(t/\tau )/\tau } , and its Fourier transform satisfies F [ g ] ( ω ) = F [ h ] ( τ ω ) {\displaystyle {\mathcal {F}}(\omega )={\mathcal {F}}(\tau \omega )} .

The Allan variance is then σ 2 ( τ ) = 1 2 ( y k y k 1 ) 2 ¯ = 1 2 ( g x ˙ f ) ( k τ ) 2 ¯ {\displaystyle \sigma ^{2}(\tau )={\frac {1}{2}}{\overline {(y_{k}-y_{k-1})^{2}}}={\frac {1}{2}}{\overline {(g\ast {\dot {x}}_{f})(k\tau )^{2}}}} , and the discrete averaging can be approximated by a continuous averaging: 1 K k = 1 K 1 2 ( y k y k 1 ) 2 1 K τ 0 K τ 1 2 ( g x ˙ f ) ( t ) 2 d t {\displaystyle {\frac {1}{K}}\sum _{k=1}^{K}{\frac {1}{2}}(y_{k}-y_{k-1})^{2}\approx {\frac {1}{K\tau }}\int _{0}^{K\tau }{\frac {1}{2}}(g\ast {\dot {x}}_{f})(t)^{2}dt} , which is the total power of the signal ( g x ˙ f ) {\displaystyle (g\ast {\dot {x}}_{f})} , or the integral of its power spectrum:

σ 2 ( 1 ) {\displaystyle \sigma ^{2}(1)} is approximately the area under the green curve; when τ {\displaystyle \tau } increases, S [ g ] ( ω ) {\displaystyle S(\omega )} shrinks on the x-axis, and the green curve shrinks on the x-axis but expands on the y-axis; when S [ x ˙ f ] ( ω ) ω α {\displaystyle S(\omega )\propto \omega ^{-\alpha }} , the combined effect of both is that σ 2 ( τ ) τ α 1 {\displaystyle \sigma ^{2}(\tau )\propto \tau ^{\alpha -1}}

σ 2 ( τ ) 0 S [ g x ˙ f ] ( ω ) d ω = 0 S [ g ] ( ω ) S [ x ˙ f ] ( ω ) d ω = 0 S [ h ] ( τ ω ) S [ x ˙ f ] ( ω ) d ω {\displaystyle \sigma ^{2}(\tau )\approx \int _{0}^{\infty }S(\omega )d\omega =\int _{0}^{\infty }S(\omega )\cdot S(\omega )d\omega =\int _{0}^{\infty }S(\tau \omega )\cdot S(\omega )d\omega } In words, the Allan variance is approximately the power of the fluctuation after bandpass filtering at ω 1 / τ {\displaystyle \omega \sim 1/\tau } with bandwidth Δ ω 1 / τ {\displaystyle \Delta \omega \sim 1/\tau } .


For 1 / f α {\displaystyle 1/f^{\alpha }} fluctuation, we have S [ x ˙ f ] ( ω ) = C / ω α {\displaystyle S(\omega )=C/\omega ^{\alpha }} for some constant C {\displaystyle C} , so σ 2 ( τ ) τ α 1 σ 2 ( 1 ) τ α 1 {\displaystyle \sigma ^{2}(\tau )\approx \tau ^{\alpha -1}\sigma ^{2}(1)\propto \tau ^{\alpha -1}} . In particular, when the fluctuating component x ˙ f {\displaystyle {\dot {x}}_{f}} is a 1/f noise, then σ 2 ( τ ) {\displaystyle \sigma ^{2}(\tau )} is independent of the averaging time τ {\displaystyle \tau } , meaning that the clock frequency does not become more stable by simply averaging for longer. This contrasts with a white noise fluctuation, in which case σ 2 ( τ ) τ 1 {\displaystyle \sigma ^{2}(\tau )\propto \tau ^{-1}} , meaning that doubling the averaging time would improve the stability of frequency by 2 {\displaystyle {\sqrt {2}}} .

The cause of the noise floor is often traced to particular electronic components (such as transistors, resistors, and capacitors) within the oscillator feedback.

Humans

In brains, pink noise has been widely observed across many temporal and physical scales from ion channel gating to EEG and MEG and LFP recordings in humans. In clinical EEG, deviations from this 1/f pink noise can be used to identify epilepsy, even in the absence of a seizure, or during the interictal state. Classic models of EEG generators suggested that dendritic inputs in gray matter were principally responsible for generating the 1/f power spectrum observed in EEG/MEG signals. However, recent computational models using cable theory have shown that action potential transduction along white matter tracts in the brain also generates a 1/f spectral density. Therefore, white matter signal transduction may also contribute to pink noise measured in scalp EEG recordings, particularly if the effects of ephaptic coupling are taken into consideration.

It has also been successfully applied to the modeling of mental states in psychology, and used to explain stylistic variations in music from different cultures and historic periods. Richard F. Voss and J. Clarke claim that almost all musical melodies, when each successive note is plotted on a scale of pitches, will tend towards a pink noise spectrum. Similarly, a generally pink distribution pattern has been observed in film shot length by researcher James E. Cutting of Cornell University, in the study of 150 popular movies released from 1935 to 2005.

Pink noise has also been found to be endemic in human response. Gilden et al. (1995) found extremely pure examples of this noise in the time series formed upon iterated production of temporal and spatial intervals. Later, Gilden (1997) and Gilden (2001) found that time series formed from reaction time measurement and from iterated two-alternative forced choice also produced pink noises.

Electronic devices

Main article: Flicker noise

The principal sources of pink noise in electronic devices are almost invariably the slow fluctuations of properties of the condensed-matter materials of the devices. In many cases the specific sources of the fluctuations are known. These include fluctuating configurations of defects in metals, fluctuating occupancies of traps in semiconductors, and fluctuating domain structures in magnetic materials. The explanation for the approximately pink spectral form turns out to be relatively trivial, usually coming from a distribution of kinetic activation energies of the fluctuating processes. Since the frequency range of the typical noise experiment (e.g., 1 Hz – 1 kHz) is low compared with typical microscopic "attempt frequencies" (e.g., 10 Hz), the exponential factors in the Arrhenius equation for the rates are large. Relatively small spreads in the activation energies appearing in these exponents then result in large spreads of characteristic rates. In the simplest toy case, a flat distribution of activation energies gives exactly a pink spectrum, because d d f ln f = 1 f . {\displaystyle \textstyle {\frac {d}{df}}\ln f={\frac {1}{f}}.}

There is no known lower bound to background pink noise in electronics. Measurements made down to 10 Hz (taking several weeks) have not shown a ceasing of pink-noise behaviour. (Kleinpenning, de Kuijper, 1988) measured the resistance in a noisy carbon-sheet resistor, and found 1/f noise behavior over the range of [ 10 5.5 H z , 10 4 H z ] {\displaystyle } , a range of 9.5 decades.

A pioneering researcher in this field was Aldert van der Ziel.

Flicker noise is commonly used for the reliability characterization of electronic devices. It is also used for gas detection in chemoresistive sensors by dedicated measurement setups.

In gravitational wave astronomy

Noise curves for a selection of gravitational-wave detectors as a function of frequency

1/f noises with α near 1 are a factor in gravitational-wave astronomy. The noise curve at very low frequencies affects pulsar timing arrays, the European Pulsar Timing Array (EPTA) and the future International Pulsar Timing Array (IPTA); at low frequencies are space-borne detectors, the formerly proposed Laser Interferometer Space Antenna (LISA) and the currently proposed evolved Laser Interferometer Space Antenna (eLISA), and at high frequencies are ground-based detectors, the initial Laser Interferometer Gravitational-Wave Observatory (LIGO) and its advanced configuration (aLIGO). The characteristic strain of potential astrophysical sources are also shown. To be detectable the characteristic strain of a signal must be above the noise curve.

Climate dynamics

Pink noise on timescales of decades has been found in climate proxy data, which may indicate amplification and coupling of processes in the climate system.

Diffusion processes

Many time-dependent stochastic processes are known to exhibit 1/f noises with α between 0 and 2. In particular Brownian motion has a power spectral density that equals 4D/f, where D is the diffusion coefficient. This type of spectrum is sometimes referred to as Brownian noise. The analysis of individual Brownian motion trajectories also show 1/f spectrum, albeit with random amplitudes. Fractional Brownian motion with Hurst exponent H also show 1/f power spectral density with α=2H+1 for subdiffusive processes (H<0.5) and α=2 for superdiffusive processes (0.5<H<1).

Origin

See also: Supersymmetry § Supersymmetry in dynamical systems

There are many theories about the origin of pink noise. Some theories attempt to be universal, while others apply to only a certain type of material, such as semiconductors. Universal theories of pink noise remain a matter of current research interest.

A hypothesis (referred to as the Tweedie hypothesis) has been proposed to explain the genesis of pink noise on the basis of a mathematical convergence theorem related to the central limit theorem of statistics. The Tweedie convergence theorem describes the convergence of certain statistical processes towards a family of statistical models known as the Tweedie distributions. These distributions are characterized by a variance to mean power law, that have been variously identified in the ecological literature as Taylor's law and in the physics literature as fluctuation scaling. When this variance to mean power law is demonstrated by the method of expanding enumerative bins this implies the presence of pink noise, and vice versa. Both of these effects can be shown to be the consequence of mathematical convergence such as how certain kinds of data will converge towards the normal distribution under the central limit theorem. This hypothesis also provides for an alternative paradigm to explain power law manifestations that have been attributed to self-organized criticality.

There are various mathematical models to create pink noise. The superposition of exponentially decaying pulses is able to generate a signal with the 1 / f {\displaystyle 1/f} -spectrum at moderate frequencies, transitioning to a constant at low frequencies and 1 / f 2 {\displaystyle 1/f^{2}} at high frequencies. In contrast, the sandpile model of self-organized criticality, which exhibits quasi-cycles of gradual stress accumulation between fast rare stress-releases, reproduces the flicker noise that corresponds to the intra-cycle dynamics. The statistical signature of self-organization is justified in It can be generated on computer, for example, by filtering white noise, inverse Fourier transform, or by multirate variants on standard white noise generation.

In supersymmetric theory of stochastics, an approximation-free theory of stochastic differential equations, 1/f noise is one of the manifestations of the spontaneous breakdown of topological supersymmetry. This supersymmetry is an intrinsic property of all stochastic differential equations and its meaning is the preservation of the continuity of the phase space by continuous time dynamics. Spontaneous breakdown of this supersymmetry is the stochastic generalization of the concept of deterministic chaos, whereas the associated emergence of the long-term dynamical memory or order, i.e., 1/f and crackling noises, the Butterfly effect etc., is the consequence of the Goldstone theorem in the application to the spontaneously broken topological supersymmetry.

Audio testing

Pink noise is commonly used to test the loudspeakers in sound reinforcement systems, with the resulting sound measured with a test microphone in the listening space connected to a spectrum analyzer or a computer running a real-time fast Fourier transform (FFT) analyzer program such as Smaart. The sound system plays pink noise while the audio engineer makes adjustments on an audio equalizer to obtain the desired results. Pink noise is predictable and repeatable, but it is annoying for a concert audience to hear. Since the late 1990s, FFT-based analysis enabled the engineer to make adjustments using pre-recorded music as the test signal, or even the music coming from the performers in real time. Pink noise is still used by audio system contractors and by computerized sound systems which incorporate an automatic equalization feature.

In manufacturing, pink noise is often used as a burn-in signal for audio amplifiers and other components, to determine whether the component will maintain performance integrity during sustained use. The process of end-users burning in their headphones with pink noise to attain higher fidelity has been called an audiophile "myth".

See also

Footnotes

  1. ^ Field, D. J. (1987). "Relations between the statistics of natural images and the response properties of cortical cells" (PDF). J. Opt. Soc. Am. A. 4 (12): 2379–2394. Bibcode:1987JOSAA...4.2379F. CiteSeerX 10.1.1.136.1345. doi:10.1364/JOSAA.4.002379. PMID 3430225.
  2. "Glossary: Pink Noise". Sound on Sound. Retrieved November 22, 2022.
  3. ^ Davis, Gary; Jones, Ralph (1987). The Sound Reinforcement Handbook. Hal Leonard. p. 107. ISBN 0-88188-900-8.
  4. Szendro, P (2001). "Pink-Noise Behaviour of Biosystems". European Biophysics Journal. 30 (3): 227–231. doi:10.1007/s002490100143. PMID 11508842. S2CID 24505215.
  5. Downey, Allen (2012). Think Complexity. O'Reilly Media. p. 79. ISBN 978-1-4493-1463-7. Visible light with this power spectrum looks pink, hence the name.
  6. Baxandall, P. J. (November 1968). "Noise in Transistor Circuits: 1 - Mainly on fundamental noise concepts" (PDF). Wireless World. pp. 388–392. Archived (PDF) from the original on 2016-04-23. Retrieved 2019-08-08.
  7. ^ Das, Abhranil (2022). Camouflage detection & signal discrimination: theory, methods & experiments (corrected) (PhD). The University of Texas at Austin. doi:10.13140/RG.2.2.10585.80487.
  8. ^ Kogan, Shulim (1996). Electronic Noise and Fluctuations in Solids. . ISBN 978-0-521-46034-7.
  9. Mandelbrot, B. B.; Van Ness, J. W. (1968). "Fractional Brownian motions, fractional noises and applications". SIAM Review. 10 (4): 422–437. Bibcode:1968SIAMR..10..422M. doi:10.1137/1010093.
  10. Mandelbrot, Benoit B.; Wallis, James R. (1969). "Computer Experiments with Fractional Gaussian Noises: Part 3, Mathematical Appendix". Water Resources Research. 5 (1): 260–267. Bibcode:1969WRR.....5..260M. doi:10.1029/WR005i001p00260.
  11. "Noise". www.sfu.ca. Retrieved 2024-02-06.
  12. ^ Press, W. H. (1978). "Flicker noises in astronomy and elsewhere". Comments in Astrophysics. 7 (4): 103–119. Bibcode:1978ComAp...7..103P.
  13. Handel, P. H.; Chung, A. L. (1993). Noise in Physical Systems and 1/"f" Fluctuations. New York: American Institute of Physics.
  14. Adak, Olgun; Rosenthal, Ethan; Meisner, Jeffery; Andrade, Erick F.; Pasupathy, Abhay N.; Nuckolls, Colin; Hybertsen, Mark S.; Venkataraman, Latha (2015-05-07). "Flicker Noise as a Probe of Electronic Interaction at Metal–Single Molecule Interfaces". Nano Letters. 15 (6): 4143–4149. Bibcode:2015NanoL..15.4143A. doi:10.1021/acs.nanolett.5b01270. ISSN 1530-6984. PMID 25942441.
  15. Bak, P.; Tang, C.; Wiesenfeld, K. (1987). "Self-Organized Criticality: An Explanation of 1/ƒ Noise". Physical Review Letters. 59 (4): 381–384. Bibcode:1987PhRvL..59..381B. doi:10.1103/PhysRevLett.59.381. PMID 10035754. S2CID 7674321.
  16. Josephson, Brian D. (1995). "A trans-human source of music?" in (P. Pylkkänen and P. Pylkkö, eds.) New Directions in Cognitive Science, Finnish Artificial Intelligence Society, Helsinki; pp. 280–285.
  17. ^ Gardner, M. (1978). "Mathematical Games—White and brown music, fractal curves and one-over-f fluctuations". Scientific American. 238 (4): 16–32. doi:10.1038/scientificamerican0478-16.
  18. Voss, R. F.; Clarke, J. (1975). "'1/f Noise' in Music and Speech". Nature. 258 (5533): 317–318. Bibcode:1975Natur.258..317V. doi:10.1038/258317a0. S2CID 4182664.
  19. ^ Voss, R. F.; Clarke, J. (1978). "1/f noise" in music: Music from 1/f noise". Journal of the Acoustical Society of America. 63 (1): 258–263. Bibcode:1978ASAJ...63..258V. doi:10.1121/1.381721.
  20. Voss, R.F. (May 1979). "1/F (Flicker) Noise: A Brief Review". 33rd Annual Symposium on Frequency Control. pp. 40–46. doi:10.1109/FREQ.1979.200297. S2CID 37302662.
  21. Schick, K. L.; Verveen, A. A. (October 1974). "1/f noise with a low frequency white noise limit". Nature. 251 (5476): 599–601. Bibcode:1974Natur.251..599S. doi:10.1038/251599a0. ISSN 1476-4687. S2CID 4200003.
  22. Hinkley, N.; Sherman, J. A.; Phillips, N. B.; Schioppo, M.; Lemke, N. D.; Beloy, K.; Pizzocaro, M.; Oates, C. W.; Ludlow, A. D. (2013-09-13). "An Atomic Clock with 10 –18 Instability". Science. 341 (6151): 1215–1218. arXiv:1305.5869. Bibcode:2013Sci...341.1215H. doi:10.1126/science.1240420. ISSN 0036-8075. PMID 23970562. S2CID 206549862.
  23. Vessot, Robert F. C. (1976-01-01), Meeks, M. L. (ed.), "5.4. Frequency and Time Standards††This work was supported in part by contract NSR 09-015-098 from the National Aeronautics and Space Administration.", Methods in Experimental Physics, Astrophysics, vol. 12, Academic Press, pp. 198–227, doi:10.1016/S0076-695X(08)60710-3, retrieved 2023-07-17
  24. Destexhe, Alain; Bédard, Claude (2020), "Local Field Potentials: LFP", in Jaeger, Dieter; Jung, Ranu (eds.), Encyclopedia of Computational Neuroscience, New York, NY: Springer, pp. 1–12, doi:10.1007/978-1-4614-7320-6_548-2, ISBN 978-1-4614-7320-6, S2CID 243735998, retrieved 2023-07-26
  25. Kerr, W.T.; et al. (2012). "Automated diagnosis of epilepsy using EEG power spectrum". Epilepsia. 53 (11): e189–e192. doi:10.1111/j.1528-1167.2012.03653.x. PMC 3447367. PMID 22967005.
  26. Douglas, PK; et al. (2019). "Reconsidering Spatial Priors in EEG Source Estimation : Does White Matter Contribute to EEG Rhythms?". 2019 7th International Winter Conference on Brain-Computer Interface (BCI). IEEE. pp. 1–12. arXiv:2111.08939. doi:10.1109/IWW-BCI.2019.8737307. ISBN 978-1-5386-8116-9. S2CID 195064621.
  27. Douglas, PK; Blair, G.; et al. (2024). "Towards a white matter ephaptic coupling model of 1/f spectra". 2024 12th International Winter Conference on Brain-Computer Interface (BCI). IEEE. pp. 1–3. doi:10.1109/BCI60775.2024.10480498.
  28. Van Orden, G.C.; Holden, J.G.; Turvey, M.T. (2003). "Self-organization of cognitive performance". Journal of Experimental Psychology: General. 132 (3): 331–350. doi:10.1037/0096-3445.132.3.331. PMID 13678372.
  29. Pareyon, G. (2011). On Musical Self-Similarity, International Semiotics Institute & University of Helsinki. "On Musical Self-Similarity" (PDF).
  30. "Noise in Man-generated Images and Sound".
  31. Anger, Natalie (March 1, 2010). "Bringing New Understanding to the Director's Cut". The New York Times. Retrieved on March 3, 2010. See also original study Archived 2013-01-24 at the Wayback Machine
  32. Gilden, David L; Thornton, T; Mallon, MW (1995). "1/ƒ Noise in Human Cognition". Science. 267 (5205): 1837–1839. Bibcode:1995Sci...267.1837G. doi:10.1126/science.7892611. ISSN 0036-8075. PMID 7892611.
  33. Gilden, D. L. (1997). "Fluctuations in the time required for elementary decisions". Psychological Science. 8 (4): 296–301. doi:10.1111/j.1467-9280.1997.tb00441.x. S2CID 145051976.
  34. Gilden, David L (2001). "Cognitive Emissions of 1/ƒ Noise". Psychological Review. 108 (1): 33–56. CiteSeerX 10.1.1.136.1992. doi:10.1037/0033-295X.108.1.33. ISSN 0033-295X. PMID 11212631.
  35. Weissman, M. B. (1988). "1/ƒ Noise and other slow non-exponential kinetics in condensed matter". Reviews of Modern Physics. 60 (2): 537–571. Bibcode:1988RvMP...60..537W. doi:10.1103/RevModPhys.60.537.
  36. Dutta, P. & Horn, P. M. (1981). "Low-frequency fluctuations in solids: 1/f noise". Reviews of Modern Physics. 53 (3): 497–516. Bibcode:1981RvMP...53..497D. doi:10.1103/RevModPhys.53.497.
  37. Kleinpenning, T. G. M. & de Kuijper, A. H. (1988). "Relation between variance and sample duration of 1/f Noise signals". Journal of Applied Physics. 63 (1): 43. Bibcode:1988JAP....63...43K. doi:10.1063/1.340460.
  38. Kleinpenning, T. G. M.; de Kuijper, A. H. (1988-01-01). "Relation between variance and sample duration of 1/ f noise signals". Journal of Applied Physics. 63 (1): 43–45. Bibcode:1988JAP....63...43K. doi:10.1063/1.340460. ISSN 0021-8979.
  39. Aldert van der Ziel, (1954), Noise, Prentice–Hall
  40. Hei Wong (2003). "Low-frequency noise study in electron devices: review and update". Microelectronics Reliability. 43 (4): 585–599. Bibcode:2003MiRe...43..585W. doi:10.1016/S0026-2714(02)00347-5.
  41. Alexander A. Balandin (2013). "Low-frequency 1/f noise in graphene devices". Nature Nanotechnology. 8 (8): 549–555. arXiv:1307.4797. Bibcode:2013NatNa...8..549B. doi:10.1038/nnano.2013.144. PMID 23912107. S2CID 16030927.
  42. Smulko, Janusz; Scandurra, Graziella; Drozdowska, Katarzyna; Kwiatkowski, Andrzej; Ciofi, Carmine; Wen, He (2024). "Flicker Noise in Resistive Gas Sensors—Measurement Setups and Applications for Enhanced Gas Sensing". Sensors. 24 (2): 405. Bibcode:2024Senso..24..405S. doi:10.3390/s24020405. PMC 10821460. PMID 38257498.
  43. Moore, Christopher; Cole, Robert; Berry, Christopher (19 July 2013). "Gravitational Wave Detectors and Sources". Retrieved 17 April 2014.
  44. Jim Shelton (2018-09-04). "Think pink for a better view of climate change". YaleNews. Retrieved 5 September 2018.
  45. Moon, Woosok; Agarwal, Sahil; Wettlaufer, J. S. (2018-09-04). "Intrinsic Pink-Noise Multidecadal Global Climate Dynamics Mode". Physical Review Letters. 121 (10): 108701. arXiv:1802.00392. Bibcode:2018PhRvL.121j8701M. doi:10.1103/PhysRevLett.121.108701. PMID 30240245. S2CID 52243763.
  46. Norton, M. P. (2003). Fundamentals of noise and vibration analysis for engineers. Karczub, D. G. (Denis G.) (2nd ed.). Cambridge, UK: Cambridge University Press. ISBN 9780511674983. OCLC 667085096.
  47. Krapf, Diego; Marinari, Enzo; Metzler, Ralf; Oshanin, Gleb; Xu, Xinran; Squarcini, Alessio (2018-02-09). "Power spectral density of a single Brownian trajectory: what one can and cannot learn from it". New Journal of Physics. 20 (2): 023029. arXiv:1801.02986. Bibcode:2018NJPh...20b3029K. doi:10.1088/1367-2630/aaa67c. ISSN 1367-2630.
  48. Krapf, Diego; Lukat, Nils; Marinari, Enzo; Metzler, Ralf; Oshanin, Gleb; Selhuber-Unkel, Christine; Squarcini, Alessio; Stadler, Lorenz; Weiss, Matthias; Xu, Xinran (2019-01-31). "Spectral Content of a Single Non-Brownian Trajectory". Physical Review X. 9 (1): 011019. arXiv:1902.00481. Bibcode:2019PhRvX...9a1019K. doi:10.1103/PhysRevX.9.011019. ISSN 2160-3308.
  49. ^ Kendal WS, Jørgensen BR (2011). "Tweedie convergence: a mathematical basis for Taylor's power law, 1/f noise and multifractality" (PDF). Phys. Rev. E. 84 (6): 066120. Bibcode:2011PhRvE..84f6120K. doi:10.1103/physreve.84.066120. PMID 22304168.
  50. Jørgensen B, Martinez JR, Tsao M (1994). "Asymptotic behaviour of the variance function". Scandinavian Journal of Statistics. 21: 223–243.
  51. Taylor LR (1961). "Aggregation, variance and the mean". Nature. 189 (4766): 732–735. Bibcode:1961Natur.189..732T. doi:10.1038/189732a0. S2CID 4263093.
  52. Eisler Z, Bartos I, Kertesz J (2008). "Fluctuation scaling in complex systems: Taylor's law and beyond". Advances in Physics. 57 (1): 89–142. arXiv:0708.2053. Bibcode:2008AdPhy..57...89E. doi:10.1080/00018730801893043. S2CID 119608542.
  53. Kendal WS (2015). "Self-organized criticality attributed to a central limit-like convergence effect". Physica A. 421: 141–150. Bibcode:2015PhyA..421..141K. doi:10.1016/j.physa.2014.11.035.
  54. Milotti, Edoardo (2002-04-12). "1/f noise: a pedagogical review". arXiv:physics/0204033.
  55. Shapoval, Alexander; Shnirman, Mikhail (2024). "Explanation of flicker noise with the Bak-Tang-Wiesenfeld model of self-organized criticality". Physical Review E. 110: 014106. arXiv:2212.14726. doi:10.1103/PhysRevE.110.014106.
  56. O'Brien, Kevin P.; Weissman, M. B. (1992-10-01). "Statistical signatures of self-organization". Physical Review A. 46 (8): R4475–R4478. Bibcode:1992PhRvA..46.4475O. doi:10.1103/PhysRevA.46.R4475. PMID 9908765.
  57. "Noise in Man-generated Images and Sound". mlab.uiah.fi. Retrieved 2015-11-14.
  58. "DSP Generation of Pink Noise". www.firstpr.com.au. Retrieved 2015-11-14.
  59. McClain, D. (May 1, 2001). "Numerical Simulation of Pink Noise" (PDF). Preprint. Archived from the original (PDF) on 2011-10-04.
  60. Timmer, J.; König, M. (1995-01-01). "On Generating Power Law Noise". Astronomy and Astrophysics. 300: 707–710. Bibcode:1995A&A...300..707T.
  61. Ovchinnikov, I.V. (2016). "Introduction to supersymmetric theory of stochastics". Entropy. 18 (4): 108. arXiv:1511.03393. Bibcode:2016Entrp..18..108O. doi:10.3390/e18040108. S2CID 2388285.
  62. Ovchinnikov, I.V.; Schwartz, R. N.; Wang, K. L. (2016). "Topological supersymmetry breaking: Definition and stochastic generalization of chaos and the limit of applicability of statistics". Modern Physics Letters B. 30 (8): 1650086. arXiv:1404.4076. Bibcode:2016MPLB...3050086O. doi:10.1142/S021798491650086X. S2CID 118174242.
  63. Loar, Josh (2019). The Sound System Design Primer. Routledge. pp. 274–276. ISBN 9781351768184.
  64. Eckstein, Matt (30 August 2018). "Sound System Commissioning - Say What?". AE Design. Retrieved November 22, 2022.
  65. Cox, Tyler. "What is Pink Noise and What Does It Do?". Yamaha Insights. Yamaha Pro Audio. Retrieved November 22, 2022.
  66. Lacanette, Kerry (1990). "Create an Accurate Noise Generator". Electronic Design. Vol. 38. Hayden. p. 108.
  67. Thomas, Christian (April 30, 2021). "Headphone burn-in isn't real". Soundguys. Retrieved November 22, 2022.
  1. Though in practice, since there are no ideal clocks, t {\displaystyle t} is actually the ticks of a much more accurate clock.

References

External links

Noise (physics and telecommunications)
General
Noise in...
Class of noise
Engineering
terms
Ratios
Related topics
Denoise
methods
General
2D (Image)
Categories: