对shape context的一种改进可图形集后直接运行

本人在作者shape context上改动，如要使用 请注明来源

function [cxcyEL]=bookstein（XYbeta_k）;
% [cxcyEL]=bookstein（XYbeta_k）;
%
% Bookstein PAMI89

N=size（X1）;
Nb=size（Y1）;

if N~=Nb
   error（‘number of landmarks must be equal‘）
end

% compute distances between left points
r2=dist2（XX）;

K=r2.*log（r2+eye（NN））; % add identity matrix to make K zero on the diagonal
P=[ones（N1） X];
L=[K  P
   P‘ zeros（33）];
V=[Y‘ zeros（23）];
if nargin>2
   % regularization
   L（1:N1:N）=L（1:N1:N）+beta_k*eye（NN）;
end
invL=inv（L）;

c=invL*V‘;
cx=c（:1）;
cy=c（:2）;

if nargout>2
   % compute bending energy （w/o regularization）
   Q=c（1:N:）‘*K*c（1:N:）;
   E=mean（diag（Q））;
end

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

     文件        659  2000-09-05 08:28  scmatlab1\bookstein.m

     文件        287  2013-05-30 19:42  scmatlab1\calculateEdit.m

     文件     648546  2013-05-29 16:52  scmatlab1\Cs.mat

     文件     648477  2013-05-29 16:52  scmatlab1\Cs1.mat

     文件    1188339  2013-06-03 17:59  scmatlab1\dismat.mat

     文件        201  2013-06-03 19:08  scmatlab1\dist.mat

     文件        820  2000-09-05 08:26  scmatlab1\dist2.m

     文件        501  2013-06-03 19:12  scmatlab1\dist_bw_sc.m

     文件       1488  2013-05-30 20:02  scmatlab1\DProgramming.m

     文件        630  2013-05-28 22:25  scmatlab1\extract.m

     文件        680  2013-05-21 20:00  scmatlab1\extract_longest.m

     文件        432  2000-09-05 08:27  scmatlab1\hist_cost_2.m

     文件       4517  2013-05-30 19:41  scmatlab1\html\calculateEdit.html

     文件      11250  2000-09-05 08:28  scmatlab1\hungarian.m

     文件        222  2004-06-09 09:46  scmatlab1\i2s.m

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_01.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_02.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_03.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_04.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_05.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_06.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_07.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_08.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_09.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_10.bmp

     文件      19326  2004-07-07 00:01  scmatlab1\mpeg7\01_11.bmp

     文件      18934  2004-07-07 00:01  scmatlab1\mpeg7\01_12.bmp

     文件      14382  2004-07-07 00:01  scmatlab1\mpeg7\01_13.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_14.bmp

     文件      20166  2004-07-07 00:01  scmatlab1\mpeg7\01_15.bmp

............此处省略1464个文件信息