资源描述：
《基于遗传算法的PID整定原理及matlab仿真程序.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、主程序：chap5_2.m%GA(GenericAlgorithm)ProgramtooptimizePIDParametersclearall;closeall;globalrinyouttimefSize=30;CodeL=3;MinX(1)=zeros(1);MaxX(1)=20*ones(1);MinX(2)=zeros(1);MaxX(2)=1.0*ones(1);MinX(3)=zeros(1);MaxX(3)=1.0*ones(1);Kpid(:,1)=MinX(1)+(MaxX(1)-MinX(1))*ran
2、d(Size,1);Kpid(:,2)=MinX(2)+(MaxX(2)-MinX(2))*rand(Size,1);Kpid(:,3)=MinX(3)+(MaxX(3)-MinX(3))*rand(Size,1);G=100;BsJ=0;%***************StartRunning***************forkg=1:1:Gtime(kg)=kg;%******Step1:EvaluateBestJ******fori=1:1:SizeKpidi=Kpid(i,:);[Kpidi,BsJ]=chap5_
3、2f(Kpidi,BsJ);BsJi(i)=BsJ;end[OderJi,IndexJi]=sort(BsJi);BestJ(kg)=OderJi(1);BJ=BestJ(kg);Ji=BsJi+1e-10;%Avoidingdevidingzerofi=1./Ji;%Cm=max(Ji);%fi=Cm-Ji;[Oderfi,Indexfi]=sort(fi);%ArrangingfismalltobiggerBestfi=Oderfi(Size);%LetBestfi=max(fi)BestS=Kpid(Indexfi(S
4、ize),:);%LetBestS=E(m),mistheIndexfibelongtomax(fi)kgBJBestS%******Step2:SelectandReproductOperation******fi_sum=sum(fi);fi_Size=(Oderfi/fi_sum)*Size;fi_S=floor(fi_Size);%SelectingBiggerfivaluer=Size-sum(fi_S);Rest=fi_Size-fi_S;[RestValue,Index]=sort(Rest);fori=Siz
5、e:-1:Size-r+1fi_S(Index(i))=fi_S(Index(i))+1;%AddingresttoequalSizeendk=1;fori=Size:-1:1%SelecttheSizethandReproducefirstlyforj=1:1:fi_S(i)TempE(k,:)=Kpid(Indexfi(i),:);%SelectandReproducek=k+1;%kisusedtoreproduceendend%************Step3:CrossoverOperation*********
6、***Pc=0.90;fori=1:2:(Size-1)temp=rand;ifPc>temp%CrossoverConditionalfa=rand;TempE(i,:)=alfa*Kpid(i+1,:)+(1-alfa)*Kpid(i,:);TempE(i+1,:)=alfa*Kpid(i,:)+(1-alfa)*Kpid(i+1,:);endendTempE(Size,:)=BestS;Kpid=TempE;%************Step4:MutationOperation**************Pm=0.1
7、0-[1:1:Size]*(0.01)/Size;%Biggerfi,smallerPmPm_rand=rand(Size,CodeL);Mean=(MaxX+MinX)/2;Dif=(MaxX-MinX);fori=1:1:Sizeforj=1:1:CodeLifPm(i)>Pm_rand(i,j)%MutationConditionTempE(i,j)=Mean(j)+Dif(j)*(rand-0.5);endendend%GuaranteeTempE(Size,:)belongtothebestindividualTe
8、mpE(Size,:)=BestS;Kpid=TempE;endBestfiBestSBest_J=BestJ(G)figure(1);plot(time,BestJ);xlabel(’Times’);ylabel(’BestJ’);figure(2);plot(timef,rin,’r’
