fftt.m

基于FFT的大地电磁信号资料处理MATLAB方法实现

function [k,f]= fftt(x,cy);
% 对信号X进行FFT变化，采样频率为cy
% 
  % o=fft(x,length(x));
   
      o=fft(x,16);
       
   f=cy*(0:round(length(x)/2))/length(x); %频率分辨率=采样频率/采样点数
   k=o(1:round(length(x)/2+1));
   
end