高斯核支持向量数据描述（SVDD）方法的源码

高斯核支持向量数据描述（SVDD）方法的源码

function [K] = kernel（kerxy）
% Calculate kernel function.   
% x: 输入样本n1×d的矩阵n1为样本个数d为样本维数
% y: 输入样本n2×d的矩阵n2为样本个数d为样本维数
% ker  核参数（结构体变量）
% the following fields:
%   type   - linear :  k（xy） = x‘*y
%            poly   :  k（xy） = （x‘*y+c）^d
%            gauss  :  k（xy） = exp（-0.5*（norm（x-y）/s）^2）  
%                                        % norm命令计算矩阵或向量（x-y）的2-范数
%            tanh   :  k（xy） = tanh（g*x‘*y+c）
%   degree - Degree d of polynomial kernel （positive scalar）.
%   offset - Offset c of polynomial and tanh kernel （scalar negative for tanh）.
%   width  - Width s of Gauss kernel （positive scalar）.
%   gamma  - Slope g of the tanh kernel （positive scalar）.
% ker = struct（‘type‘‘linear‘）;
% ker = struct（‘type‘‘ploy‘‘degree‘d‘offset‘c）;
% ker = struct（‘type‘‘gauss‘‘width‘s）;
% ker = struct（‘type‘‘tanh‘‘gamma‘g‘offset‘c）;
% K: 输出核参数n1×n2的矩阵
%-------------------------------------------------------------%

switch ker.type
    case ‘linear‘
        K = x*y‘;         % n*n的矩阵
    case ‘poly‘
        d = ker.degree;
        c = ker.offset;
        K = （x*y‘+c）.^d;
    case ‘gauss‘
        s = ker.width;
        rows = size（x1）;
        cols = size（y1）;   
%         K=zeros（rowscols）;
%         for i = 1:rows
%             for j = 1:cols
%                 tmp = norm（x（i:）-y（j:））;
%                 K（ij） = exp（-（tmp^2）/（s^2））;
%             end
%         end
      tmp = zeros（rowscols）;        
        for i = 1:rows
            for j = 1:cols
                tmp（ij） = norm（x（i:）-y（j:））;
               
            end
        end
     K = exp（-（tmp.^2）/（s.^2））;    
    case ‘tanh‘
        g = ker.gamma;
        c = ker.offset;
        K = tanh（g*x*y‘+c）;
    otherwise
        K = 0;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1871  2013-04-01 22:23  exp15_svdd\kernel.m

     文件       4030  2013-04-06 14:00  exp15_svdd\sample_banana_200.mat

     文件       1337  2013-05-10 16:35  exp15_svdd\svddSim.m

     文件       1701  2013-04-02 11:52  exp15_svdd\svddTrain.m

     文件       1277  2014-05-24 09:14  exp15_svdd\svdd_banana_gauss.m

     文件         18  2014-05-24 10:03  exp15_svdd\说明.txt

     目录          0  2014-05-24 10:03  exp15_svdd

----------- ---------  ---------- -----  ----

                10234                    7