La traducción de esta página aún no se ha actualizado a la versión más reciente. Haga clic aquí para ver la última versión en inglés.

movmedian

Mediana móvil

Sintaxis

M = movmedian(A,k)

M = movmedian(A,[kb kf])

M = movmedian(___,dim)

M = movmedian(___,nanflag)

M = movmedian(___,Name,Value)

Descripción

M = movmedian(A,k) devuelve un arreglo de valores locales de mediana de k puntos, donde se calcula cada mediana con una ventana deslizante de longitud k a través de elementos vecinos de A. Cuando k es impar, la ventana se centra sobre el elemento en la posición actual. Cuando k es par, la ventana se centra sobre los elementos actuales y anteriores. El tamaño de la ventana está truncado automáticamente en los puntos finales cuando no hay suficientes elementos para rellenar la ventana. Cuando la ventana está truncada, la mediana se toma únicamente sobre los elementos que rellenan la ventana. M es del mismo tamaño que A.

• Si A es un vector, movmedian opera en la longitud del vector A.

• Si A es un arreglo multidimensional, movmedian opera en la primera dimensión de A cuyo tamaño no es igual a 1.

M = movmedian(A,[kb kf]) calcula la mediana con una ventana de longitud kb+kf+1 que incluye el elemento en la posición actual, los elementos kb hacia atrás, y los elementos kf hacia adelante.

M = movmedian(___,dim) especifica la dimensión de A en la que operar para cualquiera de las sintaxis anteriores. Por ejemplo, si A es una matriz, movmedian(A,k,2) opera en las columnas de A, calculando la mediana deslizante de k elementos para cada fila.

M = movmedian(___,nanflag) especifica si incluir u omitir valores NaN en A. Por ejemplo, movmedian(A,k,"omitnan") ignora los valores NaN al calcular cada mediana. De forma predeterminada, movmedian incluye valores NaN.

M = movmedian(___,Name,Value) especifica parámetros adicionales para la mediana móvil utilizando uno o más argumentos de par nombre-valor. Por ejemplo, si x es un vector de valores de tiempo, movmedian(A,k,"SamplePoints",x) calcula la mediana móvil relativa a los tiempos en x.

Ejemplos

Mediana móvil centrada de un vector

Calcule la mediana móvil centrada en tres puntos de un vector fila. Cuando hay menos de tres elementos en la ventana en los puntos finales, tome la media sobre los elementos que están disponibles.

A = [4 8 6 -1 -2 -3 -1 3 4 5];
M = movmedian(A,3)

M = 1×10

    6.0000    6.0000    6.0000   -1.0000   -2.0000   -2.0000   -1.0000    3.0000    4.0000    4.5000

Mediana móvil final de un vector

Calcule la mediana móvil final en tres puntos de un vector fila. Cuando hay menos de tres elementos en la ventana en los puntos finales, tome la media sobre los elementos que están disponibles.

A = [4 8 6 -1 -2 -3 -1 3 4 5];
M = movmedian(A,[2 0])

M = 1×10

     4     6     6     6    -1    -2    -2    -1     3     4

Mediana móvil de una matriz

Calcule la mediana móvil centrada en tres puntos para cada fila de una matriz. La ventana comienza en la primera fila y se desliza horizontalmente hasta el final de la fila. Después, se mueve hacia la segunda fila, y así sucesivamente. El argumento de dimensión es dos, que desliza la ventana a través de las columnas de A.

A = [4 8 6; -1 -2 -3; -1 3 4]

A = 3×3

     4     8     6
    -1    -2    -3
    -1     3     4

M = movmedian(A,3,2)

M = 3×3

    6.0000    6.0000    7.0000
   -1.5000   -2.0000   -2.5000
    1.0000    3.0000    3.5000

Mediana móvil sin incluir valores faltantes

Cree un vector fila que contenga valores NaN.

A = [4 8 NaN -1 -2 -3 NaN 3 4 5];

Calcule la mediana móvil centrada en tres puntos de un vector fila, sin valores NaN. Para ventanas que contienen un valor NaN, movmedian se calcula con elementos que no son NaN.

M = movmedian(A,3,"omitnan")

M = 1×10

    6.0000    6.0000    3.5000   -1.5000   -2.0000   -2.5000         0    3.5000    4.0000    4.5000

Puntos de muestra para una mediana móvil

