原始程序
%主程序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.表示运行结束!