这段代码 % 变异操作 for i = 1 : population_size if rand() < mutation_rate % 对个体进行随机变异 mutated_individual = mutation(new_population(i, :)); % 将变异后的个体加入新种群中 new_population(i, :) = mutated_individual; end end如何封装为函数

可以将这段代码封装为一个名为 "mutation_population" 的函数，其输入参数为当前种群、变异率和变异操作函数，输出为变异后的新种群。以下是一个示例实现：

function new_population = mutation_population(population, mutation_rate, mutation_func)
    [population_size, genome_size] = size(population);
    new_population = population;
    for i = 1 : population_size
        if rand() < mutation_rate
            mutated_individual = mutation_func(population(i, :));
            new_population(i, :) = mutated_individual;
        end
    end
end

使用示例：

% 定义变异操作函数
function mutated_individual = my_mutation(individual)
    % 实现个体的变异操作
    % ...
    mutated_individual = individual;
end

% 调用 mutation_population 函数进行变异
old_population = ...; % 原始种群数据
mutation_rate = 0.01; % 变异率
new_population = mutation_population(old_population, mutation_rate, @my_mutation);