fftp.m

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

function [k,f]= fftp(x,y,cy);
% 对信号X进行FFT变化，采样频率为cy
% Y=1，2，3；1为功率谱，2为实部谱，3为虚部谱;K为幅值，F为频率
   o=fft(x,length(x));
      if y==1
   k=o.*conj(o)/length(x);
   elseif  y==2
       k=real(o)/length(x); 
   elseif  y==3
       k=imag(o)/length(x);
   end
       
   f=cy*(0:round(length(x)/2))/length(x); %频率分辨率=采样频率/采样点数
   plot(f,k(1:round(length(x)/2+1)));
end