Hi,
Your code is on the right track for performing an FFT on angular velocity data and analyzing the frequency spectrum of the signal. However, there are several improvements and additions you can make to enhance its accuracy and functionality. Notably, when displaying the power in decibels (dB), you need to properly scale the FFT result before converting it to dB. Additionally, identifying the noise frequency range and determining the cutoff frequency requires further steps.
Here's an example code with some improvements and additional explanations:
% Load the data
data = xlsread('ho130Gyrox.xlsx');
% FFT of data
N = length(data); % Length of data
Fs = 2000; % Sampling frequency
dt = 1/Fs; % Time per sample
Y = fft(data); % FFT
Y_abs = abs(Y/N); % Normalize FFT result to scale it properly
Y_db = 20*log10(Y_abs); % Convert amplitude spectrum to dB
% Generate frequency axis
frequencies = (0:N-1)*(Fs/N)/1000; % Frequency axis in kHz
% Plot
figure;
plot(frequencies(1:N/2), Y_db(1:N/2));
xlabel('Frequency (kHz)');
ylabel('Power (dB)');
title('Data Frequency Spectrum');
% Identify noise frequency range and determine cutoff frequency
% This step depends on your specific data and noise characteristics
% For example, if noise is identified above 300 Hz, set cutoff frequency
cutoff_freq = 300; % Example cutoff frequency in Hz
cutoff_index = cutoff_freq/(Fs/N); % Find corresponding index in FFT result
% Display cutoff frequency on plot for visualization
hold on;
line([cutoff_freq/1000, cutoff_freq/1000], ylim, 'Color', 'red', 'LineStyle', '--');
legend('Frequency Spectrum', 'Cutoff Frequency');
Improvements and Additions:
- FFT Scaling: The FFT result is divided by the length of the data N to correctly scale the amplitude. This ensures that the FFT magnitude does not depend on the length of the data.
- Power in Decibels: The amplitude spectrum is converted into decibels to better represent the dynamic range.
- Visualizing Cutoff Frequency: Adding a specific cutoff frequency to the plot helps visually identify which frequency components to filter out for noise reduction. However, the actual cutoff frequency should be selected based on the analysis of your data.
I hope this helps!
