Спектральный поток для звуковых сигналов и слуховых спектрограмм
Читайте в звуковом файле, вычислите поток с помощью параметров по умолчанию, и затем постройте результаты.
[audioIn,fs] = audioread('Counting-16-44p1-mono-15secs.wav'); flux = spectralFlux(audioIn,fs); t = linspace(0,size(audioIn,1)/fs,size(flux,1)); plot(t,flux) xlabel('Time (s)') ylabel('Flux')
Читайте в звуковом файле и затем вычислите mel спектрограмму с помощью melSpectrogram функция. Вычислите поток mel спектрограммы в зависимости от времени. Постройте график результатов.
[audioIn,fs] = audioread('Counting-16-44p1-mono-15secs.wav'); [s,cf,t] = melSpectrogram(audioIn,fs); flux = spectralFlux(s,cf); plot(t,flux) xlabel('Time (s)') ylabel('Flux')
Читайте в звуковом файле.
[audioIn,fs] = audioread('Counting-16-44p1-mono-15secs.wav');Вычислите поток спектра мощности в зависимости от времени. Вычислите поток для Окон Хэмминга на 50 мс данных с перекрытием на 25 мс. Используйте диапазон от 62,5 Гц до fs/2 для вычисления потока. Постройте график результатов.
flux = spectralFlux(audioIn,fs, ... 'Window',hamming(round(0.05*fs)), ... 'OverlapLength',round(0.025*fs), ... 'Range',[62.5,fs/2]); t = linspace(0,size(audioIn,1)/fs,size(flux,1)); plot(t,flux) xlabel('Time (s)') ylabel('Flux')
Спектральный поток измеряет изменение в последовательных спектрах. Чтобы вычислить спектральный поток для сигнала передачи потокового аудио, необходимо ввести по крайней мере две системы координат аудиоданных.
Создайте dsp.AudioFileReader возразите, чтобы читать в покадровых аудиоданных. Создайте dsp.AsyncBuffer буферизовать аудио в перекрытые системы координат. Создайте dsp.SignalSink регистрировать спектральное вычисление потока.
fileReader = dsp.AudioFileReader('Counting-16-44p1-mono-15secs.wav');
buff = dsp.AsyncBuffer;
logger = dsp.SignalSink;В цикле аудиопотока:
Читайте в системе координат аудиоданных из вашего источника.
Запишите аудиоданные в dsp.AsyncBuffer
Если система координат данных доступна от буфера, считайте систему координат и один транзитный участок данных с перекрытием, равным выборкам на систему координат. Это представляет две новых аудио системы координат.
Вычислите спектральный поток для двух новых аудио систем координат.
Регистрируйте спектральный поток для более позднего графического вывода. Поскольку поток задан текущей системой координат и предыдущей системой координат, и потому что условие перед первой системой координат неизвестно функции, спектральный поток выводит поток нуля для первой системы координат. Регистрируйте только второе значение выход от spectralFlux.
fs = fileReader.SampleRate; samplesPerFrame = round(fs*0.05); samplesOverlap = round(fs*0.025); samplesPerHop = samplesPerFrame - samplesOverlap; win = hamming(samplesPerFrame); while ~isDone(fileReader) audioIn = fileReader(); write(buff,audioIn); while buff.NumUnreadSamples >= samplesPerHop audioBuffered = read(buff,samplesPerFrame+samplesPerHop,samplesPerFrame); flux = spectralFlux(audioBuffered,fs, ... 'Window',win, ... 'OverlapLength',samplesOverlap); logger(flux(end)) end end release(fileReader)
Постройте записанные данные.
plot(logger.Buffer)
ylabel('Flux')
Спектральный поток вычисляется как описано в [1]:
где
sk является спектральным значением в интервале k.
b 1 и b 2 является ребрами полосы в интервалах, по которым можно вычислить спектральный поток.
P является типом нормы. Можно задать тип нормы с помощью NormType.
[1] Scheirer, E. и М. Слэни. "Конструкция и Оценка Устойчивого Различителя Речи/Музыки Мультифункции". Международная конференция IEEE по вопросам Акустики, Речи и Обработки сигналов. Объем 2, 1997, стр 1221–1224.