clc;
clear;
%C=imread（‘lena.bmp‘）;
%C=imread（‘stone.jpg‘）;
C=imread（‘1.bmp‘）;
if?isrgb（C）
????I=rgb2gray（C）;
else
????I=C;
end
%I=imnoise（I‘gaussian‘00.005）;
%figure
subplot（131）
imshow（I）
tic
IDX?=?otsu（I2）;
subplot（132）;
imagesc（IDX）?axis?image?off?
colormap（gray）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%確定基本粒子群算法的迭代次數(shù)和粒子個(gè)數(shù)
Iter?=?50;
global??popsize
popsize?=?20;??%%粒子個(gè)數(shù)
global?dimsize
dimsize=2;??%%粒子維數(shù)
%初始化基本粒子群優(yōu)化算法的參數(shù)
Wmax?=?0.9;?????????????%最大慣性因子
Wmin?=?0.1;?????????????%最小慣性因子
c1?=?2.0;???????????????%個(gè)體學(xué)習(xí)因子
c2?=?2.0;???????????????%群體學(xué)習(xí)因子
xmax?=?254;?????????????%粒子允許的最小位置
xmin?=?2;?
speedmax=10;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%初始化
????????for?k1=1:popsize
???????????for?i=1:popsize
??????????????%Xpop（i1:dimsize）=?unifrnd（xminxmax1dimsize）;
?????????????????%pop（i1:dimsize）=[?unifrnd（xmin85）?unifrnd（85170）unifrnd（170255）];???%%三閾值
??????????????????Xpop（i1:dimsize）=unifrnd（22541dimsize）;???%%兩閾值
???????????end
???????????for?j=1:popsize
???????????????Xpop（jdimsize+1:2*dimsize）=?unifrnd（-speedmaxspeedmax1dimsize）;
???????????end
???????????pop（k11:2*dimsize?）=?Xpop（k11:2*dimsize）;
????????end
????????
????????BestFit?=?0;
????????pBestFit=zeros（popsize1）;?%%%%%%%%
????????%初始化粒子歷史最優(yōu)位置pBest和全局最優(yōu)位置gBest
????????pBest?=?zeros（popsizedimsize）;
????????gBest?=?zeros（1dimsize）;
????????
????????%迭代更新
????????for?it?=?1:Iter
????????????for?t=1:popsize
????????????????%計(jì)算適應(yīng)度pFit并更新歷史最優(yōu)適應(yīng)度pBestFit和全局最佳適應(yīng)度gBestFit
????????????????Theta?=?pop（t1:dimsize）;

????????????????%當(dāng)前粒子的適應(yīng)度pFit
????????????????pFit?=?ThetaFit（ThetaI）;

????????????????%更新個(gè)體歷史最優(yōu)適應(yīng)度pBestFit
????????????????if?pFit?>?pBestFit（t）
????????????????????pBestFit（t）?=?pFit;
????????????????????pBest（t1:dimsize）?=?pop（t1:dimsize）;
????????????????end
????????end
????????%更新當(dāng)前時(shí)刻全局最優(yōu)位置gBest、前一時(shí)刻全局最優(yōu)適應(yīng)度PregBest
????????%和全局最佳適應(yīng)度gBestFit
????????[gBestFitind]?=?max（pBestFit）;???%%pBestFit（j:）表示第j個(gè)子群的適應(yīng)度值，ind（j）表示位置
????????gBest（1:dimsize）?=?pBest（ind1:dimsize）;??
????????%位置向量和速度向量更新
????????for?j=1:popsize
????????????????for?t?=?dimsize+1:2*dimsize
????????????????????%慣性權(quán)重更新并產(chǎn)生隨機(jī)數(shù)r1、r2
????????????????????W?=?Wmax?-?it*（Wmax?-?Wmin）/Iter;
????????????????????r1?=?rand;
????????????????????r2?=?rand;

????????????????????%速度向量Vpos更新
????????????????????pop（jt）?=?W*pop（jt）?+?c1*r1*（pBest（jt-dimsize）?-?pop（jt-dimsize））?+?...
????????????????????????????c2*r2*（gBest（t-dimsize）?-?pop（jt-dimsize））;
????????????????????%限制更新后粒子速度在[VminVmax]范圍內(nèi)
????????????????????if?pop（jt）?>?speedmax
????????????????????????pop（jt）?=?speedmax;
????????????????????end
???????????????????if?pop（jt）?????????????????????????pop（jt）?=?-speedmax;
???????????????????end

????????????????????%粒子位置更新
????????????????????pop（jt-dimsize）?=?pop（jt-dimsize）?+?pop（jt）;
???

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-10?15:39??源碼\
?????文件???????56910??2016-04-27?00:29??源碼\1.bmp
?????文件???????66614??2016-04-27?00:08??源碼\lena.bmp
?????文件????????3936??2016-04-27?00:41??源碼\main.m
?????文件????????1888??2016-04-27?12:31??源碼\main_Tsallis_2D.m
?????文件???????76885??2012-10-30?15:53??源碼\stone.jpg
?????文件??????889611??2016-04-27?12:34??源碼\畢業(yè)論文?（修復(fù)的）.docx
?????目錄???????????0??2018-10-10?15:39??演示\
?????文件?????5437006??2016-04-13?13:32??演示\9780.avi
?????文件??????????64??2018-10-10?16:34??說(shuō)明.txt
?????目錄???????????0??2018-10-10?15:38??需求\
?????文件???????13002??2016-04-04?21:27??需求\新建?Microsoft?Word?文檔.docx