Transformada de Fourier de tiempo discreto de la secuencia de pasos unitarios

10

De los libros de texto sabemos que la DTFT de viene dada poru[n]

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

Sin embargo, no he visto un libro de texto DSP que al menos pretenda dar una derivación más o menos sólida de .(1)

Proakis [1] deriva la mitad derecha del lado derecho de al establecer en la de , y dice que es válido excepto (que por supuesto es correcto). Luego declara que en el polo de la tenemos que agregar un impulso delta con un área de , pero eso me parece más una receta que cualquier otra cosa.z = e j ω Z u [ n ] ω = 2 π k Z π(1)z=ejωZu[n]ω=2πkZπ

Oppenheim y Schafer [2] mencionan en este contexto

Aunque no es del todo sencillo de mostrar, esta secuencia se puede representar mediante la siguiente transformada de Fourier:

que es seguido por una fórmula equivalente a . Desafortunadamente, no se tomaron la molestia de mostrarnos esa prueba "no completamente sencilla".(1)

Un libro que en realidad no conocía, pero que encontré al buscar una prueba de es Introducción al procesamiento de señales digitales y diseño de filtros de BA Shenoi. En la página 138 hay una "derivación" de , pero desafortunadamente está mal. Hice una pregunta de "DSP-puzzle" para que la gente muestre lo que está mal con esa prueba.](1)(1)

Entonces mi pregunta es:

¿Alguien puede proporcionar una prueba / derivación de que sea sólida o incluso rigurosa mientras sea accesible para ingenieros con inclinaciones matemáticas? No importa si solo se ha copiado de un libro. Creo que sería bueno tenerlo en este sitio de todos modos.(1)

Tenga en cuenta que incluso en matemáticas. SE puede encontrar casi nada relevante: esta pregunta no tiene respuestas, y una tiene dos respuestas, una de las cuales es incorrecta (idéntica al argumento de Shenoi), y la otra usa la "propiedad de acumulación" , con lo que estaría contento, pero luego hay que probar esa propiedad, lo que lo pone de nuevo al principio (porque ambas pruebas básicamente prueban lo mismo).

Como nota final, se me ocurrió algo así como una prueba (bueno, soy ingeniero), y también lo publicaré como respuesta algunos días a partir de ahora, pero me encantaría recopilar otras pruebas publicadas o no publicadas. que son simples y elegantes y, lo más importante, que son accesibles para los ingenieros de DSP.

PD: No dudo de la validez de , solo me gustaría ver una o varias pruebas relativamente sencillas.(1)


[1] Proakis, JG y DG Manolakis, Procesamiento de señal digital: principios, algoritmos y aplicaciones , 3a edición, Sección 4.2.8

[2] Oppenheim, AV y RW Schafer, Procesamiento de señal de tiempo discreto , 2ª edición, p. 54)



Inspirado por un comentario de Marcus Müller, me gustaría mostrar que como lo da la ecuación. cumple el requisitoU(ω)(1)

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

Si es la DTFT de , entoncesU(ω)u[n]

V(ω)=11ejω

debe ser la DTFT de

v[n]=12sign[n]

(donde definimos ), porquesign[0]=1

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

Entonces tenemos

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

de donde se deduce que

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

Con esto obtenemos

