matrix - Julia 中矩阵列的平均值

Question

我有一个大矩阵，其中包含表单的浮点条目

[ a b c d 
  e f g h
  i j k l
  m n o p ]

一些值是离群值，所以我想用相应列中最近的 k 个条目对每个条目取平均值并保留形状。换句话说，对于 k = 3 有这样的东西:

[        a                  b                   c                  d 
      (e + a)/2          (f + b)/2          (g + c)/2          (h + d)/2
    (e + a + i)/3      (f + b + j)/3      (g + c + k)/3      (h + d + l)/3
    (e + i + m)/3      (f + j + n)/3      (g + k + o)/3      (h + l + p)/3   ] 


 etc.

Answer 1

您可以使用 RollingFunctions 和 mapslices 执行此操作:

julia> a = reshape(1:16, 4, 4)
4×4 reshape(::UnitRange{Int64}, 4, 4) with eltype Int64:
 1  5   9  13
 2  6  10  14
 3  7  11  15
 4  8  12  16

julia> using RollingFunctions

julia> mapslices(x -> runmean(x, 3), a, dims = 1)
4×4 Matrix{Float64}:
 1.0  5.0   9.0  13.0
 1.5  5.5   9.5  13.5
 2.0  6.0  10.0  14.0
 3.0  7.0  11.0  15.0