Julia - 数组中值的乘法

时间:2018-04-18 15:47:36

标签: arrays julia broadcasting elementwise-operations

如果我有一个包含随机值的数组A,我想为每个长度为A的i定义一个数组B,B[i] = (A[i])²

首先,我尝试了以下代码:

using Distributions

A = rand(Uniform(1,10),1,20)

B = A

for i in 1:20
       B[i] = (A[i])^2
end

完成这些操作后,我有A = B

A
1×20 Array{Float64,2}:
 26.0478  5.36654  99.675  23.18  …  1.54846  91.3444  9.41496  2.91666

B
1×20 Array{Float64,2}:
 26.0478  5.36654  99.675  23.18  …  1.54846  91.3444  9.41496  2.91666

所以我尝试了另一种方法:

B = A^2

出现以下错误:

ERROR: DimensionMismatch("A has dimensions (1,20) but B has dimensions (1,20)")

例如,如果我这样做,B = 2*A它可以正常工作......

有些想法可以提供帮助吗?

谢谢

1 个答案:

答案 0 :(得分:5)

我认为最简单的解决方案是通过.使用广播。

julia> A = rand(Uniform(1,10),1,20)
1×20 Array{Float64,2}:
 8.84251  1.90331  8.5116  2.50216  …  1.67195  9.68628  4.05879  1.50231

julia> B = A .^ 2
1×20 Array{Float64,2}:
 78.19  3.62257  72.4473  6.26081  …  2.7954  93.8241  16.4738  2.25693

这可以满足您的期望:

julia> A
1×20 Array{Float64,2}:
 8.84251  1.90331  8.5116  2.50216  …  1.67195  9.68628  4.05879  1.50231

julia> B
1×20 Array{Float64,2}:
 78.19  3.62257  72.4473  6.26081  …  2.7954  93.8241  16.4738  2.25693

这个.运算符可以处理任何函数,它可以将所有函数广播(应用)到数组/对象的每个元素。

以下是自定义函数的示例:

julia> my_func(a) = a * 5
my_func (generic function with 1 method)

julia> my_func.(1:5)
5-element Array{Int64,1}:
  5
 10
 15
 20
 25

有关详细信息,请查看the docs