MATLAB中movsum函数用法

发布于:2025-03-09 ⋅ 阅读:(16) ⋅ 点赞:(0)

目录

语法

说明

示例

向量的中心移动和

向量的尾部移动和

矩阵的移动和

排除缺失值的移动和

基于样本点计算移动和

仅返回满窗口和


        movsum函数的功能是移动总和。

语法

M = movsum(A,k)
M = movsum(A,[kb kf])
M = movsum(___,dim)
M = movsum(___,nanflag)
M = movsum(___,Name,Value)

说明

        ​M= movsum(A,k) 返回由局部 k 个数据点总和组成的数组,其中每个总和基于 A 的相邻元素的长度为 k 的滑动窗计算得出。当 k 为奇数时,窗以当前位置的元素为中心。当 k 为偶数时,窗以当前元素及其前一个元素为中心。当没有足够的元素填满窗时,窗将自动在端点处截断。当窗口被截断时,只根据窗口内的元素计算总和。M 与 A 的大小相同。

  • 如果 A 是向量,movsum 将沿向量 A 的长度运算。

  • 如果 A 是多维数组,则 movsum 沿 A 的大小不等于 1 的第一个维度进行运算。

        ​M = movsum(A,[kb kf]) 通过长度为 kb+kf+1 的窗口计算和,其中包括当前位置的元素、前面的 kb 个元素和后面的 kf 个元素。​

        M = movsum(___,dim) 为上述任一语法指定 A 的运算维度。例如,如果 A 是矩阵,则 movsum(A,k,2) 沿 A 的列运算,计算每行的 k 个元素的移动总和。

        M = movsum(___,nanflag) 指定包含还是省略 A 中的 NaN 值。例如,movsum(A,k,"omitnan") 在计算每个和时忽略所有 NaN 值。默认情况下,movsum 包括 NaN 值。

        M = movsum(___,Name,Value) 使用一个或多个名称-值参数指定和的其他参数。例如,如果 x 是时间向量,则 movsum(A,k,"SamplePoints",x) 相对于 x 中的时间计算 A 的移动和。

示例

向量的中心移动和

         计算行向量的三点中心移动和。当端点处的窗口中少于三个元素时,将根据可用元素计算和。

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

    12    18    13     3    -6    -6    -1     6    12     9

向量的尾部移动和

        计算行向量的三点尾部移动和。当端点处的窗口中少于三个元素时,movsum 将根据可用数量的元素计算和。

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

     4    12    18    13     3    -6    -6    -1     6    12

矩阵的移动和

        计算矩阵中每行的三点中心移动和。窗从第一行开始,沿水平方向移动到该行的末尾,然后移到第二行,依此类推。维度参数为 2,即跨A的列移动窗。

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

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

M = movsum(A,3,2)
M = 3×3

    12    18    14
    -3    -6    -5
     2     6     7

排除缺失值的移动和

        创建一个包含 NaN 值的行向量。

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

        计算向量的三点中心移动和,不包括 NaN 值。对于包含任一 NaN 值的窗,movsum 使用非 NaN 元素进行计算。

M = movsum(A,3,"omitnan")
M = 1×10

    12    12     7    -3    -6    -5     0     7    12     9

基于样本点计算移动和

        根据时间向量 t,计算 A 中数据的 3 小时中心移动和。

A = [4 8 6 -1 -2 -3];
k = hours(3);
t = datetime(16,1,1,0,0,0) + hours(0:5)
t = 1x6 datetime
   01-Jan-0016 00:00:00   01-Jan-0016 01:00:00   01-Jan-0016 02:00:00   01-Jan-0016 03:00:00   01-Jan-0016 04:00:00   01-Jan-0016 05:00:00

M = movsum(A,k,"SamplePoints",t)
M = 1×6

    12    18    13     3    -6    -5

仅返回满窗口和

        计算行向量的三点中心移动和,但在输出中舍弃使用的点数少于三个的计算。也就是说,只返回从满的三元素窗口计算的和,而舍弃端点计算。

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

    18    13     3    -6    -6    -1     6    12

参数说明

A — 输入数组

        输入数组,指定为向量、矩阵或多维数组。

k — 窗长度

        窗长度,指定为数值或持续时间标量。当 k 是正整数标量时,中心和包括当前位置的元素以及周围的相邻元素。

        例如,movsum(A,3) 计算一个由局部三点和组成的数组。

movsum(A,3) computation. The elements in the sample window are 1, 3, and 5, so the resulting local sum is 9.

[kb kf] — 定向窗长度

        定向窗长度,指定为包含两个元素的数值或持续时间行向量。当 kb 和 kf 为正整数标量时,将基于 kb+kf+1 个元素进行计算。该计算包括当前位置的元素、当前位置前面的 kb 个元素以及当前位置后面的 kf 个元素。

        例如,movsum(A,[2 1]) 计算一个由局部四点和组成的数组。

movsum(A,[2 1]) computation. The elements in the sample window are 4, 1, 3, and 5, so the resulting local sum is 13.

dim — 沿其运算的维度

        沿其运算的维度,指定为正整数标量。如果不指定维度,则默认为第一个大于 1 的数组维度。

        维度 dim 表示 movsum 运算所沿的维度,即指定窗的移动方向。

        以一个 m×n 输入矩阵 A 为例:

  • movsum(A,k,1) 为 A 的每列计算包含 k 个元素的移动和,并返回一个 m×n 矩阵。

    movsum(A,k,1) column-wise operation

  • movsum(A,k,2) 为 A 的每行计算包含 k 个元素的移动和,并返回一个 m×n 矩阵。

    movsum(A,k,2) row-wise operation

nanflag — 缺失值条件

        将可选的参数对组指定为 Name1=Value1,...,NameN=ValueN,其中 Name 是参数名称,Value 是对应的值。名称-值参数必须出现在其他参数之后,但参数对组的顺序无关紧要。

Endpoints — 窗口接近端点时的处理方法

窗口接近端点时的处理方法,指定为以下选项之一:

描述
"shrink" 在窗接近输入端点时缩小窗口的大小,从而只包括现有元素。
"discard" 当窗口不与现有元素完全重叠时,不输出任何和。
"fill" 将不存在的元素替换为 NaN 值。
数值或逻辑标量 将不存在的元素替换为指定的数值或逻辑值。

SamplePoints — 用于计算和的样本点

        用于计算和的样本点,指定为向量。样本点表示数据在 A 中的位置。样本点不需要均匀采样。默认情况下,样本点向量为 [1 2 3 ... ]。

        移动窗相对于样本点定义,样本点必须已排序并且包含唯一的元素。例如,如果 t 是与输入数据对应的时间向量,则 movsum(rand(1,10),3,"SamplePoints",t) 将有一个代表 t(i)-1.5 和 t(i)+1.5 之间时间间隔的窗。

        当样本点向量的数据类型为 datetime 或 duration 时,移动窗长度的类型必须为 duration。

        如果样本点的间距不均匀并且指定了 Endpoints,则其值必须为 "shrink"。


网站公告

今日签到

点亮在社区的每一天
去签到