function oo= fftpl(x,y);
% 對(duì)EH-4的低頻采樣信號(hào)進(jìn)行FFT變化，采樣頻率為12K
% Y=1，2，3；1為功率譜，2為實(shí)部譜，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); %頻率分辨率=采樣頻率/采樣點(diǎn)數(shù)
   plot(f,k(1:round(length(x)/2+1)));
end