计算信号信噪比

学习计算信噪比时遇到的一些知识点：函数自相关、内积、正交、帕斯瓦尔能量方程等 

第一步(提取信号)：已知信号频率时，利用傅里叶级数展开 ,将信号表示出来；

​​​​​​​

第二步（加窗）：防止频谱泄露，有如下关系式：M是信号周期数，N是采样点数，fin是输入信号频率，fs是采样频率。利用窗函数对信号进行周期截断,使用汉宁窗hann(N);将信号与窗函数相乘；

第三步（除去窗函数增益） ：提取出来的信号除以窗函数增益，如下

  

第四步（计算功率）：将提取出的信号做FFT，取其模值并平方 得到信号的能量；

ssignal=(abs(fft((signal(1:N).*w)'))).^2;    % Signal PSD

第五步（获取噪声信号） ：用总的信号减去提取的信号就可以得到噪声信号；将提取出的噪声信号做FFT，取其模值并平方 得到噪声的能量；

snoise=(abs(fft((noise(1:N).*w)'))).^2;        % Noise PSD

第六步（找出带宽内的能量）：分别将带宽内的信号能量、噪声能量加起来

pwsignal=sum(ssignal(1:fB));                % Signal power
pwnoise=sum(snoise(1:fB));                    % Noise power

第七步（计算信噪比）

snr=pwsignal/pwnoise;

参考文章

傅里叶变换前传：基础知识（卷积、内积、正交）_M＆Q的博客-CSDN博客_信号的内积

什么是窗函数？

什么是泄露？

【信号与系统|吴大正】4：信号分解、傅里叶变换与信号谱(下)_kodoshinichi的博客-CSDN博客_傅里叶变换分解信号