Onde está a falha nesta derivação do DTFT da sequência de etapas unitárias


11

Esta questão está relacionada a outra questão minha, em que peço derivações da transformada de Fourier em tempo discreto (DTFT) da sequência de etapas unitárias . Durante minha busca por derivações, encontrei uma que é incrivelmente simples. Vi pela primeira vez na página 138 deste livro de BA Shenoi. Eu também me deparei com isso em mathematics.SE nesta resposta .u[n]

Como o argumento é curto e simples, repetirei aqui por conveniência.

A sequência de etapas da unidade pode ser escrita como com Obviamente, Aplicando o DTFT nos dois lados de dá que é o DTFT de . De obtemos De e obtemos o DTFT de

(1)u[n]=f[n]+12
(2)f[n]={12,n012,n<0
(3)f[n]f[n1]=δ[n]
(3)
(4)F(ω)(1ejω)=1
F(ω)f[n](4) (5)(1)u[n]U(ω)=F(ω)+πδ(ω)=1
(5)F(ω)=11ejω
(5)(1)u[n]
(6)U(ω)=F(ω)+πδ(ω)=11ejω+πδ(ω),πω<π
onde usei , .- π ω < πDTFT{1}=2πδ(ω)πω<π

Eq. para o DTFT de é sem dúvida correto. No entanto, a derivação é falha.u [ n ](6)u[n]

A questão é: encontre e explique a falha na derivação acima.

Anexe sua resposta com a etiqueta de spoiler >!.


1
o que me perturba é que é um sinal de energia finito , não um sinal de energia finita , que é o que obtemos quando somamos esses dois sinais de energia infinitos. f[n]
22818 Robert De Bristow-johnson

Além disso, ?
DTFT{x[n]=1}=2πk=+δ(ω2kπ)
precisa saber é o seguinte

Obrigado rapazes por suas respostas! Eu votei todos eles, e cada um resulta em uma boa discussão sobre aspectos não tão conhecidos da DTFT de sinais estranhos (ou seja, aqueles que não estão em ou ). Só posso aceitar uma e esperarei um pouco mais por novas respostas ou alterações nas respostas existentes. Também adicionarei minha própria resposta mais tarde. 212
Matt L.

1
Matt, decididamente não é energia finita. um número infinito de amostras cujo quadrado é não é adicionado como um número finito. 1f[n]14
22818 Robert Bristow-Johnson

1
@ robertbristow-johnson: O que você acha perturbador nisso? Se os sinais se cancelam em todos os lugares, exceto por um número finito de pontos, é isso que obtemos.
Matt L.

Respostas:


6

Existem infinitos sinais que mantêm a seguinte igualdade: A única coisa que importa é que e, em seguida, o restante dos coeficientes de pode ser determinado sob a restrição de que a Eq. estados (isto é, a subtração de amostras consecutivas deve ser 0 para n 0 ). Em outras palavras, Eq. ( 1 ) será alcançado por qualquer sinal y [ n ] tal que y [ 0 ] = y [ - 1 ]

y[n]y[n1]=δ[n](1)
y[0]y[1]=1y(1)0n0(1)y[n]
y[0]=y[1]+1y[n]=y[n1] n0
Outra maneira de ver isso é que qualquer função que é basicamenteu[n] com um deslocamento (um valor constante adicionado) será satisfatória(1) . Isso explica a afirmação derobert bristow-johnson em sua resposta: os diferenciadores destroem essas informações (como obter uma derivada em tempo contínuo destrói evidências de qualquer valor constante na função original).

Para resumir, acredito que a prova é falha porque o procedimento seguido poderia usar qualquer função da forma u[n]+C com CR , e isso levaria a muitas funções com a mesma transformação de Fourier, o que está realmente errado como a transformada de Fourier é uma bijeção. Talvez o autor tenha deliberadamente decidido ignorar qualquer coisa relacionada aos valores de DC, consciente de que, para mostrar que F(ω) é o DTFT de f[n]ele precisaria da propriedade de acumulação (cuja prova mais popular deriva do DTFT da unidade stepgo, uma prova bastante circular). A prova não está estritamente errada , pois tudo o que afirma (as fórmulas para F(ω) e U(ω) , a decomposição da etapa unitária, a equação da diferença) é verdadeiro, mas exigiria que a propriedade de acumulação mostrasse por que F(ω) não tem nenhum delta do Dirac.


Você está totalmente no caminho certo! Você tem uma idéia de como essa falha pode ser resolvida, ou seja, como fazer isso corretamente?
Matt L.

@MattL. Definir uma condição inicial para resolveria o problema e determinaria o sinal univocamente. Essa condição inicial determinaria o valor DC do sinal y [ n ] , que aparece no DTFT como uma constante multiplicadora de um impulso Dirac (de acordo com a propriedade de acumulação). Eu acho que na prova fornecida, isso funciona porque o sinal f [ n ] não tem valor DC, pois é simétrico em torno de 0 , e, portanto, o DTFT está correto nesse caso. Mas o fato de o sinal não ter CC deve ser declarado, como é fundamental, acredito. y[n]y[n]f[n]0
Tendero

Há muitas boas respostas e é difícil escolher qual delas aceitar. Mas este foi o mais apreciado pela comunidade, e também acho que indica mais claramente o erro na derivação. Obrigado a todos!
Matt L.

4

Fiquei impressionado com o número de respostas que recebi (10 respostas até agora!). Claro, todos eles receberam meu voto positivo. Foi divertido, obrigado pessoal por seus pensamentos, comentários, etc. Eu sei que agora a maioria de vocês já sabe qual é a falha, pelo menos a que eu quis dizer. As pessoas expressam as coisas de maneira diferente, e sempre há espaço para mal-entendidos, então tentarei formular claramente o que considero a falha mais importante nessa derivação. Estou ciente do fato de que nem todos concordam e tudo bem. Estou feliz por poder discutir esse tipo de tópico esotérico de DSP com mentes afiadas como vocês! Aqui vamos nós.

Minha primeira alegação é que todas as equações da minha pergunta estão corretas. No entanto, a derivação e a motivação de algumas delas são totalmente erradas e enganosas, e essa "derivação" só pode existir porque o autor sabia como deveria ser o resultado.

Eq. (3) na questão ( ) está correto para a sequência dada f [ n ] (Eq. ( 2 ) na questão), mas é claramente também correto para todas as seqüências da forma f [ n ] = u [ n ] + c com alguma constante arbitrária c . Portanto, de acordo com a derivação, o DTFT F resultante (f[n]f[n1]=δ[n]f[n](2)

(1)f[n]=u[n]+c
c deve ser o DTFT de todas as seqüências da forma ( 1 ) , independentemente do valor da constante c . Naturalmente, isso não faz sentido, porque o DTFT é único. Especificamente, usando essa mesma "prova", eu poderia "mostrar" que F ( ω ), conforme indicado na Eq. ( 5 ) da minha pergunta (ou Eq. ( 3 ) abaixo) é realmente o DTFT de u [ n ] que estamos procurando. Então, por que se preocupar em dividir u [ n ] como na Eq. ( 1 ) da pergunta?F(ω)(1)cF(ω)(5)(3)u[n]u[n](1)

No entanto, é verdade que os DTFTs de todas as seqüências satisfazem a Eq. ( 4 ) na pergunta (repetida aqui por conveniência): F ( ω ) ( 1 - e - j ω ) = 1 Mas agora vem a falha matemática real: da Eq. ( 2 ) é incorreto concluir F ( ω ) = 1(1)(4)

(2)F(ω)(1ejω)=1
(2) Eq. (3)é apenas uma das infinitas soluções possíveis de(2)e, convenientemente, é a que o autor precisa para chegar ao resultado final correto. Eq. (3)é o DTFT def[n]em(1)comc=-1
(3)F(ω)=11ejω
(3)(2)(3)f[n](1) , mas a partir da derivação dada, não há como saber isso.c=12

(2)all(1)c(2)

(4)F(ω)=11ejω+αδ(ω)
α(4)(2)
1+α(1ejω)δ(ω)=1+α(1ejω)|ω=0δ(ω)=1+0δ(ω)=1
F(ω)(4)(2)

α(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ω)

Não tenho muita certeza de que não possa haver outra função além do impulso Dirac delta (e seus derivados) que possua essa propriedade. Mas está tudo bem, sua resposta está bem escrita. Eu voto positivo. Obrigado.
AlexTP2

2

A falha segue a palavra "Obviamente", se essa deveria ser a função Dirac Delta.

Aqui está o rascunho de uma resposta para sua outra pergunta que eu nunca publiquei:

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

Eu não acho que uma prova seja possível. Pode ser um caso de uma "definição funcional" com as propriedades desejadas.

X2π(ω)=n=+x[n]ejωn
U=n=0+ejωn
U = lim N [ 1 - e - j ω N
U=limNn=0N1ejωn
U=limN[1ejωN1ejω]
U=11ejωlimN[ejωN1ejω]
ω=0 , é claro que age como um Dirac Delta. Por que o coeficiente deve ser , não sei. Pode ter a ver com a área do círculo unitário. Quando , o denominador pode ser fatorado fora do limite e o numerador apenas salta ao longo do círculo da unidade e nunca atinge um limite. Configurá-lo para zero é um ato de definição. πω0

Provar que a definição funciona de maneira desejável é uma questão diferente.

A prova da página 138 está errada (pelo menos) porque:

δ(t)=lima012a[u(t+a)u(ta)]=dudt
O que não é semelhante de maneira alguma a como eles o definem. δ(n)=u2(n)u2(n1)

Situação interessante, espero que isso ajude. Estou ansioso pelo que você tem a dizer.

Ced


Obrigado pela sua resposta! Mas observe que estamos falando de seqüências de tempo discreto aqui, não do impulso delta do Dirac. Em tempo discreto, é igual a zero em todos os lugares, exceto onde é . Portanto, nada de feio como em tempo contínuo. Então Eq. (3) na minha pergunta é realmente correto (em tempo discreto!). δ[n]n=01
Matt L.

Cedron, acho que esta equação: pode ser interessante nesta outra questão do @MattL. Talvez você deva pensar em pensar um pouco mais e postá-lo lá, se quiser.
U=11ejωlimN[ejωN1ejω]
Tendero

2

se você me permitir dividir por zero, posso provar que . quando você diz há um problema em multiplicar algo por zero (quando ) e esperando que o produto seja igual a um.1=2

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

3
moral da história: os diferenciadores destroem a informação. um diferenciador não sabe a diferença entre e . não tente dividir por zero para obter essas informações perdidas. u[n]u[n]12
22818 Robert De Bristow-johnson

então, o autor deve declarar que para corrigir a falha? w2πk
precisa saber é o seguinte

bem, o DTFT de é igual a 1 para qualquer , mesmo quando . δ[n] ωω=2kπ
22418 Robert Bristow-Johnson

por que ? ... é o domínio de que restringe o de ? δ[n]1ejwF(w)
precisa saber é o seguinte

porque quando Matt diz e diz para integrar ambos os lados para obter (que é o que Matt diz essencialmente na Eq.5), então temos um sorta "multiplique e divida pelo mesmo truque de quantidade para obter 1" . mas às vezes ele está multiplicando e dividindo por zero. multiplicar por zero destrói as informações. dividir isso por zero não recuperará as informações.
f[n]f[n1]=δ[n]
u[n]
22810 Robert Pattinson-Johnson

2

A Equação (4) deve ser escrita como Para , que não é (5). Não sei como consertar a prova sem evitar (3).

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(ω)

@AlexTP Por que a divergência desse motivo de soma para dizer que não existe? Muitos sinais são tais que seus "DTFTs de soma" divergem, mas eles têm um DTFT definido no entanto (em termos de deltas do Dirac). Não estou dizendo que o raciocínio está errado, estou apenas tentando entender (não estou muito familiarizado com a teoria da distribuição e isso). F(ω)
Tendero

Você está certo de que as somas não convergem no sentido convencional, mas, no entanto, as DTFTs existem no sentido distributivo. Afinal, o mesmo argumento seria válido para o DTFT de , que também existe nesse sentido. Portanto, há outra falha que é menos matematicamente complexa do que o problema da divergência, com o qual nos acostumamos. u[n]
Matt L.

1
@MattL. Eu não sei se é óbvio, mas eu não ver a diferença entre e porque também temos ! u[n]f[n]u[n]u[n1]=δ[n]
AlexTP18

2

Acho que descobri a melhor maneira de expressar a falha nessa prova. Então, eu vou dar outra facada.

A escolha de em (1) é arbitrária. Vamos substituí-lo por . Siga a prova e termine com: 12x

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

Não há nada na prova que restrinja a ser , ele pode assumir qualquer valor finito e a prova funciona da mesma forma. x12

Além disso, se você der o passo que dei na minha última resposta e encontrar (4), é expresso como

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

Em seguida, incluindo-o em (5) e (6), você obtém:

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

O que, como indiquei anteriormente, é inconsistente com a definição para chegar lá.

Essa prova não mostra que e parece indicar que, para qualquer x definido, resultados inconsistentes se seguirão. Portanto, volto à afirmação em minha primeira resposta de que o valor de para o coeficiente de é um ato de definição, não um truísmo matemático. x=12πδ(ω)

Talvez haja alguma outra situação que faça o valor correto, mas essa prova não o fornece. x=12

Ced


1

Isso é uma resposta aos comentários da minha primeira resposta. Por causa da camuflagem do spoiler, estou postando como uma resposta separada.

Eu ia postar minha outra resposta para a outra pergunta, mas não o fiz devido à minha falta de experiência nessa área. Eu o publiquei ontem, excluí-o, depois cancelei a exclusão e descobri como empregar tags de spoiler.

Claramente, a função definida no problema não é a função Dirac Delta. Procurei o DTFT na Wikipedia e o DTFT para a função Dirac Delta é um deles. Vou chamar o do problema . δδδp

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

Tomando o DTFT das partes esquerda e direita. Não sei se tenho as anotações corretas, mas a matemática deve ser clara. Usando a definição que está sendo provada.

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

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

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

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

Portanto, o RHS de (4) está incorreto, exceto quando . [Edit: Doh, é o Dirac Delta, então esta afirmação está errada. Eu acho que deve estar correto, exceto "indefinido" em . Análise real era a minha matemática menos favorita. Estou deixando isso em paz agora.] ω=2kπω=2kπ

Ced

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

Acompanhamento:

É claro que o DTFT de deve ser 1 quando conectado à definição de um DTFT. Portanto, como recebi uma resposta diferente ao usar a definição a ser provada, significa que a definição a ser provada não está correta (no sentido matemático). Além disso, se você levar a correção até o final da prova, chegará a uma definição diferente. Assumir que a afirmação é verdadeira é usada para provar que é falsa.δp


Você realmente mostrou que porque o termo é igual a zero. Esse é o caso, porque para qualquer função que é contínua em você tem e se (que é o caso aqui), o termo inteiro desaparece. Fp(ω)=1(1ejω)δ(ω)f(ω)ω=0f(ω)δ(ω)=f(0)δ(ω)f(0)=0
Matt L.

1

Para mim, uma primeira falha aparece entre (3) e (4): este é um exemplo da divisão integral descuidada / soma infinita clássica. São necessárias condições para permitir a equação: padrão ou podem não ser nítidas o suficiente. Isso pode estar relacionado aqui, devido à forma , ao teorema da derivação de Fubini , ou: Quando podemos trocar soma infinita e derivada discreta? As formas de locomoção podem girar em torno da monotonia ou de somas do tipo Cesaro, mas pensarei nisso por mais tempo.

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

1

então Matt,

Não sei por que você não acha problemático comparar sinais de energia com sinais de energia, mas suponha que modifiquemos ligeiramente a definição de :f[n]

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

para alguns .α>0

agora temos sinais de energia finita e todos os DTFTs devem ser comparáveis.

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

Eu me pergunto o que são os DTFTs? e então o que acontece quando deixamos ? acho que ainda há o problema de diferenciadores destruindo informações (e a correspondente destruição de informações multiplicando por 0 no domínio da frequência) que é um problema. mas talvez possamos perder o problema de comparar classes de sinal que não compartilham o mesmo espaço de Hilbert.α0

mas, infelizmente, são quase duas da manhã e eu não vou lidar com isso agora.


Essa coisa é boa e é uma opção para calcular o DTFT desses sinais não deteriorantes, tomando o limite . Experimente e tenho certeza que você terá sucesso, mas é doloroso. Existem maneiras mais fáceis de obter o mesmo resultado. A prova fornecida pode realmente ser modificada para que funcione no IMHO (veja minha resposta). α 0αα0
Matt L.
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.