Guide d’analyse en bandes d’octave : regroupement FFT (binning) vs méthode par banc de filtres

L'analyse en bandes d'octave peut être mise en œuvre de deux façons fondamentalement différentes : le regroupement de raies FFT (intégration des raies de PSD/FFT en bandes d'octave entières 1/1 et fractionnaires 1/3) et un véritable banc de filtres en bandes d'octave (filtres passe‑bande orientés normes + moyennage RMS/Leq). Dans cet article, nous comparons le fonctionnement de ces deux méthodes, les cas où leurs résultats coïncident, ceux où ils divergent (échelle, ENBW de fenêtre, pondération aux bords de bande, latence, réponse aux transitoires), et la manière dont OpenTest prend en charge les deux approches pour l'acoustique, le NVH et les mesures de conformité.

Pour une explication détaillée des concepts, lisez ceci → Analyse en bandes d'octave : fondements mathématiques et ingénierie

Bancs de filtres en bandes d'octave (véritable banc d'octave / banc de filtres CPB)

Filtres passe‑bande parallèles + détecteur d'énergie + moyennage temporel

Un analyseur par banc de filtres (véritable octave) réalise généralement :

1. Concevoir un filtre passe‑bande H_b(z) (ou H_b(s)) pour chaque fréquence centrale de bande.
2. Faire fonctionner les filtres en parallèle pour obtenir les signaux de bande y_b(t).
3. Calculer la puissance / la valeur quadratique moyenne de la bande et appliquer un moyennage temporel pour produire les niveaux de bande en sortie.

Pour être comparables entre instruments, les réponses en module des filtres doivent respecter les masques de tolérance (classe) IEC/ANSI pour le jeu de filtres spécifié. [1][3]

IIR vs FIR : pourquoi les IIR (biquads en cascade) sont courants en pratique

