File:Colored noise.png

DescriptionColored noise.png	English: Example of noise signals with different power spectra.
Date	16 juillet 2019
Source	https://twitter.com/j_bertolotti/status/1151129795008200707
Auteur	Jacopo Bertolotti
Autorisation (Réutilisation de ce fichier)	https://twitter.com/j_bertolotti/status/1030470604418428929
Autres versions	Ce fichier est la source d’une image extraite : Colored noise (cropped).png.

n = 4096;
\[Phi] = RandomReal[{0, 2 \[Pi]}, n];
\[Epsilon][i_, \[Alpha]_] := Sum[Sqrt[k^-\[Alpha] Abs[(2 \[Pi])/n]^(1 - \[Alpha])]*Cos[(2 \[Pi] i k)/n + \[Phi][[k]]], {k, 1, n/2}];
red = Table[\[Epsilon][i, 2.], {i, 1, n}];
white = Table[\[Epsilon][i, 0.], {i, 1, n}];
pink = Table[\[Epsilon][i, 1.], {i, 1, n}];
blue = Table[\[Epsilon][i, -1.], {i, 1, n}];
GraphicsGrid[{{
ListPlot[blue, Joined -> True, PlotStyle -> LightBlue, LabelStyle -> {LightBlue, Bold}, AxesLabel -> {"t", ""}, Ticks -> None, Axes -> {True, False}, PlotLabel -> "Blue noise"],
ListLogLogPlot[Evaluate[(Abs[Fourier[blue]]^2)[[2 ;; n/2]]], Joined -> True, PlotRange -> All, PlotStyle -> {LightBlue, Thick}, Axes -> False, Frame -> True, FrameLabel -> {"\[Omega]", "Power spectrum"}, LabelStyle -> {LightBlue, Bold}, GridLines -> Automatic]
    }, {
ListPlot[white, Joined -> True, PlotStyle -> White, LabelStyle -> {White, Bold}, AxesLabel -> {"t", ""}, Ticks -> None, Axes -> {True, False}, PlotLabel -> "White noise"], ListLogLogPlot[Evaluate[(Abs[Fourier[white]]^2)[[2 ;; n/2]]], Joined -> True, PlotRange -> All, PlotStyle -> {White, Thick}, Axes -> False, Frame -> True, FrameLabel -> {"\[Omega]", "Power spectrum"}, 
LabelStyle -> {White, Bold}, GridLines -> Automatic]
    }, {
ListPlot[pink, Joined -> True, PlotStyle -> Pink, LabelStyle -> {Pink, Bold}, AxesLabel -> {"t", ""}, Ticks -> None, Axes -> {True, False}, PlotLabel -> "Pink Noise"], ListLogLogPlot[Evaluate[(Abs[Fourier[pink]]^2)[[2 ;; n/2]]], Joined -> True, PlotRange -> All, PlotStyle -> {Pink, Thick}, Axes -> False, Frame -> True, FrameLabel -> {"\[Omega]", "Power spectrum"}, 
LabelStyle -> {Pink, Bold}, GridLines -> Automatic]
    }, {
ListPlot[red, Joined -> True, PlotStyle -> Red, LabelStyle -> {Pink, Bold}, AxesLabel -> {"t", ""}, Ticks -> None, Axes -> {True, False}, PlotLabel -> "Red Noise"],
ListLogLogPlot[Evaluate[(Abs[Fourier[red]]^2)[[2 ;; n/2]]], Joined -> True, PlotRange -> All, PlotStyle -> {Red, Thick}, Axes -> False, Frame -> True, FrameLabel -> {"\[Omega]", "Power spectrum"}, 
    LabelStyle -> {Red, Bold}, GridLines -> Automatic]
    }}, ImageSize -> Large, Background -> Black]

Moi, en tant que détenteur des droits d’auteur sur cette œuvre, je la publie sous la licence suivante :

Ce fichier est disponible selon les termes de la licence Creative Commons CC0 Don universel au domaine public.

La personne qui a associé une œuvre avec cet acte l’a placée dans le domaine public en renonçant mondialement à tous ses droits sur cette œuvre en vertu des lois relatives au droit d’auteur, ainsi qu’à tous les droits juridiques connexes et voisins qu’elle possédait sur l’œuvre, sans autre limite que celles imposées par la loi. Vous pouvez copier, modifier, distribuer et utiliser cette œuvre, y compris à des fins commerciales, sans qu’il soit nécessaire d’en demander la permission. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse