Faltungssatz

Gast #4165987
Lesenswert?
• ▲
▼
Ich habe eine Verständnisfrage zum Faltungstheorem. Es liegen z.B. zwei 
beliebige Folgen der selben Länge vor: a = [1 2 3 4 3 2 1] und b = [7 6 
5 4 3 2 5]. Aufgabe ist die Faltung dieser. Es ergibt sich mittels 
Matlab conv(a,b,'same') c = [60 71 72 68 56 45 36]. Jetzt das Ganze 
mittels Faltungstheorem: A = fft(a), B = fft(b), c = ifft(A.*B) = [63 65 
74 82 83 77 68]. Das ist nicht das Selbe. Vielleicht kann mich einer 
Aufklären, an welcher Stelle ich das Theorem falsch verstanden habe, und 
wie das korrekt im Beispiel funktioniert. Danke euch.
Gast #4166079
Lesenswert?
• ▲
▼
Hallo,

ich muss gestehen, dass ich keine Ahnung habe, was conv(a,b,'same') 
berechnet. Aber A = fft(a), B = fft(b), c = ifft(A.*B) ist definitiv die 
zyklische Faltung.

Die Frage ist nun, wolltest du die zyklsiche Faltung berechnen?
Oder willst du eine normale Faltung, in dem Fall muss bei der fft 
Methode noch Nullen angefügt werden.

PS. conv(a,b) gibt die standard Faltung an.

Mein Vorschlag wäre, das ganze einfach mal per Hand zu rechnen und dann 
einen kleinen Algorithmus für die zyklische Faltung zu schreiben.

NotValid
Gast #4166452
Lesenswert?
• ▲
▼
Moin,

You get what you pay for - Will sagen: Wenn du Fouriertransformierst, 
dann kriegst du was raus, was so aussieht, als wuerden sich deine 
Eingangsdaten andauernd wiederholen. Auch wenn du das nirgends explizit 
angibst.

Und wenn du 'same' als Option bei der Faltung angibst, dann bedeutet 
das, dass dich vom Ergebnis der Faltung nur irgendwas aus der Mitte 
interessiert, keinesfalls jedoch das komplette Ergebnis...

guggstu:
1
a = [1 2 3 4 3 2 1 0 0 0 0 0 0 0];
2
b = [7 6 5 4 3 2 5 0 0 0 0 0 0 0];
3
A = fft(a)
4
B = fft(b)
5
c = ifft(A.*B)
6
conv(a,b)

Da sieht sich's doch schon aehnlicher find' ich...

Gruss
WK

Antwort schreiben

Hinweis: der ursprüngliche Beitrag ist mehr als 6 Monate alt. Bitte hier nur auf die ursprüngliche Frage antworten, für neue Fragen einen neuen Beitrag erstellen.

Formatierung
  • [c]C-Code[/c]
  • [code]Code in anderen Sprachen, ASCII-Zeichnungen[/code]
  • [math]Formel in LaTeX-Syntax[/math]
  • [[Titel]] – Link zu einem Artikel
  • Verweis auf einen anderen Beitrag: Rechtsklick auf dessen Beitrags-ID (#8078628), „Adresse kopieren“, in den Text einfügen

Mehr zur Formatierung

Längeren Sourcecode als Anhang, nicht in den Text. JPG nur für Fotos und Scans; Zeichnungen und Screenshots als PNG oder GIF. Siehe Bildformate.

Schon registriert? Anmelden

Wird nicht angezeigt. Wir schicken dir einen Link; sobald du damit bestätigst, erscheint dein Beitrag.

Mit dem Abschicken legst du einen kostenlosen Account an und erkennst die Nutzungsbedingungen und die Datenschutzerklärung an.