Matlab 自带遗传算法函数 ga() 运行报错:“ fitnessfcn() 输入参数的数目不足。”

原始程序

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

 

 

©️2020 CSDN 皮肤主题: 书香水墨 设计师:CSDN官方博客 返回首页