fft matlab程序

%% fft
% N:数据点个数
% fs:采样频率
%a:最大的幅值
%p:幅值最大的频率，频率分辨率：fs/N （Hz）
%gs=True:画出频谱图


function  [a,p]=fft_temp(signal,N,fs,gs)

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);
if(gs==True)
    figure()
    title('FFT');
    plot(f(1:(N/2)),P1(1:(N/2)));
    xlabel('频率');
    ylabel('幅值')
end