fft matlab程序

%% fft
% N:数据点个数
% fs:采样频率
function  [a,p]=fft_temp(signal,N,fs)

n=0:N-1;
T=1/fs;
orig_fft=fft(signal,N);
P2=abs(orig_fft/N);
P1=P2(1:N/2+1);
P1(2:end-1)=2*P1(2:end-1);
f=fs*(0:(N))/N;
[a,n]=max(P1);
p=f(n);