原始程序
%主程序
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 放进函数文件里面,如下