核密度估计

提供二维的核密度估计，详细注解，希望能对大家有所帮助

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  第三次作业（核密度估计）\