0

这是在 Julia 中编码的遗传算法部分。代码编写如下：

popc= [individual(rand(0:1,nvar),[]) for i in 1:nc/4,j in 1:2];
for k=1:nc/4
    #select firdt parent
    i1=rand(1:npop);
    p1=pop[i1];
    #select second parent
    i2=rand(1:npop);
    if i1==i2
        i2=rand(1:npop);
    end
    p2=pop[i2]
    #apply crossover
    m=singlepointcrossover(p1.position,p2.position);
    append!(popc[k,1].position, m[1]);
    append!(popc[k,2].position, m[2]);
end

function singlepointcrossover(x1,x2)

    nvar=length(x1);

    cutpoint=rand(1:nvar-1);

    y1=append!(x1[1:cutpoint],x2[cutpoint+1:end]);

    y2=append!(x2[1:cutpoint],x1[cutpoint+1:end]);

    return y1,y2
end

但它有这个错误。你能帮帮我吗？为什么会这样？

ArgumentError: invalid index: 1.0
getindex(::Array{individual,2}, ::Float64, ::Int64) at abstractarray.jl:883
macro expansion at GA.juliarc.jl:87 [inlined]
anonymous at <missing>:?
include_string(::String, ::String) at loading.jl:522
include_string(::String, ::String, ::Int64) at eval.jl:30
include_string(::Module, ::String, ::String, ::Int64, ::Vararg{Int64,N} where N) at eval.jl:34
(::Atom.##102#107{String,Int64,String})() at eval.jl:82
withpath(::Atom.##102#107{String,Int64,String}, ::String) at utils.jl:30
withpath(::Function, ::String) at eval.jl:38
hideprompt(::Atom.##101#106{String,Int64,String}) at repl.jl:67
macro expansion at eval.jl:80 [inlined]
(::Atom.##100#105{Dict{String,Any}})() at task.jl:80
4

2 回答 2

2

问题是/运算符为整数参数提供浮点结果，而浮点结果不能用于索引Array. 您可以Array使用Integer.

/(x, y)

右除法运算符：乘以右边x的倒数y。给出整数参数的浮点结果。

for k=1:nc/4

1:nc/4将创建一个Float64范围k，aFloat64稍后用于在您的代码中的索引中append!(popc[k,1].position, m[1]);。因此，您应该制作k一个Integer.

如果nc是一个整数，你应该使用欧几里得除法div(nc, 4)或简单地nc ÷ 4，或位移运算符nc >> 2和nc >>> 2（对于欧几里得除以 2^n，你应该移位 n）。他们都会给出整数参数的整数结果。

如果nc它本身是一个浮点数，您可能应该使用@Colin T Bowers 指出的选项之一。

popc= [individual(rand(0:1,nvar),[]) for i in 1:nc/4,j in 1:2];

第一行没有错误，因为这里没有i用于索引。最好nc/4用我上面列出的选项之一替换。

于 2019-05-05T13:29:14.463 回答
1

Julia 中Float64的分数始终输出，即使答案可以精确转换为Int.

重要的是，在您的情况下，请注意Int可用于索引数组，但Float64不能。所以你需要调整：

for k=1:nc/2

至

for k=1:Int(nc/2)

这样您的索引k将是类型Int，而不是Float64。

如果nc不能保证是偶数，那么您可能需要使用floor(Int, nc/2)or ceil(Int, nc/2)，这取决于哪个更合适。

于 2019-05-05T11:48:28.387 回答