• Avantages des IIR : ordre plus faible pour une pente donnée, calcul réduit, adaptés au temps réel / systèmes embarqués ; stables lorsqu'ils sont mis en œuvre en SOS/biquads.
• Avantages des FIR : une phase linéaire est possible (utile lorsque la forme d'onde importe) ; la conception et la vérification peuvent être plus simples.

Pour les sorties en niveaux de bande, la phase n'est généralement pas la principale préoccupation, de sorte que les bancs de filtres IIR sont courants.

Traitement multirate : « l'arme secrète » des bancs de filtres CPB

Les bandes CPB de basses fréquences sont très étroites. Les implémenter à la fréquence d'échantillonnage complète est inefficace. Une stratégie courante consiste à regrouper les bandes par octave et à sous‑échantillonner pour les groupes de basses fréquences :

• Effectuer un passe‑bas puis décimer (par ex. par 2 par octave) pour les groupes de basses fréquences.
• Implémenter les filtres passe‑bande correspondants à la fréquence d'échantillonnage réduite.
• Assurer un anti‑repliement (anti‑aliasing) adéquat avant la décimation.

Moyennage temporel / pondération temporelle : les niveaux de bande sont des statistiques, pas des valeurs instantanées

Les niveaux de bande nécessitent généralement un moyennage temporel. Les options courantes incluent le RMS par blocs, le moyennage exponentiel ou le Leq (niveau d'énergie équivalente). Dans le contexte des sonomètres, la norme IEC 61672‑1 définit les pondérations temporelles Rapide / Lente (Rapide ~125 ms, Lente ~1 s). [5][6]

Conséquence pour l'ingénierie : des constantes de temps différentes produisent des lectures différentes, la pondération temporelle doit donc être précisée dans les rapports.

Comment valider qu'un banc de filtres se comporte « comme la norme »

• Balayage sinusoïdal : vérifier le comportement en bande passante et l'isolation des bandes adjacentes ; observer les effets de délai temporel.
• Bruit rose/blanc : vérifier les niveaux de bande moyens et la variance / le temps de stabilisation ; contrôler le comportement de la bande passante effective.
• Impulsion / échelon : examiner la sonnerie (ringing) et la réponse temporelle (critique pour les transitoires).
• Vérifier par recoupement avec un instrument / une implémentation de référence conforme connue.

Des définitions de bandes aux filtres numériques conformes : un flux de travail de bout en bout (conceptuel)

1. Choisir le système de bandes : base 10 / base 2, la fraction 1/b (généralement b=3), générer les fréquences centrales exactes fm et les bords de bande f1/f2.
2. Choisir l'objectif de performance : quelle édition de norme et quelle classe / quel masque de tolérance ?
3. Choisir la structure de filtre : IIR en SOS pour le temps réel ; FIR ou filtrage aller‑retour si une phase nulle / à phase linéaire est requise.
4. Concevoir chaque passe‑bande : projeter correctement f1/f2 dans le domaine numérique (par ex. pré‑déformation pour la transformée bilinéaire).
5. Mettre en œuvre le multirate si nécessaire : décimer pour les groupes de basses fréquences avec un filtrage anti‑repliement suffisant.
6. Vérifier : réponse en module par rapport au masque ; essais sur bruit pour la bande passante effective ; balayages / impulsions pour la réponse temporelle.
7. Calibrer et documenter : unités et grandeurs de référence, moyennage / pondération temporelle, détails de la méthode.

Réponse temporelle expliquée : retard de groupe, sonnerie et moyennage façonnent tous les lectures

Un analyseur de niveaux par bande est un système en domaine temporel (filtre → détecteur d'énergie → lisseur), de sorte que les lectures sont régies par plusieurs échelles de temps :

• Retard de groupe du filtre : à quel retard les événements apparaissent‑ils dans chaque bande.
• Sonnerie / décroissance du filtre : pendant combien de temps une impulsion courte « résonne » dans une bande.
• Moyennage d'énergie / pondération temporelle : compromis entre résolution temporelle et fluctuation du niveau de sortie.

Ainsi, pour les transitoires (impacts, démarrages/arrêts, balayages), différentes implémentations conformes peuvent produire des niveaux de crête et des évolutions temporelles différents — ce qui est cohérent avec la mise en garde de l'ANSI. [3]

Règle empirique : pour les contributions en régime établi, utiliser un moyennage plus long pour plus de stabilité ; pour la localisation de transitoires, raccourcir le moyennage mais accepter une plus grande variabilité et figer précisément les détails de l'algorithme.

Pièges courants en temps réel

• Oublier l'anti‑repliement dans la chaîne de décimation : les bandes de basses fréquences deviennent contaminées par l'aliasing.
• Instabilité numérique des sections IIR à fort Q et basses fréquences : utiliser des SOS/biquads et une précision suffisante.
• Moyennage en dB : toujours moyenner dans le domaine de l'énergie / de la valeur quadratique moyenne, puis convertir en dB.
• Supposer que les énergies de bande doivent se sommer exactement à l'énergie totale : les filtres normalisés ne sont pas nécessairement complémentaire en puissance ; vérifier plutôt selon des critères cohérents avec les normes.

Analyse en bancs de filtres en bandes d'octave dans OpenTest

OpenTest prend en charge l'analyse en bandes d'octave en utilisant une approche par banc de filtres :
1) Connectez le dispositif, par exemple SonoDAQ Pro
2) Sélectionnez les voies et ajustez les paramètres. Pour un microphone externe, activez l'IEPE et passez en mode de mesure de signal acoustique.
3) Dans la section Octave-Band Analysis sous Mode de mesure, choisissez l'algorithme conforme à l'IEC 61260‑1. Il prend en charge l'analyse en temps réel, le moyennage linéaire, le moyennage exponentiel et la fonction de maintien de crête (peak hold).
4) Après avoir configuré les paramètres, cliquez sur le bouton Test pour démarrer la mesure.
5) Un enregistrement unique peut être analysé simultanément en bandes 1/1 d'octave, 1/3 d'octave, 1/6 d'octave, 1/12 d'octave, 1/24 d'octave et 1/24 d'octave.

Figure 1 : analyse en bancs de filtres en bandes d'octave dans OpenTest

Regroupement FFT (FFT binning) et synthèse FFT

Regroupement FFT : convertir un spectre étroit en intégrales de bandes CPB

1. Estimer le spectre (FFT unique, PSD de Welch ou STFT).
2. Intégrer / sommer à l'intérieur de chaque bande d'octave / de fraction d'octave pour obtenir la puissance de bande.

