Subversion Repositories Projects

function VibTestSpec(file)

% Read logfile

log=csvread(file);

fs=11111; % Sample freq

for col=1:10

    % Take one of the measurements

    inSignal=log(:,col)';

    inSignal = inSignal - mean(inSignal);

    n=length(inSignal);

    % Take FFT

    FFT = fft(inSignal);

    %spec = FFT.*conj(FFT)/n

    spec = abs(FFT)/n;

    spec = spec(1:n/2);

    specs(:,col) = spec;

end

f = int32(fs/n*(0:n/2-1)); % Freq

plot(f,specs);

title('Spectrum');

xlabel('Hz');

axis([1 400 0 15])