基于FFT的频率估计matlab程序

function pinlv（fNSNR）
%利用FFT对加了高斯白噪声的单频正弦信号进行频率估计，f:信号频率;N:FFT点数；SNR:信噪比


n=0:N-1;
x1=2*sin（2*pi*f*n）;
%a=x1.^2;
%s1=sum（a）/N
x=awgn（x1SNR‘measured‘‘dB‘）;
%b=x.^2;
%s2=sum（b）/N
figure（1）
plot（nx）
X=fft（xN）;
X=abs（X）;
figure（2）
stem（nX）
n1=0:round（N/2）-1;
X1=X（1:round（N/2））;
figure（3）
stem（n1X1）;
m1=max（X1）;%最大峰值
k=find（X1==m1）;%最大峰值对应的点
f1=1/N*（k-1）
a=X1（k）; %保存最大值
X1（k）=0