联合分布适配JDA代码

function [ZA] = JDA（XsXtYsYt0options）

% Transfer Feature Learning with Joint Distribution Adaptation.  
% M. Long J. Wang G. Ding J. Sun and P.S. Yu.
% IEEE International Conference on Computer Vision （ICCV） 2013.

% Contact: Mingsheng Long （longmingsheng@gmail.com）

% Load algorithm options
addpath（genpath（‘../liblinear/matlab‘））;

if nargin < 5
    error（‘Algorithm parameters should be set!‘）;
end
if ~isfield（options‘k‘）
    options.k = 100;
end
if ~isfield（options‘lambda‘）
    options.lambda = 0.1;
end
if ~isfield（options‘ker‘）
    options.ker = ‘primal‘;
end
if ~isfield（options‘gamma‘）
    options.gamma = 1.0;
end
if ~isfield（options‘data‘）
    options.data = ‘default‘;
end
k = options.k;
lambda = options.lambda;
ker = options.ker;
gamma = options.gamma;
data = options.data;

fprintf（‘JDA:  data=%s  k=%d  lambda=%f
‘dataklambda）;

% Set predefined variables
X = [XsXt];
X = X*diag（sparse（1./sqrt（sum（X.^2））））;
[mn] = size（X）;
ns = size（Xs2）;
nt = size（Xt2）;
C = length（unique（Ys））;

% Construct MMD matrix
e = [1/ns*ones（ns1）;-1/nt*ones（nt1）];
M = e*e‘*C;
if ~isempty（Yt0） && length（Yt0）==nt
    for c = reshape（unique（Ys）1C）
        e = zeros（n1）;
        e（Ys==c） = 1/length（find（Ys==c））;
        e（ns+find（Yt0==c）） = -1/length（find（Yt0==c））;
        e（isinf（e）） = 0;
        M = M + e*e‘;
    end
end
M = M/norm（M‘fro‘）;

% Construct centering matrix
H = eye（n）-1/（n）*ones（nn）;

% Joint Distribution Adaptation: JDA
if strcmp（ker‘primal‘）
    [A~] = eigs（X*M*X‘+lambda*eye（m）X*H*X‘k‘SM‘）;
    Z = A‘*X;
else
    K = kernel（kerX[]gamma）;
    [A~] = eigs（K*M*K‘+lambda*eye（n）K*H*K‘k‘SM‘）;
    Z = A‘*K;
end

fprintf（‘Algorithm JDA terminated!!!

‘）;

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-11-09 22:26  CodeData
     文件       12292  2013-11-09 22:26  CodeData.DS_Store
     目录           0  2013-11-09 22:25  CodeDatacode
     文件        6148  2013-11-09 22:05  CodeDatacode.DS_Store
     文件        1792  2013-11-09 22:07  CodeDatacodeJDA.m
     目录           0  2013-11-09 22:26  __MACOSX
     目录           0  2013-11-09 22:26  __MACOSXCodeData
     目录           0  2013-11-09 22:26  __MACOSXCodeDatacode
     文件         120  2013-11-09 22:07  __MACOSXCodeDatacode._JDA.m
     文件        1297  2013-11-09 22:07  CodeDatacodekernel.m
     文件         164  2013-11-09 22:07  __MACOSXCodeDatacode._kernel.m
     文件        1998  2013-11-09 22:07  CodeDatacode
un_office.m
     文件         120  2013-11-09 22:07  __MACOSXCodeDatacode._run_office.m
     文件        1633  2013-11-09 22:07  CodeDatacode
un_other.m
     文件         120  2013-11-09 22:07  __MACOSXCodeDatacode._run_other.m
     文件        1999  2013-11-09 22:18  CodeDatacode
un_pie.m
     文件         120  2013-11-09 22:18  __MACOSXCodeDatacode._run_pie.m
     文件         164  2013-11-09 22:25  __MACOSXCodeData._code
     目录           0  2013-11-09 22:04  CodeDatadata
     文件        6148  2013-04-14 15:39  CodeDatadata.DS_Store
     文件      104970  2013-01-09 11:14  CodeDatadataamazon_SURF_L10.mat
     文件        3692  2013-01-09 11:14  CodeDatadataamazon_SURF_L10_imgs.mat
     文件      150751  2013-01-09 11:14  CodeDatadataCaltech10_SURF_L10.mat
     文件        4251  2013-01-09 11:14  CodeDatadataCaltech10_SURF_L10_imgs.mat
     文件     3016148  2013-04-05 22:12  CodeDatadataCOIL_1.mat
     目录           0  2013-11-09 22:26  __MACOSXCodeDatadata
     文件         164  2013-04-05 22:12  __MACOSXCodeDatadata._COIL_1.mat
     文件     3016146  2013-04-05 22:12  CodeDatadataCOIL_2.mat
     文件         164  2013-04-05 22:12  __MACOSXCodeDatadata._COIL_2.mat
     文件       18348  2013-01-09 11:14  CodeDatadatadslr_SURF_L10.mat
     文件         893  2013-01-09 11:14  CodeDatadatadslr_SURF_L10_imgs.mat
............此处省略16个文件信息