function?[ctrainctest]?=?centerTrainTestKernels（trkerneltskernel）

%?[ctrainctest]?=?centerTrainTestKernels（trkerneltskernel）?:?centers
%?train/test?kernels.


l?=?size（trkernel1）;
j?=?ones（l1）;
ctrain?=?trkernel?-?（j*j‘*trkernel）/l?-?（trkernel*j*j‘）/l?+?（（j‘*trkernel*j）*j*j‘）/（l^2）;

if（?nargin?>?1?）
????tk?=??（1/l）*sum（trkernel1）;?%?（1?x?l）
????tl?=?ones（size（tskernel1）1）;?%?（n?x?1）
????ctest?=?tskernel?-?（?tl?*?tk）;?%?（?n?x?l?）
????tk?=?（1/（size（tskernel2）））*sum（ctest2）;?%?（?n?x?1?）???
????ctest?=?ctest?-?（tk?*?j‘）;?%?（?n?x?l?）
end