DFT e equivalência multiplicação / convolução

7

Existe uma explicação simples ou potencialmente intuitiva para, com a DFT, a multiplicação de vetores em um domínio ser equivalente à convolução circular das transformações dos vetores no outro domínio?

Como uma DFT é apenas a multiplicação por uma matriz quadrada (especial), o que acontece com essa matriz e multiplicação de matrizes que permite a dualidade acima?

hotpaw2
fonte

Respostas:

5

Como você diz corretamente, o DFT pode ser representado por uma multiplicação de matrizes, a matriz de Fourier F. Por outro lado, o DFT "transforma" uma convolução cíclica em uma multiplicação (como todas as variantes de transformada de Fourier como DFT, DTFT, FT têm uma propriedade semelhante de transformar convolução em multiplicação) e vice-versa.

Para entender isso na imagem da matriz, observe que também a convolução (circular) com uma determinada sequência pode ser representada por uma multiplicação da matriz. Mais especificamente, essa é uma matriz circulante, um tipo especial de matriz de Toeplitz.

tão y=cx com a convolução cíclica pode ser escrita como y=C(c)xcom denotando a matriz circulante formada a partir das entradas do vetor .Cc

Se "transformarmos" essa equação com a DFT (isto é, multiplicação por ), obteremosF

y^=FC(c)FHx^

com y^=Fy e x^=Fx os respectivos DFTs (observação FH representa o IDFT).

O ponto é agora que FC(c)FHé sempre uma matriz diagonal, porque todas as matrizes circulantes são diagonalizadas pela matriz de Fourier. Isso significa que os vetores próprios das matrizes circulantes são dados apenas pelas linhas da matriz de Fourier.

É claro que isso é consistente com a imagem da convolução, porque o DFT transforma a convolução em uma multiplicação em todo o elemento. Além disso, os elementos diagonais dessa matriz são apenas o DFT decou, eqivalentemente, os autovalores da matriz circulante formados a partir de c.

Andreas H.
fonte
Claro, acabei de inserir a matriz de identidade entre C(c) e x. Observe queFHF=Eu. Isto é para que na equação também a transformada de Fourier dex parece. FHF=Eué porque a matriz de Fourier é uma matriz unitária, assim como a DFT é uma transformação unitária.
Andreas H.
Lado esquerdo de x. Começar comy=C(c)x, multiplique por F Da esquerda: y^=FC(c)x depois insira Eu=FHF e pegue y^=FC(c)FHFx qual é y^=FC(c)FHx^
Andreas H.
Não, FCFH é o que resulta da derivação, FHCFestaria incorreto.
Andreas H.
2

Aliás, a DFT é a única transformação linear bijetiva que troca convolução e multiplicação por termos (até a permutação dos coeficientes, obviamente). Isso não é difícil de provar, mas não encontrei nenhuma referência sobre esse resultado antes de explicá -lo em Music Through Fourier Space , Thm. 1.11 (Springer 2016). É mais confuso no caso contínuo, porque é preciso escolher bem os espaços funcionais envolvidos.

Talvez esse recíproco também possa ser provado usando matrizes circulantes e diagonalização simultânea.

E. Amiot
fonte