我们有时候做算法分析,需要自己给信号添加噪声来验证自己算法的鲁棒性。
以语音信号为例,给信号添加指定信噪比的高斯白噪声程序如下:
data=xx(10000:26000); %读入数据
%添加噪声
snr = 10; %设定信噪比,单位db
%输出参数data是带噪语音,noise为加在信号上的白噪声
[data,noise] = Gnoisegen(data,snr);
添加噪声的函数Gnoisegen:
function [y,noise] = Gnoisegen(x,snr)
noise=randn(size(x)); % 用randn函数产生高斯白噪声
Nx=length(x); % 求出信号x长
signal_power = 1/Nx*sum(x.*x); % 求出信号的平均能量
noise_power=1/Nx*sum(noise.*noise);% 求出噪声的能量
noise_variance = signal_power / ( 10^(snr/10) ); % 计算出噪声设定的方差值
noise=sqrt(noise_variance/noise_power)*noise; % 按噪声的平均能量构成相应的白噪声
y=x+noise; % 合成带噪语音