12π(UU)(ω)=12π[(πδ(ω)+V(ω))(πδ(ω)+V(ω))]=12π[π2δ(ω)+2πV(ω)+(VV)(ω)]=π2δ(ω)+V(ω)+π2δ(ω)=U(ω)q.e.d.
Matt L.
fuente
waaah No rompas mi mundo La duda en esa fórmula introduce un reino de caos. Por ejemplo, , y por lo tanto (con un prefactor de definición de FT cont. Que depende de 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: No hay duda sobre esa fórmula, es correcta. La pregunta es cómo mostrarlo de una manera que un ingeniero de mente simple pueda entender. Y funciona para el DTFT dado, no hay problema. u2[n]=u[n]
Matt L.
Me considero muy simple, y eso significa que me preocupo cuando las cosas no se sienten "seguras" cuando no puedo ver cómo se derivan.
Marcus Müller
1
Veo que lo que busca no es probar si la ecuación es correcta o no, sino más bien derivar de manera rigurosa y directa de los primeros principios y la definición de DTFT. Entonces, cada vez que se quiera hacer una prueba rigurosa que implique impulsos , supongo que deberíamos referirnos mejor a los libros citados de la teoría de la función generalizada: Lighthill-1958 se cita en Opp & Schafer para una discusión sobre la función del impulso y su uso en las transformadas de Fourier. Todas las demás pruebas dependerán inevitablemente de las pruebas realizadas en esas referencias y serán insuficientes para reemplazar una prueba rigurosa. U(w)
Fat32
1
@ Fat32: Ese es un punto de vista válido. Sin embargo, creo que es posible una derivación razonablemente sólida si aceptamos transformaciones básicas como , y si estamos contentos de definir integrales por su valor principal de Cauchy. DTFT{1}=2πδ(ω)
Matt L.

Respuestas:

3

Cedron Dawg publicó un punto inicial interesante en esta respuesta . Comienza con estos pasos:

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

Resulta que el término dentro del límite se puede expandir de la siguiente manera :

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

El factor común fuera de los corchetes se puede expresar como :

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

La parte real dentro de los corchetes también es igual a :

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

Por otro lado, la parte imaginaria se puede reescribir como :

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

Reescribiendo el término original obtenemos 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)

donde utilicé y el límite tampoco se ve afectado como .M=N1M

De acuerdo con la séptima definición en este sitio :

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

Hasta ahora tenemos eso:

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

Si pudiéramos demostrar que el segundo término a la derecha de la igualdad es en algún sentido, entonces hemos terminado. Lo pregunté en matemáticas. SE y, de hecho, esa secuencia de funciones tiende a la distribución cero. Entonces, tenemos eso:0

U(ω)=11ejωlimN[ejωN1ejω]=11ejω+πδ(ω)+jlimMcos[ω(M+1/2)]2sin(ω/2)=11ejω+πδ(ω)
Tendero
fuente
¡Esto esta muy bien! Lo revisé y todo parece ser correcto, por lo que esa parte imaginaria debe tender a cero en algún sentido. Lo pensaré por un momento.
Matt L.
@MattL. ¡Avíseme si puede progresar!
Tendero
@MattL. ¡La prueba finalmente está completa!
Tendero
¡Buen trabajo! Había descubierto que el término coseno tendería a cero debido al lema de Riemann-Lebesgue, pero mi problema era el caso . Porque la primera fórmula se basa en la suma geométrica, que solo es válida para . Después de todo, todo funciona de alguna manera, pero sigue siendo un defecto menor. Tengo otra derivación que no divide el término , en el que el caso se maneja con un poco más de cuidado, pero sigue siendo una "prueba de ingeniero" . Podría publicarlo cuando tenga más tiempo. ω=0ω01/(1ejω)ω=0
Matt L.
2

Proporcionaré dos pruebas relativamente simples que no requieren ningún conocimiento de la teoría de distribución. Para obtener una prueba que calcule la DTFT mediante un proceso de límite utilizando los resultados de la teoría de distribución, vea esta respuesta de Tendero .

Solo mencionaré (y no detallaré) la primera prueba aquí, porque la publiqué como respuesta a esta pregunta , cuyo propósito era mostrar que cierta prueba publicada es defectuosa.

u[n]

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

(1)

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

u[n]

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

u[n]U(ω)UR(ω)U(ω)h[n]=h[n]u[n]h[n]H(ω)=12π(HU)(ω)u[n]U(ω)

UI(ω)=Im{U(ω)}u[n]

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

(4)

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

(3)(5)

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

(6)

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

UI(ω)ωu[n]c=0(3)(7)

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