Correction/filtrage semi-automatisé (DStudio v5)

[thierryvalk]

Oui mais non, tu a quoté en modifiant le contenu.

Et tout cela ne m'éclaire pas sur ce que je dois faire.

[alka]

Oui mais non, tu a quoté en modifiant le contenu.

j'ai édité pour la précision des propos

Et tout cela ne m'éclaire pas sur ce que je dois faire.

tu dois faire une fenetre pardi !

le plus simple serait de faire une fenetre comme REW, tukey 25.

[thierryvalk]

Évidemment la fenêtre Tukey est la plus complexe.
Pas grave, j'ai trouvé fonction qui va bien.

Si j'ai bien compris.
Je crée un tableau de la longueur de la fenêtre.
Je lui applique pour chaque valeur la fonction du fenêtrage. On devrait avoir un tableau dont les valeurs vont de 0 à 1 pour terminer a nouveau à 0.

Je crée un nouveau tableau de longueur égale à l'impulse contenant tous des 0.
Je prend le nombre de samples qui détermine le début de la fenêtre . A partir de ce nombre je copie le tableau de la fenêtre dans ce tableau.
Je multiplie chaque valeur de ce tableau par celui de l'impulse et devrait obtenir une impulse fenêtrée.

L'astuce est que tout comme la position du t0, on ne touche pas à l'impulse de départ, faudra donc voir si l'on calcule a chaque fois le fenêtrage ou si l'on conserve l'impulse fenêtrée.

[alka]

oui c'est ça. Ca revient a multiplier les valeurs de l'IR dans la partie fenetrée par un coefficient qui va de 0 à 1. Toute l'idée est de rendre cette IR fenetrée sans rupture brutale, allant doucement de 0 à 0 sans toutefois altérer l'information utile.

un algorithme plus direct :
Créer une nouvelle table IR appelée IRF (F=fenetrée)
Permet d'éviter de toucher a la première table pour pouvoir y revenir si on veut changer le fenetrage.

puis

Code : Tout sélectionner

for t=0 à Tf                 // Tf = longueur de fenetre
   IRF(t) = K(t) * IR(to+t)  // t0 = début de la fenetre
next t

