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