clc;clear;
close?all;
%-----確定實際問題參數集------
wmin=0;wmax=0.9;????????%慣性權重
itmax=50;???????????????%最大循環次數
c1=2;c2=2;??????????????%自身因子和全局因子
for?iter=1:itmax
????W（iter）=wmax-（（wmax-wmin）/itmax）*iter;%自適應權重
end
a=0;b=1;???????????????%變量取值范圍
N=30;??????????????????%粒子個數
D=1;???????????????????%維數
step=0.01;?????????????%步長
f=‘exp（-2*log（2）*（（（x-0.1）./0.8）.^2））.*（（sin（5*pi*x））.^6）‘;
x=[a:step:b];
xf=x;
y=eval（f）;
x=[];
vmfit=[];
x=a+（b-a）*rand（N1）;?????????????????%位置初始化
V=wmin+（wmin-wmax）*rand（N1）;????????%向量初始化
for?i=1:N
????F（i1）=exp（-2*log（2）*（（（x（i1）-0.1）/0.8）.^2））*（（sin（5*pi*x（i1）））.^6）;%第一代個體適應度計算
end
media=mean（F（:1））;
vmfit=[vmfit?media];
[CI]=max（abs（F（:1）））;??????????????
B（11）=C;
gbest（11）=x（I1）;??????????????????%群體最優位置
for?p=1:N
????G（p1）=gbest（11）;??????????????%更新群體最優位置記錄
end
Fbest（11）=exp（-2*log（2）*（（（G（11）-0.1）/0.8）.^2））*（（sin（5*pi*G（11）））.^6）;%群體最優位置適應度
for?i=1:N
????pbest（i1）=x（i1）;??????????????%個體最優位置，第一代群體設為當前位置
end
Fb（11）=exp（-2*log（2）*（（（pbest（11）-0.1）/0.8）.^2））*（（sin（5*pi*pbest（11）））.^6）;%個體最優位置適應度

%------程序主循環----------
for?j=2:itmax
????V（:j）=W（j-1）*V（:j-1）+c1*rand*（pbest（:j-1）-x（:j-1））;%+c2*rand*（G（:j-1）-x（:j-1））;%向量更新
????x（:j）=x（:j-1）+V（:j）;?????????%位置更新
????for?xx=1:N
????????if?x（xxj）????????????x（x