Transformée de Fourier à temps discret de la séquence de pas unitaire

10

D'après les manuels, nous savons que la DTFT de est donnée paru[n]

(1)U(ω)=πδ(ω)+11ejω,πω<π

Cependant, je n'ai pas vu de manuel DSP qui prétend au moins donner une dérivation plus ou moins saine de .(1)

Proakis [1] dérive la moitié droite du côté droit de en définissant dans la transformation de , et dit qu'il est valide sauf pour (ce qui est bien sûr correct). Il déclare ensuite qu'au pôle de la transformation nous devons ajouter une impulsion delta avec une zone de , mais cela me semble plus une recette qu'autre chose.(1)z=ejωZu[n]ω=2πkZπ

Oppenheim et Schafer [2] mentionnent dans ce contexte

Bien qu'il ne soit pas complètement simple à montrer, cette séquence peut être représentée par la transformée de Fourier suivante:

qui est suivie d'une formule équivalente à . Malheureusement, ils n'ont pas pris la peine de nous montrer cette preuve "pas tout à fait simple".(1)

Un livre que je ne connaissais pas, mais que j'ai trouvé en recherchant une preuve de est Introduction to Digital Signal Processing and Filter Design par BA Shenoi. À la page 138, il y a une "dérivation" de , mais malheureusement c'est faux. J'ai posé une question "DSP-puzzle" pour que les gens montrent ce qui ne va pas avec cette preuve.]( 1 )(1)(1)

Donc , ma question est:

Quelqu'un peut-il fournir une preuve / dérivation de qui soit solide ou même rigoureuse tout en étant accessible aux ingénieurs mathématiquement enclins? Peu importe qu'il soit simplement copié d'un livre. Je pense que ce serait bien de l'avoir sur ce site de toute façon.(1)

