此代码是对非支配点进行排序的代码的一部分。 F
是一个保留pareto前沿的单元格数组。但是朱莉娅给出了这个错误:
BoundsError: attempt to access 0-element Array{Any,1} at index [1]
setindex!(::Array{Any,1}, ::Any, ::Int64) at essentials.jl:183
macro expansion at Nondominated sorting.juliarc.jl:41 [inlined]
请你帮我一下。我该如何解决?
#=pop is a structure including
struct individual
position
cost
Rank
Dominationset
Dominatedcount
Crowdingdostance
end
this pop is used for this code as input
5-element Array{individual,1}:
individual(Any[-3, 4, -2], Any[1.0, 1.0], Any[], Any[], Any[], Any[])
individual(Any[4, 4, 1], Any[1.0, 1.0], Any[], Any[], Any[], Any[])
individual(Any[1, 4, 4], Any[1.0, 1.0], Any[], Any[], Any[], Any[])
individual(Any[4, 4, -2], Any[1.0, 1.0], Any[], Any[], Any[], Any[])
individual(Any[1, 4, 2], Any[0.999999, 1.0], Any[], Any[], Any[], Any[])=#
using JuMP,DataStructures
npop=length(pop);
l=0;
F=[]; # F{1}=[] in matlab it is correct
for i=1:npop
for j=i+1:npop
p=pop[i];
q=pop[j];
if (dominates(p.cost,q.cost))
append!(p.Dominationset,j);
append!(q.Dominatedcount,l+1)
end
if (dominates(q.cost,p.cost))
append!(q.Dominationset,i);
append!(p.Dominatedcount,l+1)
end
pop[i]=p;
pop[j]=q;
end
if pop[i].Dominatedcount==[]
F[1]=[F;i]; #F{1}=[F{1} i] in matlab is used
append!(pop[i].Rank,1);
end
end
主导功能如下:
function dominates(x,y)
b=all(x.<=y) && any(x.<y);
return(b);
end
答案 0 :(得分:2)
您的问题在这里:
F[1]=[F;i]
F
是一个空向量,因此您不能为元素1
分配任何内容。这与Matlab有所不同,后者会自动(并且危险地是imho)动态扩展阵列。而是使用:push!
:
push!(F, ??)
我添加问号是因为我不明白您实际上要在F中输入什么。将F
本身放入向量中,然后 then 将其分配给F
的第一个元素?很抱歉,这对我来说毫无意义。 (编辑:您的Matlab代码F{1}=[F{1} i]
在这里也没有意义。您是否有意写F = [F, i]
?)
此外,您确定需要在所有地方使用Any
吗?无类型的struct
和Any
向量将使您的代码非常慢。