function dom_vector = checkDomination(fitness) Np = size(fitness,1); dom_vector = zeros(Np,1); all_perm = nchoosek(1:Np,2); % Possible permutations 可能的排列 all_perm = [all_perm; [all_perm(:,2) all_perm(:,1)]]; d = dominates(fitness(all_perm(:,1),:),fitness(all_perm(:,2),:)); dominated_particles = unique(all_perm(d==1,2)); dom_vector(dominated_particles) = 1; end

这段代码是一个用于检查支配关系（domination）的函数实现，下面是代码的解释：

1. 首先，获取输入参数 fitness 的行数 Np，也就是个体的数量。

2. 创建一个大小为 Np 的零向量 dom_vector，用于标记被支配的个体。

3. 使用 nchoosek 函数生成了一个包含所有可能排列的矩阵 all_perm，其中每一行表示一对个体的索引。

4. 为了考虑支配关系的双向性，将 all_perm 复制一份，并交换每一对个体的索引，得到一个扩展的 all_perm。

5. 调用 dominates 函数，将 all_perm 中的第一列对应的个体与第二列对应的个体的适应度值传入。dominates 函数用于判断一个个体是否支配另一个个体。

6. 根据 dominates 函数的返回值 d，找到所有被支配的个体的索引，并使用 unique 函数去除重复的索引。

7. 将被支配的个体在 dom_vector 中对应位置设置为 1，表示这些个体被其他个体所支配。

综上所述，该函数通过比较个体间的适应度值来确定支配关系，并将被支配的个体在 dom_vector 中标记为 1。这可以用于多目标优化问题中的非支配排序和粒子群优化等算法。