循环矩阵
充分利用循环矩阵及其特性的是核相关滤波跟踪算法的另一个重要特征,它不仅涉及到目标采样,而且巧妙的将目标特征的频域空间与岭回归相结合,实现了目标特征的快速学习与检测。
首先考虑一维样本的情况,设
x=[x0,x1,x2,⋯,xN−1]
表示一行图像像素样本,
x
的循环矩阵表示为:
X=C(x)=⎡⎣⎢⎢⎢⎢⎢⎢⎢x0xN−1xN−2⋮x1x1x0xN−1⋮x2x2x1x0⋮x3⋯⋯⋯⋱⋯xN−1xN−2xN−3⋮x0⎤⎦⎥⎥⎥⎥⎥⎥⎥(1)
其中, x 又称为 X 的生成向量。
循环矩阵有一个重要的性质为:任意循环矩阵均可以被离散傅立叶变换矩阵对角化,即:
X=C(x)=F⋅diag(F(x))⋅FH(2)
其中, F(⋅) 表示离散傅立叶变换(Discrete Fourier Transform, DFT), F 表示DFT矩阵,具有如下形式:
F=1n√⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢111⋮11ww2⋮wN−11w2w4⋮w2(N−1)⋯⋯⋯⋱⋯1wN−1w2(N−1)⋮w(N−1)2⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥(3)
其中, w=e−j2π/N 为旋转因子,且 F 是一个酉矩阵,满足 FHF=FFH=I 。式(1)的证明过程见文献[49]. 通过式(2.7),以及 F 的性质,可以继续推导出循环矩阵具有以下特性:
已知 A,B 都为循环矩阵,其生成向量分别是 a,b ,且具有相同的向量长度,则有以下公式成立:
A+B=C(a+b)(4)
A⋅B=C(F−1(F(a)⊙F(b)))(5)
其中 F−1(⋅) 表离散傅立叶反变换(IDFT), ⊙ 表对位相乘运算。证明:
A+B=F⋅diag(F(a))⋅FH+F⋅diag(F(b))⋅FH=F⋅diag(F(a)+F(b))⋅FH=C(F−1(F(a)+F(b)))=C(a+b)
A⋅B=F⋅diag(F(a))⋅FH⋅F⋅diag(F(b))⋅FH=F⋅diag(F(a))⋅diag(F(b))⋅FH=F⋅diag(F(a)⊙F(b))⋅FH=C(F−1(F(a)⊙F(b)))
证毕已知 X 为循环矩阵,其生成向量为 x ,且向量元素为实数,则有以下公式成立:
XT=F⋅diag(F∗(x))⋅FH(6)
X−1=F⋅diag(F(x))−1⋅FH(7)证明:由于已知 X 为实矩阵,即 X=X∗ ,且 F 为对称酉矩阵,则:
XT=(F⋅diag(F(x))⋅FH)T=(FH)T⋅diag(F(x))⋅FT=F∗⋅diag(F(x))⋅FT=(F∗⋅diag(F(x))⋅FT)∗=F⋅diag(F∗(x))⋅FH
X−1=(F⋅diag(F(x))⋅FH)−1=(FH)−1⋅diag(F(x))−1⋅F−1=F⋅diag(F(x))−1⋅FH
证毕已知 X 为循环矩阵,其生成向量为 x , y 为与 x 具有相同向量长度的向量,则有以下公式成立:
F(Xy)=F∗(x)⊙F(y)(8)证明:
Xy=C(x)y=⎡⎣⎢⎢⎢⎢⎢⎢⎢x0xN−1xN−2⋮x1x1x0xN−1⋮x2x2x1x0⋮x3⋯⋯⋯⋱⋯xN−1xN−2xN−3⋮x0⎤⎦⎥⎥⎥⎥⎥⎥⎥⎡⎣⎢⎢⎢⎢⎢⎢⎢y0y1y2⋮yN−1⎤⎦⎥⎥⎥⎥⎥⎥⎥
上式可以表示为向量 x~=[x0,xN−1,xN−2,⋯,x1]T 与向量 y 的循环卷积,即:
Xy=C(x)y=x~⊗y
其中 ⊗ 表示循环卷积。
令 Fx~y=F(x~⊗y) , Fx~=F(x~) , Fy=F(y) ,根据卷积定理,有Fx~y(n)=Fx~(n)⋅Fy(n)
即
Fx~y=Fx~⊙Fy=F(x~)⊙F(y)
根据DFT的定义,
Fx~(u)=∑n=0N−1x~(n)wnu=x0w0+xN−1wu+xN−2w2u+⋯+x1w(N−1)u=x0w(N−N)u+x1w(N−1)u+x2w(N−2)u+⋯+xN−1wu=∑n=0N−1xnw(N−n)u
利用旋转因子 w 的周期性,
wN=e−j2πNN=1;w(N−n)u=wNu⋅w−nu=(wN)u⋅w−nu=w−nu
可得:
Fx~(u)=∑n=0N−1xnw−nu
而向量 x 的DFT表达式为:
Fx(u)=∑n=0N−1x(n)wnu=∑n=0N−1xnwnu
可以发现, Fx~ 与 Fx 是互为共轭的关系,由此可得:
F(Xy)=F(x~)⊙F(y)=F∗(x)⊙F(y)
证毕
以上几条性质,对相关滤波器的证明非常重要。