Cette approche est courante dans les traitements logiciels / hors ligne, car une seule FFT fournit un spectre à haute résolution qui peut être re‑regroupé dans n'importe quel système de bandes (1/1, 1/3, 1/12, …).

Défi clé n°1 : mise à l'échelle FFT et corrections de fenêtre

Après une FFT, la mise à l'échelle dépend de vos définitions : normalisation 1/N, amplitude vs puissance vs PSD, spectre simple face vs double face, et fenêtrage. Pour les mesures de bruit, l'ENBW est cruciale ; l'ignorer peut introduire des décalages systématiques. [7]

Une normalisation PSD pratique (forme périodogramme)

# convertir en PSD simple face : multiplier par 2 sauf DC (et Nyquist si présent)

Cela donne une PSD en unités (unité d'entrée)²/Hz et permet des vérifications de cohérence énergétique en intégrant la PSD sur la fréquence.

Deux auto‑contrôles rapides pour la mise à l'échelle

• Test sur bruit blanc : générer un bruit de variance connue σ² ; intégrer la PSD simple face sur 0..fs/2 et retrouver ≈σ² (en tenant compte du facteur ×2).
• Test sur tonalité pure : générer un sinus d'amplitude A (RMS=A/√2) ; l'intégration de l'énergie spectrale doit retrouver ≈A²/2 (sous réserve de la fuite et du choix de la fenêtre).

Si les deux tests sont satisfaits, votre mise à l'échelle FFT est probablement correcte ; le pondération partielle des raies et le regroupement en bandes d'octave deviennent alors pertinents.

Défi clé n°2 : les bords de bande coïncident rarement avec les raies → pondération partielle des raies

Des décisions tout‑ou‑rien aux bords de bande provoquent des erreurs en escalier, particulièrement aux basses fréquences où les bandes sont étroites. Utiliser une pondération basée sur le chevauchement (section 4.2.4) pour les raies situées aux frontières de bande.

La mise à zéro (zero‑padding) résout‑elle le problème de désalignement des bords ? (idée reçue courante)

Le zero‑padding interpole le spectre affiché mais n'améliore pas la vraie résolution fréquentielle (qui est fixée par la longueur de la fenêtre d'origine). Il peut réduire l'effet « en escalier » visuel mais ne peut pas transformer des bandes de basses fréquences ne contenant qu'1–2 raies en estimations fiables de niveaux de bande. Les solutions fondamentales sont des fenêtres plus longues ou le traitement multirate / les bancs de filtres.

Défi clé n°3 : compromis temps–fréquence (la longueur de fenêtre fixe la précision en basses fréquences et le délai)

La résolution FFT est Δf = fs/N. Les bandes 1/3 d'octave de basses fréquences peuvent ne faire que quelques hertz de large, de sorte qu'obtenir suffisamment de raies par bande nécessite un N très grand, ce qui augmente la latence et lissage des transitoires.

Cause fondamentale : le 1/3 d'octave est à Q constant, alors que la STFT utilise des raies de pas fréquentiel constant Δf

En CPB, la largeur de bande est proportionnelle à la fréquence (Δf_bande ∝ f, Q constant). En STFT, l'espacement des raies est constant (Δf_raie constant). Par conséquent, aux basses fréquences CPB nécessite un Δf_raie extrêmement fin (fenêtres longues), tandis qu'aux hautes fréquences la résolution est surabondante.

Pistes de solution : STFT à fenêtre longue vs STFT multirate vs CQT / ondelettes

• STFT à fenêtre longue : la plus simple, mais avec une latence élevée et un lissage important des transitoires.
• STFT multirate : sous‑échantillonner le contenu de basses fréquences et appliquer la FFT à une fréquence fs plus faible, de manière analogue aux bancs de filtres multirate.
• Transformée à Q constant (CQT) / ondelettes : résolution naturellement logarithmique, mais l'adaptation aux masques IEC/ANSI exige une calibration / validation supplémentaire. [4]

Pour les mesures de conformité, les bancs de filtres orientés normes sont préférés ; pour la recherche / l'extraction de descripteurs, la CQT / les ondelettes peuvent être attractives.