Calcule una mediana móvil centrada en 3 horas de los datos en A de acuerdo con el vector de tiempo t.

A = [4 8 6 -1 -2 -3];
k = hours(3);
t = datetime(2016,1,1,0,0,0) + hours(0:5)

t = 1x6 datetime
   01-Jan-2016 00:00:00   01-Jan-2016 01:00:00   01-Jan-2016 02:00:00   01-Jan-2016 03:00:00   01-Jan-2016 04:00:00   01-Jan-2016 05:00:00

M = movmedian(A,k,"SamplePoints",t)

M = 1×6

    6.0000    6.0000    6.0000   -1.0000   -2.0000   -2.5000

Devolver únicamente medianas de ventana completa

Calcule la mediana móvil centrada en tres puntos de un vector fila, pero descarte cualquier cálculo que utilice menos de tres puntos desde la salida. En otras palabras, devuelva únicamente las medianas calculadas de una ventana completa de tres elementos, descartando los cálculos de punto final.

A = [4 8 6 -1 -2 -3 -1 3 4 5];
M = movmedian(A,3,"Endpoints","discard")

M = 1×8

     6     6    -1    -2    -2    -1     3     4

Argumentos de entrada

A — Arreglo de entrada

vector | matriz | arreglo multidimensional

Arreglo de entrada, especificado como vector, matriz o arreglo multidimensional.

Tipos de datos: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | logical

k — Longitud de la ventana

escalar numérico o de duración

Longitud de la ventana, especificada como escalar numérico o de duración. Cuando k es un escalar entero positivo, la mediana centrada incluye el elemento en la posición actual, además de los vecinos que lo rodean.

Por ejemplo, movmedian(A,3) calcula un arreglo de valores de mediana locales de tres puntos.

[kb kf] — Longitud de ventana direccional

vector fila numérico o de duración que contiene dos elementos

Longitud de ventana direccional, especificada como vector fila numérico o de duración que contiene dos elementos. Cuando kb y kf son escalares enteros positivos, el cálculo es sobre los elementos kb+kf+1. El cálculo incluye el elemento en la posición actual, los elementos kb anteriores a la posición actual, y los elementos kf posteriores a la posición actual.

Por ejemplo, movmedian(A,[2 1]) calcula un arreglo de valores de mediana locales de cuatro puntos.

dim — Dimensión en la que operar

escalar entero positivo

Dimensión en la que operar, especificada como escalar entero positivo. Si no especifica la dimensión, el valor predeterminado es la primera dimensión del arreglo cuyo tamaño no es igual a 1.

La dimensión dim es la que indica sobre cuál opera movmedian, es decir, la dirección en la que se desliza la ventana especificada.

Considere una matriz de entrada de m por n, A:

• movmedian(A,k,1) calcula la mediana deslizante de k elementos para cada columna de A y devuelve una matriz de m por n.

  

• movmedian(A,k,2) calcula la mediana deslizante de k elementos para cada fila de A y devuelve una matriz de m por n.

  

nanflag — Condición de valor faltante

"includemissing" (predeterminado) | "includenan" | "omitmissing" | "omitnan"

Condición de valor faltante, especificada como uno de estos valores:

• "includemissing" o "includenan": incluyen los valores NaN en A al calcular cada mediana. Si un elemento de la ventana es NaN, el elemento correspondiente en M es NaN. "includemissing" y "includenan" presentan el mismo comportamiento.

• "omitmissing" u "omitnan": ignoran todos los valores NaN en A y calculan cada mediana en menos puntos. Si todos los elementos de la ventana son NaN, el elemento correspondiente en M es NaN. "omitmissing" y "omitnan" presentan el mismo comportamiento.

Argumentos de par nombre-valor

Especifique pares opcionales de argumentos como Name1=Value1,...,NameN=ValueN, donde Name es el nombre del argumento y Value es el valor correspondiente. Los argumentos nombre-valor deben aparecer después de otros argumentos, pero el orden de los pares no importa.

En las versiones anteriores a R2021a, utilice comas para separar cada nombre y valor, y encierre Name entre comillas.

Ejemplo: M = movmedian(A,k,"Endpoints","fill")

Endpoints — Método para tratar las ventanas cercanas a los puntos finales

"shrink" (predeterminado) | "discard" | "fill" | escalar numérico o lógico

