¿Dónde está la falla en esta derivación de la DTFT de la secuencia de pasos unitarios


11

Esta pregunta está relacionada con esta otra pregunta mía donde solicito derivaciones de la transformada de Fourier de tiempo discreto (DTFT) de la secuencia de pasos unitarios u[n] . Durante mi búsqueda de derivaciones encontré una que es asombrosamente simple. Lo vi por primera vez en la página 138 de este libro de BA Shenoi. También lo encontré en matemáticas . SE en esta respuesta .

Como el argumento es corto y simple, lo repetiré aquí por conveniencia.

La secuencia de pasos unitarios se puede escribir como

(1)u[n]=f[n]+12
con
(2)f[n]={12,n012,n<0
Obviamente,
(3)f[n]f[n1]=δ[n]
Aplicando el DTFT en ambos lados de(3)da
(4)F(ω)(1ejω)=1
dondeF(ω)es la DTFT def[n]. De(4)obtenemos
(5)F(ω)=11ejω
De(5)y(1)obtenemos la DTFT deu[n]
(6)U(ω)=F(ω)+πδ(ω)=11ejω+πδ(ω),πω<π
donde he usadoDTFT{1}=2πδ(ω),πω<π.

Eq. (6) para la DTFT de u[n] es sin duda correcta. Sin embargo, la derivación es defectuosa.

La pregunta es: encontrar y explicar la falla en la derivación anterior.

Anteponga su respuesta con la etiqueta de spoiler >!.


1
Lo que me perturba es que es una señal de potencia finita , no una señal de energía finita , que es lo que obtenemos cuando sumamos estas dos señales de energía infinita. f[n]
robert bristow-johnson

Además, ¿no es ?
DTFT{x[n]=1}=2πk=+δ(ω2kπ)
robert bristow-johnson

Gracias a todos por sus respuestas! Los he votado a todos, y cada uno de ellos resulta en una agradable discusión sobre aspectos no tan conocidos del DTFT de señales extrañas (es decir, los que no están en o 2 ). Solo puedo aceptar una, y esperaré un poco más por nuevas respuestas o cambios en las respuestas existentes. También agregaré mi propia respuesta más adelante. 12
Matt L.

1
Matt, definitivamente no es energía finita. un número infinito de muestras que cuadran para ser 1f[n] no suman para ser un número finito. 14
Robert Bristow-Johnson

1
@ robertbristow-johnson: ¿Qué es lo que te preocupa de eso? Si las señales se cancelan entre sí en todas partes, excepto por un número finito de puntos, entonces eso es lo que obtenemos.
Matt L.

Respuestas:


7

Hay infinitas señales que hacen que se cumpla la siguiente igualdad:

y[n]y[n1]=δ[n](1)
Lo único que importa es quey[0]y[1]=1 , y luego el resto de los coeficientes dey pueden determinarse bajo la restricción de que Eq. (1) estados (es decir, la sustracción de muestras consecutivas debe ser0 paran0 ). En otras palabras, la ecuación. (1) se logrará mediante cualquier señaly[n] tal que
y[0]=y[1]+1y[n]=y[n1] n0
Otra forma de ver esto es que cualquier función que sea básicamenteu[n] con un desplazamiento (un valor constante agregado) satisfará(1) . Esto explica la declaración hecha porRobert Bristow-Johnson en su respuesta : los diferenciadores destruyen esta información (como tomar una derivada en tiempo continuo destruye la evidencia de cualquier valor constante en la función original).

En resumen, creo que la prueba es defectuosa porque el procedimiento seguido podría usar cualquier función de la forma u[n]+C con CR , y esto conduciría a que muchas funciones tengan la misma transformada de Fourier, lo que de hecho es incorrecto como la transformada de Fourier es una biyección. Tal vez el autor decidió deliberadamente ignorar cualquier cosa relacionada con los valores DC, consciente de que para mostrar que F(ω) es la DTFT de f[n]necesitaría la propiedad de acumulación (cuya prueba más popular se deriva de la DTFT del paso unitario, ergo, una prueba bastante circular). La prueba no es estrictamente incorrecta , ya que todo lo que dice (las fórmulas para F(ω) y U(ω) , la descomposición del paso unitario, la ecuación de diferencia) es cierto, pero requeriría la propiedad de acumulación para mostrar por qué F(ω) no tiene ningún delta de Dirac.


Estás totalmente en el camino correcto! ¿Tiene una idea de cómo se podría resolver esta falla, es decir, cómo hacerlo correctamente?
Matt L.

@MattL. Establecer una condición inicial para haría el truco y determinaría la señal de forma unívoca. Esa condición inicial determinaría el valor de CC de la señal y [ n ] , que aparece en la DTFT como una constante multiplicando un impulso de Dirac (de acuerdo con la propiedad de acumulación). Creo que en la prueba dada, esto funciona porque la señal f [ n ] no tiene valor de CC ya que es simétrica alrededor de 0 , por lo que el DTFT es correcto en ese caso. Pero el hecho de que la señal no tiene DC debería ser declarado, ya que es fundamental, creo. y[n]y[n]f[n]0
Tendero

Hay muchas buenas respuestas y es difícil elegir cuál aceptar. Pero esta fue muy apreciada por la comunidad, y también creo que señala más claramente el error en la derivación. ¡Gracias a todos!
Matt L.

4

Estaba abrumado por la cantidad de respuestas que recibí (¡10 respuestas hasta ahora!). Por supuesto, todos ellos obtuvieron mi voto positivo. Esto fue divertido, gracias a todos por sus pensamientos, comentarios, etc. Sé que a estas alturas la mayoría de ustedes saben cuál es la falla, al menos a la que me refería. Las personas expresan las cosas de manera diferente, y siempre hay lugar para malentendidos, por lo que intentaré formular claramente lo que creo que es la falla más importante en esa derivación. Soy consciente del hecho de que no todos estarán de acuerdo y eso está bien. ¡Estoy feliz de poder discutir este tipo de temas DSP esotéricos con mentes tan agudas como ustedes! Aquí vamos.

Mi primer reclamo es que todas y cada una de las ecuaciones de mi pregunta son correctas. Sin embargo, la derivación y motivación de algunos de ellos es totalmente errónea y engañosa, y esa "derivación" solo puede existir porque el autor sabía cómo se suponía que debía ser el resultado.

Eq. (3) en la pregunta ( ) es correcta para la secuencia dada f [ n ] (Ec. ( 2 ) en la pregunta), pero claramente también correcto para todas las secuencias de la forma f [ n ] = u [ n ] + c con alguna constante arbitraria c . Entonces, de acuerdo con la derivación, la DTFT resultante Ff[n]f[n1]=δ[n]f[n](2)

(1)f[n]=u[n]+c
cde la pregunta? debe ser la DTFT de todas las secuencias de la forma ( 1 ) , independientemente del valor de la constante c . Por supuesto, eso no tiene sentido porque el DTFT es único. Específicamente, usando esa misma "prueba" podría "mostrar" que F ( ω ) como se indica en la ecuación. ( 5 ) de mi pregunta (o la ecuación ( 3 ) a continuación) es en realidad la DTFT de u [ n ] que estamos buscando. Entonces, ¿por qué molestarse en dividir u [ n ] como en la ecuación? ( 1 )F(ω)(1)cF(ω)(5)(3)u[n]u[n](1)

Sin embargo, es cierto que las DTFT de todas las secuencias satisfacen la ecuación. ( 4 ) en la pregunta (repetida aquí por conveniencia): F ( ω ) ( 1 - e - j ω ) = 1 Pero ahora viene la falla matemática real: de la ecuación. ( 2 ) es incorrecto concluir F ( ω ) = 1(1)(4)

(2)F(ω)(1ejω)=1
(2) Eq. (3)es solo una de las infinitas soluciones posibles de(2), y resulta conveniente que sea la que necesita el autor para llegar al resultado final correcto. Eq. (3)es la DTFT def[n]en(1)conc=-1
(3)F(ω)=11ejω
(3)(2)(3)f[n](1) , pero a partir de la derivación dada no hay forma de saberlo.c=12

Entonces, ¿cómo podemos evitar ese error matemático y usar para derivar las DTFT de a l l secuencias ( 1 ) , con cualquier constante c ? La conclusión correcta de ( 2 ) es F ( ω ) = 1(2)all(1)c(2)con alguna constanteαaún indeterminada. Al enchufar(4)en el lado izquierdo de(2) seobtiene1+α(1-e-jω)δ(ω)=1+α(1-e-jω)| ω=0δ

(4)F(ω)=11ejω+αδ(ω)
α(4)(2) Entonces, todas las funciones F ( ω ) dadas por ( 4 ) satisfacen ( 2 )
1+α(1ejω)δ(ω)=1+α(1ejω)|ω=0δ(ω)=1+0δ(ω)=1
F(ω)(4)(2) , según sea necesario.

La constante en ( 4 ) se puede determinar a partir del valor de f [ n ] en n = 0 : f [ 0 ] = 1 + c = 1α(4)f[n]n=0

(6)f[0]=1+c=12πππF(ω)dω=12πππdω1ejω+α2π
(6)
(7)PVππdω1ejω=π
(6)(7)
(8)α=π(1+2c)
c=12α=0f[n]c=0f[n]=u[n]α=πu[n]
(9)U(ω)=11ejω+πδ(ω)


F(ω)F(ω)

(4)(2)ω=0ω=0(1ejω)(1ejω)

No estoy realmente seguro de que no pueda haber otra función que no sea el impulso delta de Dirac (y sus derivados) que tenga esta propiedad. Pero está bien, su respuesta está bien escrita. Voto a favor. Gracias.
AlexTP

2

La falla sigue a la palabra "Obviamente", si se supone que esa es la función Delta de Dirac.

Aquí está el borrador de una respuesta para su otra pregunta que nunca publiqué:

-------------------------------------------------- -------------

No creo que sea posible una prueba. Este puede ser un caso de una "definición funcional" que tiene las propiedades deseadas.

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

Probar que la definición funciona de manera deseable es una cuestión diferente.

La prueba de la página 138 es incorrecta (al menos) porque:

δ(t)=lima012a[u(t+a)u(ta)]=dudt
δ(n)=u2(n)u2(n1)

Situación interesante, espero que esto ayude. Estoy ansioso por lo que tienes que decir.

Sección de la economía


δ[n]n=01

U=11ejωlimN[ejωN1ejω]

2

1=2

F(ω)(1ejω)=1
ω=2kπ for kZ

3
u[n]u[n]12

w2πk

δ[n] ωω=2kπ

δ[n]1ejwF(w)

f[n]f[n1]=δ[n]
u[n]

2

limNn=Nn=Nf[n]ejωn(1ejωn)+(ejωNf[N]+ejωNf[N])ejω=1
f[n]=u[n]
limNn=Nn=Nf[n]ejωn(1ejωn)+ejωNejω=1

F(ω)ω2kπk

F(ω)ω=2kπn=1sin(ωn)ωω=2kπU(ω)

F(ω)

u[n]

1
u[n]f[n]u[n]u[n1]=δ[n]

2

Creo que he descubierto la mejor manera de expresar la falla en esta prueba. Así que voy a darle otra puñalada.

12x

U(ω)=11ejω+2πxδ(ω)

x12

Además, si toma el paso que hice en mi última respuesta y encuentra (4) se expresa como

F(ω)(1ejω)=1+2πx(1ejω)δ(ω)

Seguido de incluirlo en (5) y (6) obtienes:

U(ω)=11ejω+4πxδ(ω)

Lo cual, como señalé anteriormente, es inconsistente con la definición de llegar allí.

x=12πδ(ω)

x=12

Sección de la economía


1

Esto es en respuesta a los comentarios en mi primera respuesta. Debido al encubrimiento del spoiler, lo publico como una respuesta separada.

Iba a publicar mi otra respuesta a la otra pregunta, pero no lo hice debido a mi falta de experiencia en esta área. Lo publiqué ayer, lo eliminé, luego lo borré, luego descubrí cómo emplear etiquetas de spoiler.

δδδp

δp[n]=f[n]f[n1]=u[n]u[n1]

Tomando el DTFT de las partes izquierda y derecha. No estoy seguro de tener las anotaciones correctas, pero las matemáticas deberían ser claras. Usando la definición que se está probando.

Fp(ω)=Fu(ω)Fu(ω)ejω

Fp(ω)=[11ejω+πδ(ω)][ejω1ejω+(πejω)δ(ω)]

Fp(ω)=1ejω1ejω+π(1ejω)δ(ω)

Fp(ω)=1+π(1ejω)δ(ω)1

ω=2kπω=2kπ

Sección de la economía

==============================

Seguimiento:

δp


Fp(ω)=1(1ejω)δ(ω)f(ω)ω=0f(ω)δ(ω)=f(0)δ(ω)f(0)=0

1

(a[n]b[n])cω[n]=a[n]cω[n]b[n]cω[n]
12f[n]f[n1]

1

entonces Matt

f[n]

f[n]{ 12eαnn012eαnn<0

α>0

ahora tenemos señales de energía finita y las DTFT deberían ser comparables.

f[n]f[n1]={12(eαneα(n1))n>012(1+eα)n=012(eαneα(n1))n<0={12(1eα)eαnn>012(1+eα)n=012(eα1)eαnn<0

α0

pero, por desgracia, son casi las 2 am y no voy a lidiar con eso ahora.


αα0
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.