Amplitudenspektrum so richtig?

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Illmatic

Forum-Anfänger


	Beiträge: 17

	Anmeldedatum: 19.06.12

	Wohnort: Hamburg

	Version: R2011b

Verfasst am: 28.06.2013, 13:40 Titel: Amplitudenspektrum so richtig?

Hi,
kann jemand mit mehr Ahnung als ich mir kurz sagen ob meine Funktion zum Erzeugen eines Amplitudenspektrums so OK ist?

Code:

function [freq,X,N] = simplefft(measured_signal,sampling_rate)
L=length(measured_signal);
window = hann(length(measured_signal));
measured_signal=measured_signal.*window;
N=2^nextpow2(L);
X = fft(measured_signal,N)/L;
freq = sampling_rate/2*linspace(0,1,N/2+1);
cutOff = ceil(N/2);
X = abs(X(1:cutOff));
freq = freq(1:cutOff);
end

Funktion ohne Link?

Vielen Dank!

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 13.08.2013, 13:14 Titel:

Ich würde sagen, dass das eigentlich passt. Hätte nur die Frequenz so erzeugt:

Code:

freq = linspace(0,sampling_rate/2,N/2);

Funktion ohne Link?

macht aber ja keinen großen unterschied Wink

DSP

Forum-Meister


	Beiträge: 2.117

	Anmeldedatum: 28.02.11

	Wohnort: ---

	Version: R2014b

Verfasst am: 13.08.2013, 16:48 Titel:

Die Skalierung stimmt bei folgenden Teilen nicht

- window
- X darf nicht nur durch L sondern muss durch N geteilt werden, wenn N>L
- nur der Gleichsignalanteil, also X(1) darf durch L geteilt werden...alle anderen nur durch L/2 (wenn du die Nyquistfrequenz weglässt).

Siehe hier...

http://www.gomatlab.de/fft-umfassendes-beispiel-t777.html

Hier noch ein weiteres Skript zu dem Thema:

Code:

function [mag, mag_dB, fv] = FFT_betragsspektrum( signal, nfft, fa, scale) ;
% Input:
% Signal im Zeitbereich
% nfft = Anzahl Messwerte für fft
% wenn nfft > length(sig) -> fft(sig,nfft) führt Zeropadding durch
% fa = Abtastfreq.
% scale: 0 = keine Impulsantwort als Eingang, 1 = Impulsantwort
% Output:
% Magnitude des Spektrums linear und dB skaliert
% Frequenzvektor fv in [Hz] von 0...fa/2

% un-,gerade Anzahl Messwerte?
if mod(nfft,2) == 0;
k = (nfft/2) + 1;
else
nfft = nfft + 1;
k = (nfft/2) + 1;
end

fn = fa/2; % Nyquistfreq.
df = fa/nfft; % Frequenzauflösung des Spektrums
% Frequenzvektor: Darstellung bis Nyquistfreq.
fv = 0: df : fn;

sig = signal(:);

% Signal transformieren
Fy = fft(sig,nfft);

% Betrag - nur positives Freq.spektrum
if scale == 0
Fy_pos0 = abs(Fy(1:k));

% Skalierung
mag = [Fy_pos0(1)/nfft ;Fy_pos0(2:k-1)/(nfft/2);Fy_pos0(k)/nfft];
else % nicht durch nfft teilen bei Impulsantwort
mag = abs(Fy(1:k));
end

% Umrechnung in dB
mag_dB = 20*log10(mag + eps); % eps = kleine Konstante zur Vermeidung von log(0)

Funktion ohne Link?


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.