提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档
前言
在通信系统仿真中,经常需要评估通信模型的噪声性能,而最常用到的噪声就是加性高斯白噪声。MATLAB中的awgn函数用于在输入信号上增加高斯白噪声,可以很方便地用于高斯白噪声信道的仿真。下面使用awgn函数在信号上叠加高斯白噪声,噪声的功率以信号功率为参照,使叠加信号的信噪比等于设定的SNR值。
一、仿真代码
代码如下:
f = 1; % 信号频率,单位Hz
fs = 1000; % 采样率,单位samples per second
T = 10; % 信号时长,单位s
t = 0:1/fs:T;
x = sin(2*pi*f*t);
snr = 20; % 单位dB
y = awgn(x, snr,'measured','db'); % 首先测量输入信号x的功率,然后加入噪声使输出信号y的信噪比=snr(dB)
figure()
subplot(2,1,1);
plot(t,x);
ylim([-2 2]);
grid on;
title('正弦信号x')
subplot(2,1,2);
plot(t,y);
ylim([-2 2]);
grid on;
title('叠加了加性高斯白噪声后的正弦信号')
% 计算噪声功率(也即噪声方差)
n = y - x;
var(n) % 噪声功率=0.005,是信号功率的1%
二、关于仿真代码的几点说明
(1)awgn使用的噪声相对功率的衡量方式是SNR(信号噪声功率比),默认取dB形式。
(2)awgn函数的第三个参数为信号功率,当设置为’measured’时,awgn首先测量测量输入信号x的功率,然后根据设定的信噪比(SNR)加入具有特定功率的噪声,使输出信号y的信噪比=snr(dB)。关于信号功率的设置,推荐使用’measured’方法,因为这种方法简单,并且SNR容易控制。
(3)awgn函数的第三个参数设置为一个数值时,如0,表示的是信号功率为0 dBW。
(4)awgn函数还有其他用法,可以通过在命令行输入help awgn得到。
三、仿真结果
仿真结果如下:
该结果与上一篇文章中使用randn的方法一致,randn方法中可以更加清楚的看到信噪比是如何设置的。上一篇文章:
MATLAB使用randn仿真加性高斯白噪声信道
总结
下一篇文章将给出一个通信系统示例,说明如何在一个QPSK通信系统中使用加性高斯白噪声(Additive Gaussian white noise)。