findmaxima.m

来自「Matlab时频分析工具箱,希望能对大家有所帮助啊」· M 代码 · 共 48 行

function [maximaout,yout] = findmaxima(x,F)
%FINDMAXIMA  Find location of local maxima
%  Usage: [maximaout,yout] = findmaxima(x,F);
%
%  F is number of maxima to find.
%  

% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
% 
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
% 
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <http://www.gnu.org/licenses/>.

%  Google for local maxima matlab

% Unwrap to vector
x = x(:);
% Identify whether signal is rising or falling
upordown = sign(diff(x));
% Find points where signal is rising before, falling after
maxflags = [upordown(1)<0; diff(upordown)<0; upordown(end)>0];
maxima   = find(maxflags);

% Sort them with the largest one first.
values=x(maxima);
[Y,ii]=sort(values);
ii=flipud(ii);

% tlen handles the case where fewer than F local maxima are found.
tlen=min(length(maxima),F);

% Keep the largest tlen local maxima, sort them by index.
ikeep=sort(maxima(ii(1:tlen)));

maximaout=zeros(F,1);
yout=zeros(F,1);

% Handle the case where fewer than F local maxima are found.
maximaout(1:tlen)=ikeep;
yout(1:tlen)=x(ikeep);