function [k,f]= fftt(x,cy);
% 對(duì)信號(hào)X進(jìn)行FFT變化，采樣頻率為cy
% 
  % o=fft(x,length(x));
   
      o=fft(x,16);
       
   f=cy*(0:round(length(x)/2))/length(x); %頻率分辨率=采樣頻率/采樣點(diǎn)數(shù)
   k=o(1:round(length(x)/2+1));
   
end