以上程序运行显示错误使用 gamultiobj (line 248) Fitness function must be a function handle.

抱歉，我之前的回答有误。您的问题确实是适应度函数没有正确指定函数句柄。请将以下行:

problem.objective = @objectives;

改为:

problem.fitnessfcn = @objectives;

即可解决该问题。完整的代码如下：

% 定义目标函数
function [f1, f2] = objectives(x)
    f1 = x(1);
    f2 = (1 + x(2)) / x(1);
end

% 定义优化问题
problem.fitnessfcn = @objectives;
problem.x0 = [1, 1];
problem.lb = [0, 0];
problem.ub = [Inf, Inf];

% 使用NSGA-II算法求解Pareto前沿
options = optimoptions(@gamultiobj, 'PopulationSize', 100, 'MaxGenerations', 200);
[result, fval] = gamultiobj(problem, 2, [], [], [], [], problem.lb, problem.ub, options);

% 绘制Pareto前沿
figure;
plot(fval(:,1), fval(:,2), 'o');
xlabel('f1');
ylabel('f2');
title('Pareto Front');

请注意，我们将原来的`problem.objective`改为了`problem.fitnessfcn`来指定适应度函数的函数句柄。