Transformee de fourier dans les applications audio-numeriques

Soyez le premier à donner votre avis sur cette source.

Snippet vu 28 334 fois - Téléchargée 32 fois

Commenter

Contenu du snippet

Le module Fourier comprend les fonctions utiles dans les applications destinees à manipuler le son numerique. Les commentaires fournissent les definitions de base et expliquent le principe. La transformee de Fourier sert a analyser le spectre de fréquences d'un echantillon sonore. La transformee inverse permet de calculer le son apres filtrage des magnitudes.

Source / Exemple : 


Module ModuleFourier
    
    ' Utilitaires pour calculs audio-numériques

    ' FFT (Fast Fourier Transform) : TRANSFORMEE RAPIDE DE FOURIER

    ' PCM (Pulse Code Modulation) : modulation par impulsions
    ' Sampling-Rate : Cadencement (Fréquence des impulsions) 

    ' La transformée de Fourier permet de calculer le spectre de fréquences 
    ' correspondant à l'échantillon PCM de longueur 2^P2 .
    ' Exemple pour un bloc de 8192 impulsions (2^13),dans un échantillon cadencé à 44100 hz
    ' on aura une fréquence de base de 44100 / 8192 =   5.38330078125 Hz
    ' C'est à dire des magnitudes et angles correspondants aux fréquences :
    ' 10.76 ,   16.14 , 21.53 ,   26.91 , 32.29 ,   37.68 , 43.06 ,   48.44 , 53.83   ...
    ' ... 430.66 ,  436.04 , 441.43 , 446.81 ... jusqu'à  22050 Hz (fréquence de Nyquist)
    ' REMARQUE : On ne connaitra jamais la magnitude exacte du " La 440 " par exemple (diapason)
    ' Pour celà , il faudrait un échantillon cadencé à : (440 / 82) * 8192 = 43957 Hz

    Public Structure NombreComplex
        Public Réel As Double
        Public Imaginaire As Double
    End Structure

    Public Function Fourier(ByVal Suite() As NombreComplex) As NombreComplex()

        'Calcule la transformée d'une suite de nombres complexes

        Dim Nbr As Integer ' Nombre de complexes dans la suite
        Dim P2 As Byte ' Puissance de 2 du nombre entier Nbr
        Dim ISuit, ITrans As Integer ' Compteurs de boucles
        Dim Fin, Taille As Integer ' Butées
        Dim dAngl, dAr As Double
        Dim Alpha, Beta As Double
        Dim k, I, J, L As Integer
        Dim AR, TR, Ti, AI As Double

        ' Nbr est la taille de la suite à transformer
        Nbr = UBound(Suite) + 1 ' Obligatoirement: une puissance de 2
        P2 = Puissance(Nbr)

        Dim Transformée(Nbr - 1) As NombreComplex

        For ISuit = 0 To Nbr - 1
            ITrans = Ivrs(ISuit, P2)
            Transformée(ITrans).Réel = Suite(ISuit).Réel
            Transformée(ITrans).Imaginaire = Suite(ISuit).Imaginaire
        Next ISuit

        Fin = 1
        Taille = 2
        Do While (Taille <= Nbr)
            dAngl = 2 * Math.PI / Taille
            Alpha = 2 * ((System.Math.Sin(0.5 * dAngl)) ^ 2)
            Beta = System.Math.Sin(dAngl)

            I = 0
            Do While (I < Nbr)
                AR = 1 ' Cosinus(0)
                AI = 0 ' Sinus(0)

                J = I
                For L = 1 To Fin
                    k = J + Fin
                    TR = AR * Transformée(k).Réel - AI * Transformée(k).Imaginaire
                    Ti = AI * Transformée(k).Réel + AR * Transformée(k).Imaginaire
                    Transformée(k).Réel = Transformée(J).Réel - TR
                    Transformée(k).Imaginaire = Transformée(J).Imaginaire - Ti
                    Transformée(J).Réel = Transformée(J).Réel + TR
                    Transformée(J).Imaginaire = Transformée(J).Imaginaire + Ti
                    dAr = Alpha * AR + Beta * AI
                    AI = AI - (Alpha * AI - Beta * AR)
                    AR = AR - dAr
                    J = J + 1
                Next L
                I = I + Taille
            Loop
            Fin = Taille
            Taille = Taille * 2
        Loop
        Return Transformée
    End Function

    Public Function Inverse(ByVal Transformée() As NombreComplex) As NombreComplex()

        'Calcule la transformée inverse d'une suite de nombres complexes

        Dim Nbr As Integer ' Nombre de complexes dans la suite
        Dim P2 As Byte ' Puissance de 2 du nombre entier Nbr
        Dim ISuit, ITrans As Integer ' Compteurs de boucles
        Dim Fin, Taille As Integer ' Butées
        Dim dAngl, dAr As Double
        Dim Alpha, Beta As Double
        Dim k, I, J, L As Integer
        Dim AR, TR, Ti, AI As Double

        ' Nbr est la taille de la Transformée à transformer
        Nbr = UBound(Transformée) + 1 ' Obligatoirement: une puissance de 2
        P2 = Puissance(Nbr)

        Dim Suite(Nbr - 1) As NombreComplex

        For ISuit = 0 To Nbr - 1
            ITrans = Ivrs(ISuit, P2)
            Suite(ITrans).Réel = Transformée(ISuit).Réel
            Suite(ITrans).Imaginaire = Transformée(ISuit).Imaginaire
        Next ISuit

        Fin = 1
        Taille = 2
        Do While (Taille <= Nbr)
            dAngl = -2 * Math.PI / Taille
            Alpha = 2 * ((System.Math.Sin(0.5 * dAngl)) ^ 2)
            Beta = System.Math.Sin(dAngl)

            I = 0
            Do While (I < Nbr)
                AR = 1 ' Cosinus(0)
                AI = 0 ' Sinus(0)

                J = I
                For L = 1 To Fin
                    k = J + Fin
                    TR = AR * Suite(k).Réel - AI * Suite(k).Imaginaire
                    Ti = AI * Suite(k).Réel + AR * Suite(k).Imaginaire
                    Suite(k).Réel = Suite(J).Réel - TR
                    Suite(k).Imaginaire = Suite(J).Imaginaire - Ti
                    Suite(J).Réel = Suite(J).Réel + TR
                    Suite(J).Imaginaire = Suite(J).Imaginaire + Ti
                    dAr = Alpha * AR + Beta * AI
                    AI = AI - (Alpha * AI - Beta * AR)
                    AR = AR - dAr
                    J = J + 1
                Next L
                I = I + Taille
            Loop
            Fin = Taille
            Taille = Taille * 2
        Loop
        For I = 0 To UBound(Suite)
            Suite(I).Réel = 2 * Suite(I).Réel \ Nbr
            Suite(I).Imaginaire = 2 * Suite(I).Imaginaire \ Nbr
        Next I
        Return Suite
    End Function

    Public Function Magnitude(ByVal Transformée() As NombreComplex) As Double()
        Dim Id As Integer
        Id = UBound(Transformée) \ 2
        Dim M(Id) As Double
        For Id = 0 To UBound(M)
            M(Id) = Math.Sqrt(Math.Pow(Transformée(Id).Réel, 2) + Math.Pow(Transformée(Id).Imaginaire, 2))
        Next
        Return M
    End Function

    Public Function Ivrs(ByVal IdPuls As Integer, ByVal P2 As Byte) As Integer
        Dim I As Short
        Ivrs = 0

        For I = 1 To P2
            Ivrs = (Ivrs * 2) Or (IdPuls And 1)
            IdPuls = IdPuls \ 2
        Next I

    End Function

    Public Function Puissance(ByVal Entier As Integer) As Byte
        Return Int(Math.Log(Entier) / Math.Log(2))
    End Function

