Auto Power Spektrum

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

bet05031

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 02.12.08

	Wohnort: Ingolstadt

	Version: R2008a

Verfasst am: 15.01.2009, 10:18 Titel: Auto Power Spektrum

Hallo zusammen,

ich bin seit Tagen am rumversuchen ein Autopowerspektrum eines Schwingsignals zu berechnen. Aber ich bekomm das einfach nicht gebacken.

Der Hintergrund ist das ich in einer Aufnahmesofware die Daten verarbeiten und ausgeben kann und nu versuche das identische Signal in MatLab zu berechnen! Ich bekomme also die Zeitrohdaten (Beschleunigung über der Zeit) in MatLab rein und muss über FFT und APS und zweifache Integration ein Schwingweg über der Frequenz bekommen.

Falls jemand sowas schonmal gemacht hat oder sonst eine idee hat wäre ich sehr dankbar!

Mfg

Pat

Maddy

Ehrenmitglied


	Beiträge: 494

	Anmeldedatum: 02.10.08

	Wohnort: Greifswald

	Version: ---

Verfasst am: 15.01.2009, 10:46 Titel:

Wie sehen denn deine bisherigen Versuche aus?

Fragen die mir in Bezug auf das Problem einfallen:
Hast du die FFT normiert und reskaliert bzw. ist dies notwendig?
_________________

>> why
The computer did it.

bet05031

Themenstarter

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 02.12.08

	Wohnort: Ingolstadt

	Version: R2008a

Verfasst am: 15.01.2009, 10:57 Titel:

Ich hab die FFT nach einem Skriptbeispiel aus dem Forum berechnt!
Ich poste mal! Hoffe es ist einigermaße übersichtlich:

Code:

% Zeitbereich
% ----------------------------------

fa = 8000; % Abtastfrequenz
fn = fa/2; % Nyquistfrequenz
N = 8000; % gewünschte FFT-Länge (N=2^x, sonst wird der DFT-Algorithmus verwendet!)
df = fa/N; % Frequenzauflösung

% Fensterfunktion
% ----------------------

win = hann(N)';
%y_win = y.*win; % Fensterung ohne Amplitudenkorrektur
y_win = y.*win*N/sum(win); % Fensterung mit Amplitudenkorrektur
max_y = max(abs(y_win))*1.1;

fig = figure(fig+1);
plot(y_win)
axis([0 N -max_y max_y])
title('Datensatz nach Fensterung mit Hann-Fenster')
ylabel('Amplitude')
xlabel('N Stützstellen')
grid

% Berechnung der FFT
% ------------------
H = fft(y_win, N);
% Berechnung des Amplitudengangs aus dem komplexen Frequenzvektor H:
amplH = abs(H);
% Amplitudenskalierung (Normierung auf N) und verschieben der Elemente des
% Amplitudenvektors, so dass die Darstellung des Amplitudengangs von -fn...0...fn
% erfolgen kann:
amplitudengang = fftshift(amplH/N);

Funktion ohne Link?

Die Daten kann ich wegen Geheimhaltung leider nicht mit posten! Sorry


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.