%%?清空環(huán)境
clc
clear

%%?參數(shù)初始化
%粒子群算法中的兩個參數(shù)
c1?=?1.49445;
c2?=?1.49445;
wmax=0.9;%慣性因子最大值
wmin=0.4;%慣性因子最小值
D=21;%粒子維數(shù)
maxgen=1000;???%?進化次數(shù)??
sizepop=50;???%種群規(guī)模
u=3;%混沌系數(shù)
Vmax=0.3;
Vmin=-0.3;
popmax=1;
popmin=-1;

%%?產(chǎn)生初始粒子和速度
for?i=1:sizepop
????%隨機產(chǎn)生一個種群
????pop（i:）=rands（1D）;????%初始種群
????V（i:）=rands（1D）;??%初始化速度
????%計算適應度
????fitness（i）=fun5（pop（i:））;???%粒子的適應值
end

%%?個體極值和群體極值
[bestfitness?bestindex]=min（fitness）;
zbest=pop（bestindex:）;???%全局最佳
gbest=pop;????%個體最佳
fitnessgbest=fitness;???%個體最佳適應度值
fitnesszbest=bestfitness;???%全局最佳適應度值

%%?迭代尋優(yōu)
for?i=1:maxgen
????????w=wmax-i*（wmax-wmin）./maxgen;
????for?j=1:sizepop
????????
????????%速度更新
????????V（j:）?=?w*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）=fun5（pop（j:））;?
???
????end
????
????for?j=1:sizepop
????????
????????%個體最優(yōu)更新
????????if?fitness（j）?????????????gbest（j:）?=?pop（j:）;
????????????fitnessgbest（j）?=?fitness（j）;
????????end
????????
????????%群體最優(yōu)更新
????????if?fitness（j）?????????????zbest?=?pop（j:）;
????????????fitnesszbest?=?fitness（j）;
????????end
????end?
??
????%%對粒子群最優(yōu)位置進行混沌優(yōu)化
??????y（1:）=（zbest-popmin）/（popmax-popmin）;%將最優(yōu)位置映射到Logistic方程的定義域[01]
??????fitness（1）=fun5（y（1:））;?
????????for?t=1:sizepop-1?%通過Logistic方程進行M次迭代，得到混沌序列
????????????for?e=1:D
????????y（t+1e）=u*y（te）*（1-y（te））;?
????????????end
????????y（t+1:）=popmin+（popmax-popmin）*y（t+1:）;%將混沌序列逆射到原解空間
????????fitness（t+1）=fun5（y（t+1:））;?%計算混沌變量可行解序列的適應度值
????????end
[ybestfitness?ybestindex]=min（fitness）;%尋找最優(yōu)混沌可行解矢量
??ybest=y（ybestindex:）;
????????ran=1+fix（rand（）*sizepop）;%產(chǎn)生一隨機數(shù)1~sizepop之間
????????pop（ran:）=ybest;
????yy（i）=fitnesszbest??
????????
end
%%?結果分析
plot（yy‘m‘）
title（‘最優(yōu)個體適應度曲線‘‘fontsize‘15）;
xlabel（‘迭代次數(shù)‘‘fontsize‘15）;ylabel（‘適應度值‘‘fontsize‘15）;
hold?on
display（‘混沌粒子群算法輸出結果‘）;
zbest
minbest=min（yy）
meanbest=mean（yy）
stdbest=std（yy）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3931??2012-05-17?19:25??FIR濾波器設計\autid1.fig

?????文件??????78330??2012-05-17?19:25??FIR濾波器設計\autid1.jpg

?????文件???????2522??2012-05-08?21:17??FIR濾波器設計\CPSO.m

?????文件???????8061??2012-05-08?21:30??FIR濾波器設計\FIR.fig

?????文件??????51106??2012-05-08?21:30??FIR濾波器設計\FIR.jpg

?????文件???????4328??2012-05-08?21:43??FIR濾波器設計\FIR2.fig

?????文件??????65332??2012-05-08?21:43??FIR濾波器設計\FIR2.jpg

?????文件???????4930??2012-05-17?19:19??FIR濾波器設計\FIR3.fig

?????文件??????81586??2012-05-17?19:19??FIR濾波器設計\FIR3.jpg

?????文件????????422??2012-05-09?16:31??FIR濾波器設計\firplot.m

?????文件???????1017??2012-05-09?16:38??FIR濾波器設計\firplot2.m

?????文件????????439??2012-05-09?16:36??FIR濾波器設計\firplot3.m

?????文件????????364??2012-04-26?21:06??FIR濾波器設計\fun5.m

?????文件???????3943??2012-05-08?21:44??FIR濾波器設計\GAPSO.m

?????文件????????460??2012-05-09?15:05??FIR濾波器設計\plot1.m

?????文件??????30214??2012-05-17?19:28??FIR濾波器設計\single1.fig

?????文件?????120172??2012-05-17?19:28??FIR濾波器設計\single1.jpg

?????文件??????86382??2012-05-17?18:55??FIR濾波器設計\single2.jpg

?????目錄??????????0??2012-10-06?19:42??FIR濾波器設計

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

???????????????543539????????????????????19