How to get the FFT of a PRC file

2 views (last 30 days)
I have a code to plot information by extracting a file
% %%
clear
[filename, pathname] = uigetfile('*.raw;*.prc', 'Pick raw or processed data file');
N=str2double(filename(5:6));
% load processed file
fid = fopen([pathname filename],'r','b');
A= fread(fid,inf,'*single')';
prcdata=reshape(A,N,[])';
plot((1:size(prcdata,1))./1.16E6, prcdata (:,6))
title(strrep(filename,'_','-'))
fclose(fid);
I want to plot the FFT (fast fourier transform) of the data I plotted in the above code. I tried the following code
% %%
clear
plot(fft(1:size(prcdata,1))./1.16E6, prcdata (:,6))
and also
% %%
y= (prcdata (:,6));
z= fft(y);
plot z;
But none seem to be working I want to exchange the plot I obtained with the plot for the FFT for the same

Accepted Answer

Roger
Roger on 28 Oct 2015
y= prcdata(:,6); x=1:size(prcdata,1); z= fft(y); plot(x,z);
  2 Comments
Annonymous User
Annonymous User on 28 Oct 2015
thank you for the code got the warning when I ran the code Warning: Imaginary parts of complex X and/or Y arguments ignored > In Myfile (line 12) It seems not to be working I need to line 12 was plot (x,z); i was looking for x and y axis to be the same values as my code so tried y=prcdata(:,6)/1.16E6; x=(1:size(prcdata,1)/1.16e6); z=fft(y) plot (x,z); so got the error saying it should be of the same values
Annonymous User
Annonymous User on 28 Oct 2015
thank you the code is working now :) cheers

Sign in to comment.

More Answers (0)

Categories

Find more on Fourier Analysis and Filtering in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!