Quelle est l'explication la plus lucide et intuitive pour les différents FT - CFT, DFT, DTFT et la série Fourier?


30

Même après les avoir étudiés pendant un certain temps, j'ai tendance à oublier [si je suis déconnecté pendant un certain temps] comment ils sont liés les uns aux autres et ce que chacun représente [car ils ont des noms qui sonnent tellement similaires]. J'espère que vous arriverez à une explication si intuitive et si belle mathématiquement qu'ils seront intégrés dans ma mémoire pour toujours et ce fil servira de rafraîchissement super rapide chaque fois que j'en aurai besoin [ou quelqu'un d'autre].


2
Devrait probablement commencer par la série de Fourier
endolith

Connaissez-vous la dualité Pontryagin?
Lorem Ipsum,

@yoda - Non. Pourriez-vous s'il vous plaît développer ou me signaler quelques bonnes références? [Je vais bien sûr le
chercher

1
"Steve on Image Processing": les transformées de Fourier répondent exactement à cette question.
nobar

Je ne sais pas quand réécrire une réponse ici (sauf si requis). Pourtant, une réponse possible est donnée dans Puis-je étudier la transformée de Fourier à temps continu et traiter le reste comme des cas particuliers en suivant la piste de dualité de Pontryagin proposée par @LoremIpsum
Laurent Duval

Réponses:


24

J'ai écrit ce document en complément d' Oppenheim et Willsky . Veuillez consulter le tableau 4.1 à la page 14, reproduit ci-dessous. (Cliquez pour agrandir l'image.) J'ai écrit ce tableau spécifiquement pour répondre à des questions telles que la vôtre.

Comparaison des séries de Fourier et de la transformée de Fourier.

Notez les similitudes et les différences entre les quatre opérations:

  1. "Série": périodique dans le temps, discrète en fréquence
  2. "Transform": apériodique dans le temps, continu en fréquence
  3. "Continuous Time": continu dans le temps, apériodique en fréquence
  4. "Discrete Time": discret dans le temps, périodique en fréquence

J'espère que ces notes vous seront utiles! N'hésitez pas à distribuer comme vous le souhaitez.


1
Bon résumé. Notez que la "série de Fourier à temps discret" référencée dans le tableau ci-dessus est généralement appelée la transformée de Fourier discrète (DFT).
Jason R

Pour taper un peu, cette réponse est en effet un bon résumé comme le dit Jason R, et quelque chose qui vaut la peine d'être en permanence sur dsp.SE pour que tout le monde puisse s'y connecter pour référence future, mais elle ne répond pas vraiment à la question qui s'est posée pour une explication intuitive de ces problèmes (la lucidité étant vraisemblablement un bonus supplémentaire et non absolument requise car elle est mentionnée dans le titre mais pas dans le texte de la question).
Dilip Sarwate,

2
Une excellente réponse Steve - Je pense que c'est ce que recherche le PO. Court, doux et précis.
Spacey

S'agit-il d'une erreur d'impression au bas de la page 2 de votre document? Il est dit: . Cela ne voulait-il pas dire ? - x ( t ) b ( t - t 0 ) d t = x ( t 0 )x(t)b(tt0)=x(t0)b(tt0)x(t)b(tt0)dt=x(t0)
mbaitoff

1
Pas une faute de frappe. Vos deux affirmations sont vraies, mais j'avais l'intention d'écrire la première car cette section du guide décrit les définitions axiomatiques de base de l'impulsion unitaire. La deuxième instruction est alors dérivée de ces définitions: . X(t)δ(t-t0)t=X(t0)δ(t-t0)t=X(t0)δ(t-t0)t=X(t0)
Steve Tjoa

9

Pour une explication lucide et correcte de ces concepts, vous devez parcourir certains des manuels standard (Oppenheim-Schafer, Proakis-Manolakis ou "Understanding Digital Signal Processing" de Richard Lyons qui est un très bon livre mais relativement moins populaire) . Mais en supposant une discussion à la table basse, je ferai des déclarations extrêmement vagues dans ce qui suit. :)

