I have a data of velocity (get from accelerometer) and force data (get from hammer). Now, I need to get mobility from them. I already know that mobility=velocity/force. But I don't know how code it. I just begin learning about signal processing for a week.
function [X, freqs] = easy_fft(x, fs) L = length(x); NFFT = 2^nextpow2(length(x)); X = fft(x, NFFT); X = 2/L*X(1:NFFT/2 + 1); % 2/L is the required scaling freqs = fs/2*linspace(0,1,NFFT/2+1); % Frequency vector end
Then do
fs = 1/y(2, 1); % Convert time step to sample rate [A, freq] = easy_fft(y(:, 2), fs); % Acceleration [F, ~] = easy_fft(y(:, 3), fs); % Force
H = A./F;
omega = 2*pi*freq; % Convert to rad/s mobility = H./(1i*omega); % Acceleration = i*omega*Velocity for a periodic steady-state response
To replicate the plot you showed, take the absolute value of mobility and put it on a log scale vs. frequency.
semilogy(freq, abs(mobility));
Apologies if there are any bugs in what I wrote, I didn't test it on anything. Hopefully you get the idea.
Matlabsolutions.com provides guaranteed satisfaction with a
commitment to complete the work within time. Combined with our meticulous work ethics and extensive domain
experience, We are the ideal partner for all your homework/assignment needs. We pledge to provide 24*7 support
to dissolve all your academic doubts. We are composed of 300+ esteemed Matlab and other experts who have been
empanelled after extensive research and quality check.
Matlabsolutions.com provides undivided attention to each Matlab
assignment order with a methodical approach to solution. Our network span is not restricted to US, UK and Australia rather extends to countries like Singapore, Canada and UAE. Our Matlab assignment help services
include Image Processing Assignments, Electrical Engineering Assignments, Matlab homework help, Matlab Research Paper help, Matlab Simulink help. Get your work
done at the best price in industry.