PSD (méthode de Welch - par moyennage).

Espace de nom : tsd::fourier

Prototype

tuple<Vecf, Vecf> psd_welch(const Veccf &x, entier N, cstring fen="hn")

Paramètres

xSignal à analyser.
NDimension des sous-fenêtres (autrement dit, la résolution de la PSD résultante).
fenChoix du fenêtrage.

Retourne

Un tuple de deux éléments : vecteur de fréquences (normalisées, entre -0.5 et 0.5), et le spectre.

Description

Cette fonction est une estimation de la PSD par la technique du moyennage : plusieur PSD sont calculées sur des fenêtres (pondérées et avec du recouvrement), puis moyennées.

Exemple

  soit n = 8 * 1024;
  soit x = sigexp(0.02, n) + 0.1 * randn(n);
  soit [fr1, S1] = psd(x);
  soit [fr2, S2] = psd_welch(x, 1024);
  soit [fr3, S3] = psd_welch(x, 256);

  Figures f;
  f.subplot().plot(x);
  f.gcf().titres("Signal");
  f.subplot().plot(fr1, S1);
  f.gcf().titres("PSD - sans moyennage", "Frequence (Hz)", "Magnitude (dB)");
  f.subplot().plot(fr2, S2);
  f.gcf().titres("PSD - moyennage par blocs de 1024 pts", "Frequence (Hz)", "Magnitude (dB)");
  f.subplot().plot(fr3, S3);
  f.gcf().titres("PSD - moyennage par blocs de 256 pts", "Frequence (Hz)", "Magnitude (dB)");

PSD et moyennage :

Voir aussi

psd(), psd_sousesp()