該程序為標準遺傳算法優化主函數，染色體為整數編碼，供初學者學習。

clear;

clc;

data=xlsread('data.xlsx');

maxgen=400;

sizepop=60;

pcross=0.7;

pmutation=0.2;

lenchrom=5;

bound=[1 2;1 5;1 3;1 3;1 4];

individuals=struct('fitness',zeros(1,sizepop),'chrom',[]);

for i=1:sizepop

    individuals.chrom(i,:)=ceil(rand(5,1).*(bound(:,2)-bound(:,1))+bound(:,1));

    x=individuals.chrom(i,:);

    individuals.fitness(i)=my_func(x,data);

end

[bestfitness,bestindex]=max(individuals.fitness);

bestchrom=individuals.chrom(bestindex,:);

trace=[];

for i=1:maxgen

    

    individuals=Select(individuals,sizepop);

    individuals.chrom=Cross(pcross,lenchrom,individuals.chrom,sizepop);

    individuals.chrom=Mutation(pmutation,lenchrom,individuals.chrom,sizepop,bound);

    for j=1:sizepop

        x=individuals.chrom(j,:);

        individuals.fitness(j)=my_func(x,data);

    end

    

    [newbestfitness,newbestindex]=max(individuals.fitness);

    

    if bestfitness<newbestfitness

        bestfitness=newbestfitness;

        bestchrom=individuals.chrom(newbestindex,:);

    end

    trace=[trace; bestfitness];

    disp(['**********第',num2str(i),'代尋優**********']);

end

disp(['最優結果:',num2str(bestchrom)]);

disp(['最優適應度:',num2str(bestfitness)]);

figure;

plot(trace);

grid on;

xlabel('迭代次數');

ylabel('適應度');

title('迭代結果');