function?[gbestxmfvMWn]=yueshu_APSO_FCG_norm（testNc1c2MwmaxwminxmaxxminymaxyminNmaxNmin）
%有約束的基于適應(yīng)度和連續(xù)世代策略的自適應(yīng)粒子群算法：調(diào)節(jié)w、N初始化為范圍內(nèi)隨機(jī)
%test函數(shù)為增廣目標(biāo)函數(shù)
%停止條件為：最大迭代次數(shù)和連續(xù)9次最優(yōu)適應(yīng)度沒(méi)有更新
%全局最優(yōu)增廣適應(yīng)度：gbest
%目標(biāo)函數(shù)的最小值的自變量值：xm
%目標(biāo)函數(shù)的最小值：fv
%最大迭代次數(shù)：M
%粒子數(shù)目動(dòng)態(tài)軌跡：n
%test增廣適應(yīng)度函數(shù)
%粒子數(shù)目初始值：N
%學(xué)習(xí)因子1：c1
%學(xué)習(xí)因子2：c2
%最大權(quán)重：wmax
%最小權(quán)重：wmin
%優(yōu)化變量范圍：wmax、wmin、ymaxymin

format?long;
%%?在粒子范圍內(nèi)隨機(jī)初始化
x1=10*rand（N1）;
x2=10*rand（N1）;
x=[x1x2];
v=rand（N2）;
W=zeros（MN）;%存放慣性因子，方便繪圖
n=zeros（M1）;%存放種群規(guī)模
gbest=zeros（1M）;%預(yù)定義全局最優(yōu)適應(yīng)度
%%?求解有約束的初始化粒子的歷史最優(yōu)適應(yīng)度和位置
for?i=1:N
????f（i）=test（x（i:）0xmaxxminymaxymin）;%求解增廣適應(yīng)度
????p（i:）=x（i:）;%歷史最優(yōu)
end
%比較適應(yīng)度
pg=x（N:）;?????%pg為全局最優(yōu)
a=N;
for?i=1:（N-1）
????if?f（i）????????pg=x（i:）;
????????v1=v（i:）;%存放歷史最優(yōu)位置的速度
????????a=i;
????end
end

%%?粒子群算法更新
for?t=1:M?
????n（t）=N;
????for?j=1:N?
????????fv（j）=test（x（j:）txmaxxminymaxymin）;??%計(jì)算增廣目標(biāo)函數(shù)
????end??????????????
%%?自適應(yīng)慣性因子???????????
????fvag=sum（fv）/N;?????????%適應(yīng)度平均值
????fmin=min（fv）;???????????%適應(yīng)度最小值
????for?i=1:N
????????if?fv（i）????????????w=wmin+（fv（i）-fmin）*（wmax-wmin）/（fvag-fmin）;%減少w
????????else
????????????w=wmax;%適應(yīng)度差，增強(qiáng)全局
????????end
????????W（ti）=w;%慣性因子
%%?更新位置和速度?????????????????
????????v（i:）=w*v（i:）+c1*rand*（p（i:）-x（i:））+c2*rand*（pg-x（i:））;
????????x（i:）=x（i:）+v（i:）;??
%%?更新有約束的歷史最優(yōu)和全局最優(yōu)????????
????????if?test（x（i:）txmaxxminymaxymin）????????????f（i）=test（x（i:）txmaxxminymaxymin）;%歷史最優(yōu)適應(yīng)度
????????????p（i:）=x（i:）;%歷史最優(yōu)位置
????????????a=i;
????????end
????????if?f（i）????????????pg=p（i:）;
????????????v1=v（i:）;%全局最優(yōu)位置的速度
????????end?????????
????end
????gbest（t）=test（pgtxmaxxminymaxymin）;%全局最優(yōu)適應(yīng)度
%%?連續(xù)世代調(diào)整種群規(guī)模
????if?t>=4?????%更新4代及其以上時(shí)候使用連續(xù)世代策略（因?yàn)樵O(shè)定連續(xù)三代，所以從第四代開(kāi)始統(tǒng)計(jì)
????????if?gbest（t）==gbest（t-1）&&gbest（t）