ToolFreq2Bark.m

% ======================================================================
%> @brief converts frequency to bark
%>
%> @param fInHz: frequency
%> @param cModel: 'Schroeder','Terhardt', 'Zwicker', or 'Traunmuller'
%>
%> @retval bark value
% ======================================================================
function [mel] = ToolFreq2Bark(fInHz, cModel)

    if (nargin < 2)
        cModel  = 'Schroeder';
    end

    % set function handle
    hPitchFunc  = str2func (['aca' cModel]);
    
    mel         = hPitchFunc(fInHz);
end

function [mel] = acaSchroeder(f)
    mel         = 7 * asinh(f/650);
end

function [mel] = acaTerhardt(f)
    mel         = 13.3 * atan(0.75 * f/1000);
end

function [mel] = acaZwicker(f)
    mel         = 13 * atan(0.76 * f/1000) + 3.5 * atan(f/7500);
end

function [mel] = acaTraunmuller(f)
    mel         = 26.81/(1+1960./f) - 0.53;
end