tomwalters@0: % support file for 'aim-mat'
tomwalters@0: %
tomwalters@0: % (c) 2011, University of Southampton
bleeck@3: % Maintained by Stefan Bleeck (bleeck@gmail.com)
bleeck@3: % download of current version is on the soundsoftware site: 
bleeck@3: % http://code.soundsoftware.ac.uk/projects/aimmat
bleeck@3: % documentation and everything is on http://www.acousticscale.org
bleeck@3: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: %
tomwalters@0: %    Estimate fr from fpeak 
tomwalters@0: %    Toshio IRINO
tomwalters@0: %    10 June 98
tomwalters@0: %
tomwalters@0: % function [fpeak, ERBw] = Fr2Fpeak(n,b,c,fr)
tomwalters@0: %    INPUT:  n,b,c : gammachirp param.
tomwalters@0: %            fr    : fr
tomwalters@0: %    OUTPUT: fpeak : peak freq.
tomwalters@0: %            ERBw  : ERBw(fr)
tomwalters@0: %
tomwalters@0: function [fpeak, ERBw] = Fr2Fpeak(n,b,c,fr)
tomwalters@0: 
tomwalters@0: if nargin < 4, help Fr2Fpeak; end;
tomwalters@0: 
tomwalters@0: n = n(:);
tomwalters@0: b = b(:);
tomwalters@0: c = c(:);
tomwalters@0: fr = fr(:);
tomwalters@0: 
tomwalters@0: [dummy ERBw] = Freq2ERB(fr);
tomwalters@0: fpeak = fr + c.*ERBw.*b./n;