Сгенерируйте 100 секунд розового шума со скоростью дискретизации 44,1 кГц.

fs = 44.1e3;
duration = 100;

y = pinknoise(duration*fs);

Постройте график средней степени спектральной плотности (PSD) сгенерированного розового шума.

[~,freqVec,~,psd] = spectrogram(y,round(0.05*fs),[],[],fs);
meanPSD = mean(psd,2);

semilogx(freqVec,db(meanPSD,"power"))
xlabel('Frequency (Hz)')
ylabel('PSD (dB/Hz)')
title('Power Spectral Density of Pink Noise (Averaged)')
grid on

Сгенерируйте 500 секунд розового шума со скоростью дискретизации 16 кГц.

fs = 16e3;
duration = 500;

y = pinknoise(duration*fs);

Постройте график относительной вероятности амплитуды розового шума. Амплитуда всегда ограничена между $-$1 и 1.

histogram(y,"Normalization","probability","EdgeColor","none")
xlabel("Amplitude")
ylabel("Probability")
title("Relative Probability of Pink Noise Amplitude")
grid on

Создайте 5-секундный стерео-розовый шумовой сигнал со скоростью дискретизации 48 кГц.

fs = 48e3;
duration = 5;
numChan = 2;

pn = pinknoise(duration*fs,numChan);

Слушайте стерео розовый сигнал шума.

sound(pn,fs)

Каналы функции розового шума генерируются независимо. Обратите внимание, что внедорожные коэффициенты корреляции близки к нулю (некоррелированы).

R = corrcoef(pn(:,1),pn(:,2))

R = 2×2

    1.0000   -0.0190
   -0.0190    1.0000

Коррелированный и некоррелированный розовый шум имеет различные психоакустические эффекты. Когда шум коррелирует, звук менее эмбиентный и более централизованный. Чтобы прослушать коррелированный розовый шум, отправьте один канал сигнала розового шума на свое стерео устройство. Эффект наиболее выражен при использовании наушников.

sound([pn(:,1),pn(:,1)],fs)

Чтение в аудио файла.

[audioIn,fs] = audioread("MainStreetOne-16-16-mono-12secs.wav");

Создайте сигнал розового шума того же размера и типа данных, что и audioIn.

noise = pinknoise(size(audioIn),'like',audioIn);

Добавьте розовый шум в аудиосигнал и затем прослушайте первые 5 секунд.

noisyMainStreet = noise + audioIn;
sound(noisyMainStreet(1:fs*5,:),fs)

The pinknoise функция генерирует аппроксимацию $-$29,5 дБ сигнала, который близок к степени аудиосигнала.

noisePower = sum(noise.^2,1)/size(noise,1);
signalPower = sum(audioIn.^2,1)/size(audioIn,1);
snr = 10*log10(signalPower./noisePower)

snr = 1.9791

noisePowerdB = 10*log10(noisePower)

noisePowerdB = -29.6665

signalPowerdB = 10*log10(signalPower)

signalPowerdB = -27.6874

Смешайте входной аудио с сгенерированным розовым шумом при 8 дБ ОСШ.

desiredSNR = 8;
scaleFactor = sqrt(signalPower./(noisePower*(10^(desiredSNR/10))));

noise = noise.*scaleFactor;

Проверьте результирующий ОСШ 8 дБ, а затем прослушайте первые 5 секунд.

noisePower = sum(noise.^2,1)/size(noise,1);
snr = 10*log10(signalPower./noisePower)

snr = 8.0000

noisyMainStreet = noise + audioIn;
sound(noisyMainStreet(1:fs*5,:),fs)