%讀取所有訓練數(shù)據(jù)與測試數(shù)據(jù)
clear?all;
C=cell（21）;
for?i=1:21
????C{i}=imread（strcat（‘frame_0‘num2str（i-1‘%02d‘）‘0.jpg‘））;
end

%核密度值估計
h=80;
medpr=ones（288384）;
finpr=zeros（288384）;
tic;%核函數(shù)值運算記時
for?n=1:20
????for?i=1:3
????????A=double（C{21}（::i））;
????????B=double（C{n}（::i））;
????????sinpr=ones（288384）-（1/h^2）*（A-B）.*（A-B）;
????????sinpr（sinpr<0）=0;%滿足核函數(shù)條件；
????????medpr=medpr.*sinpr;
????end
????finpr=finpr+medpr;
????medpr=ones（288384）;
end
pr=（15/（8*pi*20*h^3））*finpr;
toc;

%繪制核密度函數(shù)估計值圖像
mesh（pr）;

%繪制二值化圖像
pr（pr<2*10^（-7））=0;
pr（pr>2*10^（-7））=1;
pr=imcomplement（pr）;
figure;?
imshow（pr）;

????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????51009??2010-10-29?12:18??0000.jpg
?????文件???????50900??2010-10-29?12:20??0010.jpg
?????文件???????50806??2010-10-29?12:20??0020.jpg
?????文件???????51003??2010-10-29?12:20??0030.jpg
?????文件???????50636??2010-10-29?12:20??0040.jpg
?????文件???????50620??2010-10-29?12:21??0050.jpg
?????文件???????50545??2010-10-29?12:21??0060.jpg
?????文件???????50446??2010-10-29?12:21??0070.jpg
?????文件???????50397??2010-10-29?12:21??0080.jpg
?????文件???????50421??2010-10-29?12:21??0090.jpg
?????文件???????50472??2010-10-29?12:22??0100.jpg
?????文件???????50355??2010-10-29?12:22??0110.jpg
?????文件???????50318??2010-10-29?12:22??0120.jpg
?????文件???????50285??2010-10-29?12:22??0130.jpg
?????文件???????50216??2010-10-29?12:22??0140.jpg
?????文件???????50307??2010-10-29?12:23??0150.jpg
?????文件???????50265??2010-10-29?12:23??0160.jpg
?????文件???????50276??2010-10-29?12:23??0170.jpg
?????文件???????50327??2010-10-29?12:23??0180.jpg
?????文件???????50256??2010-10-29?12:23??0190.jpg
?????文件???????50176??2010-10-29?12:24??0200.jpg
?????文件?????????700??2018-11-10?19:36??KED_cpu.m