KPCA-故障检测代码已优化

基于核主元分析（KPCA）的工业过程故障检测，代码已优化，运行效率高，有详细的注释，附有训练数据和测试数据。

function [K] = computeKernelMatrix（kerxy）
% Calculate kernel function.   
%
% x: 输入样本d×n1的矩阵n1为样本个数d为样本维数
% y: 输入样本d×n2的矩阵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）
%            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‘;
        
    case ‘ploy‘
        d = ker.degree;
        c = ker.offset;
        K = （x‘*y+c）.^d;
        
    case ‘gauss‘
        s = ker.width;
        sx = sum（x.^22）;
        sy = sum（y.^22）;
        K = exp（0.5*（bsxfun（@minusbsxfun（@minus2*x*y‘sx）sy‘））/s^2）;
    case ‘tanh‘
        g = ker.gamma;
        c = ker.offset;
        K = tanh（g*x‘*y+c）;
    otherwise
        K = 0;
end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1451  2018-03-17 09:23  computeKernelMatrix.m
     文件        3220  2018-03-17 09:22  main.m
     文件       27697  2018-03-17 08:47  testData.mat
     文件       27448  2018-03-17 08:48  trainData.mat