Matlab FFT / IFFT

Gast #646045
Hallo!
Ich spiele gerade etwas mit Matlab rum und habe mir eine Zeitfunktion 
bestehend aus drei Sinunsfunktionen unterschiedlicher Frequenz 
generiert. Diese Wandele ich ins Spektrum und wieder in den Zeitbereich 
und wieder ins Spektrum. Nun wundere ich mich, warum die Amplitudenhöhen 
nicht übereinstimmen! ...richtig, ich wandele ins 
Leistungsdichtespektrum (also keine Phaseninfo)! Also ist mein 
Zeitsignal aus der Rücktransformation natürlich Phasenfalsch - was aber 
doch auf mein daraus entstehendes Spektrum keinen Einfluss hat.
Hier skizziere ich mal den Code:

%Wahl des Testsignals (oben auch ändern)
ZeitSig=10*sin(2*pi*Frequenz*t)+5*sin(2*pi*Frequenz2*t)+8*sin(2*pi*Frequ 
enz3*t);

%FFT vom Signal x
Y= fft(ZeitSig,8192);
% Leistungsdichte Pyy
Pyy=Y.*conj(Y)/8192;

ZeitSigRueck = ifft(Y,8192);
ZeitSigAusPyy = real(ifft(Pyy,8192));

%FFT vom Signal ZeitSigAusPyy
Y= fft(ZeitSigAusPyy,8192);
Pyyy=Y.*conj(Y)/8192;

Ergebnis: Pyy ungleich Pyyy!!! Warum?? ...1000 Dank!
Gast #646097
Die Phase schmeiße ich weg (oder besser: In meiner späteren Anwendung 
habe ich sie gar nicht, da ein Spektrumanalyser sie nicht messen kann) - 
möchte quasi nur ein Zeitsignal, was ebenfalls das selbe Spektrum hat 
wie das Spektrum der original-Zeitfunktion!
Bei einem Zeitsignal, welches aus mehreren Sinusfunktionen gleicher 
Amplitude besteht klappt das! ...doch wehe sie haben unterschiedliche 
Amplituden! ...wie kommt das? Danke für weitere Tipps
Gast #744139
Matthias hat Recht. Mit dem Imaginärteil schmeißt man NICHT (NUR) die 
Phase weg! Das sind komplexe Werte, und keine Polarkoordinaten! Also 
statt real einfach abs benutzen, und natürlich sqrt nicht vergessen.

Antwort schreiben

Eine eigene Frage? Dieser Thread ist über 19 Jahre alt. Stell sie in einem eigenen Thread – dort lesen mehr Leute mit, und sie wird öfter und schneller beantwortet.

Neuen Thread starten

Eine Lösung, Ergänzung oder Rückmeldung zu diesem Thread? Dann antworte hier.

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.