Synthèse FFT : construction d'un filtrage par bande dans le domaine fréquentiel

La synthèse FFT rapproche l'approche FFT d'un banc de filtres :

• Définir un poids fréquentiel W_b[k] par bande (mur de brique ou lisse / inspiré du masque).
• Calculer Y_b[k] = X[k]·W_b[k] puis appliquer l'IFFT pour obtenir y_b[n].
• Calculer le RMS / les moyennes de bande à partir de y_b[n].

Cette approche permet facilement d'implémenter un filtrage à phase nulle (non causal). Pour une conformité stricte aux normes IEC/ANSI, W_b et la normalisation doivent être soigneusement conçus et validés.

Rendre la synthèse FFT pseudo‑continue : OLA, double fenêtrage et normalisation d'amplitude

Pour produire des signaux temporels continus par bande, utilisez l'overlap‑add (OLA) : segmentation, fenêtrage, FFT, application de W_b, IFFT, fenêtre de synthèse puis OLA. Choisir des fenêtres d'analyse / de synthèse satisfaisant les conditions COLA (constant overlap‑add), par ex. Hann avec recouvrement de 50 %, afin d'éviter une modulation périodique du niveau.

Si l'objectif est de reproduire les filtres normalisés, comment choisir W_b ?

W_b[k] dépend de ce que vous cherchez à reproduire :

• Reproduire l'intégration à bords abrupts : W_b vaut strictement 0/1 à l'intérieur de [f1,f2].
• Reproduire le comportement des filtres IEC/ANSI : |W_b(f)| approxime le masque normalisé et la bande passante effective (fait correspondre ∫|W_b|²).
• Reproduire la complémentarité énergétique pour la reconstruction : concevoir Σ_b |W_b(f)|² ≈ 1 (section 7.6).

Vous ne pouvez généralement pas satisfaire parfaitement ces trois critères simultanément ; définissez votre priorité (conformité vs décomposition / reconstruction) dès le départ.

Bancs de filtres fréquentiels conservant l'énergie : pourquoi Σ|W_b|² est important

Si vous souhaitez que les énergies de bande se somment à l'énergie totale (à l'erreur numérique près), une conception courante vise une complémentarité de puissance approchée :

Les masques IEC/ANSI n'imposent pas nécessairement une stricte complémentarité, il ne faut donc pas supposer une additivité exacte en contexte de conformité.

Stratégies de Welch / de moyennage : comment rendre les niveaux de bande FFT stables

• Utiliser le moyennage de Welch (segmenter, fenêtrer, recouvrir, moyenner les spectres de puissance).
• Moyenner dans le domaine de la puissance (|X|² ou PSD), puis convertir en dB.
• Pour les signaux non stationnaires, envisager la STFT afin d'obtenir des matrices temps–bande.
• Indiquer dans le rapport le type de fenêtre, le recouvrement, le nombre de moyennes et le traitement ENBW/CG.

Analyse par regroupement FFT dans OpenTest

OpenTest prend en charge l'analyse en bandes d'octave basée sur le regroupement FFT :
1) Connectez le dispositif, par exempleSonoDAQ Pro
2) Sélectionnez les voies et ajustez les paramètres. Pour un microphone externe, activez l'IEPE et passez en mode de mesure de signal acoustique.
3) Dans la section Octave-Band Analysis sous Mode de mesure, choisissez l'algorithme basé sur la FFT.
4) Un enregistrement unique peut être analysé simultanément en bandes 1/1 d'octave, 1/3 d'octave, 1/6 d'octave, 1/12 d'octave et 1/24 d'octave.

Figure 2 : analyse en bandes d'octave par regroupement FFT dans OpenTest

Banc de filtres vs FFT / synthèse FFT : différences, conditions d'équivalence et compromis

Tableau comparatif

Dans quels cas les deux méthodes donnent‑elles (presque) les mêmes résultats ?

Les résultats moyennés par bande coïncident généralement de près lorsque :

