阿里云网站域名备案,网站域名自己做,东莞企业做网站,前端网站建设和维护在信号处理中经常需要把噪声叠加到信号上去#xff0c;在叠加噪声时往往需要满足一定的信噪比#xff0c;这样产生二个问题#xff0c;其一噪声是否按指定的信噪比叠加#xff0c;其二怎么样检验带噪信号中信噪比满足指定的信噪比。在MATLAB中可以用randn产生均值为0方差为…在信号处理中经常需要把噪声叠加到信号上去在叠加噪声时往往需要满足一定的信噪比这样产生二个问题其一噪声是否按指定的信噪比叠加其二怎么样检验带噪信号中信噪比满足指定的信噪比。在MATLAB中可以用randn产生均值为0方差为1的正态分布白噪声但在任意长度下xrandn(1,N)x不一定是均值为0方差为1(有些小小的偏差)这样对后续的计算会产生影响。在这里提供3个函数用于按一定的信噪比把噪声叠加到信号上去同时可检验带噪信号中信噪比。把白噪声叠加到信号上去function [Y,NOISE] noisegen(X,SNR)% noisegen add white Gaussian noise to a signal.% [Y, NOISE] NOISEGEN(X,SNR) adds white Gaussian NOISE to X. The SNR is in dB.NOISErandn(size(X));NOISENOISE-mean(NOISE);signal_power 1/length(X)*sum(X.*X);noise_variance signal_power / ( 10^(SNR/10) );NOISEsqrt(noise_variance)/std(NOISE)*NOISE;YXNOISE;其中X是纯信号SNR是要求的信噪比Y是带噪信号NOISE是叠加在信号上的噪声。检验带噪信号的信噪比信噪比的定义为信号能量 (纯信号)^2SNR-------------------------------------------噪声能量 (带噪信号-纯信号)^2function snrSNR_singlech(I,In)% 计算信噪比函数% I :original signal% In:noisy signal(ie. original signal noise signal)snr0;Pssum(sum((I-mean(mean(I))).^2));%signal powerPnsum(sum((I-In).^2)); %noise powersnr10*log10(Ps/Pn);其中I是纯信号In是带噪信号snr是信噪比
