陆振波的svm的matlab代码的解释,陆振波SVM的MATLAB代码解释

陆振波的SVM的MATLAB代码的解释

%构造训练样本

n = 50;

randn('state',6);

x1 = randn(2,n); %2行N列矩阵 y1 = ones(1,n); %1*N个1 x2 = 5+randn(2,n); %2*N矩阵 y2 = -ones(1,n); %1*N个-1

figure;

plot(x1(1,:),x1(2,:),'bx',x2(1,:),x2(2,:),'k.');

%x1(1,:)为x1的第一行，x1(2,:)为x1的第二行

axis([-3 8 -3 8]); title('C-SVC') hold on;

X = [x1,x2]; %训练样本d*n矩阵，n为样本个数，d为特征向量个数 Y = [y1,y2]; %训练目标1*n矩阵，n为样本个数，值为+1或-1

%训练支持向量机

function svm = svmTrain(svmType,X,Y,ker,p1,p2)

options = optimset; % Options是用来控制算法的选项参数的向量 http://doc.docsou.comrgeScale = 'off'; options.Display = 'off';

switch svmType case 'svc_c',

C = p1;

n = length(Y);

H = (Y'*Y).*kernel(ker,X,X);

f = -ones(n,1); %f为1*n个-1,f相当于Quadprog函数中的c A = []; b = [];

Aeq = Y; %相当于Quadprog函数中的A1,b1 beq = 0;

lb = zeros(n,1); %相当于Quadprog函数中的LB，UB ub = C*ones(n,1);

a0 = zeros(n,1); % a0是解的初始近似值

[a,fval,eXitflag,output,lambda]

quadprog(H,f,A,b,Aeq,beq,lb,ub,a0,options); %a是输出变量，它是问题的解

=

% Fval是目标函数在解a 处的值 % Exitflag>0,则程序收敛于解x

Exitflag=0,则函数的计算达到了最大次数 Exitflag<0,则问题无可行解，或程序运行失败

% Output 输出程序运行的某些信息