How to get the FFT of a PRC file

3 次查看(过去 30 天)
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

采纳的回答

Roger
Roger 2015-10-28
y= prcdata(:,6); x=1:size(prcdata,1); z= fft(y); plot(x,z);
  2 个评论
Annonymous User
Annonymous User 2015-10-28
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 2015-10-28
thank you the code is working now :) cheers

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Fourier Analysis and Filtering 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by