2

你能帮帮我吗?我想在 Julia 中生成一个结构体。在 matalab 中,我们可以使用 repmat 作为结构体。但在朱莉娅这是不可能的。Julia中有没有像repmat这样的功能?

例如:

npop=20;
struct individual
    position
    cost
end

individual1=individual([],[])

repmat(individual1,npop,1)
repeat(individual1,npop,1)

repmat(individual1,npop,1)
ERROR: MethodError: no method matching repmat(::individual, ::Int64, ::Int64)
Closest candidates are:
  repmat(::Union{AbstractArray{T,1}, AbstractArray{T,2}} where T, ::Int64, ::Int64) at abstractarray
math.jl:311
  repmat(::Union{AbstractArray{T,1}, AbstractArray{T,2}} where T, ::Integer, ::Integer) at abstracta
rraymath.jl:334
  repmat(::AbstractArray{T,1} where T, ::Int64) at abstractarraymath.jl:325
  ...
Stacktrace:
 [1] macro expansion at C:\Users\admin\AppData\Local\JuliaPro-0.6.4.1\pkgs-0.6.4.1\v0.6\Atom\src\rep
l.jl:118 [inlined]
 [2] anonymous at .\<missing>:?

julia>

julia> repmat(individual1,npop,1)
ERROR: MethodError: no method matching repmat(::individual, ::Int64, ::Int64)
Closest candidates are:
  repmat(::Union{AbstractArray{T,1}, AbstractArray{T,2}} where T, ::Int64, ::Int64) at abstractarray
math.jl:311
  repmat(::Union{AbstractArray{T,1}, AbstractArray{T,2}} where T, ::Integer, ::Integer) at abstracta
rraymath.jl:334
  repmat(::AbstractArray{T,1} where T, ::Int64) at abstractarraymath.jl:325
  ...
Stacktrace:
 [1] macro expansion at C:\Users\admin\AppData\Local\JuliaPro-0.6.4.1\pkgs-0.6.4.1\v0.6\Atom\src\rep
l.jl:118 [inlined]
 [2] anonymous at .\<missing>:?

如果使用循环,则有可能

pop=individual([],[])
for i=2:npop

    pop=[pop;individual1]

end

julia> pop
20-element Array{individual,1}:
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])
 individual(Any[], Any[])

但我在 Julia 中找不到像 repmat 这样的函数。

4

1 回答 1

2

如果我理解你想要正确实现的目标(假设你想做一种 ABM),那么这就是实现它的方法:

[individual([],[]) for i in 1:npop]

这样,每个人都将被单独分配(这可能是您想要的)。

position附带说明一下,cost出于性能原因,最好添加类型。这是一个示例,假设position它将是包含位置的两个元素元组的cost向量,并且是浮点向量:

struct individual
    position::Vector{Tuple{Float64,Float64}}
    cost::Vector{Float64}
end

然后你应该写:

[individual(Tuple{Float64,Float64}[],Float64[]) for i in 1:npop]

编辑

请注意,struct在 Julia 中是不可变的。这意味着您不能重新绑定positioncost指向的值。因此,您有两种选择。

选项 1. 将结构定义为可变的,如下所示:

mutable struct individual
    position
    cost
end

然后您编写的代码将起作用。

选项 2 更新而不重新绑定

在不重新绑定的情况下更新不可变结构中的向量。例如(使用您的原始定义)

pop = [individual([],[]) for i in 1:npop]
for i=1:npop
    append!(pop[i].position, rand(0:1,10))
end

或根据需要简单地对其进行初始化:

pop = [individual(rand(0:1,10),[]) for i in 1:npop]

作为附加评论,请注意,在 Julia 中,您通常不需要像在(rand(0:1,10))'. 对向量进行操作通常足够好。并且在选项 2 中(正常情况struct下,您将无法将 2 维对象附加到[]1 维)。

于 2019-05-03T19:30:13.357 回答