matlab中的信噪比

matlab中的信噪比

  (2011-03-22 19:52:49)

以高斯噪声为例：若有用信号s(n)的最大幅度am，要求得到的信噪比为p，则p=10log10[(am^2)/b^2]，用这个公式反推出高斯噪声的方差b^2，若s(n)单通道实信号，则Matlab程序就是x=s+b*randn(size(s))；若s(n)是正交双通道信号，则Matlab程序就是x=s+b/sqrt(2)*randn(size(s))。 

 如果s(n)是一个N行、2列的复信号，前后两列各表示实部和虚部，则b/sqrt(2)*randn(size(s))产生的也是N×2的高斯分布噪声，实部和虚部的方差均为b/sqrt(2)。实部和虚部分别产生也可以，但不能用b*randn(size(s))。第一，如果这样产生噪声，那么最终信号的信噪比应该用p=10log10[(am^2)/(2*b^2)]；第二，不能用size(s)，应该用size(c)，c为s(n)的实（虚）部列矢量的长度。

Matlab中计算信噪比方式：

%===========================Happy===================================%
function snr=SNR(I,In)
% 计算信号噪声比函数
% by Qulei
% I :original signal
% In:noisy signal(ie. Original signal + noise signal)
% snr=10*log10(sigma2(I2)/sigma2(I2-I1))

[row,col,nchannel]=size(I);

snr=0;
if nchannel==1%gray image
Ps=sum(sum((I-mean(mean(I))).^2));%signal power
Pn=sum(sum((I-In).^2));%noise power
snr=10*log10(Ps/Pn);
elseif nchannel==3%color image
for i=1:3
Ps=sum(sum((I(:,:,i)-mean(mean(I(:,:,i)))).^2));%signal power
Pn=sum(sum((I(:,:,i)-In(:,:,i)).^2));%noise power
snr=snr+10*log10(Ps/Pn);
end
snr=snr/3;
end
%===================================================================%

转：

      符号的平均功率           每比特信号能量/符号周期                
SNR信噪比  = -------------------------  =  ----------------------------------------- =   
                            噪声的平均功率         噪声单边功率普密度*符号带宽                  
            
                            Eb/T     在限带情况下 T*B=1         Eb
                       = ----------  ==================    ---------
                             N0*B                                                 N0