function?[d]?=?kde2（?data?）
subplot（121）;
[hdxy]=kdde（data）;
surf（xyd‘linestyle‘‘none‘）;
subplot（122）;
plot（data（:1）data（:2）‘r.‘）;
end
function?[hdxy]=kdde（data）

if?nargin<2
????n=2^8;
end
n=2^ceil（log2（n））;?%?round?up?n?to?the?next?power?of?2;
global?N;
N=size（data1）;
if?nargin<3
????MAX=max（data[]1）;?MIN=min（data[]1）;?Range=MAX-MIN;
????MAX_XY=MAX+Range/4;?MIN_XY=MIN-Range/4;
end
scaling=MAX_XY-MIN_XY;
if?N<=size（data2）
????error（‘data?has?to?be?an?N?by?2?array?where?each?row?represents?a?two?dimensional?observation‘）
end
global?A2?I
transformed_data=（data-repmat（MIN_XYN1））./repmat（scalingN1）;
%bin?the?data?uniformly?using?regular?grid;
initial_data=ndhist（transformed_datan）;
%?discrete?cosine?transform?of?initial?data
a=?dct2d（in

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????14402??2012-11-20?16:10??第三次作業（核密度估計）\Point2-1.txt
?????文件????????3630??2012-12-25?15:02??第三次作業（核密度估計）\kde2.m
?????文件??????579946??2012-12-25?12:34??第三次作業（核密度估計）\核密度估計.pptx
?????目錄???????????0??2012-12-25?20:43??第三次作業（核密度估計）\