• Vous comparez des niveaux de bande moyennés (et non des courbes de crête transitoires).
• Le signal est approximativement stationnaire et la durée d'observation est suffisamment longue.
• La résolution FFT est suffisamment fine pour que chaque bande contienne assez de raies (en particulier pour la bande la plus basse).
• La mise à l'échelle FFT est correcte (traitement du spectre simple face, Δf, normalisation de fenêtre U, ENBW/CG si nécessaire).
• Une pondération partielle des raies est utilisée aux bords de bande.

Pourquoi les différences augmentent pour les transitoires et les événements courts

Les différences sont induites par des échelles de temps non concordantes : les bancs de filtres ont un retard de groupe et une sonnerie dépendant de la bande mais une sortie continue ; la STFT utilise une fenêtre fixe qui définit à la fois la résolution fréquentielle et le lissage temporel. Si la durée de l'événement est comparable à la longueur de la fenêtre ou à la réponse impulsionnelle du filtre, les résultats dépendent fortement des détails d'implémentation.

Bilan d'erreur : d'où viennent généralement les divergences (et comment les localiser rapidement)

1. Mauvais moyennage / combinaison en dB : il faut toujours moyenner et sommer dans le domaine de l'énergie.
2. Mise à l'échelle FFT incohérente : conventions 1/N, spectre simple face vs double face, Δf, normalisation de fenêtre U.
3. Oublier les corrections de fenêtre : ENBW pour le bruit ; gain cohérent / fuite pour les tonalités.
4. Utilisation des fréquences nominales pour calculer les bords de bande au lieu des définitions exactes.
5. Absence de pondération partielle des raies aux frontières de bande (particulièrement préjudiciable aux basses fréquences).
6. Problèmes de multirate / anti‑repliement dans les bancs de filtres.
7. Constantes de temps / fenêtres de moyennage différentes entre méthodes.
8. Vraies différences de méthode : l'intégration à bords abrupts (brick‑wall) par regroupement de raies vs les jupes / pentes progressives des filtres normalisés induisent des décalages systématiques.

Une forte approche de débogage : commencer par faire correspondre la valeur quadratique moyenne totale en utilisant du bruit blanc (mise à l'échelle / ENBW / pondération partielle des raies), puis valider les fréquences centrales de bande et l'isolation des bandes adjacentes à l'aide de balayages sinusoïdaux ou de tonalités.

Liste de contrôle d'ingénierie : rendre l'analyse 1/3 d'octave correcte, stable et reproductible

Choisir une méthode : conformité → banc de filtres ; statistiques hors ligne → regroupement FFT

• Pour les réglementations / essais de type / comparabilité des instruments : privilégier les bancs de filtres conformes IEC/ANSI et indiquer l'édition de la norme et la classe. [1][3]
• Pour le traitement hors ligne, les grands jeux de données ou des définitions de bandes flexibles : le regroupement FFT peut être efficace, mais la mise à l'échelle et la pondération aux frontières doivent être rigoureuses.
• Si vous avez besoin de signaux temporels par bande (modulation, enveloppe, etc.) : envisager la synthèse FFT ou des bancs de filtres explicites.

Choix des paramètres FFT à partir de la bande la plus basse (exemple)

Exemple : fs=48 kHz, bande de plus basse fréquence d'intérêt : 20 Hz (1/3 d'octave). Sa largeur de bande n'est que de quelques hertz. Si vous voulez au moins M=10 raies par bande, vous pouvez avoir besoin de Δf_raie ≤ largeur de bande/10, ce qui implique un N très grand (par ex. ~100 k points ; 2^17=131072). Cela illustre pourquoi, en temps réel et pour la conformité, on privilégie souvent les bancs de filtres.

Erreurs typiques qui empêchent la concordance des résultats

• Sommer la magnitude |X| au lieu de la puissance |X|² ou de la PSD.
• Moyenner en dB au lieu de moyenner en puissance / valeur quadratique moyenne linéaire.
• Ignorer l'ENBW / la mise à l'échelle de la fenêtre pour le bruit. [7]
• Calculer les bords de bande à partir des fréquences nominales.
• Ne pas préciser les conventions de pondération / de moyennage temporel (Rapide / Lent / Leq). [5][6]

