function g = myParzen( x , xi , h )
%本程序完成Parzen窗估計和估計分類
%  x: 測試樣本點
%  xi: 訓練樣本點
%  h: 窗寬參數
%  g: 判決函數
[lx,ly] = size( xi );
Vn = 4/3 * pi * (h/2)^lx;

g = 0;
for i = 1:ly
    g = g + exp(-(x-xi(:,i))'*(x-xi(:,i))/(2*h^2));
end
g = g / (Vn*ly);