clcclearclose allAlpha=0.1;X=load('d:\data.txt');for k=1:21 eval(sprintf('syms a%d',k));enda=[a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,a15,a16,a17,a18,a19,a20,a21];[n,p]=size(X);x=zeros(n,p);Xjmax=max(X);Xjmin=min(X);for i=1:n x(i,:)=(X(i,:)-Xjmin)./(Xjmax-Xjmin);end%构造投影指标for i=1:n Z(i)=sum(a.*x(i,:));end%计算投影指标函数meanZ=mean(Z);Sa=0;for k=1:n sa=(Z(i)-meanZ).^2; Sa=Sa+sa; Sa=sqrt(Sa/n);endR=Alpha*Sa;Da=0;for k=1:n rik=abs(Z(i)-Z(k));%if R>rik Da=Da+rik;%endendQa=Sa*Da; %选择操作function ret =select(individuals,sizepop)individuals.fitness=1./(individuals.fitness);sumfitness=sum(individuals.fitness);sumf=individuals./fitness;index=[];for i =1:sizepop pick=rand;while pick==0 pick=rand;endfor j =1:sizepop pick =pick-sumf(j);if pick<0 index=[index j];break;endendendindivlduals.chrom=indivuduals.chrom(index,:);individuals.fitness=individuals.fitness(index);ret=individuals;%交叉操作function ret=Cross(pcross,lenchrom,chrom,sizepop,bound)for i=1:sizepop pick = rand(1,2);while prod(pick)==0 pick=rand(1,2);end index=ceil(pick.*sizepop); pick=rand;while pick==0 pick=rand;endif pick>pcrosscontinue;end flag=0;while flag==0 pick=rand;while pick==0 pick=rand;end pos=ceil(pick.*sum(lenchrom)); pick=rand; v1=chrom(index(1),pos); v2=chrom(index(2),pos); chrom(index(1),pos)=pick*v2+(1-pick)*v1; chrom(index(2),pos)=pick*v1+(1-pick)*v2; flag1=test(lenchrom,bound,chrom(index(1),:),fcode); flag2=test(lenchrom,bound,chrom(index(2),:),fcode);if flag1*flag2==0 flag=0;else flag=1;endendendret= chrom;% 变异操作function ret=Mutation(pmutation,lenchrom,chrom,sizepoppop,bound)pick=rand;while pick==0 pick=rand;endindex=ceil(pick*sizepop);pick=rand;if pick>pmutationcontinue;endflag=0;while flag==0 pick =rand;while pick==0 pick=rand;end pos=ceil(pick*sum(lenchrom)); v=chrom(i,pos); v1=v-bound(pos,1); v2=bound(pos,2)-v; pick=rand;if pick>0.5 delta=v2*(1-pick^((1-pop(1)/pop(2))^2)); chrom(i,pos)=v+delta...
