fftt.m

来自「基于FFT的大地电磁信号资料处理MATLAB方法实现」· M 代码 · 共 11 行

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