%維數n=3
%設置參數邊界
bounds =[1 4;1 7;1 7];
%遺傳算法優化
%二進制編碼
options=[1e-6 2];
%initPop=initializega(10,bounds,'f553',[],options);
%[m endPop] = ga(bounds,'f553',[],initPop,[1e-6 1 1],'maxGenTerm',100,'normGeomSelect',…
[0.08],['arithXover'],[2], 'nonUnifMutation',[2 5 3]);

[p,endPop,bestSols,trace]=ga(bounds,'f553');

%性能跟蹤
plot(trace(:,1),trace(:,3),'b-');
hold on
plot(trace(:,1),trace(:,2),'r-');
xlabel('Generation');
ylabel('Fittness');
legend('解的變化','種群平均值的變化');