function?[EdgeNumberofaera]=pcnn（X）
%?功能：采用PCNN算法進行邊緣檢測
%?輸入：X—輸入的灰度圖像
%?輸出：Edge—檢測到的???????Numberofaera—表明了在各次迭代時激活的塊區(qū)域
figure（1）;
imshow（X）;?
X=double（X）;
%?設定權值
Weight=[0.07?0.1?0.07;0.1?0?0.1;0.07?0.1?0.07];?
WeightLI2=[-0.03?-0.03?-0.03;-0.03?0?-0.03;-0.03?-0.03?-0.03];
d=1/（1+sum（sum（WeightLI2）））;
%%%%%%測試權值%%%%%%
WeightLI=[-0.03?-0.03?-0.03;-0.03?0.5?-0.03;-0.03?-0.03?-0.03];
d1=1/（sum（sum（WeightLI）））;
%%%%%%%%%%%%%%%%%%
Beta=0.4;??
Yuzhi=245;
%衰減系數(shù)
Decay=0.3;???
[ab]=size（X）;
V_T=0.2;
%門限值
Threshold=zeros（ab）;??
S=zeros（a+2b+2）;
Y=zeros（ab）;
%點火頻率
Firate=zeros（ab）;?
n=1;
%統(tǒng)計循環(huán)次數(shù)
count=0;?
Tempu1=zeros（ab）;?
Tempu2=zeros（a+2b+2）;?
%%%%%%圖像增強部分%%%%%%
Out=zeros（ab）;
Out=uint8（Out）;
for?i=1:a
for?j=1:b
?if（i==1|j==1|i==a|j==b）
??Out（ij）=X（ij）;
?else??
??H=[X（i-1j-1）??X（i-1j）?X（i-1j+1）;
?????X（ij-1）???X（ij）???X（ij+1）;
????X（i+1j-1）?X（i+1j）?X（i+1j+1）];?
?temp=d1*sum（sum（H.*WeightLI））;
?Out（ij）=temp;
?end
?end
end
figure（2）;
imshow（Out）;?
%%%%%%%%%%%%%%%%%%%
for?count=1:30?
?for?i0=2:a+1
????for?i1=2:b+1
?????????V=[S（i0-1i1-1）??S（i0-1i1）?S（i0-1i1+1）;
?????????????S（i0i1-1）???S（i0i1）???S（i0i1+1）;
?????????????S（i0+1i1-1）?S（i0+1i1）?S（i0+1i1+1）];
???????????L=sum（sum（V.*Weight））;
???????????V2=[Tempu2（i0-1i1-1）??Tempu2（i0-1i1）?Tempu2（i0-1i1+1）;
???????????????Tempu2（i0i1-1）???Tempu2（i0i1）???Tempu2（i0i1+1）;
???????????????Tempu2（i0+1i1-1）?Tempu2（i0+1i1）?Tempu2（i0+1i1+1）];????????F=X（i0-1i1-1）+sum（sum（V2.*WeightLI2））;
%保證側(cè)抑制圖像無能量損失
F=d*F;?
U=double（F）*（1+Beta*double（L））;??????????????????????????
Tempu1（i0-1i1-1）=U;
????if?U>=Threshold（i0-1i1-1）|Threshold（i0-1i1-1）<60
??????T（i0-1i1-1）=1;
??????Threshold（i0-1i1-1）=Yuzhi;
???????%點火后一直置為1
Y（i0-1i1-1）=1;????
?????else
????????T（i0-1i1-1）=0;
????????Y（i0-1i1-1）=0;
?????????????????end
????????????end
?????????end
???Threshold=exp（-Decay）*Threshold+V_T*Y;
???%被激活過的像素不再參與迭代過程
?????if?n==1
????????S=zeros（a+2b+2）;
????????else
????????S=Bianhuan（T）;
?????end
?????n=n+1;
?????count=count+1;?
?????Firate=Firate+Y;
????figure（3）;
????imshow（Y）;
????Tempu2=Bianhuan（Tempu1）;
end
???Firate（find（Firate<10））=0;
???Firate（find（Firate>=10））=10;
???figure（4）;
???imshow（Firate）;
%%%%%%子函數(shù)?%%%%%%%
function?Y=Jiabian（X）
[mn]=size（X）;
Y=zeros（m+2n+2）;
for?i=1:m+2
????for?j=1:n+2
????????if?i==1&j~=1&j~=n+2
????????????Y（ij）=X（1j-1）;
????????elseif?j==1&i~=1&i~=m+2
????????????Y（ij）=X（i-11）;
????????elseif?i~=1&j==n+2&i~=m+2
????????????Y（ij）=X（i-1n）;
????????elseif?i==m+2&j~=1&j~=n+2
????????????Y（ij）=X（mj-1）;
????????elseif?i==1&j==1
????????????Y（ij）=X（ij）;
????????elseif?i==1&j==n+2
????????????Y（ij）=X（1n）;
????????elseif?i==（m+2）&j==1
????????????Y（ij）=X（m1）;
????????elseif?i==m+2&j==n+2
????????????Y（ij）=X（mn）;
????????else
????????????Y（ij）=X（i-1j-1）;
????????end
????end
end
%%%%%%子函數(shù)%%%%%%
function?Y=Bianhuan（X）
[mn]=size（X

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

?????文件??????14059??2011-06-06?21:07??3.25?基于脈沖耦合神經(jīng)網(wǎng)絡的圖像分割\lena.JPG

?????文件???????3695??2012-02-15?23:28??3.25?基于脈沖耦合神經(jīng)網(wǎng)絡的圖像分割\pcnn.m

?????文件??????24064??2011-09-26?22:47??3.25?基于脈沖耦合神經(jīng)網(wǎng)絡的圖像分割\程序運行說明.doc

?????目錄??????????0??2019-06-18?09:26??3.25?基于脈沖耦合神經(jīng)網(wǎng)絡的圖像分割

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

????????????????41818????????????????????4