fftph.m

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

function [k]= fftph(x,y);
  % 对EH-4的高频采样信号进行FFT变化，采样频率为192K
% Y=1，2，3；1为功率谱，2为实部谱，3为虚部谱
   o=fft(x,length(x));
   if y==1
  % k=log(o.*conj(o)/length(x));
  k=abs(o)/length(x);
   elseif  y==2
       k=real(o)/length(x); 
   elseif  y==3
       k=imag(o)/length(x);
   end
       
   f=192000*(0:round(length(x)/2))/length(x)
   plot(f,k(1:round(length(x)/2+1)));
end