I want to add noise (10dB,20dB) to my speech samples. I have used this code to add noise but i do not know how much dB of noise is been added because this is random noise
5 views (last 30 days)
Show older comments
[s, fs] = audioread(file);
rng default
noise = 0.01 * rand(size(s));
noisySig = s + noise;
SNR4545= snr(s,noise);
% calc SNR using equation
sigPow_check = 10 * log10(mean(s.^2)); % signal power
noisePow_check = 10 * log10(mean(noise.^2));% noise power
SNR_check = sigPow_check - noisePow_check;
fprintf('\n calculation with equation: \n');
fprintf('SNR: %.2f dB \n', SNR_check) ;
fprintf('noise power: %.2f dB \n', noisePow_check) ;
fprintf('signal power: %.2f dB \n', sigPow_check) ;
Results: calculation with equation: SNR: 8.91 dB noise power: -44.78 dB signal power: -35.87 dB
2 Comments
Souarv De
on 5 Mar 2022
SNR = 20;
Noise = randn(1,length(s)); % Generate initial noise;
Signal_Power = sum(abs(s).*abs(s))/length(s); % Signal Power
Noise_Power = sum(abs(Noise).*abs(Noise))/length(s); % Noise Power
K = (Signal_Power/Noise_Power)*10^(-SNR/10); % Scale factor xdB SNR
New_Noise = sqrt(K)*Noise; % Change Noise vector
x = S + New_Noise;
studentmatlaber
on 6 Mar 2022
Here you have fixed it as SNR = 20. But shouldn't you write as SNR = Signal_Power/Noise_Power?
Answers (0)
See Also
Categories
Find more on Fourier Analysis and Filtering in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!