Mercurial > hg > aimmat
annotate aim-mat/tools/getmaximumrightof.m @ 4:537f939baef0 tip
various bug fixes and changed copyright message
author | Stefan Bleeck <bleeck@gmail.com> |
---|---|
date | Tue, 16 Aug 2011 14:37:17 +0100 |
parents | 20ada0af3d7d |
children |
rev | line source |
---|---|
tomwalters@0 | 1 % tool |
tomwalters@0 | 2 % |
tomwalters@0 | 3 % INPUT VALUES: |
tomwalters@0 | 4 % |
tomwalters@0 | 5 % RETURN VALUE: |
tomwalters@0 | 6 % |
tomwalters@0 | 7 % |
bleeck@3 | 8 % (c) 2011, University of Southampton |
bleeck@3 | 9 % Maintained by Stefan Bleeck (bleeck@gmail.com) |
bleeck@3 | 10 % download of current version is on the soundsoftware site: |
bleeck@3 | 11 % http://code.soundsoftware.ac.uk/projects/aimmat |
bleeck@3 | 12 % documentation and everything is on http://www.acousticscale.org |
bleeck@3 | 13 |
tomwalters@0 | 14 |
tomwalters@0 | 15 function [pos_max,maxval,index]=getmaximumrightof(where,maxpos,minpos,maxs,mins) |
tomwalters@0 | 16 % usage: [pos_max,maxval]=getmaximumrightof(where,womins,minvals,womaxs,maxvals) |
tomwalters@0 | 17 % returns the maximum that is right of point "where" |
tomwalters@0 | 18 % mins and maxs must be complete sets of minimums and maximums |
tomwalters@0 | 19 % interwoven |
tomwalters@0 | 20 |
tomwalters@0 | 21 pos_max=[]; |
tomwalters@0 | 22 maxval=[]; |
tomwalters@0 | 23 index=-1; |
tomwalters@0 | 24 |
tomwalters@0 | 25 if where > maxpos(end) |
tomwalters@0 | 26 return; |
tomwalters@0 | 27 error('getmaximumrightof:: no maximum right of point'); |
tomwalters@0 | 28 end |
tomwalters@0 | 29 |
tomwalters@0 | 30 |
tomwalters@0 | 31 nr =length(where); |
tomwalters@0 | 32 for j=1:nr |
tomwalters@0 | 33 cwhere=where(j); |
tomwalters@0 | 34 nr_maxs=length(maxpos); |
tomwalters@0 | 35 for i=1:nr_maxs |
tomwalters@0 | 36 if fround(maxpos(i),5) > fround(cwhere,5) |
tomwalters@0 | 37 pos_max(j)=maxpos(i); |
tomwalters@0 | 38 maxval(j)=maxs(i); |
tomwalters@0 | 39 index=j; |
tomwalters@0 | 40 break; |
tomwalters@0 | 41 end |
tomwalters@0 | 42 end |
tomwalters@0 | 43 end |
tomwalters@0 | 44 |
tomwalters@0 | 45 % nr_maxs=length(womaxs); |
tomwalters@0 | 46 % for i=nr_maxs:-1:1 |
tomwalters@0 | 47 % if womaxs(i) < where |
tomwalters@0 | 48 % pos_max=womaxs(i); |
tomwalters@0 | 49 % maxval=maxvals(i); |
tomwalters@0 | 50 % return; |
tomwalters@0 | 51 % end |
tomwalters@0 | 52 % end |
tomwalters@0 | 53 |