博客

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)

%改为复信号
X= cos(0:pi/1000000:6*pi)+i*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(X)         %求出噪声功率
SNR_10=10*log10(sigPower/noisePower)            %由信噪比定义求出信噪比，单位为db

%我的方法
power_y = (X*X')/(length(X)-1);
power_yn = (Y*Y')/(length(Y)-1);
power_noise2 = ((Y-X)*(Y-X)')/(length(Y)-1);
10*log10(power_y/power_noise2)

%总结：
%（1）我的方法没错，只是求snr时，10*log10(power_y/power_noise2)错写成了10*log(power_y/power_noise2)