Calcul DSP à partir de la FFT
Utilisez ce calculateur premium pour convertir des amplitudes de FFT en densité spectrale de puissance (DSP, ou PSD en anglais) en W/Hz et dB/Hz. L’outil gère l’échantillonnage, la taille FFT, la correction liée à la fenêtre et le choix d’un spectre mono-côté ou bi-côté.
Calculateur DSP
Entrez les magnitudes FFT séparées par des virgules. Pour un spectre mono-côté, fournissez généralement les bins de 0 à N/2 inclus. Pour un spectre bi-côté, fournissez tous les bins.
Résultats
Saisissez vos paramètres puis cliquez sur “Calculer la DSP”.
Comprendre le calcul de la DSP à partir de la FFT
Le calcul de la densité spectrale de puissance à partir de la FFT est une opération fondamentale en traitement du signal, en vibration, en acoustique, en électronique, en radar, en communication numérique et en instrumentation. En pratique, de nombreux ingénieurs obtiennent une FFT d’un signal mesuré, mais veulent aller plus loin qu’un simple tracé d’amplitude. Ils cherchent à exprimer l’énergie du signal par unité de bande passante, ce qui correspond à la DSP, notée souvent PSD en anglais pour Power Spectral Density. Cette grandeur s’exprime typiquement en W/Hz, V²/Hz ou encore en dB/Hz selon le contexte de mesure.
La FFT, à elle seule, fournit une représentation discrète du contenu fréquentiel. Cependant, pour passer d’une magnitude FFT brute à une DSP exploitable, il faut tenir compte de plusieurs facteurs: la fréquence d’échantillonnage, la longueur de la FFT, la largeur de bande d’un bin, le type de fenêtre utilisé et la convention de représentation mono-côté ou bi-côté. Sans ces corrections, les comparaisons entre acquisitions, fenêtres ou tailles FFT peuvent être trompeuses.
Pourquoi la DSP est plus utile qu’une simple amplitude FFT
Une FFT en amplitude est très utile pour repérer des pics, des harmoniques ou des composantes parasites, mais elle dépend fortement du nombre d’échantillons et du fenêtrage. La DSP normalise l’énergie par Hertz. Cela signifie qu’un bruit large bande mesuré avec une FFT de 1024 points ou de 8192 points peut être comparé de manière plus rigoureuse, car la densité tient compte de la bande fréquentielle représentée par chaque bin. C’est particulièrement important dans les analyses de bruit électronique, de bruit vibratoire et de bruit de phase.
- Elle rend les niveaux de bruit comparables entre différentes résolutions fréquentielles.
- Elle permet l’intégration sur une bande pour retrouver une puissance totale.
- Elle améliore l’interprétation physique d’un spectre mesuré.
- Elle est adaptée aux signaux aléatoires, où une simple amplitude ne suffit pas.
Formule générale utilisée dans le calculateur
Pour un signal discret de longueur N, échantillonné à Fs et transformé par FFT non normalisée, une estimation courante de la DSP bi-côté s’écrit:
Sxx[k] = |X[k]|² / (Fs × N × U)
où |X[k]| est la magnitude du bin k, et U représente le facteur de puissance moyen de la fenêtre. Pour une fenêtre rectangulaire, ce facteur vaut 1. Pour d’autres fenêtres, il devient inférieur à 1 parce qu’une partie de l’énergie temporelle est atténuée par le fenêtrage. Dans une représentation mono-côté pour signal réel, les bins internes sont multipliés par 2 afin de conserver la puissance totale lorsque l’on replie la partie négative sur la partie positive.
Rôle de la fréquence d’échantillonnage et de la résolution fréquentielle
La fréquence d’échantillonnage détermine la fréquence maximale observable sans repliement, généralement jusqu’à Fs/2 pour un signal réel. La largeur de bande d’un bin FFT est égale à Δf = Fs / N. Cette quantité est centrale pour comprendre la relation entre puissance et densité. Une puissance de bin devient une densité lorsqu’elle est rapportée à cette largeur de bande. C’est d’ailleurs l’une des raisons pour lesquelles une simple augmentation de N modifie l’apparence d’un bruit dans le spectre brut, alors que la DSP reste conceptuellement stable.
| Paramètre | Formule | Exemple avec Fs = 1024 Hz | Interprétation |
|---|---|---|---|
| Résolution fréquentielle avec N = 1024 | Δf = Fs / N | 1 Hz | Chaque bin couvre 1 Hz de bande |
| Résolution fréquentielle avec N = 4096 | Δf = Fs / N | 0,25 Hz | Meilleure finesse de séparation fréquentielle |
| Fréquence de Nyquist | Fs / 2 | 512 Hz | Limite haute du spectre utile pour un signal réel |
| Durée d’observation avec N = 4096 | N / Fs | 4 s | Observation plus longue, meilleure résolution en fréquence |
Impact du choix de fenêtre sur la DSP
Lorsqu’un signal ne contient pas un nombre entier de périodes dans la fenêtre d’analyse, la fuite spectrale apparaît. C’est pour limiter cette fuite que l’on applique une fenêtre de type Hann, Hamming ou Blackman. En contrepartie, la fenêtre modifie le niveau de puissance et la largeur du lobe principal. Une bonne estimation DSP doit donc intégrer une correction énergétique liée à la fenêtre. C’est précisément ce que fait ce calculateur grâce au facteur U.
Les valeurs suivantes sont très utilisées en pratique pour la puissance moyenne de fenêtre:
| Fenêtre | Facteur de puissance U | ENBW approximative en bins | Usage courant |
|---|---|---|---|
| Rectangulaire | 1,0000 | 1,00 | Bonne résolution, fuite spectrale élevée |
| Hann | 0,3750 | 1,50 | Compromis très populaire en analyse générale |
| Hamming | 0,3974 | 1,36 | Bonne réduction des lobes secondaires |
| Blackman | 0,3046 | 1,73 | Excellente réduction de fuite, lobe principal plus large |
Mono-côté ou bi-côté: quelle convention choisir?
Pour un signal réel, le spectre FFT bi-côté est symétrique. Dans ce cas, il est souvent plus pratique de représenter uniquement les fréquences positives, donc un spectre mono-côté. Cependant, si l’on fait cela, il faut veiller à conserver la puissance totale. C’est pourquoi les bins internes sont multipliés par 2, à l’exception du bin DC et du bin de Nyquist lorsqu’il existe. Si vous travaillez sur des signaux complexes I/Q, il est généralement préférable de rester en bi-côté.
- Choisissez bi-côté pour les signaux complexes ou pour une analyse mathématique complète.
- Choisissez mono-côté pour les signaux réels et une lecture plus intuitive du spectre.
- Appliquez la correction de facteur 2 uniquement aux bins internes du spectre mono-côté.
- Vérifiez toujours la cohérence entre le nombre de bins saisis et le mode choisi.
Exemple concret de calcul DSP à partir d’une FFT
Supposons un signal échantillonné à 1024 Hz, avec une FFT de longueur 8, sous fenêtre de Hann. Vous obtenez les magnitudes mono-côté suivantes: 0,15, 1,20, 3,80, 1,10, 0,30. La résolution fréquentielle vaut alors 128 Hz. Le facteur énergétique de Hann est 0,375. Pour chaque bin, le calculateur élève la magnitude au carré, divise par Fs × N × U, puis applique la correction mono-côté sur les bins internes. Le résultat produit une DSP par fréquence positive, ainsi qu’un affichage en dB/Hz pour faciliter les comparaisons logarithmiques.
Ce procédé est particulièrement utile lorsqu’on veut identifier si un pic dominant correspond à une raie quasi déterministe ou s’il s’inscrit dans un bruit large bande. En intégrant la DSP sur l’ensemble des bins, on retrouve aussi une estimation de la puissance totale associée au spectre saisi. Cela aide à contrôler la cohérence énergétique de l’analyse.
Erreurs fréquentes à éviter
- Confondre amplitude FFT, puissance de bin et densité spectrale de puissance.
- Oublier de corriger l’effet de la fenêtre sur l’énergie.
- Comparer des spectres ayant des résolutions fréquentielles différentes sans passer par la DSP.
- Multiplier tous les bins par 2 dans un spectre mono-côté, y compris DC et Nyquist.
- Utiliser une FFT de signal complexe comme si elle était un signal réel.
- Travailler avec des unités incohérentes, par exemple V au lieu de V²/Hz.
Comment interpréter un résultat en dB/Hz
L’échelle en dB/Hz s’obtient via 10 × log10(DSP). Elle comprime les écarts dynamiques importants et permet une lecture plus naturelle lorsque plusieurs composantes diffèrent de plusieurs ordres de grandeur. Si la DSP est exprimée en V²/Hz, le résultat en dB/Hz reste relatif à cette même base. Dans des contextes RF, on rencontre aussi des variantes comme dBm/Hz, mais cela suppose une référence de puissance explicite, souvent 1 mW sur une impédance donnée.
Cas d’usage industriels et scientifiques
En vibration mécanique, la DSP aide à caractériser l’énergie vibratoire répartie selon la fréquence, ce qui permet de comparer un environnement aléatoire à une spécification d’essai. En électronique, elle sert à quantifier le bruit de fond d’un capteur, d’un amplificateur ou d’un convertisseur. En acoustique, elle facilite l’analyse de bruit large bande et la distinction entre tonalités et composantes aléatoires. En télécommunications, elle intervient dans l’étude du bruit de canal, des interférences et des signaux modulés.
Références fiables pour approfondir
Si vous souhaitez aller plus loin, consultez des sources institutionnelles et universitaires reconnues. Les documents pédagogiques de la Rice University expliquent clairement la relation entre puissance et densité spectrale. La National Institute of Standards and Technology met à disposition des ressources de référence sur les mesures fréquentielles et les bonnes pratiques métrologiques. Pour les fondements du traitement numérique du signal, les supports de cours du MIT OpenCourseWare restent également très pertinents.
Bonnes pratiques pour obtenir une DSP fiable
- Choisissez une fréquence d’échantillonnage adaptée à la bande fréquentielle utile.
- Assurez-vous qu’un éventuel filtrage anti-repliement est correctement appliqué en amont.
- Sélectionnez une fenêtre cohérente avec votre objectif: résolution fine ou faible fuite spectrale.
- Documentez toujours la convention de normalisation FFT utilisée.
- Travaillez avec des moyennes de plusieurs FFT pour stabiliser l’estimation sur les signaux bruités.
- Vérifiez la cohérence entre puissance intégrée et énergie attendue dans le domaine temporel.
Conclusion
Le calcul de la DSP à partir de la FFT n’est pas seulement une transformation mathématique de plus. C’est la passerelle entre une visualisation fréquentielle brute et une mesure physiquement interprétable. En tenant compte de la taille FFT, de la fréquence d’échantillonnage, du choix de fenêtre et du mode mono-côté ou bi-côté, vous pouvez convertir des magnitudes FFT en densités spectrales cohérentes, comparables et exploitables. Ce calculateur a été conçu pour vous fournir un résultat immédiat, un graphique lisible et une méthode de calcul transparente afin de sécuriser vos analyses de bruit, de vibration ou de spectre.