%GA?Program?tooptimize?PID?Parameters
%基于GA算法的PID參數(shù)優(yōu)化
clc?%?清屏
clear?all;?%?刪除workplace變量
close?all;?%?關(guān)掉顯示圖形窗口

size=30;
codel=3;
minx（1）=zeros（1）;
maxx（1）=50*ones（1）;
minx（2）=zeros（1）;
maxx（2）=50*ones（1）;
minx（3）=zeros（1）;
maxx（3）=50*ones（1）;
kpid（:1）=minx（1）+（maxx（1）-minx（1））*rand（size1）;
kpid（:2）=minx（2）+（maxx（2）-minx（2））*rand（size1）;
kpid（:3）=minx（3）+（maxx（3）-minx（3））*rand（size1）;
G=100;
BsJ=0;
for?kg=1:1:G
????time（kg）=kg;
????for?i=1:1:size
????????kpidi=kpid（i:）;
????????BsJ=pid_GA（kpidi）;
????????BsJi（i）=BsJ;
????end
[OderJiIndexJi]=sort（BsJi）;
BestJ（kg）=OderJi（1）;
BJ=BestJ（kg）;
Ji=BsJi+1e-10;
fi=1./Ji;
[OderfiIndexfi]=sort（fi）;
Bestfi=Oderfi（size）;
BestS=kpid（Indexfi（size）:）;
%?select?and?reproduct?operation
fi_sum=sum（fi）;
fi_size=（Oderfi/fi_sum）*size;
fi_s=floor（fi_size）;
r=size-sum（fi_s）;
Rest=fi_size-fi_s;
[RestvalueIndex]=sort（Rest）;
for?i=size:-1:size-r+1
????fi_s（Index（i））=fi_s（Index（i））+1;
end
k=1;
for?i=size:-1:1
????for?j=1:1:fi_s（i）
????????TempE（k:）=kpid（Indexfi（i）:）;
????????k=k+1;
????end
end
%?crossover?operation
Pc=0.90;
for?i=1:2:（size-1）
????temp=rand;
????if?Pc>temp
????????alfa=rand;
????????TempE（i:）=alfa*kpid（i+1:）+（1-alfa）*kpid（i:）;
????????TempE（i+1:）=alfa*kpid（i:）+（1-alfa）*kpid（i+1:）;
????end
end
TempE（size:）=BestS;
kpid=TempE;
%Mution?operation
Pm=0.1-[1:1:size]*（0.01）/size;
Pm_rand=rand（sizecodel）;
Mean=（maxx+minx）/2;
Dif=（maxx-minx）;
for?i=1:1:size?
????for?j=1:1:codel
????????if?Pm（i）>Pm_rand（ij）
????????????TempE（ij）=Mean（j）+Dif（j）*（rand-0.5）;
????????end
????end
end
%guarantee?TempE（size:）belong?tu?the?best?individual
TempE（size:）=BestS;
kpid=TempE;
end
Bestfi;
BestS;
Best_J=BestJ（G）;
figure（1）;
plot（timeBestJ‘LineWidth‘3）;
title（‘最優(yōu)個(gè)體適應(yīng)值‘‘fontsize‘10）;
xlabel（‘迭代次數(shù)‘）;ylabel（‘適應(yīng)值‘）;
grid?on


????????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1929??2014-05-04?12:50??GA\GA.m

?????文件???????1134??2014-05-04?12:50??GA\pid_GA.m

?????文件????????531??2014-05-04?14:52??GA\shuchuquxianGA.m

?????文件??????16565??2014-05-04?14:56??GA\ysw_PID1.slx

?????目錄??????????0??2014-06-22?22:04??GA

-----------?---------??----------?-----??----

????????????????20159????????????????????5