(éventuellement remplir IRF avec des 0 au delà de Tf s'il le faut. me rappelle plus si la fft a besoin d'une longueur d' impulse fixe.)

K(t) = coefficient dont la valeur dépend du type de fenetrage. K(t) a des valeurs de 0 à 1.
- pour rectangulaire, K(t) = constant = 1
- pour les autres, K dépend de t selon les jolies formules de wikipedia

[thierryvalk]

Oui, la FFT doit toujours avoir une longueur de 32K.
C'est un choix qui simplifie grandement la gestion de l'ensemble.Et elle retourne toujours une réponse ne fréquence en 1/48 de longueur fixe aussi.

[alka]

ah oui. maintenant que tu le dis ca me revient. J'oublie trop de choses entre deux pages du sujet

mais ou est passé mon magnésium déjà ?

[thierryvalk]

Ça ressemble à ceci du tukey 0.25 ?

Pour le moment c'est bidon, fenêtre théorique de -50 à +450 samples ce qui ne semble pas juste, mais c'est du détail.
Edit : correction de l'axe X, image changée.

Question 2, faut' il afficher l'impulse non fenêtrée, fenêtrée ou les 2 ?

[alka]

la fenetre tukey de REW ne ressemble pas à ça. La tienne a deux poles alors que pour REW il n'y en a qu'un seul. Tu l'as trouvé où ton algorithme tukey ?

J'ai pas rew sur ce pc et peux pas faire copie d'écran. SDans REW, afficher l'Impulse d'une mesure quelconque, et décocher tout en ne laissant que window. Puis clic sur IR window et essayer les différentes formes. REW définit sa fenetre comme -t1 et +t2 p/r à t0. Son t0 est dans la fenetre, pas l'un des bords. En pricnipe, t0 est sur le pic, t1 est le début un peu avant le pic et t2 est la fin, assez loin du pic. Tukey 0.25 change de courbure selon la distance de t0 au bord t1 ou t2. C'est plutot adapté.

Sur ton exemple, on voit que cette forme de fenetrage touche au pic de l'impulse, pas bien.
C'est pas bete l'histoire de REW de faire une forme différente au début et a la fin.

Pour l'affichage, je dirai montrer au choix l'une ou l'autre. Les deux simultanément c'est trop brouillon et on ne voit rien.
Montrer en fond d'écran la courbe de la fenetre comme tu fais là c'est bien aussi, ca éduque les foules

[thierryvalk]

Dans ma version de REW, on peut choisir le type de fenêtrage pour Left et Right et par défaut tous les 2 sont sur Tukey 0.25.
Et en jouant comme tu le décris, c'est juste que la pente change en se rapprochant de 0.

Par contre, selon moi et pas spécialiste du tout, ce n'est plus vraiment du Tukey 0.25 mais avec une valeur de 0.25 qui s'adapte en fonction du temps.

[thierryvalk]

voici du 0.25 au-dessus et 0.5 en-dessous.

Edit : je crois comprendre.
REW fait d'une certaine manière 2 fenêtrages, l'un en Left et l'autre ne Right alors que moi un seul.
Et je présume que la pente est fonction de la largeur de la fenêtre.

[alka]

Edit : je crois comprendre.
REW fait d'une certaine manière 2 fenêtrages, l'un en Left et l'autre ne Right alors que moi un seul.
Et je présume que la pente est fonction de la largeur de la fenêtre.

deux fenetres oui. Largeur variable du tukey probable car il différentie dans sa liste "Tukey 10ms" et "Tukey 0.25" sans unité. Fort a parier que le 0.25 représente une fraction de la largeur de la 1/2 fenetre concernée. 25% peut etre.

Cela dit il y a autre chose : même en grossisant, la forme n'est pas la meme entre REW et DS5.

edit : image de tukey dans REW.

[alka]

tukey n'a qu'une courbure chez REW (un seul pole) alors que tu en a deux bien visibles.
0.25 a l'air de correspondre a 25% du bord pour le départ de la courbure chez REW. démarre a 750ms pour la 1/2 fenetre entre t0=0 et T=1000

[thierryvalk]

Voici par rapport au t0 :

Pour ce qui de la forme, désolé, mais si tu fais une recherche sur Google tu verras que c'est REW qui a tort.
Je dis pas que c'est pas bon mais pas du tukey.

[thierryvalk]

Et mon REW donne ceci :

Ce qui semble très rapprochant avec mon résultat.

Edit : tu as mis 1000ms, c'est p'têtre une peu trop.

[alka]

c'est la différence d'affichage dans REW : %Fs ou dBFs qui change la forme de la courbe. j'étais en dB FS
OK pour l'allure de la fenetre qui a l'air bonne.
Reste que le détail du cadrage gauche/droite a régler pour bien prendre l'impulse comme il faut.

[thierryvalk]

c'est la différence d'affichage dans REW : %Fs ou dBFs qui change la forme de la courbe. j'étais en dB FS

En effet l'axe Y est soit en log soit en lin, étrange tout de même de ne pas avoir alors une portion de cercle dans le bas.

Reste que le détail du cadrage gauche/droite a régler pour bien prendre l'impulse comme il faut.

Quoi ? Comment ? elle n'est pas correcte ma dernière courbe ?

[alka]

j'étais resté sur l'image de REW !
oui parfait. maintenant la fenetre ne rabote plus le pic de l'impulse.


t0 n'est pas éditable. Tu mets d'office t0 sur le pic ?

[thierryvalk]

Faut suivre.
t0 peut être positionné au peak, au début de l'impulse ou manuellement via les buttons + et - qui sont à 1 sample.