有关snr计算程序

TJ球

我有一组adc点数值
在matlab上进行fft运算
同时为了突显谐波
我加了hann视窗
再依图表上显示主波及谐波db数值
进行手动运算
主波=2175.3246
一次谐波=0.1865
二次谐波=1.0568
snr=20*log10(2175.3246/sqrt(0.1865*0.1865+1.0568*1.0568))=66.137 db
想请问各位的是：
1.我计算snr的程序是否有误      
2.我加了hann视窗是否会影响snr数值
以下是我的程式码
x=load('cc.txt');
M=length(x);
x=x.*hann(M);
N=4095;
fs=96000;
y=fft(x,N);
fn=0:N-1;
f=fs*fn/N;
plot( f(1:N/2),(abs( y(1:N/2))));
axis([0 10000 -5 2800]); grid;

谢谢