Amostragem de uma função contínua: o delta de Kronecker ou Dirac?

12

Eu tenho lido alguns artigos em processamento de sinais e estou muito confuso sobre o assunto no título da minha pergunta. Considere-se uma função contínua de tempo , f ( t ) , que amostra a vezes irregulares t k , onde k = 1 , 2 , . . . , N . Para mim, faz sentido que a função amostrada seja: f s ( t ) = N k = 1 δ t , t k f ( t )tf(t)tkk=1,2,...,N onde δ t , t k éodeltade Kronecker(igual a 1 quando t = t k , zero em outro lugar). Entretanto,neste artigo, o autor define o sinal amostrado como: f s ( t ) = 1

fs(t)=k=1Nδt,tkf(t),               (1)
δt,tk1t=tk ondeδ(t-tk)é a função delta de Dirac e eu realmente não entendo por que o1/Naparece aqui (o autor afirma que a função de amostragem é na verdade uma soma ponderada das funções delta s(t)=C N k = 1 wkδ(
fs(t)=1Nk=1Nf(t)δ(ttk),   (2)
δ(ttk)1/N e aqui ele escolheC=wk=1. Eu realmente não entendi o porquê). Esta última afirmação não faz muito sentido para mim: o sinal amostrado teria amplitude infinita emt=tk!
s(t)=Ck=1Nwkδ(ttk)k=1Nwk,
C=wk=1t=tk

fs(t)(2)f(t)(1)f(t)

Néstor
fonte

Respostas:

9

Modelar o processo de amostragem via multiplicação de um sinal de tempo contínuo por um trem de impulsos de Dirac é a interpretação mais comum em minha experiência. Se você se aprofundar o suficiente, encontrará algumas divergências sobre a precisão matemática dessa abordagem *, mas eu não me preocuparia; é apenas um modelo conveniente para o processo. Não há geradores de impulso dentro do ADC do seu telefone celular, gerando raios periódicos que multiplicam suas entradas analógicas.

Como você observou, não é possível calcular a transformação Fourier em tempo contínuo da função delta Kronecker, pois seu domínio não é contínuo (está limitado aos números inteiros). A função delta Dirac, ao contrário, possui uma simples transformada de Fourier, e é fácil mostrar o efeito de multiplicar um sinal por um trem de impulsos de Dirac devido à sua propriedade de peneirar.

*: Por exemplo, se você for matematicamente preciso, diria que o delta do Dirac não é uma função, mas uma distribuição . Mas no nível da engenharia, essas questões são realmente apenas semânticas.

Edit: Vou abordar o comentário abaixo. Você forneceu seu modelo mental do processo de amostragem como:

fs(t)=k=1Ntk-ϵktk+ϵkf(t)δ(t-tk)dt.

fs(t)tϵk>0 0

fs(t)=k=1Nf(tk),

o que não está correto. Em vez disso, o modelo para o sinal amostrado é:

fs(t)=k=-f(t)δ(t-kT)

O que é muito semelhante ao anterior, exceto a generalização de um trem de impulso infinitamente longo ao longo do eixo do tempo e a suposição de que os dados são amostrados uniformemente em instantes de tempo tk=kT. A transformada de Fourier do sinal resultante é:

Fs(ω)=fs(t)ejωtdt=k=f(t)δ(tkT)ejωtdt=k=f(t)δ(tkT)ejωtdt=k=f(kT)ejωkT

If we define the discrete, sampled version of the signal f(t) to be x[n]=f(nT), then you are left with:

Fs(ω)=n=x[n]ejωn

which is exactly the definition of the discrete-time Fourier transform.

Jason R
fonte
How would you address the fact that the amplitude is "infinite"? What I have usually thought is that you actually don't "sample" the signal at a discrete time tk, but rather you integrate the signal for a given time Δtk. However, this interpretation would violate any form of calculating the Fourier Transform by the same reason as the Kronecker delta. Also...why does the author of the paper in the link I give divides the Dirac comb by N? That doesn't make any sense to me.
Néstor
1
In practice, you're right. There is always some effective "integration time" of the analog signal owing to the finite bandwidth of any ADC's analog front end. However, the theoretical construct isn't limited by such concerns. Roughly speaking, the "infinite height" of the impulse is balanced by its "zero width", such that it integrates to unity. If you apply your short-time integration interpretation in this case (multiplying by an impulse, integrating for an infinitestimally-short time period), then by the sifting property, you would get x[n]=x(nT), as is usually presented.
Jason R
Yeah, but my concern is not with the interpretation, but rather with taking the fourier transform of it. Suppose I write the sampling process that we are talking about as:
fs(t)=k=1Ntkϵktk+ϵkf(t)δ(ttk)dt,
how would you take the fourier transform of that? I know the trick when ϵk0, but that doesn't make much sense to me (and it's even harder to do the FT!). Even supposing that is the way, I don't get the same window function as in the paper of Roberts et al. that I cited. And I insist...that 1/N doesn't make any sense to me.
Néstor
Ok with your edit and with my error in my comment above. However, I still can't make up my mind about the fact that the impulse has infinite amplitude at t=tk thanks to Dirac's Delta, i.e., f(t=tk), while what we actually observe (and want to model) is f(t=tk)=f(tk)
Néstor