correlation.m

本人编的一个程序

function [fx, dfx] = correlation(x,g,samples,aaa,param,constant)

bbb = sum(x.^2);
kvec = exp(-param.kernelparam*(aaa + bbb - 2*(samples*x)));
gkvec = g'*kvec;
fx = -gkvec^2;
dfx = -4*gkvec*param.kernelparam*((g.*kvec)'*samples - gkvec*x')';