Flux de validation recommandé (quelle que soit l'implémentation)

• Test de tonalité au centre (ou balayage) : vérifier que l'énergie culmine dans la bonne bande et que le rejet des bandes adjacentes est conforme aux attentes.
• Bruit blanc / rose : vérifier la forme spectrale attendue en niveaux de bande et évaluer la stabilité / le temps de moyennage.
• Comparaison inter‑implémentations : comparer votre implémentation à une référence connue sur des signaux identiques ; isoler les différences de mise à l'échelle, de définition et de jupe de filtre.
• Consigner et figer les paramètres (définition de bande, fenêtrage, moyennage) dans le rapport d'essai.

Liste de contrôle pour la reproductibilité : inclure ces éléments dans les rapports afin que d'autres puissent recalculer vos niveaux

• Définition des bandes : base 10 ou base 2 ? b dans 1/b ? valeurs exactes ou nominales utilisées pour le calcul ? fréquence de référence fr ?
• Implémentation : banc de filtres normalisé (IIR/FIR, multirate) vs regroupement / synthèse FFT ; versions des logiciels / bibliothèques.
• Échantillonnage / prétraitement : fs, suppression de tendance / suppression de la composante continue, filtrage anti‑repliement, rééchantillonnage.
• Moyennage temporel : Leq / RMS par blocs / exponentiel ; constantes de temps, taille de bloc, recouvrement, nombre de trames moyennées ; contexte Rapide / Lent si pertinent.
• Détails FFT (le cas échéant) : type de fenêtre, N, pas de trame (hop), zero‑padding, normalisation de la PSD, traitement simple face, ENBW / gain cohérent, pondération partielle des raies.
• Calibration / unités : unités d'entrée et grandeurs de référence (par ex. 20 µPa), facteurs de calibration des capteurs et dates.
• Définition de la sortie : RMS vs crête vs puissance de bande ; conventions 10log vs 20log ; éventuelles étapes d'agrégation de bandes.

Si vous ne devez retenir qu'une phrase : documentez « définition de bande + moyennage temporel + traitement de la mise à l'échelle / du fenêtrage FFT (le cas échéant) ». La plupart des désaccords disparaissent.

Formules rapides et exemple numérique (prêts pour le code / le rapport)

Constantes pour les bandes un tiers d'octave en base 10

G = 10^(3/10) ≈ 1.995262
r = 10^(1/10) ≈ 1.258925 # adjacent center-frequency ratio
k = 10^(1/20) ≈ 1.122018 # edge multiplier about center
f1 = fm / k
f2 = fm * k

Exemple : la bande un tiers d'octave à 1 kHz

fm = 1000 Hz
f1 = 1000 / 1.122018 ≈ 891.25 Hz
f2 = 1000 * 1.122018 ≈ 1122.02 Hz
Δf ≈ 230.77 Hz
Q ≈ 4.33

OpenTest intègre les deux méthodes. Téléchargez‑le et commencez dès maintenant -> ou remplissez le formulaire ci‑dessous ↓ pour planifier une démonstration en direct.

Découvrez plus de fonctionnalités et de cas d'application sur www.opentest.com.

Références

[1] Extrait PDF IEC 61260‑1:2014 (iTeh) : https://cdn.standards.iteh.ai/samples/13383/3c4ae3e762b540cc8111744cb8f0ae8e/IEC-61260-1-2014.pdf

[3] Aperçu PDF ANSI S1.11‑2004 (ASA/ANSI) : https://webstore.ansi.org/preview-pages/ASA/preview_ANSI%2BS1.11-2004.pdf

[4] Note d'application HEAD acoustics : FFT – Analyse 1/n d'octave – Ondelette (description de banc de filtres) : https://cdn.head-acoustics.com/fileadmin/data/global/Application-Notes/SVP/FFT-nthOctave-Wavelet_e.pdf

[5] IEC 61672‑1:2013 (page IEC) : https://webstore.iec.ch/en/publication/5708

[6] NTi Audio Know‑how : pondération temporelle Rapide / Lente (contexte IEC 61672‑1) : https://www.nti-audio.com/en/support/know-how/fast-slow-impulse-time-weighting-what-do-they-mean

[7] MathWorks : exemple de définition de l'ENBW : https://www.mathworks.com/help/signal/ref/enbw.html