fftpl.m

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

function oo= fftpl(x,y);
% 对EH-4的低频采样信号进行FFT变化，采样频率为12K
% Y=1，2，3；1为功率谱，2为实部谱，3为虚部谱
   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=12000*(0:round(length(x)/2))/length(x); %频率分辨率=采样频率/采样点数
   plot(f,k(1:round(length(x)/2+1)));
end