Pour un signal temporel continu général, vous ne vous attendriez pas à ce qu'une fréquence particulière soit absente, donc sa transformée de Fourier (ou la transformée de Fourier continue) serait une courbe continue avec un support éventuellement -inf à + inf.

Pour un signal continu périodique (période T), Fourier a exprimé le signal comme une combinaison de sinus et cosinus ayant la même période (T, T / 2, T / 3, T / 4, ...). En effet, le spectre de ce signal est une série de pointes aux emplacements 1 / T, 2 / T, 3 / T, 4 / T, ... C'est ce qu'on appelle la représentation de la série de Fourier. Il existe un théorème qui dit que la représentation en série de Fourier de tout signal temporel continu périodique converge vers le signal lorsque vous incluez de plus en plus de sinus et cosinus (ou exponentielles complexes) dans le sens carré moyen.

Moralité jusqu'à présent: périodicité dans le temps => spectre épineux

Passons au temps discret ... Que se passe-t-il si vous échantillonnez un signal horaire continu? Il doit être clair que pour un signal suffisamment élevé, vous ne pourrez pas reconstruire le signal. Si vous ne faites aucune hypothèse sur les fréquences dans le signal, étant donné le signal échantillonné, il n'y a aucun moyen de dire quel est le vrai signal. En d'autres termes, différentes fréquences sont représentées de manière équivalente dans le signal à temps discret. En parcourant quelques mathématiques, vous pouvez obtenir le spectre du signal échantillonné à partir du signal continu d'origine. Comment? Vous décalez le spectre du signal temporel continu de quantités + -1 / T, + -2 / T, ... et ajoutez toutes les copies décalées (avec une certaine mise à l'échelle). Cela vous donne un spectre continu qui est périodique avec la période 1 / T. (Remarque: le spectre est périodique en raison de l'échantillonnage dans le temps, le signal temporel ne t doivent être périodiques) Puisque le spectre est continu, vous pouvez aussi bien le représenter avec une seule de ses périodes. Il s'agit de la DTFT (Transformée de Fourier à "Temps Discret"). Dans le cas où votre signal de temps continu d'origine a des fréquences ne dépassant pas + -1 / 2T, les copies décalées du spectre ne se chevauchent pas et, par conséquent, vous pouvez récupérer le signal de temps continu d'origine en sélectionnant une période du spectre ( le théorème d'échantillonnage de Nyquist).

Une autre façon de se souvenir: signal temporel hérissé => périodicité du spectre

Que se passe-t-il si vous échantillonnez un signal périodique à temps continu avec une période d'échantillonnage T / k pour certains k? Eh bien, le spectre du signal en temps continu était épineux d'être avec, et l'échantillonner par un diviseur de T signifie que les pointes dans les copies décalées tombent exactement sur des multiples de 1 / T, donc le spectre résultant est un spectre périodique épineux . signal temporel périodique épineux <=> spectre périodique épineux (en supposant que la période et la fréquence d'échantillonnage sont "bien liées" comme ci-dessus.) C'est ce que l'on appelle la DFT (Discrete Fourier Transform). La FFT (Fast Fourier Transform) est une classe d'algorithmes pour calculer efficacement la DFT.

La façon dont la DFT est invoquée est la suivante: Supposons que vous souhaitiez analyser une séquence de N échantillons dans le temps. Vous pouvez prendre DTFT et traiter l'une de ses périodes, mais si vous supposez que votre signal est périodique avec la période N, alors DTFT se réduit à DFT et vous n'avez que N échantillons d'une période de DTFT qui caractérisent complètement le signal. Vous pouvez mettre à zéro le signal dans le temps pour obtenir un échantillonnage plus fin du spectre et (beaucoup plus de ces propriétés).

Tout ce qui précède n'est utile que s'il est accompagné d'une étude du DSP. Ce qui précède ne sont que quelques lignes directrices très approximatives.


7

Soit une fonction bornée de période T , c'est-à-dire pour tous les nombres réels t , x ( t + T ) = x ( t ) . À titre d'exemple particulier, cos ( 2 π t / T ) est une telle fonction. On veut trouver la "meilleure" approximation a n cos ( 2 π n t / T ) pour cette fonction où l'on souhaite choisir le coefficient a nX(t)TtX(t+T)=X(t)cos(2πt/T)unencos(2πnt/T)unende sorte que l'erreur quadratiqueest aussi petite que possible. En développant l'intégrande, nous avons une erreur quadratique = T 0 x 2 ( t )

0T(X(t)-unencos(2πnt/T))2t,
L'intégrale la plus à gauche est l'énergie E fournie par une période de x ( t ) tandis que l'intégrale la plus à droite a la valeur T / 2 et nous voyons donc que l' erreur quadratique = E - 2 a n T 0 x ( t ) cos ( 2 π n t / T )
erreur quadratique=0TX2(t)t-2unen0TX(t)cos(2πnt/T)t+(unen)20Tcos2(2πnt/T)t.
EX(t)T/2 À présent. poura>0, la fonction quadratiqueaz2+bz+ca un minimum àz=-b/2a (à mi-chemin entre les racines(-b/2a)±
erreur quadratique=E-2unen0TX(t)cos(2πnt/T)t+(unen)2T2.
une>0unez2+bz+cz=-b/2une!!) et ainsi, puisque nous avons exprimé l'erreur quadratique en fonction quadratique d'unn, le choix d'unnqui minimise l'erreur quadratique est unn=2(-b/2une)±b2-4unec/2uneunenunen De même, choisir b n comme b n = 2
unen=2T0TX(t)cos(2πnt/T)t.
bn minimise l'erreur quadratique entre x ( t ) et b n sin ( 2 π n t / T ) . Ainsi, nous voyons que la série de Fourier n'est rien d'autre qu'une astuce bon marché pour trouver l'approximation d'erreur quadratique minimale d'une fonction périodique x ( t ) en termes de signaux sinus et cosinus de la même période et de leurs harmoniques.
bn=2T0TX(t)péché(2πnt/T)t
X(t)bnpéché(2πnt/T)X(t)

4

Endolith a raison, si vous commencez réellement avec la série de Fourier et voyez comment elle est étendue à la transformée de Fourier, alors les choses commencent à avoir beaucoup de sens. Je donne une brève explication à ce sujet dans la première moitié de cette réponse .

Une bonne (peut-être pas simple) façon de regarder la famille de la transformée de Fourier (par laquelle je veux dire les 4 que vous avez énumérées ci-dessus), est à travers les lunettes de dualité Pontryagin . Cela vous permet de mémoriser les différentes transformations par les domaines d'origine et transformés.

RRRR

nZ/nZnZ/nZ

TZZT

Cette réponse n'est pas complète et je vais peut-être m'appuyer sur cette réponse pour clarifier quelques points lorsque j'en ai le temps, mais d'ici là, cela pourrait être quelque chose à mâcher jusqu'à ce que vous obteniez une explication plus intuitive de quelqu'un d'autre. Essayez également de lire des variantes de l'analyse de Fourier sur Wikipédia.


3

Je pense que la chose la plus importante est de comprendre fondamentalement pourquoi avons-nous besoin de transformations de Fourier. Ce sont l'une des nombreuses transformations de signaux possibles, mais aussi l'une des plus utiles. Une transformation transforme fondamentalement un signal dans un autre domaine qui peut nous donner un aperçu du signal dans ce domaine, ou il se peut que le domaine soit mathématiquement facile à travailler. Une fois que nous avons fini de travailler dans ce domaine, nous pouvons prendre la transformation inverse pour obtenir plus facilement le résultat souhaité.

Les blocs de construction les plus élémentaires de la théorie des fouriers sont les monotones (sinus et cosinus). Nous pouvons décomposer un signal en ses composantes de fréquence (monotones) en utilisant les mathématiques de Fourier. Ainsi, la transformée de Fourier transforme fondamentalement un signal du domaine temporel en domaine fréquentiel. Le coefficient de chacun des monotones de la série de Fourier nous renseigne sur la force de cette composante de fréquence dans le signal. Les transformées de Fourier (CFT, DFT) nous donnent explicitement une vue du domaine fréquentiel du signal. Dans la nature, les sinus et les cosinus sont les formes d'onde les plus importantes. Les signaux synthétiques comme les ondes carrées ou les signaux présentant de fortes fluctuations sont moins susceptibles de se produire naturellement et ne s'étonnent pas, de façon surprenante, de la gamme infinie de fréquences, comme l'expliquent très clairement les transformées de Fourier. Les gens doutaient qu'un signal puisse être transformé en somme de sinus / cosinus. Fourier a montré que la forme d'onde carrée (qui est loin des sinus / cosinus) peut en effet l'être. Le bruit blanc contient toutes les fréquences de force égale.

De plus, si vous travaillez avec des séries de Fourier, les coefficients ainsi que le terme de phase peuvent être considérés comme ceux requis pour superposer correctement les formes d'onde sinosoïdales constitutives de sorte que la superposition est en effet le signal requis dont vous prenez la transformation. Lorsque vous travaillez avec des transformées de Fourier, les nombres complexes ont implicitement les termes de phase et la magnitude requise de chacun des monotones. (l'intégration est à peu près comme la sommation. continue => intégration, discrète => sommation)

