ethe

function [ freq ] = HPS(audiofile)
    CORRECTFACTOR=1;
    
    [sample, fs] = audioread(audiofile);
    fftn = fs;
    sample = sample .* hann(length(sample));
    
    fftr = abs(fft(sample,fftn));
    fftr = fftr(1:fftn/2);