翻译自Euclidean Distance Matrices: Essential theory, algorithms, and applications
EDMs是点之间的平均距离矩阵。该文的目标是介绍EMD在信号处理领域的应用,展示EDM如何被用来设计算法--对距离数据进行修复和去噪。同时,介绍了其在麦克风的位置校准(microphone position calibration)、超声波断层扫描(ultrosound tomography)、room reconstruction from echoes以及相位恢复的应用。
引言
假设你有一份瑞士火车时刻表但是没有其地图,但是这足够去重建一份阿尔卑斯的粗略的地图,即便火车时间不能体现距离甚至部分时间点是不知道的。
我们经常处理距离数据因为它们是容易测量和评估的。比如说,在无线传感器网络,传感器节点可以测量由其它节点传送过来的信号包的强度或由它们的邻节点发送的脉冲到达时间(TOA),这属于自我定位(self-localization)。
有时候数据不是矩阵的,但是我们可以像心理测验学一样寻找一个矩阵表示。实际上,心理测验学是很多跟EDM有关的方法形成的起源,包括多维标度分析(MDS)--用多维空间中的点表示不同刺激物之间的感知或心理测量关系(百度)。
EDM是对点集的一个有用的描述和算法设计的好起点。一个经典的任务就是恢复原始点设置:它仅仅需要对一个对称矩阵进行特征值分解(EVD)。事实上,大部分欧几里德距离问题重建点集伴随着以下问题:
1)距离数据有噪声
2)一些距离数据缺失
3)距离无法识别(unlabelled)
距离几何有两个基本问题:一是给定一个矩阵,判断它是否是EDM;二是给定一些不完整的距离数据,在已知嵌入维数下--形成点的最小仿射空间的维度--判断是否存在一种点配置(configuration of points)
FroM PoIntS to EdMs And BAcK
X为d*n矩阵,
xi和xj之间的平均距离为,并定义
||.||表示欧几里德范数,即
1表示n*1的列向量,即[1 1 1...1 ]的转置,1的转置即1*n的行向量;edm(X)实际就是
类似地,令 ,则
公式3和4揭示了一个重要的属性:X的秩最多为d,的秩也最多为d,公式3的另外两个加数秩为1,由此可得:一个EDM的秩最多为d+2 (T1)
这个定理说明EDM的秩和其点的数目无关,在许多应用中,d是3或者更少然而n可以是上千。根据定理1,这样的矩阵秩最多为5。定理1中最重要的是点集的仿射维数。任何仿射子空间都是一个线性子空间的转变,即重要的唯一性(essential uniqueness)。如下图所示,从一个仿射子空间的所有点中取出任意一点可唯一地描述一个包含零向量的平行线性子空间。
ESSENTIAL UNIQUENESS
当处理一个逆问题时(inverse problem),我们需要明白哪些是可重获的哪些是已丢失的。重建点配置通常会增加其规模(size),成对距离的数量远比坐标描述的规模大,即。很明显刚性变换不改变固定点的距离,就像我们从公式3和4中得知的:edm(X)=,故用代数方法(翻转、反射)不会改变距离。因此对于旋转点集(Q为d*d的正交矩阵,即):
通过d*1的列向量b得到平移矩阵
因此
这说明我们无法只通过距离数据得到点的具体位置,不同的重建步骤会得到不同的点配置,如下图所示:
RECONSTRUCTION THE POINT SET FROM DISTANCE
公式3给了我们一种从距离矩阵中计算点配置的方法。
假设点x1是原点,D的第一列包含了点向量的平方范数,d1是D的第一列,即
由于的对角线元素恰好是,则
由于G是对称半正定矩阵(PSD),,其中
所以,即为重建点配置,。
于是,我们有定理2:当且仅当对任何满足条件的s
是一个PSD时,D是一个EDM。
很容易得知,s=e1。平移点集使得x1被平移到原点,通过对edm(x)左乘矩阵右乘矩阵,我们将最终得到,重建点配置将在原点有第一个点(have the first point at the origin)。
另一方面,s=(1/n)1时,坐标系的原点将是点集的质心,因此矩阵被称为几何中心矩阵
为说明这一点,我们定义X的质心为所有点集的平均值
然后从点集中减去这一向量,即