Je pense qu'une fois que vous avez compris le thème d'un concept, tout le reste n'est que des détails que vous devrez vous-même comprendre en lisant des livres. La lecture de l'application des transformées de Fourier à divers domaines vous donnera une meilleure perception.


2

Un DFT est une transformation d'un vecteur de paires de nombres d'un espace orthogonal à un autre. Très couramment fait comme un calcul numérique. Pour une raison quelconque, lorsque vous prenez un groupe de chiffres du monde réel, le 2e groupe de chiffres s'avère souvent assez proche de quelque chose de très utile.

Je me souviens de l' efficacité déraisonnable des mathématiques dans les sciences naturelles , en particulier en ce qui concerne l'application de la DFT à de nombreux systèmes qui semblent être approximés par différents types d'équation différentielle du 2e degré, même le son de la cuillère à café que je viens de déposer.

Les 3 autres XYZ-FT font des hypothèses sur l'existence de certaines entités infinies mythiques pour aider les solutions symboliques à s'intégrer sur le tableau blanc avant que le café ne refroidisse trop. Ce sont les "vaches sphériques" du traitement du signal. Les séries DTFT et Fourier prétendent qu'un vecteur peut être étendu à l'infini au prix d'une densité infinie de l'autre entité. La série de Fourier prétend que les deux entités peuvent être des fonctions continues infinies.

Prenez suffisamment de cours de mathématiques et on pourrait même déterminer toutes les définitions et hypothèses nécessaires pour rendre ces entités fictives exactes et complètes en deux sens.


Qu'entend-on par «espace orthogonal» dans votre première phrase? Quel est l'espace orthogonal à , ou quelle propriété spéciale possède l'espace que vous le distinguez des autres espaces ordinaires en lui conférant l'adjectif «orthogonal»?
Dilip Sarwate

Peut-être "orthonormal" le terme le plus correct pour les espaces vectoriels?
hotpaw2

XyX,y=0UNEUNEUNETUNEUNETles vecteurs dans l'espace sont orthogonaux entre eux ou sont orthogonaux et ont aussi une longueur unitaire? Si oui, pouvez-vous donner un exemple d'un tel espace?
Dilip Sarwate

Le produit scalaire entre tous les sinus ou cosinus qui sont exactement périodiques dans une longueur d'ouverture DFT est zéro, à l'exception des fonctions de fréquence identiques. Même si N est supérieur au nombre de grains de café dans le sac. Faites-leur une amplitude unitaire pour l'orthonormal.
hotpaw2

NNN NN
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.