function [k,f]= fftp(x,y,cy);
% 對信號X進行FFT變化，采樣頻率為cy
% Y=1，2，3；1為功率譜，2為實部譜，3為虛部譜;K為幅值，F為頻率
   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=cy*(0:round(length(x)/2))/length(x); %頻率分辨率=采樣頻率/采樣點數
   plot(f,k(1:round(length(x)/2+1)));
end