function y=kernel_bmin(b,xdat)

%求核密度函數中窗寬的目標函數
xmin=min(xdat);
xmax=max(xdat);

x=linspace(xmin,xmax,10);
yy=kernel_density(x,xdat,b,'epan');

y=hist_xzh1(xdat);  %核密度函數擬合頻率直方圖
y=sum((yy-y).^2);