%%?I.?清空環境
clc
clear
%%?II.?繪制目標函數曲線
figure
[xy]?=?meshgrid（-2:0.1:2-2:0.1:2）;
z?=?x*cos（2*pi*y）?+?y*sin（2*pi*x）;
mesh（xyz）
hold?on
%%?III.?參數初始化
c1?=?0.49445;
c2?=?0.49445;
maxgen?=?100;?%?進化次數
sizepop?=?100;?%種群規模
Vmax?=?0.5;
Vmin?=?-0.5;
popmax%%?IV.?產生初始粒子和速度
for?i?=?1:sizepop
????%?隨機產生一個種群
????pop（i:）?=?2*rands（12）;?%初始種群
????V（i:）?=?rands（12）;?%初始化速度
????%?計算適應度
????fitness（i）?=?fun2（pop（i:））;?%染色體的適應度
end
%%?V.?個體極值和群體極值
[bestfitness?bestindex]?=?max（fitness）;
zbest?=?pop（bestindex:）;?%全局最佳
gbest?=?pop;?%個體最佳
fitnessgbest?=?fitness;?%個體最佳適應度值
fitnesszbest?=?bestfitness;?%全局最佳適應度值
%%?VI.?迭代尋優
for?i?=?1:maxgen
????for?j?=?1:sizepop
????????%?速度更新
????????V（j:）?=?V（j:）?+?c1*rand*（gbest（j:）?-?pop（j:））?+?c2*rand*（zbest?-?pop（j:））;
????????V（jfind（V（j:）>Vmax））?=?Vmax;
????????V（jfind（V（j:）????????%?種群更新
????????pop（j:）?=?pop（j:）?+?V（j:）;
????????pop（jfind（pop（j:）>popmax））?=?popmax;
????????pop（jfind（pop（j:）????????%?適應度值更新
????????fitness（j）?=?fun2（pop（j:））;
????end
????for?j?=?1:sizepop
????????%?個體最優更新
????????if?fitness（j）?>?fitnessgbest（j）
????????????gbest（j:）?=?pop（j:）;
????????????fitnessgbest（j）?=?fitness（j）;
????????end
????????%?群體最優更新
????????if?fitness（j）?>?fitnesszbest
????????????zbest?=?pop（j:）;
????????????fitnesszbest?=?fitness（j）;
????????end
????end
????yy（i）?=?fitnesszbest;
end?=?2;
popmin?=?-2;

%%?VII.輸出結果
[fitnesszbest?zbest]
plot3（zbest（1）?zbest（2）?fitnesszbest‘bo‘‘linewidth‘1.5）
figure
plot（yy）
title（‘最優個體適應度‘‘fontsize‘12）;
xlabel（‘進化代數‘‘fontsize‘12）;ylabel（‘適應度‘‘fontsize‘12）;
function?y?=?fun2（x）
%函數用于計算粒子適應度值
%x???????????input???????????輸入粒子
%y???????????output??????????粒子適應度值
y?=?x（1）*cos（2*pi*x（2））?+?x（2）*sin（2*pi*x（1））;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-03-14?16:04??粒子群\
?????文件????????2000??2019-03-13?15:53??粒子群\liangge.m
?????文件????????2184??2019-01-15?11:04??粒子群\qiujie.m
?????文件???????34304??2019-03-14?16:04??粒子群\粒子群.doc