Whew code is going really well again today, made a first test version of a simple band identification routine. Really simple, yet it works very well, almost as well as Bandleader does... I've put up a snapshot. The algorithm is from a small part of Bandleader, the algorithm there is:

% Function to find all local peaks in a vector. A local peak isdefined % as a value x(n) satisfying x(n-2) < x(n-1) < x(n) > x(n+1) > x(n+2), % and for which the peak amplitude is greater than twice the local mean, % where the local mean is taken over +/- m samples.