parzen_index.m

来自「模式识别中Parzen窗函数的实现」· M 代码 · 共 28 行

%window function: exponential distribution func.
%
% (C)2008 TangSheng
%
%---------指数分布窗函数------------%
function parzen_index(XX,YY,ZZ,Xi,n,nn,h1)
hn = h1/sqrt(n);
Vn = hn*hn;
for i = 1:sqrt(nn)
    for j = 1:sqrt(nn)
        ZZ(i,j) = 0;
    end
end
        
for i = 1:sqrt(nn)
    for j = 1:sqrt(nn)
          sum = 0;
        for k =1:n
            u = sqrt((XX(i,j)-Xi(k,1))^2+(YY(i,j)-Xi(k,2))^2);
            sum =sum+0.5*exp(-abs(u));
        end
        ZZ(i,j) = sum/n;
    end
end
mesh(XX,YY,ZZ);
name = ['指数窗估计：','N=',num2str(n),', h1=',num2str(h1)];
title(name);
end