How do I take the fft of a signal after I using the blackman window?
17 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
SCuri
il 6 Mar 2019
Risposto: Enthusiastic Student
il 8 Mar 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?
Risposta accettata
Enthusiastic Student
il 8 Mar 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 Commenti
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Blackman in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!