objvalue.m

一个遗传规划的源程序

function objvalue = objvalue(popu,S,X,Y);
% 计算每一个体的适应值(函数逼近真实值的程度)
% x为数据序列的号码
% y为瓦斯预测的真实值
% Q为函数计算得到的模型预测值

x = X;
y = Y;

% disp(S);
Q = eval(S);

% 采用最小二乘计算误差平房和最小确定较优模型结构
Error = sum((Q-y).^2);
objvalue = Error;

% 如果目标值不是数字或者是无穷大或者它不是实数而是带有虚部的复数时
if (isnan(objvalue))|(isinf(objvalue))|(~isreal(objvalue)),
%     disp('所生成的函数不是有效函数表达式！');
    objvalue = inf;
end
% disp(objvalue);
% disp('.............................................................');

return