Как нормализовать массив частот - PullRequest
Новая
       9

Как нормализовать массив частот

0 голосов
/ 26 сентября 2018

У меня есть два аудиофайла, взятые из одного файла, который я разделил на файл «сигнала» и файл «шума» - фон. Мне нужно знать доминирующие частоты, схему распределения или частоты этой записи, чтобыуметь сравнивать разные звуки, издаваемые разными животными.

Я выполнил fft для каждого файла, а затем вычел фоновые шумы из сигнала.

Мне все равно, что происходит нижеДля меня 20 кГц и выше 100 кГц - это шум, от которого нужно отказаться.

Амплитуда - это то, что я не могу контролировать, поэтому каждая запись должна быть нормализована.

Каков наилучший способ нормализовать этоданные и сделать сравнения между различными записями статистически жизнеспособными? 

function bindel=binset(raw_data_val,signal,noise)

    %in case all the recording is only noise
    if isempty(signal)
        bindel=nan;
        return
    end

    %frequancy of sampling
    %Fs= 250000;

    %extract the signal parts and noise parts
    %"signal" is an index array of all the elemnts of the 
    %"raw data" array that contain a signal
    signal_data=raw_data_val(signal);
    noise_data=raw_data_val(noise);

    %determine the size of the signal array
    L= size(signal_data,1);

    NFFT = 2^nextpow2(L(1,1));


    Y1 = fft(signal_data,NFFT)/L(1,1);

    del1=smooth(2*abs(Y1(1:NFFT/2+1)));

    Y2 = fft(noise_data,NFFT)/L(1,1);

    del2=smooth(2*abs(Y2(1:NFFT/2+1)));

    del=del1-del2;

    %combine the data into 125 bindels
    binsum=size(del)/125;
    bindel=zeros(1,125);
    for j=1:125,
        bindel(j)= sum(del((j-1)*floor(binsum(1,1))+1:j*floor(binsum(1,1))));
    end

    %%%deleting low freuqencies- testing filter set to change
    %%%everything bellow 20 khz to zero
    %%%normalizing between 1 to 0
    bindel(1:20)=0;
    bindel(100:end)=0;
    norm_bin=(bindel - min(bindel)) / ( max(bindel) - min(bindel) );


    bindel=norm_bin;

end

1 Ответ

0 голосов
/ 26 сентября 2018

Я не думаю, что есть лучший способ нормализовать спектральные данные (зависит от вопроса, на который вы пытаетесь ответить), но, учитывая, что вас не интересует абсолютная амплитуда, а распределение доминантных частот, на которые я бы положилсяПлотность и нормализация по сумме вашего спектра: 

norm_bin = bindel / sum(bindel)

Я предполагаю, что ваш NFFT одинаков для всех записей, которые вы сравниваете, если это не так, нормализуйте с учетом NFFT: 

norm_bin = bindel / mean(bindel)
...