function?pinlv（fNSNR）
%利用FFT對(duì)加了高斯白噪聲的單頻正弦信號(hào)進(jìn)行頻率估計(jì)，f:信號(hào)頻率;N:FFT點(diǎn)數(shù)；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）;%最大峰值對(duì)應(yīng)的點(diǎn)
f1=1/N*（k-1）
a=X1（k）;?%保存最大值
X1（k）=0