function [k]= fftph(x,y);
  % 對EH-4的高頻采樣信號(hào)進(jìn)行FFT變化，采樣頻率為192K
% Y=1，2，3；1為功率譜，2為實(shí)部譜，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