matlab 求信噪比

 信噪比是衡量信号质量的重要指标, 尤其是在比较滤波算法的过程中 有着重要的参考价值。

先看看维基百科的定义：

clear all;
clc ;
X= sqrt(2)*sin(0:pi/1000000:6*pi);              %产生正弦信号
Y = awgn(X,14,'measured');                      %加入信噪比为10db的噪声，加入前预估信号的功率（强度）
sigPower = sum(abs(X).^2)/length(X)             %求出信号功率
noisePower=sum(abs(Y-X).^2)/length(Y-X)         %求出噪声功率
SNR_10=10*log10(sigPower/noisePower)            %由信噪比定义求出信噪比，单位为db
b=snr(X,Y-X)                             % snr(a,b) : a是原始信号，b是噪声信号
snr(X)

输出：
SNR_10 =

   14.0033

b =

   14.0033

问题：

  MIT-BIH 那些数据吧，请问你知道原始ECG数据的信噪比是如何求出来的嘛？

参考：http://blog.sina.com.cn/s/blog_b09328370102vr8g.html

           https://zhuanlan.zhihu.com/p/33033117

         例子2：

clear all;clc; close all;
%%%  calucate the signal snr 
X = sqrt(2)*sin(0:pi/1000000:6*pi);               %产生正弦信号
%% select this mode 
Y = awgn(X,20,'measured');                          %加入信噪比为10db的噪声，加入前预估信号的功率（强度）
z=awgn( X,20);

figure
subplot(311);plot(X);
subplot(312);plot(Y)
subplot(313);plot(z)

sigPower = sum(abs(X).^2)/length(X)           %求出信号功率

noisePower=sum(abs(Y-X).^2)/length(Y-X)  %求出噪声功率

SNR=10*log10(sigPower/noisePower)          %由信噪比定义求出信噪比，单位为db\

by=snr(X,Y-X)                             % snr(a,b) : a是原始信号，b是噪声信号
bz=snr(X,z-X)

%%%%caclute the  MSE
mx=(sum(X)-sum(Y)).^2/length(X);
mz=(sum(X)-sum(z)).^2/length(X);