Confusion with FFT for a specific part of the whole signal

7 次查看（过去 30 天）

显示更早的评论

Ancalagon8 2017-2-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/323365-confusion-with-fft-for-a-specific-part-of-the-whole-signal

Locked: Rena Berman 2025-1-13

采纳的回答： Star Strider

As a newbie i am working on a signal from an accelerometer and i am trying to produce an fft for a specific part of the whole signal. I am a little bit confused about which code for the fft is the appropriate to use for my case.

here is my code:

clear all

close all

load signal

tvec=noact(:,1);

acc=noact(:,2);

%% ALL SIGNAL

L=length(signal); % Length of signal

t_max=max(signal(:,1)); % Tmax

Fs=L/t_max; % Sampling Rate

T = 1/Fs; % Sample time

%% PART OF SIGNAL

trace1=signal(Fs*25:Fs*32,:); % select part of signal

matrice1=trace1(:,2); % accelerations column

L1=length(matrice1); % Length of trace (number of samples)

t1 = (0:L1-1)*T; % Time vector

% Frequency vector

f1=(0:Fs/length(matrice1):(length(matrice1)-1)*Fs/length(matrice1));

FrecSign=abs(fft(matrice1))*T;

figure

plot(FrecSign);

xlabel Frequency(Hz)

ylabel Amplitude

title 'FFT'

NFFT = 2^nextpow2(L); % Next power of 2 from length of matrice1

Y = fft(matrice1,NFFT)/L;

f = Fs/2*linspace(0,1,NFFT/2+1); % create frequency vector

figure % 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)|')

Thanks in advance

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rena Berman 2025-1-13

(Answers Dev) Restored edit

This question is locked.

采纳的回答

Star Strider 2017-2-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/323365-confusion-with-fft-for-a-specific-part-of-the-whole-signal#answer_253338

I cannot follow our code.

However, you cannot ‘produce an fft for a specific part of the whole signal’, since you have to take the fft (link) of the entire signal if you want valid data.

If you want to view only part of the signal, use the axis function to limit the frequency (here x) axis to the range you want.

If you want to eliminate specific frequencies from your accelerometer signals (for example a constant offset, low-frequency baseline variations, high-frequency noise), use a discrete filter. There are several ways to to this, with the designfilt function being the easiest to use (in my opinion).