Find local maxima

specifies additional parameters for finding local maxima using one or more
name-value pair arguments. For example,
`TF`

= islocalmax(___,`Name,Value`

)`islocalmax(A,'SamplePoints',t)`

finds local maxima of
`A`

with respect to the time stamps contained in the time
vector `t`

.

`[`

also returns the prominence corresponding to each element of `TF`

,`P`

] = islocalmax(___)`A`

for any
of the previous syntaxes.

Compute and plot the local maxima of a vector of data.

```
x = 1:100;
A = (1-cos(2*pi*0.01*x)).*sin(2*pi*0.15*x);
TF = islocalmax(A);
plot(x,A,x(TF),A(TF),'r*')
```

Create a matrix of data, and compute the local maxima for each row.

A = 25*diag(ones(5,1)) + rand(5,5); TF = islocalmax(A,2)

`TF = `*5x5 logical array*
0 0 1 0 0
0 1 0 0 0
0 0 1 0 0
0 1 0 1 0
0 1 0 0 0

Compute the local maxima of a vector of data relative to the time stamps in the vector `t`

. Use the `'MinSeparation'`

parameter to compute maxima that are at least 45 minutes apart.

t = hours(linspace(0,3,15)); A = [2 4 6 4 3 7 5 6 5 10 4 -1 -3 -2 0]; TF = islocalmax(A,'MinSeparation',minutes(45),'SamplePoints',t); plot(t,A,t(TF),A(TF),'r*')

Specify a method for indicating consecutive maxima values.

Compute the local maxima of data that contains consecutive maxima values. Indicate the maximum of each flat region based on the first occurrence of that value.

x = 0:0.1:5; A = min(0.75, sin(pi*x)); TF1 = islocalmax(A, 'FlatSelection', 'first'); plot(x,A,x(TF1),A(TF1),'r*')

Indicate the maximum of each flat region with all occurrences of that value.

TF2 = islocalmax(A, 'FlatSelection', 'all'); plot(x,A,x(TF2),A(TF2),'r*')

Select maxima based on their prominence.

Compute the local maxima of a vector of data and their prominence, and then plot them with the data.

x = 1:100; A = peaks(100); A = A(50,:); [TF1,P] = islocalmax(A); P(TF1)

`ans = `*1×2*
1.7703 3.5548

plot(x,A,x(TF1),A(TF1),'r*') axis tight

Compute only the most prominent maximum in the data by specifying a minimum prominence requirement.

TF2 = islocalmax(A,'MinProminence',2); plot(x,A,x(TF2),A(TF2),'r*') axis tight

`A`

— Input datavector | matrix | multidimensional array | table | timetable

Input data, specified as a vector, matrix, multidimensional array, table, or timetable.

**Data Types: **`single`

| `double`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

| `logical`

| `table`

| `timetable`

`dim`

— Operating dimensionpositive integer scalar

Operating dimension, specified as a positive integer scalar. By default,
`islocalmax`

operates along the first dimension whose
size does not equal 1.

For example, if `A`

is a matrix, then
`islocalmax(A,1)`

operates along the rows of
`A`

, computing local maxima for each column.

`islocalmax(A,2)`

operates along the columns of
`A`

, computing local maxima for each row.

**Data Types: **`single`

| `double`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

Specify optional
comma-separated pairs of `Name,Value`

arguments. `Name`

is
the argument name and `Value`

is the corresponding value.
`Name`

must appear inside quotes. You can specify several name and value
pair arguments in any order as
`Name1,Value1,...,NameN,ValueN`

.

`TF = islocalmax(A,'MinProminence',2)`

`'MinProminence'`

— Minimum prominence0 (default) | nonnegative scalar

Minimum prominence, specified as the comma-separated pair consisting
of `'MinProminence'`

and a nonnegative scalar.
`islocalmax`

returns only local maxima whose
prominence is at least the value specified.

**Data Types: **`single`

| `double`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

`'ProminenceWindow'`

— Prominence windowpositive integer scalar | two-element vector of positive integers | positive duration scalar | two-element vector of positive durations

Prominence window, specified as the comma-separated pair consisting of
`'ProminenceWindow'`

and a positive integer scalar,
a two-element vector of positive integers, a positive duration scalar,
or a two-element vector of positive durations. The value defines a
window of neighboring points for which to compute the prominence for
each local maximum.

When the window value is a positive integer scalar
*k*, then the window is centered about each local
maximum and contains *k-1* neighboring elements. If
*k* is even, then the window is centered about the
current and previous elements. If a local maximum is within a flat
region, then `islocalmax`

treats the entire flat region
as the center point of the window.

When the value is a two-element vector `[b f]`

of
positive integers, then the window contains the local maximum,
`b`

elements backward, and `f`

elements forward. If a local maximum is within a flat region, then the
window starts `b`

elements before the first point of
the region and ends `f`

elements after the last point
of the region.

When the input data is a timetable or when
`'SamplePoints'`

is specified as a
`datetime`

or `duration`

vector,
the window value must be of type `duration`

, and the
window is computed relative to the sample points.

**Data Types: **`single`

| `double`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

| `duration`

`'FlatSelection'`

— Flat region indicator`'center'`

(default) | `'first'`

| `'last'`

| `'all'`

Flat region indicator for when a local maximum value is repeated
consecutively, specified as the comma-separated pair consisting of
`'FlatSelection'`

