原始程序

%主程序
m = [1  2];
R = [1.5   1];  
epsilon = [2^R(1)-1 2^R(2) - 1];
A = [-1 epsilon(1); 0 -1];
b = [0 0]';
Aeq = [1 1];
Beq = 1;%matlab自带函数包
[x_optimal, fval] = ga(@fitnessfcn, 2, A, b, Aeq, Beq, [0;0], [1;1], []);  --------------------------------------------------------------------------------
%自定义适应度函数
%变量x为二维向量
function y = fitnessfcn(x, epsilon, m, R)eta_1_ga = epsilon(1)/(x(1) - epsilon(1)*x(2));eta_2_ga = max([eta_1_ga,epsilon(2)/x(2)]);P_1_ga = gammainc(eta_1_ga,m(1)); P_2_ga = gammainc(eta_2_ga,m(2));y = R(1) * P_1_ga + R(2) * P_2_ga;   end

运行报错

fitnessfcn() 输入参数的数目不足。

 

解决方法

在 ga() 中适应度函数 fitnessfcn() 是以函数句柄被调用的，变量 nvars = 2，也即输入参数只能为两个变量，这里即二维向量 x，所以参数 epsilon, m, R 不能一起放在这里当作输入变量。 可以把 epsilon, m, R 放进函数文件里面，如下

%自定义适应度函数
%变量x为二维向量
function y = fitnessfcn(x)  %===== 删除变量 epsilon, m, R======= 增添以下3个参数赋值 ======m = [1  2];  R = [1.5   1];  epsilon = [2^R(1)-1 2^R(2) - 1];=================================eta_1_ga = epsilon(1)/(x(1) - epsilon(1)*x(2));eta_2_ga = max([eta_1_ga,epsilon(2)/x(2)]);P_1_ga = gammainc(eta_1_ga,m(1)); P_2_ga = gammainc(eta_2_ga,m(2));y = R(1) * P_1_ga + R(2) * P_2_ga;   end

运行结果

命令窗口显示：Optimization terminated: minimum fitness limit reached. 表示运行结束！