FFT-Maximum von Signal mit mehreren Fenstern driftet

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

bordeur

Forum-Newbie


	Beiträge: 1

	Anmeldedatum: 17.09.18

	Wohnort: ---

	Version: ---

Verfasst am: 17.09.2018, 22:12 Titel: FFT-Maximum von Signal mit mehreren Fenstern driftet

Hallo,

der Titel klingt seltsam, aber das hier habe ich gemacht:
ähnlich einer STFT habe ich ein Signal in mehrere Fenster aufgeteilt, eine Fensterfunktion angewendet und jeweils das Maximum des Spektrums gespeichert.

Problem: wenn ich die Anzahl der Fenster verändere, so ändert sich auch, ob die Maxima eher größer, oder kleiner werden.
Bei den meisten Fensterfunktionen zeigt sich dieses Problem, nur beim Flat-top nicht. Allerdings ist hier gar kein Drift zu beobachten, aber ich vermute das Signal driftete tatsächlich. Um das zu überprüfen, habe ich einfach mal integriert und cumtrapz zeigt tatsächlich unabhängig von der Fensteranzahl einen Drift. Bild füge ich ein.

Ist das ein Problem in meinem Programm, oder habt ihr eine andere Idee? Ich bin euch schon jetzt sehr dankbar!

Code:

load date;
number=length(date)/1000; %Anzahl der Fenster (bei 1000 knapp 2 Sekunden)
Fs=60000; %Samplefrequenz

i=1;
j=1;
while i<(length(date)-number)
fenst=date(i:(i+number-1));

%Überprüfung des Drifts mittels Integration
flaeche(j)=max(cumtrapz(abs(detrend(fenst)).*hann(length(fenst))));
flaeche2(j)=max(cumtrapz(abs(detrend(fenst))));

freqspek=abs(fft((detrend(fenst)).*hann(length(fenst))));
freqspek1=abs(fft(detrend(fenst))); %zum Vergleich ohne Fensterfkt

%zur Frequenzkontrolle
[n,m]=size(freqspek);
N=0:n-1;
Hertz=N*Fs/(n);

[WWert,Index]=max(freqspek);
Mwert(j)=WWert(1); %Maximale Amplitude in Spektrum
Fwert(j)=Hertz(Index(1)); %zugehörige Frequenz

[Wert1,Index1]=max(freqspek1); %ohne Fensterfkt
Mwert1(j)=Wert1(1); %Maximale Amplitude in Spektrum
Fwert1(j)=Hertz(Index1(1)); %zugehörige Frequenz

j=j+1;
i=i+number;
end

Funktion ohne Link?

500 ohne hann.jpg

Beschreibung:

500 Fenster (4 Sek) ohne Hanning

Download

Dateiname:

500 ohne hann.jpg

Dateigröße:

24.11 KB

Heruntergeladen:

208 mal

500 mit hann.jpg

Beschreibung:

500 Fenster (4 Sek) mit Hanning

Download

Dateiname:

500 mit hann.jpg

Dateigröße:

29.38 KB

Heruntergeladen:

203 mal

1000 integriert.jpg

Beschreibung:

1000 Fenster (2 Sek) integriert

Download

Dateiname:

1000 integriert.jpg

Dateigröße:

39.74 KB

Heruntergeladen:

193 mal

1000 ohne hann.jpg

Beschreibung:

1000 Fenster (2 Sek) ohne Hanning

Download

Dateiname:

1000 ohne hann.jpg

Dateigröße:

32.41 KB

Heruntergeladen:

204 mal

1000 mit hann.jpg

Beschreibung:

1000 Fenster (2 Sek) mit Hanning

Download

Dateiname:

1000 mit hann.jpg

Dateigröße:

44.28 KB

Heruntergeladen:

215 mal


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.