clc;clear?all;close?all;
t1=clock;
I=imread（‘a(chǎn).jpg‘）;
%J=rgb2gray（I）;%如果是彩色圖像要加此句，但去掉下面一句
?J=I;
[ab]=size（J）;
figure（1）
imshow（J）
[px]=imhist（J256）;
figure（2）
plot（p）;
title（‘灰度直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘像素點個數(shù)‘）
L=x‘;
LP=p‘/（a*b）;

n=256
c1=2;
c2=2;
wmax=0.9;
wmin=0.4;
G=100;
M=10;
X=min（L）+fix（（max（L）-min（L））*rand（1M））;
V=min（L）+（max（L）-min（L））*rand（1M）;


m=0;
for?i=1:1:n
????m=m+L（i）*LP（i）;
end
pbest=zeros（M2）;
gbest1=0;
gbest2=0;
GG=0;
t2=clock;
for?k=1:1:G
????w（k）=wmax-（wmax-wmin）*k/G;
????for?i=1:1:M
????????t=length（find（X（i）>=L））;
????????r=0;
????????s=0;
????????for?j=1:1:t
????????????r=r+LP（j）;
????????????s=s+L（j）*LP（j）;
????????end
????????W0（i）=r;
????????W1（i）=1-r;
????????U0（i）=s/r;
????????U1（i）=（m-s）/（1-r）;
????end
????for?i0=1:1:M
????????BB（i0）=W0（i0）*W1（i0）*（（U1（i0）-U0（i0））^2）;
????end
??
???for?i=1:1:M
???????if?pbest（i2）???????????pbest（i2）=BB（i）;????????????
???????????pbest（i1）=X（i）;

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

?????文件??????21034??2012-11-30?10:51??圖像分割1\a.jpg

?????文件???????4170??2012-05-18?11:14??圖像分割1\pic1.fig

?????文件??????65155??2012-05-18?11:14??圖像分割1\pic1.jpg

?????文件???????4284??2012-05-18?11:20??圖像分割1\pic2.fig

?????文件??????74224??2012-05-18?11:20??圖像分割1\pic2.jpg

?????文件???????4497??2012-05-18?11:30??圖像分割1\pic3.fig

?????文件??????62023??2012-05-18?11:30??圖像分割1\pic3.jpg

?????文件???????9096??2012-05-18?10:47??圖像分割1\picture1.jpg

?????文件?????114592??2012-05-18?10:49??圖像分割1\picture2.jpg

?????文件?????165880??2012-05-18?10:59??圖像分割1\picture3.jpg

?????文件???????1593??2012-11-30?14:49??圖像分割1\pso.m

?????目錄??????????0??2012-12-10?21:15??圖像分割1

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

???????????????526548????????????????????12