function objvalue = objvalue(popu,S,X,Y);
% 計算每一個體的適應值(函數逼近真實值的程度)

x = X;
y = Y;

Q = eval(S);

% 采用最小二乘計算誤差平方和最小確定較優模型結構
Error = sum((Q-y).^2);
objvalue = Error;

% 如果目標值不是數字或者是無窮大或者它不是實數而是帶有虛部的復數時
if (isnan(objvalue))|(isinf(objvalue))|(~isreal(objvalue)),
    objvalue = inf;
end
return