Bonjour,
je cherche à ajouter du bruit à un signal connaissant sa densité spectrale de puissance. Ce bruit peut être additif ou multiplicatif.
Dans le cas d'un bruit additif j'ai vu que, pour retirer un bruit additif on pouvait faire :
FFTSignalBruite = FFTSignal*(1-DSPbruit/(DSPSignal+DSPBruit) )
J'essaye de faire des tests sous matlab mais sans succès...
Est ce que qulqu'un pourrait m'aider avec des idées ou un petit bout de code ?
Merci
Bonjour,
Si par bruit additif, on entend qu'au signal s(t), on ajoute un bruit b(t), on obtient : srésultant(t) = spur(t) + b(t).
Donc pour la FFT : Srésultant(f) = Spur(f) + B(f) (les FFT s'ajoutent)
Bruit multiplicatif : srésultant(t) = spur(t) x b(t) donc pour la FFT : Srésultant(f) = Spur(f)*B(f) avec * signifiant produit de convolution.
"Quand les gens sont de mon avis, il me semble que je dois avoir tort."O.Wilde
Salut NicoEnac et merci pour ta réponse.
Le seul problème c'est que je ne connais pas la FFT de mon bruit. Je ne connais que sa DSP. (DSP = |FFT|²). Donc je perds l'information de phase de mon spectre.
Un début de réponse pourrait être dans :
http://publication.lal.in2p3.fr/2002...L/node258.html
où on génère un bruit par une méthode autorégressive à partir d'un bruit blanc et de la DSP désirée. J'ai l'impression que c'est compliqué et mon chef se demande s'il n'y a pas plus simple... Je suis toujours à la recherche d'un coup de pouce svp si vous avez une idée.
Il y a quelquechose que je ne comprends pas dans la méthode autorégressive.
On se fixe un ordre p (puissance de 2 en général) et on fait un calcul des coefficients par diverses méthodes (algorithme de Durbin Levinson, moindre carrés comme dans http://www.mathworks.com/access/help...ient=firefox-a )
Supposons que ma DSP soit sur 1024 échantillons. Lorsque je prends la TF inverse, j'ai une fonction d'autocorrelation sur 1024 échantillons.
J'ai vu que l'odre choisi pour l'AR est 4 voir 8 donc on résouds les équations de Yule-Walker (cf. lien précédent) en considérant que 4 ou 8 coefficient d'autocorrelation ?
J'ai l'impression qu'on perd plein d'infos...
Quelqu'un peut m'éclaircir ca svp je suis perdu...
Personne ?
