13.07.2015 Views

Modélisation de sons bruités par la Synth`ese Granulaire

Modélisation de sons bruités par la Synth`ese Granulaire

Modélisation de sons bruités par la Synth`ese Granulaire

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

du signal analysé. En réalité, l’information <strong>de</strong> phase du spectre n’est pasvraiment primordiale dans notre cas s’attachant plus à détecter <strong>de</strong>s formesspectrales associables à un jeu <strong>de</strong> <strong>par</strong>amètres. Nous avons donc choisi <strong>de</strong>travailler exclusivement sur le module du spectre |X(t, f)|. Nous supposeronsdonc que |X(t, f)| peut être décomposé selon :|X(t, f)| =K∑α k |G k (t, f)|k=1Il est important <strong>de</strong> voir que l’on travaille maintenant avec <strong>de</strong>s matriceset non plus <strong>de</strong>s vecteurs temporels. Pour simplifier les notations et <strong>la</strong> structure<strong>de</strong> <strong>la</strong> métho<strong>de</strong>s, nous avons fait en sorte <strong>de</strong> manipuler le module <strong>de</strong>sspectrogrammes sous <strong>la</strong> forme <strong>de</strong> vecteurs X[p t,f ] :X[p t,f ] = |N∑w[n]x[t ∗ d + n]e 2jπf |n=0avec l’indicep t,f = tN + f/Feavec N taille <strong>de</strong> <strong>la</strong> TFCT et d pas d’avancement pour le calcul <strong>de</strong>s TFCTrég<strong>la</strong>nt notamment leur recouvrement.Comme préce<strong>de</strong>mment, on veux donc à chaque étape i <strong>de</strong> l’algorithmetrouver dans le dictionnaire D l’atome G m(i) qui, <strong>par</strong> projection du résidu R iminimise <strong>la</strong> norme-2 <strong>de</strong> R i+1 <strong>de</strong>fini <strong>par</strong> :R i+1 = R i − α i G m(i)avec m(i) <strong>de</strong>signant l’in<strong>de</strong>x <strong>de</strong> l’atome choisi dans le dictionnaire à l’étapei. L’algorithme commence avec R 1 = X. Le but est <strong>de</strong> trouver G m(i) selon :Par ailleurs :G m(i) = arg min Gm(i)∈D‖R i+1 ‖ 2〈R i+1 , G m(i) 〉 = 〈R i , G m(i) 〉 − α i 〈G m(i) , G m(i) 〉Par le théorème <strong>de</strong> <strong>la</strong> projection othogonale on a :〈R i+1 , G m(i) 〉 = 〈R i − α i G m(i) , G m(i) 〉 = 0D’où :α i = 〈R i, G m(i) 〉〈G m(i) , G m(i) 〉 = 〈R i, G m(i) 〉47

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!