X=(x1,x2,...,xm)∈Rd×m , Y=(y1,y2,...,yl)∈Rd×l , l≤m . xi↔yi , i=1,…,l . I 为 l×l 单位矩阵。
已知对应点信息的保持:
∑li=1(xi−yi)2=tr((XY)d×(m+l)⎛⎝⎜⎜I0−I000−I0I⎞⎠⎟⎟(m+l)×(m+l)(XTYT)(m+l)×d)
其中 [I000] 为 m×m 矩阵, [I0] 为 m×l 矩阵, [I0] 为 l×m 矩阵
或者 ∑li=1(xi−yi)2=tr((X1,...,lY)d×(l+l)(I−I−II)(l+l)×(l+l)(XT1,...,lYT)(l+l)×d)
非对应点尽量远离,第
i
个点
∑li=1∑j≠i(xi−yj)2=∑li=1∑lj=1(xi−yj)2−∑li=1(xi−yi)2
其中 ∑li=1∑lj=1(xi−yj)2 = ∑mi=1∑lj=1(xi−yj)2WijX,Y (其中 WijX,Y=1,i=1,...,l,WijX,Y=0,i=l+1,...,m
=trace(
(XY)(Ω1−WY,X−WX,YΩ2)(XTYT)
)
其中
Ω1∈Rm×m
对角矩阵,
Ωii1=∑jWijX,Y={l,i=1,...,l0,i=l+1,...,m
,
Ω1=[lIl×l000]
Ω2∈Rl×l
对角矩阵,
Ωii2=∑jWijX,Y=l
,
Ω2=lIl×l
所以,原式=trace(
(XY)⎛⎝⎜⎜lIl0−El000−El0lIl⎞⎠⎟⎟(XTYT)
),其中
El
为
l×l
的元素全为1的矩阵。
保持流形
X
的局部结构:
LX=DX−WX , DiiX=∑jWijX
两个流形 X1∈Rd×m1 , X2∈Rd×m2 间的对应关系的保持:
∑m1i=1∑m2j=1(xi1−xj2)2Wi,j1,2
=trace( (X1X2)(Ω1−W2,1−W1,2Ω2)(XT1XT2) )
其中
Ω1∈Rm1×m1
对角矩阵,
Ωii1=∑jWij1,2
Ω2∈Rm2×m2
对角矩阵,
Ωii2=∑jWij2,1
保持流形 X1 , X2 的局部结构,并保持流形 X1 与 X2 之间的对应关系:
12μ2∑m1i=1∑m1j=1(xi1−xj1)2Wij1 + 12μ2∑m2i=1∑m2j=1(xi2−xj2)2Wij2 + μ1∑m1i=1∑m2j=1(xi1−xj2)2Wi,j1,2
=
μ2
trace(
X1L1XT1
)+
μ2
trace(
X2L2XT2
)+
μ1∑m1i=1∑m2j=1(xi1−xj2)2Wi,j1,2
=
μ2
trace(
(X1X2)(L1000)(XT1XT2)
)
+
μ2
trace(
(X1X2)(000L2)(XT1XT2)
)
+
μ1
trace(
(X1X2)(Ω1−W2,1−W1,2Ω2)(XT1XT2)
)
=trace( [X1X2][μ1Ω1+μ2L1−μ1W2,1−μ1W1,2μ1Ω2+μ2L2][XT1XT2] )
假设有两个流形 X1∈Rd×m1 , X2∈Rd×m2 :
∑m1i=1∑m2j=1(xi1−xj2)2
= ∑m1i=1∑m2j=1(xi1−xj2)2Wij1,2 (其中 Wij1,2=1 , W1,2 为 m1×m2 矩阵,其每个元素都为1)
=trace( (X1X2)(Ω1−W2,1−W1,2Ω2)(XT1XT2) )
其中
Ω1∈Rm1×m1
对角矩阵,
Ωii1=∑jWij1,2=m2
,
Ω1=m2Im1×m1
Ω2∈Rm2×m2
对角矩阵,
Ωii2=∑jWij2,1=m1
,
Ω2=m1Im2×m2
xi
是
p×1
列向量
||∑ni=1xi||2=(∑ni=1xi)T(∑ni=1xi)=∑ni=1∑nj=1xTixj=tr(XIXT)
,其中
I
是每个元素都为1的
xi
是
p×1
列向量,
yj
是
p×1
列向量
(∑ni=1xi)T(∑mj=1yj)=∑ni=1∑mj=1xTiyj=tr(XIYT)
,其中
I
是每个元素都为1的