How do I take the fft of a signal after I using the blackman window?
17 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
SCuri
el 6 de Mzo. de 2019
Respondida: Enthusiastic Student
el 8 de Mzo. de 2019
I have written the following code for putting the blackman window on the data :
window = blackman(length(strain_seg)); %the length of the strain_seg is 256^2
windowed_strain = strain_seg*window;
How do I perform an fft of the windowed strain? Do I need to use the rfft function?
Respuesta aceptada
Enthusiastic Student
el 8 de Mzo. de 2019
I think you might have issue with dimensions early on.
Have you tried to see what
size(strain_seg)
length(strain_seg)
actually output? I'm guessing the strain_seg has a size of [65536 1] and the blackman has a size of [1 64], which is why you have to add zeroes. I would suggest transposing one of them using the ' operator:
M = length(strain_seg);
w = blackman(M);
xw = C.*w';
X = fft(xw)
With regards to plotting the fft of the windowed strain segment you would need to know the time the segment spans. Lets say the segment you are examining was taken over 10 ms:
Ts = 10e-3; %samples taken in a 10 ms long period of time /[s]
dt = 1/Ts; %/[Hz]
freq = dt:dt:(dt*length(xw)+dt); %Frequencies
figure
loglog(freq,abs(X));
xlabel('Frequency /[Hz]')
ylabel('Magnitude of FFT')
0 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Blackman en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!