求信噪比计算公式

syy2468

function y=snr(x1,x2);%x1是原始信号，x2是降噪后信号
N=length(x1);
y1=sum(x2.^2); %去噪后信号的能量
y2=sum((x1-x2).^2);
y=10*log10(y1/y2);

信噪比是去噪后信号的能量比上取出的噪声能量。

syy2468

function y=snr(x1,x2);%x1是原始信号，x2是降噪后信号
N=length(x1);
y1=sum(x2.^2); %去噪后信号的能量
y2=sum((x1-x2).^2);
y=10*log10(y1/y2);

信噪比是去噪后信号的能量比上去除的噪声能量。

syy2468

这个程序求出的只是输入信噪比，输出信噪比怎么求呢

jytan

原帖由 zhangnan3509 于 2007-8-2 21:22 发表
这应该和eight的精华贴 相得益彰啊
http://www.chinavib.com/forum/vi ... p%3Bfilter%3Ddigest
上面是仿真信号的信噪比总结
而本贴则是针对实际数据处理降噪后，来计算信噪比的程序。使得精华贴更完善了
functi ...

偶然翻到这个老帖子，发表点自己的看法，
大家应该都注意到了，输入信噪比用“消噪后的信号能量/消去的噪声能量”来估计，
但是sum(x1-x2)^2能表示消去的噪声能量吗？
如果消噪方法不产生相移的话，
我觉得此方法似乎可行，
但有些消噪方法并不能保证消噪后的信号还跟输入信号保持相同相位。
所以用sum(x1.^2-x2.^2)来估计消去的噪声能量会不会更正确一些？

xiaowei1982

我觉得 jytan 讲的很有道理！

james551304

实际使用中，可能依据环境情况使用合适的方法。比如在超声信号处理中，信噪比是这样计算的，首先设定一个时间窗口，计算

时间窗口内染噪信号的均方根Root-mean-square，就是该信号对零的偏差（记得信号要校正为零）， 接着

计算该时间窗口内

染噪信号的最大幅值，因是射频显示，所以信噪比是SNR=20*log10(Amax/RMS)。

这是美国一个公司采取的算法，具体的原因大家讨论讨论，相互学习。

[ 本帖最后由 james551304 于 2009-10-28 22:20 编辑 ]

chenyuluoyan

如果我要比较两种去噪效果的优劣，能不能用纯净信号来做呢
假设一纯净信号x ，加噪声得到设定的信噪比SNR1 ，去噪后信号y。
那么信噪比SNR2可否这样计算 ：

Ps=sum(abs(x).^2);%纯净信号能量
Po=sum((abs(x)-abs(y)).^2);%去噪信号含有的噪声能量
SNR2=10*log10(Ps/Po);

比较不同去噪方法的SNR2的大小，从而比较其去噪效果优劣。

请大家指导一下，这样是否可行？

鲤鱼

我认为所得的结果应该是原始信号的信噪比吧。

goodstudent

信噪比可以是负值的呀

zxceric

我觉得首先应该搞清楚信噪比的定义

tsfei

回复 9 # zhangnan3509 的帖子

原始信号去噪，就说明原始信号含噪声。信噪比，就是针对信号与噪声而言的，应该是去噪之后的信号与噪声的信噪比，如果原始信号与噪声之间的信噪比，那就是带噪声信号与噪声比了，与信噪比的定义相违背，我觉得你的说是对的。

ufo158

回复 37 # chenyuluoyan 的帖子

个人理解为: 有用信号(或相对较为纯净的信号)  与 噪声的比值关系.

1. 
   
2. 
   
3. function mysnr = myFunc_SP_CalSNR(I,In)
   
4. %%
   
5. %函数功能: 计算信号信噪比.
   
6. %           理解为: 有用信号(或相对较为纯净的信号) 与 噪声的比值关系.
   
7. %输入参数:
   
8. %       I:      N*1,原始信号(仿真信号,不含噪声的信号) 或 滤波降噪***"后"***的信号.
   
9. %                   可理解为:相对纯净的信号,噪声水平相对较低.
   
10. %       In:     N*1,含有噪声的信号(仿真信号+噪声) 或 滤波降噪***"前"***的信号.
    
11. %                   可理解为:噪声水平相对较高的信号.
    
12. %       noise:  (不输入)N*1,噪声信号, noise = In-I;
    
13. %输出参数:
    
14. %       mysnr:  1*1,信号信噪比值(dB).
    
15. %       对于既是仿真信号,又对仿真信号进行降噪的时候,我做如下理解:
    
16. %              仿真信号加入噪声后的信噪比: I:原始纯净信号, In:纯净信号中加入噪声的信号之和.
    
17. %              仿真信号滤波降噪后的信噪比: I:滤波后的信号, In:滤波前的信号(滤波后的信号+滤除的噪声).
    
18. 
    
19. %
    
20. noise = In-I;%噪声
    
21. Ps=sum(sum(I.^2));%signal power,相对纯净信号的功率.
    
22. Pn=sum(sum((In-I).^2));%noise power,噪声信号功率.
    
23. mysnr=10*log10(Ps/Pn);%注意log10对数.

复制代码

按照37L的做法, 计算降噪后的SNR, 但是I使用 纯净信号, 相当于计算纯净信号的信噪比,
不够恰当.

dingdingysu

回复 10 # xuyuling 的帖子

你好。我也总是求的信噪比为负值。我想请问您是怎么解决的？

valleysilence

新来的，连eight大神的帖子都看不了

liuxianming456

长见识了，我觉得这个编程被大家完善了。
如果这个snr公式要是被用于随机共振中SNR定义，不知道如何编程，请赐教！！！