FFT-Analyse für Oberschwingungsgehalt

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

ElBarto

Forum-Anfänger


	Beiträge: 27

	Anmeldedatum: 27.09.10

	Wohnort: Bamberg

	Version: R2009b

Verfasst am: 20.12.2010, 15:24 Titel: FFT-Analyse für Oberschwingungsgehalt

Hallo zusammen,

habe ein wahrscheinlich total triviales Problem, welches ich aber bisher nicht lösen konnte.

Ich will verschiedene Signale mit Matlab auf ihren Oberschwingungsgehalt (Harmonische von 50 Hz) untersuchen.

Ich habe ein Signal, welches in Matrizenform <2001x4 double> vorliegt. Dabei ist die Matrix so aufgebaut: [t' u_L1' u_L2' u_L3']
Also eine Spalte mit Zeitwerten, drei Spalten mit Spannungswerten. Die Spalten sind logischerweise alle gleich lang.

Matlab führt erfolgreich die FFT durch, liefert z.B. für

Code:

t=linspace(0,0.2,1000)
sig1=100*sin(2*pi*50*t)+20*sin(2*pi*150*t)

Funktion ohne Link?

plausible Ergebnisse.

Nun möchte ich aber die Berechnung so gestalten, dass nur diejenigen Frequenzanteile angezeigt werden, die bei Vielfachen von 50Hz auftreten, und das möglichst genau.

Im Moment bekomme ich bei sehr verzerrten Signalen (siehe Anlage) einfach keine brauchbaren plots.

Zur Übersicht hier mein Programmcode:

Code:

t=linspace(0,0.2,1000);
sig1=100*sin(2*pi*50*t)+20*sin(2*pi*150*t);

figure(1)
plot(t,sig1);

L = length(sig1);
T = t(end);
Fs = length(sig1);
NFFT=L;
% NFFT=2^nextpow2(L);

Y = fft(sig1,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);

figure(2)
bar(f/10,2*abs(Y(1:NFFT/2+1)))
set(gca,'XLim',[0 50]);
set(gca,'XTick',(0:1:50));
grid

Funktion ohne Link?

Bin für jede Hilfe dankbar!

Gruß
Michael

PS: Nur zum Verständnis - Die x-Achse ist mit der Ordnung der Harmonischen beschriftet (50Hz = 1, 100Hz = 2, usw.)

untitled.png

Beschreibung:

Download

Dateiname:

untitled.png

Dateigröße:

20.91 KB

Heruntergeladen:

1100 mal

Scriptor

Forum-Century


	Beiträge: 217

	Anmeldedatum: 22.02.08

	Wohnort: ---

	Version: ---

Verfasst am: 21.12.2010, 00:57 Titel:

Du hast die Fensterbreite von 0.2 heißt das deine Spektralsamples um df = 1/0.2 voneinander entfernt liegen. df = 5;

Du hast also 1000 Samples. Das heißt wiederum dass du ein Spektrum von
f = [0 : 999]*df;

hast. (Um ein Symetrisches Spektrum zu erhalten nimm fftshift was korrekter ist. So erhieltest du ein Spektrum von f = [-500:499]*5).

Gehen wir aber vom ersterens aus, heißt du hast ein Spektrum von 0 : 4995.

Vll hilft dir das als Ansatz...

ElBarto

Themenstarter

Forum-Anfänger


	Beiträge: 27

	Anmeldedatum: 27.09.10

	Wohnort: Bamberg

	Version: R2009b

Verfasst am: 21.12.2010, 09:17 Titel:

Hallo Scriptor,

danke für deine Bemühungen, aber es besteht kein Handlungsbedarf mehr. Hatte versehentlich den Einschaltvorgang mit in meinem Signal betrachtet, der verrauscht den Frequenzganz natürlich recht ordentlich.

Betrachtet man das Signal im "eingeschwungenen" Zustand, kommt was brauchbares raus.

Trotzdem danke nochmal, dass du dich mit meinem Problem beschäftigt hast!

Gruß
Michael

ElBarto

Themenstarter

Forum-Anfänger


	Beiträge: 27

	Anmeldedatum: 27.09.10

	Wohnort: Bamberg

	Version: R2009b

Verfasst am: 22.12.2010, 10:25 Titel:

Frequenzganz? Naja wohl eher Frequenzgang...


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.