Método para tratar las ventanas cercanas a los puntos finales, especificado como una de estas opciones:

Valor	Descripción
"shrink"	Reduce el tamaño de la ventana cerca de los puntos finales de la entrada para incluir únicamente elementos existentes.
"discard"	No muestra valores de mediana de salida cuando la ventana no se solapa por completo con elementos existentes.
"fill"	Reemplaza los elementos no existentes con NaN.
escalar númerico lógico	Reemplaza los elementos no existentes con el valor numérico o lógico especificado.

SamplePoints — Puntos de muestra para calcular medianas

vector

Puntos de muestra para calcular medianas, especificados como vector. Los puntos de muestra representan las ubicaciones de los datos en A. Los puntos de muestra no tienen que ser muestreados de manera uniforme. De forma predeterminada, el vector de puntos de muestra es [1 2 3 ... ].

Las ventanas móviles se definen de forma relativa a los puntos de muestra, que se deben ordenar y contienen elementos únicos. Por ejemplo, si t es un vector de tiempos correspondiente a los datos de entrada, movmedian(rand(1,10),3,"SamplePoints",t) tiene una ventana que representa el intervalo de tiempo entre t(i)-1.5 y t(i)+1.5.

Cuando el vector de puntos de muestra tiene el tipo de datos datetime o duration, la longitud de la ventana móvil debe tener el tipo duration.

Si los puntos de muestra no están espaciados de manera uniforme y se especifica Endpoints, su valor debe ser "shrink".

Capacidades ampliadas

Arreglos altos

Realice cálculos con arreglos que tienen más filas de las que caben en la memoria.

Esta función es compatible con arreglos altos, con las siguientes limitaciones:

• El argumento de par nombre-valor SamplePoints no es compatible.

Para obtener más información, consulte Arreglos altos.

Generación de código C/C++
Genere código C y C++ mediante MATLAB® Coder™.

Entorno basado en subprocesos

Ejecute código en segundo plano con MATLAB® backgroundPool o acelere código con Parallel Computing Toolbox™ ThreadPool.

Esta función es totalmente compatible con entornos basados en subprocesos. Para obtener más información, consulte Ejecutar funciones de MATLAB en entornos basados en subprocesos.

Arreglos GPU

Acelere código mediante la ejecución en una unidad de procesamiento gráfico (GPU) mediante Parallel Computing Toolbox™.

La función movmedian es compatible con entradas de arreglos de GPU con estas notas y limitaciones de uso:

• La longitud máxima de ventana es 31 para los datos de gpuArray.

• El argumento de par nombre-valor SamplePoints no es compatible.

Para obtener más información, consulte Run MATLAB Functions on a GPU (Parallel Computing Toolbox).

Arreglos distribuidos

Realice particiones de arreglos grandes por toda la memoria combinada de su cluster mediante Parallel Computing Toolbox™.

Esta función es totalmente compatible con los arreglos distribuidos. Para obtener más información, consulte Run MATLAB Functions with Distributed Arrays (Parallel Computing Toolbox).

Historial de versiones

Introducido en R2016a

R2023a: Especificar una condición de valor faltante

Incluya u omita valores faltantes en el arreglo de entrada cuando calcule cada mediana usando las opciones "includemissing" u "omitmissing". Estas opciones tienen el mismo comportamiento que las opciones "includenan" y "omitnan", respectivamente.

R2023a: Rendimiento mejorado al calcular sobre una matriz con puntos de muestra

La función movmedian muestra una mejora en su rendimiento al realizar el cálculo sobre una matriz cuando hay puntos de muestra.

Por ejemplo, este código calcula la mediana móvil de una matriz de 300 por 300 con puntos de muestra correspondientes. El código es aproximadamente 1,7 veces más rápido que en la versión anterior.

function timingMovmedian
A = randn(300);
t = sort(rand(300,1));
tic
for k = 1:400
    movmedian(A,0.1,"SamplePoints",t);
end
toc
end

Los tiempos de ejecución aproximados son:

R2022b: 1,14 s

R2023a: 0.67

El código se cronometró en un sistema de prueba Windows^® 10, Intel^® Xeon^® CPU E5-1650 v4 con 3,60 GHz llamando a la función timingMovmedian.

Consulte también

Funciones

• movmean | movmad | movvar | movstd | median

Temas

• Suavizado de datos y detección de valores atípicos