Notez que même sur math.SE, il n'y a presque rien de pertinent: cette question n'a pas de réponses, et l' une a deux réponses, dont l'une est fausse (identique à l'argument de Shenoi), et l'autre utilise la "propriété d'accumulation" , ce qui me ferait plaisir, mais il faut ensuite prouver cette propriété, ce qui vous ramène au début (car les deux preuves prouvent essentiellement la même chose).

Pour finir, j'ai trouvé quelque chose comme une preuve (enfin, je suis ingénieur), et je vais également la poster comme réponse dans quelques jours, mais je serais heureux de collecter d'autres preuves publiées ou non. qui sont simples et élégants, et, surtout, qui sont accessibles aux ingénieurs DSP.

PS: Je ne doute pas de la validité de , je voudrais juste voir une ou plusieurs preuves relativement simples.(1)


[1] Proakis, JG et DG Manolakis, Traitement numérique du signal: principes, algorithmes et applications , 3e édition, section 4.2.8

[2] Oppenheim, AV et RW Schafer, Discrete-Time Signal Processing , 2e édition, p. 54.



Inspiré par un commentaire de Marcus Müller, je voudrais montrer que donné par Eq. satisfait à l'exigenceU(ω)(1)

u[n]=u2[n]U(ω)=12π(UU)(ω)

Si est la DTFT de , alorsU(ω)u[n]

V(ω)=11-e-jω

doit être le DTFT de

v[n]=12signe[n]

(où nous définissons ), carsigne[0]=1

V(ω)=U(ω)-πδ(ω)u[n]-12=12signe[n]

Donc nous avons

12π(VV)(ω)(12sign[n])2=14

d'où il résulte que

12π(VV)(ω)=DTFT{14}=π2δ(ω)

Avec cela, nous obtenons

12π(UU)(ω)=12π[(πδ(ω)+V(ω))(πδ(ω)+V(ω))]=12π[π2δ(ω)+2πV(ω)+(VV)(ω)]=π2δ(ω)+V(ω)+π2δ(ω)=U(ω)q.e.d.
Matt L.
la source
waaah. Ne brise pas mon monde. Le doute dans cette formule introduit un royaume de chaos. Par exemple, , et donc (avec un préfacteur de définition FT suivant la constante ),c DTFT ( u 2 ) ( ω )u2(t)=u(t)c
DTFT(u2)(ω)=cU(ω)U(ω)=cπU(ω)+c1+ejωU(ω)=cπ(πδ(ω)+11ejω)+cπ1+ejω+c11+ejω11+ejω=cπ2δ(ω)+2cπ1+ejω+c11+ejω11+ejω=magic?U
Marcus Müller
@ MarcusMüller: Aucun doute sur cette formule, c'est correct. La question est simplement de savoir comment le montrer d'une manière qu'un ingénieur simple d'esprit peut comprendre. Et fonctionne pour le DTFT donné, pas de problème. u2[n]=u[n]
Matt L.
Je me considère très simple d'esprit, et cela signifie que je m'inquiète lorsque les choses ne me semblent pas "en sécurité" lorsque je ne vois pas comment elles sont dérivées.
Marcus Müller
1
Je vois que ce que vous recherchez n'est pas de prouver si l'équation est correcte ou non, mais plutôt de dériver rigoureusement et directement des premiers principes et de la définition de DTFT. Ensuite, chaque fois que l'on veut faire une preuve rigoureuse impliquant des impulsions, je suppose qu'il vaut mieux se référer aux livres cités de la théorie de la fonction généralisée: Lighthill-1958 est cité dans Opp & Schafer pour une discussion sur la fonction d'impulsion et son utilisation dans les transformées de Fourier. Toutes les autres preuves s'appuieront inévitablement sur les preuves faites sur ces références et seront insuffisantes pour remplacer une preuve rigoureuse. U(w)
Fat32
1
@ Fat32: C'est un point de vue valide. Je pense, cependant, qu'une dérivation raisonnablement saine est possible si nous acceptons les transformations de base telles que , et si nous nous contentons de définir des intégrales par leur valeur principale de Cauchy. DTFT{1}=2πδ(ω)
Matt L.

Réponses:

3

Cedron Dawg a publié un point initial intéressant dans cette réponse . Cela commence par ces étapes:

U(ω)=n=0+ejωn=limNn=0N1ejωn=limN[1ejωN1ejω]=11ejωlimN[ejωN1ejω]

Il s'avère que le terme à l'intérieur de la limite peut être développé comme suit :

ejωN1ejω=1sin2(ω)+(1cos(ω))2[cos(ω)cos(Nω)+cos(Nω)sin(ω)sin(Nω)+j(sin(ω)cos(Nω)+cos(ω)sin(ω)sin(Nω))]

Le facteur commun en dehors des parenthèses peut être exprimé comme suit :

1sin2(ω)+(1cos(ω))2=14sin2(ω/2)

La partie réelle à l'intérieur des crochets est également égale à :

cos(ω)cos(Nω)+cos(Nω)sin(ω)sin(Nω)=2sin(ω/2)sin[ω(N+1/2)]

En revanche, la partie imaginaire peut être réécrite comme :

sin(ω)cos(Nω)+cos(ω)sin(ω)sin(Nω)=2sin(ω/2)cos[ω(N+1/2)]

En réécrivant le terme d'origine, nous obtenons que:

ejωN1ejω=2sin(ω2)4sin2(ω2)(sin[ω(N+1/2)]jcos[ω(N+1/2)])=sin[ω(M+1/2)]2sin(ω2)jcos[ω(M+1/2)]2sin(ω2)

où j'ai utilisé et la limite reste également inchangée comme .M M=N1M

Selon la 7ème définition de ce site :

limM12sin(ω/2)sin[ω(M+1/2)]=πδ(ω)

Jusqu'à présent, nous avons cela:

limMejω(M+1)1ejω=πδ(ω)jlimMcos[ω(M+1/2)]2sin(ω/2)

Si nous pouvions prouver que le deuxième terme à droite de l'égalité est dans un certain sens, alors nous avons terminé. Je l'ai demandé à math.SE et, en effet, cette séquence de fonctions tend vers la distribution zéro. Donc, nous avons cela:0

U(ω)=11ejωlimN[ejωN1ejω]=11ejω+πδ(ω)+jlimMcos[ω(M+1/2)]2sin(ω/2)=11ejω+πδ(ω)
Tendero
la source
C'est très gentil! Je l'ai vérifié et tout semble être correct, de sorte que la partie imaginaire doit tendre à zéro dans un certain sens. J'y penserai un peu.
Matt
@MattL. Faites-moi savoir si vous êtes en mesure de progresser!
Tendero
@MattL. La preuve est enfin complète!
Tendero
Bon travail! J'avais compris que le terme cosinus tendrait à zéro en raison du lemme de Riemann-Lebesgue, mais mon problème était le cas . Parce que la toute première formule est basée sur la somme géométrique, qui n'est valable que pour . Tout fonctionne d'une manière ou d'une autre après tout, mais c'est toujours un défaut mineur. J'ai une autre dérivation qui ne divise pas le terme , dans laquelle le cas est traité avec un peu plus de soin, mais c'est toujours une "preuve d'ingénieur" . Je pourrais l'afficher quand j'aurai plus de temps. ω=0ω01/(1ejω)ω=0
Matt L.
2

Je fournirai deux preuves relativement simples qui ne nécessitent aucune connaissance de la théorie de la distribution. Pour une preuve qui calcule la DTFT par un processus limite en utilisant les résultats de la théorie de la distribution, voir cette réponse de Tendero .

Je ne mentionnerai ici (et ne développerai pas) la première preuve ici, car je l'ai postée comme réponse à cette question , dont le but était de montrer qu'une certaine preuve publiée est défectueuse.

L'autre preuve va comme suit. Écrivons d'abord la partie paire de la séquence de pas unitaire :u[n]

(1)ue[n]=12(u[n]+u[n])=12+12δ[n]

Le DTFT de est(1)

(2)DTFT{ue[n]}=πδ(ω)+12

ce qui équivaut à la partie réelle de la DTFT de :u[n]

(3)UR(ω)=Re{U(ω)}=πδ(ω)+12

Puisque est une séquence à valeurs réelles, nous avons terminé parce que les parties réelle et imaginaire de sont liées via la transformée de Hilbert et, par conséquent, détermine uniquement . Cependant, dans la plupart des textes DSP, ces relations de transformée de Hilbert sont dérivées de l'équation (qui est valable pour toute séquence causale ), dont il résulte que . Donc, pour montrer la relation de transformation de Hilbert entre les parties réelles et imaginaires du DTFT, nous avons besoin du DTFT deu[n]U(ω)UR(ω)U(ω)h[n]=h[n]u[n]h[n]H(ω)=12π(HU)(ω)u[n], que nous voulons réellement dériver ici. La preuve devient donc circulaire. C'est pourquoi nous choisirons une manière différente de dériver la partie imaginaire de .U(ω)

Pour dériver nous écrivons la partie impaire de comme suit:UI(ω)=Im{U(ω)}u[n]

(4)uo[n]=12(u[n]u[n])=u[n1]12+12δ[n]

Prendre la DTFT de donne(4)

jUI(ω)=ejωU(ω)πδ(ω)+12=ejω(UR(ω)+jUI(ω))πδ(ω)+12=ejω(πδ(ω)+12)+ejωjUI(ω)πδ(ω)+12(5)=12(1+ejω)+ejωjUI(ω)

où j'ai utilisé . Eq. peut s'écrire(3)(5)

(6)jUI(ω)(1ejω)=12(1+ejω)

La conclusion correcte de est (voir cette réponse pour plus de détails)(6)

(7)jUI(ω)=121+ejω1ejω+cδ(ω)

Mais comme nous savons que doit être une fonction impaire de (parce que une valeur réelle), nous pouvons immédiatement conclure que . Par conséquent, à partir de et nous obtenons finalementUI(ω)ωu[n]c=0(3)(7)

U(ω)=UR(ω)+jUI(ω)=πδ(ω)+12+121+ejω1ejω=πδ(ω)+12(1+1+ejω1ejω)(8)=πδ(ω)+11ejω
Matt L.
la source