Discussion :

[oliAudio]

Bonjour à tous,

Je souhaiterais obtenir la réponse impulsionnelle d'un signal à partir d'une réponse en amplitude et en phase donnée X(f).

lorsque je fais la FFT inverse x = ifft(X), j'obtiens un signal comme ceci : \___/

J'applique donc fftshit(x) pour obtenir x0 = _/\_.

Mais, FFT(x0) et X on une réponse en phase différente! Leur réponse en amplitude sont identiques.

Si quelqu'un à une idée comment je peux obtenir une réponse impulsionnelle comme ça _/\_ avec une FFT égale à X je suis plus que preneur,

Merci, @+

Oli

[phryte]

Salut.
As-tu un bout de code à nous montrer ?

[oliAudio]

Salut phryte,

En exemple de code on peut prendre :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
% A : module du spectre S. A est paire par rapport à Fe.
% P : phase du spectre S. P est imparie par rapport à Fe.
 
S = A.*exp(i*P);
 
s = real(ifft(S)); % Par les ptés de A et P => s reel.
                       % On met real pour les erreurs d'arrondis.

Le problème est que s(t) ne ressemble pas à une réponse impulsionnelle. Il faut passer par fftshift, mais angle(fft(fftshift(s))) ≠ P !!!

(on a bien evidemment angle(fft(s)) = P)

Je souhaite un reponse impulsionnelle impulse tel que angle(fft(impulse)) = P, mais pas moyen...

++

[oliAudio]

La fft est bijective, donc modif de s(t) => modif de fft(s). La méthde ifft(spectre) pour obtenir le filtre FIR n'est pas bonne.

Je vais essayer la fonction cfirpm() pour générer le filtre FIR non linéaire.

++