LinFIR - Mesures et égalisation FIR/IIR

[pvrx]

Bonjour Arnaud,

Cette mise à jour est la bienvenue, notamment pour moi qui suis en train de me faire les dents avec HFD...
Petites suggestions pour les prochaines mises à jour :
1. Enregistrement : pour l'instant, il n'existe qu'une commande "Save". Si on fait juste quelques modifs sur un fichier LinFIR, la fonction "Save" écrase la précédente. Or, on peut avoir simplement envie de créer une version différente de l'existante afin de produire des variantes et de les conserver. Une fonction "Save as" serait la bienvenue.
2. La documentation a été mise à jour, bravo et merci. Un point de détail concernant l'export HFD : sur l'exemple que tu donnes, le tweeter est indiqué en canal 1. Ce ne sera jamais le cas. Pour les modules Hypex 3 voies, le canal 3 est toujours le moins puissant et est donc réservé à l'aigu. Mais c'est vraiment du détail...

Pascal

[Arnwald]

Bonjour Pascal,

Pour le point 1, cela est déjà possible en passant par "Export/LinFIR project".

Bien vu pour l'exemple, je vais intervertir. :-)

[Arnwald]

Nouvelle mise à jour: 1.3.1

J'ai ajouté un mécanisme pour afficher la courbe cible lors de l'utilisation de l'auto EQ. J'ai aussi amélioré l'algorithme d'auto EQ: des filtres larges (à faible Q) pouvaient entrainer un masquage de l'erreur résultante pour éviter l'empilement de filtre, au détriment de la correction globale. Ceci a été corrigé avec un nouveau système de pénalités.

[pvrx]

Bonjour Arnaud,

Comme tu pourras le voir sur le fil dédié à mon projet (viewtopic.php?p=17181#p17181),
je me pose une question de faisabilité :
Je vais effectuer avec REW plusieurs mesures du même HP, dans l'axe et hors axe, et en faire une moyenne (RMS Average) afin d'établir une courbe de référence "bonne pour égalisation".
De ce fait, je perds ma référence de phase, et je voudrais injecter dans LinFIR la réponse en fréquence moyennée et la réponse en phase mesurée dans l'axe.
Est-ce possible ? Et comment faire, ou comment contourner l'obstacle ?

Merci de ta réponse !

Pascal

[Arnwald]

Bonjour Pascal,

Je ne pense pas que ce soit une très bonne idée, cela n’a pas vraiment de sens physique.

Dans les outils avancés il existe déjà la possibilité d’égaliser en utilisant la listening window.

[pvrx]

Bonjour Arnaud,

Merci de ton retour, mais je n'ai pas encore trouvé les outils avancés...

Une petite question : dans la documentation, il serait peut-être utile de lever l'ambiguïté du texte ci-dessous, et aussi de préciser le mode opératoire.

Lorsque que le HP émettant le chirp de référence est plus loin du micro que le HP sous test (cas le plus probable), le manual offset est-il positif ou négatif ? Et comment le paramétrer ?

Pascal

[Arnwald]

Bonjour Pascal,

Je vais voir comment améliorer la documentation.

Quand le haut parleur produisant le chirp de référence est plus proche, il faut reculer l’impulsion, donc ajouter un offset positif.

[jlo]

Bonjour Arnaud et tous les autres,
le soft en version Windows est-il toujours dispo ? Le store Microsoft propose un installer de 1.1M dont le lancement ramène à ce même store avec une erreur.
Ce qui fait que je n'ai pas pu tester mais après lecture (trop partielle ?) de la doc, j'aimerai demander certaines confirmations sur le calcul de la Room Correction :
- après mesures, on obtient des IR dont au moins une seconde est utilisée : le calcul de Room Correction est-il basé sur une FFT globale de cet IR sans aucun fenêtrage ?
- en FIR (et évidemment en IIR), la correction est-elle à phase minimale, càd juste une transformée de Hilbert de l'amplitude ?
- comme on est à phase minimale, je suppose que l'IR de correction est assymétrique avec le pic près du temps zéro et donc une latence pratiquement nulle, quel que soit le nombre de taps de cette correction.

[Arnwald]

Bonjour Arnaud et tous les autres,
le soft en version Windows est-il toujours dispo ? Le store Microsoft propose un installer de 1.1M dont le lancement ramène à ce même store avec une erreur.
Ce qui fait que je n'ai pas pu tester mais après lecture (trop partielle ?) de la doc, j'aimerai demander certaines confirmations sur le calcul de la Room Correction :
- après mesures, on obtient des IR dont au moins une seconde est utilisée : le calcul de Room Correction est-il basé sur une FFT globale de cet IR sans aucun fenêtrage ?
- en FIR (et évidemment en IIR), la correction est-elle à phase minimale, càd juste une transformée de Hilbert de l'amplitude ?
- comme on est à phase minimale, je suppose que l'IR de correction est assymétrique avec le pic près du temps zéro et donc une latence pratiquement nulle, quel que soit le nombre de taps de cette correction.

Salut @jlo

En effet, la FFT est calculée sur l'IR complète sans fenêtrage pour le mode Room Calibration.

La correction de magnitude FIR est phase minimum par défaut, mais ce comportement peut être modifié avec le curseur de causalité. Une causalité de 0 équivaut à un filtre phase linéaire, et 1 donne un filtre phase minimum.
Le filtrage IIR est toujours phase minimum, par définition.
Tant que l'on n'introduit pas de filtre phase linéaire (par exemple un crossover FIR) ou de correction de phase, l'énergie de l'IR du FIR est effectivement concentrée près de t=0.

Les différents filtres FIR que l'on peut appliquer sont convolués les uns avec les autres pour former un filtre FIR unique. J'ai mis en place un algorithme de centrage de l'IR obtenue pour que le filtre soit tronqué de façon optimale selon le nombre de taps alloués:
- si tous les filtres sont phase linéaire (sans correction de phase), l'IR est simplement centrée,
- si tous les filtres sont phase minimum (sans correction de phase également), l'IR est concentrée vers t=0,
- un mélange de filtres phase linéaires et phase minimum donnera un entre deux.

[jlo]

Merci pour ces confirmations.
Les transformées IR > FFT(IR) > iFFT(FFT(IR)) donnent une IR à phase nulle avec l'énergie centrée au milieu de la longueur de la FFT.
Par curiosité, pour retrouver une IR à phase minimale, tu utilises quoi comme algorithme : Hilbert, cepstre, minrceps,... ?

[Arnwald]

J’utilise la méthode du filtrage homomorphique en passant par la transformation cepstrale. La proportion de quéfrences positives et négatives donne un comportement plus ou moins causal. :-)

