投影寻踪遗传算法matlab程序.docxVIP

clc clear close all Alpha=0.1; X=load( d:\data.txt ); for k=1:21 eval(spri ntf( syms a%d ,k)); end a=[a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,a15,a16,a17,a18,a19,a 20,a21]; [n ,p]=size(X); x=zeros( n,p); Xjmax=max(X); Xjmi n=mi n( X); for i=1: n x(i,:)=(X(i,:)-Xjmi n)./(Xjmax-Xjmi n); end %构造投影指标 for i=1: n Z(i)=sum(a.*x(i,:)); end % 计算投影指标函数 mean Z=mea n( Z); Sa=0; for k=1: n sa=(Z(i)-mea nZ)42; Sa=Sa+sa; Sa=sqrt(Sa/n); end R=Alpha*Sa; Da=0; for k=1: n rik=abs(Z(i)-Z(k)); %if Rrik Da=Da+rik; %e nd end Qa=Sa*Da; %选择操作 fun cti on ret =select(i ndividuals,sizepop) in dividuals.fit ness=1./(i ndividuals.fit ness); sumfit ness=sum(i ndividuals.fit ness); sumf= in dividuals./fit ness; in dex=[]; for i =1:sizepop pick=ra nd; while pick==0 pick=ra nd; end for j =1:sizepop pick =pick-sumf(j); if pick0 in dex=[i ndex j]; break ; end end end in divlduals.chrom=in divuduals.chrom(i ndex,:); in dividuals.fit ness=in dividuals.fit ness(i ndex); ret=i ndividuals; %交叉操作 fun cti on ret=Cross(pcross,le nchrom,chrom,sizepop,bou nd) for i=1:sizepop pick = ran d(1,2); while prod(pick)==0 pick=ra nd(1,2); end in dex=ceil(pick.*sizepop); pick=ra nd; while pick==0 pick=ra nd; end if pickpcross con ti nue ; end flag=0; while flag==O pick=ra nd; while pick==0 pick=ra nd; end pos=ceil(pick.*sum(le nchrom)); pick=ra nd; v1=chrom(i ndex(1),pos); v2=chrom(i ndex(2),pos); chrom(i ndex(1),pos)=pick*v2+(1-pick)*v1; chrom(i ndex(2),pos)=pick*v1+(1-pick)*v2; flag仁test(le nchrom,bou nd,chrom(i ndex(1),:),fcode); flag2=test(le nchrom,bou nd,chrom(i ndex(2),:),fcode); if flag1*flag2==0 flag=0; else flag=1; end end end ret= chrom; %变异操作 fun cti on ret=Mutati on( pmutati on ,le nchrom,chrom,sizepoppop,bou nd) pick=ra nd; while pick==0 pick=ra nd; end in dex=ceil(pick*sizepop); pick=ra nd; if pickpmutati on con ti nue ; end flag=0; while flag==0 pick =rand; while pick==0 pick=ra nd; end pos=ceil(pick*sum(le nchrom)); v=chrom(i,pos); v1=v-bou nd(pos,1); v2=bou nd(pos,2)-v; pick=ra nd; if pick0.5 d