and one of the following:

`'center'`

— Indicate only the center element of a flat region as the local maximum. The element of`TF`

corresponding to the center of the flat is 1, and is 0 for the remaining flat elements.`'first'`

— Indicate only the first element of a flat region as the local maximum. The element of`TF`

corresponding to the start of the flat is 1, and is 0 for the remaining flat elements.`'last'`

— Indicate only the last element of a flat region as the local maximum. The element of`TF`

corresponding to the end of the flat is 1, and is 0 for the remaining flat elements.`'all'`

— Indicate all the elements of a flat region as the local maxima. The elements of`TF`

corresponding to all parts of the flat are 1.

When using the `'MinSeparation'`

or
`'MaxNumExtrema'`

name-value pairs, flat region
points are jointly considered a single maximum point.

`'MinSeparation'`

— Minimum separation0 (default) | nonnegative scalar

Minimum separation between local maxima, specified as the
comma-separated pair consisting of `'MinSeparation'`

and a nonnegative scalar. The separation value is defined in the same
units as the sample points vector, which is ```
[1 2 3
...]
```

by default. When the separation value is greater than
0, `islocalmax`

selects the largest local maximum and
ignores all other local maxima within the specified separation. This
process is repeated until there are no more local maxima
detected.

When the sample points vector has type `datetime`

,
the separation value must have type `duration`

.

**Data Types: **`single`

| `double`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

| `duration`

`'MaxNumExtrema'`

— Maximum number of maximapositive integer scalar

Maximum number of maxima detected, specified as the comma-separated
pair consisting of `'MaxNumExtrema'`

and a positive
integer scalar. `islocalmax`

finds no more than the
specified number of most prominent maxima, which is the length of the
operating dimension by default.

**Data Types: **`single`

| `double`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

`'SamplePoints'`

— Sample pointsvector

Sample points, specified as the comma-separated pair consisting of
`'SamplePoints'`

and a vector. The sample points
represent the location of the data in `A`

. Sample
points do not need to be uniformly sampled, but must be sorted with
unique elements. By default, the sample points vector is ```
[1 2 3
...]
```

.

`islocalmax`

does not support this name-value pair
when the input data is a timetable.

**Data Types: **`double`

| `single`

| `int8`

| `int16`

| `int32`

| `int64`

| `uint8`

| `uint16`

| `uint32`

| `uint64`

| `datetime`

| `duration`

`'DataVariables'`

— Table variablesvariable name | cell array of variable names | numeric vector | logical vector | function handle

Table variables, specified as the comma-separated pair consisting of
`'DataVariables'`

and a variable name, a cell array
of variable names, a numeric vector, a logical vector, or a function
handle. The `'DataVariables'`

value indicates which
columns of an input table or timetable to operate on. This value can be
one of the following:

A character vector specifying a single table variable name

A cell array of character vectors where each element is a table variable name

A vector of table variable indices

A logical vector whose elements each correspond to a table variable, where

`true`

includes the corresponding variable and`false`

excludes itA function handle that takes a table variable as input and returns a logical scalar

The specified table variables must have numeric or
`logical`

type.

**Example: **`'Age'`

**Example: **`{'Height','Weight'}`

**Example: **`@isnumeric`

**Data Types: **`char`

| `cell`

| `double`

| `single`

| `logical`

| `function_handle`

`TF`

— Local maxima indicatorvector | matrix | multidimensional array

Local maxima indicator, returned as a vector, matrix, or multidimensional
array. `TF`

is the same size as
`A`

.

**Data Types: **`logical`

`P`

— Prominencevector | matrix | multidimensional array

Prominence, returned as a vector, matrix, or multidimensional
array. `P`

is the same size as `A`

.

If the input data has a signed or unsigned integer type, then
`P`

is an unsigned integer.

The prominence of a local maximum (or peak) measures how the peak stands out with respect to its height and location relative to other peaks.

To measure the prominence of a peak, first extend a horizontal line from the peak to the left and to the right of the peak. Find where the line intersects the data on the left and on the right, which will either be another peak or the end of the data. Mark these locations as the outer endpoints of the left and right intervals. Next, find the lowest valley in both intervals. Take the larger of these two valleys, and measure the vertical distance from that valley to the peak. This distance is the prominence.

For a vector `x`

, the largest prominence is at
most `max(x)-min(x)`

.

Calculate with arrays that have more rows than fit in memory.

Usage notes and limitations:

Tall timetables are not supported.

You must specify a value for the

`'ProminenceWindow'`

name-value pair.The

`'MaxNumExtrema'`

,`'MinSeparation'`

, and`'SamplePoints'`

name-value pairs are not supported.The value of

`'DataVariables'`

cannot be a function handle.

For more information, see Tall Arrays.

Generate C and C++ code using MATLAB® Coder™.

Usage notes and limitations:

You must enable support for variable-size arrays.

The

`'ProminenceWindow'`

name-value pair is not supported.Tables are not supported.

Find Local
Extrema | `ischange`

| `islocalmin`

| `isoutlier`

| `max`

A modified version of this example exists on your system. Do you want to open this version instead?

You clicked a link that corresponds to this MATLAB command:

Run the command by entering it in the MATLAB Command Window. Web browsers do not support MATLAB commands.

Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .

Select web siteYou can also select a web site from the following list:

Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.

- América Latina (Español)
- Canada (English)
- United States (English)

- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)

- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)