[Dagda]

Bah vous voyez, rien que sur ce message, j'ai 2 mots que je ne connais pas et la seconde phrase (fréquence plutôt que quéfrences j'imagine) j'ai pas capté le concept ...

[jlo]

J’utilise la méthode du filtrage homomorphique en passant par la transformation cepstrale. La proportion de quéfrences positives et négatives donne un comportement plus ou moins causal. :-)

Merci pour cette info, c'est intéressant.

rien que sur ce message, j'ai 2 mots que je ne connais pas et la seconde phrase (fréquence plutôt que quéfrences j'imagine) j'ai pas capté le concept ...

Bienvenu dans le monde cepstral

le cepstre est le résultat de la transformée inverse du logarithme du spectre du signal.............spectre devient cepstre, fréquence devient quéfrence, harmonique devient rahmonique, phase devient saphe, filtrage devient liftrage

L'analyse cepstrale est surtout utilisée en sismologie. En audio, c'est quelquefois utilisé pour l'analyse de l'influence des réflexions sur la courbe de réponse.

[Arnwald]

Bah vous voyez, rien que sur ce message, j'ai 2 mots que je ne connais pas et la seconde phrase (fréquence plutôt que quéfrences j'imagine) j'ai pas capté le concept ...

Nan nan, pas de faute de frappe, on parle bien du domaine quéfrentiel

En analyse cepstrale, on déroule le cercle unité complexe sur une ligne, ce qui permet de séparer les composantes causales et anticausales d'un signal.

Il y a une explication un peu "intuitive" ici:
https://www.katjaas.nl/minimumphase/minimumphase.html

D'ailleurs, je devrais parler d'homomorphic liftering et pas filtering

[Dagda]

On en apprend tout les jours

D.

[Elgrosso]

Ouch ! C’est du louchébem d’ingé-son ?

[xnwrx]

Comme avec le cepstre on fait un TF inverse du log de la TF d'un signal, on est en mesure de dissocier dans le domaine temporel les composantes sources et filtre. C'est le log qui permet ça puisque log a.b = log a + log b or a.b en fréquentiel est une convolution en temporel (un filtre)

[Dagda]

[pvrx]

Bonjour Arnaud,

Petites questions concernant l'export en Wav de mesures effectuées avec REW :
Comment faudrait-il idéalement paramétrer la fenêtre d'exportation REW ci-dessous ?
Est-ce que ce type d'exportation contient les données de distorsion ?
Est-il utile de lisser avant d'exporter ?
Si on importe successivement deux mesures effectuées à position de micro identique, par exemple le mid et l'aigu, comment préserver leur décalage temporel lors de l'importation dans LinFIR ?

Merci de ton retour !

Pascal