我有一个大矩阵,其中包含表单的浮点条目
[ 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.
最佳答案
您可以使用 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
关于matrix - Julia 中矩阵列的平均值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/72433632/