噪声叠加到信号上

把指定的噪声叠加到信号上去，有标准噪声库NOISEX-92，其中带有白噪声、办公室噪声、工厂噪声、汽车噪声、坦克噪声等等，在信号处理中往往需要把库中的噪声叠加到信号中去，而噪声的采样频率与纯信号的采样频率往往不一致，需要采样频率的校准。

% 2，把指定的噪声叠加到信号上去
% 有标准噪声库NOISEX-92，其中带有白噪声、办公室噪声、工厂噪声、汽车噪声、坦克噪声等等，在信号处理中往往需要把库中的噪声叠加到信号中去，而噪声的采样频率与纯信号的采样频率往往不一致，需要采样频率的校准。
function [YNOISE] = add_noisem（Xfilepath_nameSNRfs）
% add_noisem add determinated noise to a signal.
% X is signal and its sample frequency is fs;
% filepath_name is NOISE‘s path and name and the SNR is signal to noise ratio in dB.
[wavinfs1nbits]=wavread（filepath_name）;
if fs1~=fs
wavin1=resample（wavinfsfs1）;
end
nx=size（X1）;
%{
noisySize=size（wavin11）;
NOISE=zeros（nx1）;
if nx>noisySize
    for i=1:noisySize:nx
        
    end
end
%}
NOISE=wavin1（1:nx）;
NOISE=NOISE-mean（NOISE）;
signal_power = 1/nx*sum（X.*X）;
noise_variance = signal_power / （ 10^（SNR/10） ）;
NOISE=sqrt（noise_variance）/std（NOISE）*NOISE;
Y=X+NOISE;
% 其中X是纯信号，filepath_name是指定噪声文件（.wav）的路径和文件名，SNR是要求的信噪比，fs是信号X的采样频率，Y是带噪信号，NOISE是叠加在信号上的噪声。

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1046  2012-03-03 16:34  add_noisem.m

----------- ---------  ---------- -----  ----

                 1046                    1