Filter löschen
Filter löschen

Power spectrum ploting matlab

1 Ansicht (letzte 30 Tage)
aditi
aditi am 12 Jul. 2013
I have to draw a plot between frequency(x axis) and power(y axis).
All I have is a file with first colum time values, second column magnitude

Antworten (1)

David Sanchez
David Sanchez am 12 Jul. 2013
You need to perform a Fast Fourier Trasnsform (fft). The following code, extracted from matlab doc:
doc fft
will help you. Run it, and then, adapt it to your needs:
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid
x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
y = x + 2*randn(size(t)); % Sinusoids plus noise
plot(Fs*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
plot(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
  1 Kommentar
aditi
aditi am 12 Jul. 2013
sorry...but this is not giving the correct result.. :(

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Fourier Analysis and Filtering finden Sie in Help Center und File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by