nvldtcfun.m

数据挖掘的新方法－支持向量机书中算法例子

function L=NvlDtcFun(a);

Dim=size(a);
if Dim(1)<Dim(2),   % Conditioning the vector A as a column vector
    a=a';
end

fid=fopen('a.dat','r');
    n=fread(fid,1,'float');
    x1=fread(fid,n,'float');
    x2=fread(fid,n,'float');
fclose(fid);

% Linear Kernel Function K(x,y)=x*y;
X=[x1,x2]; % a n by 2 matrix;
L=sum(sum(diag(a)*( (X*X') )*diag(a))) - sum((X').^2)*a ;