MATLAB中awgn()加噪后自定义代码计算SNR偏差问题求助

你的计算结果偏差巨大，核心是采样率计算错误、FFT dB转换错误、功率计算逻辑混乱这三个关键问题，逐个修正后就能得到符合预期的SNR结果。

1. 采样率与频率轴的致命错误 #

原代码中fs = (0:len-1)/t;完全错误，t是时间向量，采样率应为相邻时间点的倒数：fs = 1/(t(2)-t(1));。错误的fs导致后续频率轴TF完全失真，所有频域分析失去意义。

正确的频率轴计算：

fs = 1/(t(2)-t(1)); % 正确采样率：10kHz
TF = (-len/2:len/2-1)*(fs/len); % 对称频率轴

2. FFT的dB转换错误 #

MATLAB中log()是自然对数，计算dB必须用log10()，原代码的10*log(...)会导致dB值被放大约2.3倍，直接影响功率计算结果。同时要先对FFT结果做线性域归一化，再转dB：

Y_linear = abs(fftshift(fft(y)))/len; % 归一化到单频点功率
Y_dB = 10*log10(Y_linear + eps); % 加eps避免log10(0)报错

3. 信号与噪声功率计算逻辑错误 #

原代码在dB域直接平均峰值、用负峰值找噪声底的方法极不严谨，建议在线性域计算功率后再转dB，或者直接在时域计算（最准确）：

• 时域信号功率：P_signal = mean(x.^2);
• 时域总功率：P_total = mean(y.^2);
• 噪声功率：P_noise = P_total - P_signal;
• SNR（dB）：SNR_dB = 10*log10(P_signal/P_noise);

如果非要用频域计算，需注意将频域功率密度转换为总功率（功率密度×带宽）。

修正后的完整代码 #

f = 1000;
t = 0:.0001:1;
x = cos(2*pi*f*t);
len = length(x);
fs = 1/(t(2)-t(1)); % 修正：正确采样率

% 原始信号时域图
plot(t,x)
xlim([0 .01])
title('Signal "x"')
pause

% FFT of x（修正频率轴）
X_linear = abs(fftshift(fft(x)))/len;
TF = (-len/2:len/2-1)*(fs/len);
plot(TF, 10*log10(X_linear + eps))
title('FFT of "x" (dB)')
pause

% 添加10dB SNR噪声
y = awgn(x,10);

% 带噪信号时域图
plot(t,y)
xlim([0 .01])
title('Signal "y"')
pause

% FFT of y（修正dB转换）
Y_linear = abs(fftshift(fft(y)))/len;
Y_dB = 10*log10(Y_linear + eps);
plot(TF,Y_dB)
grid
title('FFT of "y" (dB)')
xlim([500 1500])
pause

% --- 方法1：时域计算SNR（最准确）---
P_signal = mean(x.^2);
P_total = mean(y.^2);
P_noise = P_total - P_signal;
SNR_dB_time = 10*log10(P_signal/P_noise);
fprintf('时域计算SNR：%.2f dB\n', SNR_dB_time);

% --- 方法2：频域计算SNR ---
% 信号频段：900-1100Hz
sig_band = (TF >= 900) & (TF <= 1100);
P_signal_freq = sum(Y_linear(sig_band)) * fs/len; % 转换为总功率
% 噪声频段：1200-9000Hz（避开DC和信号）
noise_band = (TF >= 1200) & (TF <= 9000);
P_noise_density = mean(Y_linear(noise_band));
P_noise_freq = P_noise_density * (fs - 2200); % 噪声总带宽
SNR_dB_freq = 10*log10(P_signal_freq/P_noise_freq);
fprintf('频域计算SNR：%.2f dB\n', SNR_dB_freq);

结果说明 #

修正后，时域计算的SNR会非常接近10dB（误差来自awgn的随机噪声），频域计算也会在10dB左右波动，完全符合预期。

内容的提问来源于stack exchange，提问作者Digit_Fish