09-22-2018, 02:00 AM
(Modification du message : 09-22-2018, 02:18 AM par LeChacal619.)
Bon allez que je me couche l'esprit tranquille quand même... je m'en voudrais de lâcher l'affaire aussi vite. J'ai généré un signal avec un code matlab : un impulsion parfaite (phase linéaire, ou plus intelligemment délai de groupe constant), suivi d'une réverbération artificielle générée avec un sinus à 1000hz dont j'ai bruité fortement la phase pour étaler l'énergie de la réverbération sur tout le spectre, j'ai modulé cette réverbération par une courte rampe montante exponentielle et ensuite une décroissance exponentielle (atténuation en dBSPL linéaire sur l'échelle logarithmique donc).
Le fichier .wav généré se trouve ici. Je vous laisse juger ce qui est le mieux : la phase, le GD de REW, ou le "délai de groupe" affiché sur le spectrogramme (rappel : la réverbération à 1000Hz est environ 114ms après l'impulsion si mes calculs sont corrects...).
Si vous voulez corrigez la phase (vous ne trouvez pas que quelque chose cloche avec cette phase ? )
Phase unwrap
ou le GD de REW
GD
libre à vous de le faire. Moi, je ne le ferai pas : le délai de groupe estimé par le spectrogramme me va très bien...
Spectrogramme
Promis, si vous me dites que tout va bien et que je suis un fou, j'arrêterai là, et je retournerai à mes livres de traitement du signal, pour voir si j'arrive à comprendre quelque chose au calcul de cette phase...
PS: il y a un bug avec ce forum pour le linkage d'images ou d'url... c'est assez embêttant
PS2: la phase est à l'ouest, mais imaginez bien que la réverbération est parfaitement distincte de l'impulsion et à -42dBX (pas SPL ou FS je me ferai taper dessus) de celle-ci......
Cdlt,
Jean
Code :
% 5 secondes à 96khz
t=linspace(0,5,96000*5);
s=zeros(1,480000);
% Dirac à -6dB à 2.5secondes
s(1,240000)=10^(-6/20);
% Réverbération centrée autour de 1000hz à -20dB, décroissance
% exponentielle (= atténuation logarithmique linéaire)
s2=10^(-30/20)*cos(2*pi*1000*t+rand(1,480000)*5).*exp(-20*t);
% Rampe de montée
up=flip(exp(-1000*t(1:10000)));
s2(1:10000)=s2(1:10000).*up;
% Dirac + résonance 1000hz, pic de résonance à 240961+10000=250961 soit
% 2.614sec, soit 114ms après l'impulsion
s(240961:480000)=s(240961:480000)+s2(1:239040);
audiowrite('test_phase.wav',s,96000)
Le fichier .wav généré se trouve ici. Je vous laisse juger ce qui est le mieux : la phase, le GD de REW, ou le "délai de groupe" affiché sur le spectrogramme (rappel : la réverbération à 1000Hz est environ 114ms après l'impulsion si mes calculs sont corrects...).
Si vous voulez corrigez la phase (vous ne trouvez pas que quelque chose cloche avec cette phase ? )
Phase unwrap
ou le GD de REW
GD
libre à vous de le faire. Moi, je ne le ferai pas : le délai de groupe estimé par le spectrogramme me va très bien...
Spectrogramme
Promis, si vous me dites que tout va bien et que je suis un fou, j'arrêterai là, et je retournerai à mes livres de traitement du signal, pour voir si j'arrive à comprendre quelque chose au calcul de cette phase...
PS: il y a un bug avec ce forum pour le linkage d'images ou d'url... c'est assez embêttant
PS2: la phase est à l'ouest, mais imaginez bien que la réverbération est parfaitement distincte de l'impulsion et à -42dBX (pas SPL ou FS je me ferai taper dessus) de celle-ci......
Cdlt,
Jean