End Module

Conclusion :


Remerciements: Don Cross et Jean Debord

A voir également

Ajouter un commentaire

Commentaires

Signaler
cs_Patrice99
Messages postés
1222
Date d'inscription
jeudi 23 août 2001
Statut
Membre
Dernière intervention
9 septembre 2018
 -
Evite de coder PI ainsi : tu perds un peu de précision, utilise plutot Math.PI
Voir aussi :
http://patrice.dargenton.free.fr/vbwavecomp/index.html
http://patrice.dargenton.free.fr/vbwavecomp/VBWaveCompSrc.html
Signaler
cs_AbriBus
Messages postés
492
Date d'inscription
jeudi 28 août 2003
Statut
Membre
Dernière intervention
25 avril 2007
 5 -
Salut,
Tout ca me semble un peu obscur dans son application finale... apres quelques recherche sur le net, j'ai cru comprendre que fourier pouvait etre utilisé pour le cadage de Mp3 (?)
Serait-il possible d'avoir plus d'informations quand à l'application du calcul de la transformee de fourier ?

AbriBus
Signaler
jazzd2
Messages postés
1
Date d'inscription
lundi 26 mai 2003
Statut
Membre
Dernière intervention
22 mars 2006
 -
La transformée de fourrier, sans rentrer dans les détails, sert à pas mal de choses dans l'analyse du signal. Elle est utilisé pour faire apparaitre les fréquences d'un signal. Elle est utilisé pour le son, les signaux acquis dans l'industrie, le traitement d'image.... en gros chaque fois que tu veux analyser un signal pour trouver certaines infos c'est une des méthodes utilisées.
Signaler
cs_AbriBus
Messages postés
492
Date d'inscription
jeudi 28 août 2003
Statut
Membre
Dernière intervention
25 avril 2007
 5 -
Salut,
tout d'abord, merc jazzd2 pour l'info... mais, toujours sans rentrer dans des détails indescents, quel est le genre d'info qu'on recupere par ces methode ? la qualité, la puissance, l'integrité ...?
En fait la question, c'est quel est l'objectif final d'une telle analyse (je ne doute pas ici qu'il y en a une...)?
Signaler
cs_Patrice99
Messages postés
1222
Date d'inscription
jeudi 23 août 2001
Statut
Membre
Dernière intervention
9 septembre 2018
 -
La TF donne le spectre du signal, c'est-à-dire la distribution (répartition) des fréquences contenues dans le signal, ce qui est utile par exemple pour l'audio, pour savoir l'énergie du signal qui risque d'être perdu selon le type d'écouteur.
Voir par exemple :
http://patrice.dargenton.free.fr/vbwavecomp/index.html
Afficher les 39 commentaires
S'inscrire maintenant

Vous n'êtes pas encore membre ?

inscrivez-vous, c'est gratuit et ça prend moins d'une minute !

Les membres obtiennent plus de réponses que les utilisateurs anonymes.

Le fait d'être membre vous permet d'avoir un suivi détaillé de vos demandes et codes sources.

Le fait d'être membre vous permet d'avoir des options supplémentaires.

Du même auteur (comemich)