Multiplikation von zirkulierenden Matrizen mit einer Diagonalmatrix


8

Sei , B i eine Folge von zirkulierenden Matrizen der Größe n × n .AiBin×n

Wir wissen, dass in quadratischer Zeit berechnet werden kann (verwenden Sie FFT, um die Diagonalmatrizen zu diagonalisieren und zu addieren und IFFT anzuwenden).i=1nAiBi

Angenommen, ist eine beliebige diagonale Matrix (der Einfachheit halber sei r die n- te Wurzel der Einheit und betrachte die diagonalen Elemente als alle unterschiedlichen Potenzen kleiner als n von r ).Drnnr

Was ist die Komplexität von ? Ich vermute, dass es quadratisch ist, da ich in jedem Term die gleiche Diagonalmatrix ( O ( n ) Terme) einbinde.i=1nAiDBiO(n)

Betrachten Sie eine zirkulierende Matrix der Größe n × n mit einer ersten Reihe aus unterschiedlichen Potenzen von weniger als n von r . Sei X i und Y i für i = 1 n Diagonalmatrizen mit vollem Rang.Rn×nnrXiYii=1n

Was ist die Komplexität von ? Wieder vermute ich, dass dies quadratisch ist.i=1nXiRYi

Die in Bezug auf r definierten Matrizen und R sind künstlich. Ich suche für den Fall der allgemeinen Diagonale D und allgemeinen Voll Rang zirkulante R .DRrD R

Antworten:


3

i=1nXiRYiXiYiAXiiBYiiABR. Wenn alle Einträge vonRungleich Null sind (auch wenn es nicht der volle Rang sein sollte), kannABaus n i = 1 XiRYirekonstruiert werden. Daher ist die rechnerische Komplexität dieses Problems identisch mit der rechnerischen Komplexität der Matrixmultiplikation.i=1nXiRYiRABi=1nXiRYi

Für den Fall , wobei A i und B i zirkulierende Matrizen sind, müssen wir zuerst FFT auf A i , B i und D anwenden , um zum Fall n i = 1 zurückzukehren X i R Y i . Dazu benötigen wir 2 n + 1 FFTs der Größe n, wodurch die Komplexität um n 2 log n erhöht wird.i=1nAiDBiAiBiAiBiDi=1nXiRYi2n+1nn2logn
DDrirnDD

Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.