linregfun.m

来自「数据挖掘的新方法－支持向量机书中算法例子」· M 代码 · 共 30 行

%
%

function a1=LinRegFUN(a)

fid=fopen('a.dat','r');
    n=fread(fid,1,'float');
    e=fread(fid,1,'float');
    X=fread(fid,n,'float');
    Y=fread(fid,n,'float');
fclose(fid);
    
DmX=size(a);
DmX=size(X);
DmY=size(Y);
if DmX(1)<DmX(2),
    a=a';
end
if DmX(1)<DmX(2),
    X=X';
end
if DmY(1)<DmY(2),
    Y=Y';
end

FUN=inline('-0.5*sum(sum(   diag(a(1:n)-a(n+1:2*n))*  ((X*transpose(X)).^1)  *diag(a(1:n)-a(n+1:2*n))  )) - e*sum(a) + transpose(Y)*(a(1:n)-a(n+1:2*n))','a','X','Y','n','e');

a1=-FUN(a,X,Y,n,e);