parzen.m

来自「这是一个parzen窗口分类算法」· M 代码 · 共 16 行

function p=Parzen(xi,x,h1,f)
%xi为样本，x为概率密度函数的自变量的取值，
%h1为样本数为1时的窗宽，f为窗函数 
%返回x对应的概率密度函数值

if isempty(f)    %若没有指定窗的类型，就使用正态窗函数句柄
     
     f=@(u)(1/sqrt(2*pi))*exp(-0.5*u.^2);
end;
N=size(xi,2);
hn=h1/sqrt(N);
[X Xi]=meshgrid(x,xi);
p=sum(f((X-Xi)/hn)/hn)/N;