How to get the FFT of a PRC file

2015 10 28
1 个回答
回答已采纳
更新时间：2015 10 28
4 次查看（30 天）

0 个投票

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

1 个投票

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 个）

类别

在 帮助中心 和 File 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