0001 function h = chi2_kernel(x1,x2,gamma) 0002 ind = ((x1+x2)>0); 0003 %x1=x1(ind); x2=x2(ind); 0004 h = exp(-gamma^2*sum((x1(ind)/sum(x1)-x2(ind)/sum(x2)).^2./(x1(ind)+x2(ind)))); 0005 0006