关于2个空间信号的协方差矩阵的问题


9

每当我认为我了解了协方差矩阵时,都会有人提出不同的表述。

我目前正在阅读本文:

J. Benesty,“用于无源声源定位的自适应特征值分解算法”,J。Acoust。Soc。上午。第107卷,第1期,第384-391页(2000)

我遇到了一个我不太了解的表述。在这里,作者正在构建两个信号x1和之间的协方差矩阵。这两个信号来自不同的传感器。x2

对于一个信号的协方差矩阵,我知道我们可以通过计算回归矩阵,然后将其乘以同一矩阵的Hermitian,然后除以(原始向量的长度)来获得。这里的协方差矩阵的尺寸可以是任意的,具有最大尺寸为。NN×N

对于两个空间信号的协方差矩阵,如果将第一个信号放在矩阵的第一行中,将第二个信号放在矩阵的第二行中,然后乘以其Hermitian,再除以,则得到两个空间信号的协方差矩阵。N2×2

然而,在本文中,笔者单位计算什么样子4个 matricies,,和,然后将它们放入一个超级矩阵和呼叫协方差矩阵。R11,R12,R21R22

为什么会这样呢?这是文本的图片:

在此处输入图片说明

Answers:


6

如果每个N个元素都有两个信号向量x 2 [ n ],那么我们可以考虑两个不同的事物。x1[n]x2[n]N

  1. 如何量比较?特别是,当所述信号是嘈杂和的噪声可以被认为是共同固定的(或共同广义平稳),这些量可用于估计噪声在两个信号方差以及噪声的协方差在任何固定的采样时间。 这就是您从2 × 2得到的n=1Nxi[n]xj[n], i,j{1,2}2×2 协方差矩阵 在噪声X 1 [ Ñ ]具有方差 σ 2 1 = - [R 1 1,这可能是不同的从- [R 2 2 = σ 2 2,在噪声的方差X 2 [ Ñ ]。但是噪声与协方差R 1.2相关

    [R2×2=[σ1个2CCσ22]
    X1个[ñ]σ1个2=[R1个1个[R22=σ22X2[ñ]。现在,如果我们打算只处理在 n处发生的事情,而忽略在 n 1 n + 1 等处可能发生的事情,那么这就是我们所需要的全部信息。[R1.2=[R21个=Cññ-1个ñ+1个
  2. 除非已知该噪声为(或假定为)白噪声,以便来自不同采样时刻的噪声样本是独立的(因此是不相关的),或者我们只是假设不相关的噪声样本,否则我们将忽略一些信息,而无需考虑相关性在x 1 [ m ]之间,来自同一过程的采样在不同的时间或位置,以及x 1 [ n ]x 2 [ m ]之间的相关性X1个[ñ]X1个[]X1个[ñ]X2[],这是两个过程在不同时间或位置的样本。这些附加信息可能会导致更好的估计/解决方案。现在我们总共有噪声样本,因此需要考虑2 N × 2 N的协方差矩阵。如果我们按照作者的方式安排问题,我们有 R full = E [ X X T ],其中 X = x 1 [ 1 ] x 1 [ 2 ] x 1 [2ñ2ñ×2ñ[R充分=Ë[XXŤ] ,因此 R full = [ R x 1x 1 R x 1x 2 R x 2x 1 R x 2

    X=X1个[1个]X1个[2]X1个[ñ]X2[1个]X2[2]X2[ñ]Ť=X1个X2Ť
    ,其中Rxixj=E[xixTj]。需要注意的是- [RXXĴ是,在本质上, 互相关的函数X[1]X[2]...X[Ñ]XĴ[1
    [R充分=[[RX1个X1个[RX1个X2[RX2X1个[RX2X2]
    [RX一世XĴ=Ë[X一世XĴŤ][RX一世XĴX一世[1个]X一世[2]X一世[ñ]如果 i j,则 自相关函数,如果 i = j。如果噪声过程是白色的,不相关的,除了当 Ñ = ,然后 ř ř 简单 = [ σ 2 1ç Ç σ 2 2 ] 其中XĴ[1个]XĴ[2]XĴ[ñ]一世Ĵ一世=Ĵñ=
    [R充分[R简单=[σ1个2一世C一世C一世σ22一世]
    Ñ × Ñ同一性矩阵, σ 2 1σ 2 2 Ç是如第1项定义的上方。最终用户需要确定此噪声模型的逼真度。如果模型现实的,则通过查看 2 N × 2 N矩阵 R full不会获得任何结果, 因为所有信息都在 2 × 2矩阵 R 2 × 2中一世ñ×ñσ1个2σ22C2ñ×2ñ[R充分2×2[R2×2以上第1项。如果模型不切实际,但我们不打算(或无法)使用全部矩阵R full中的所有信息,则同上 ;我们将用只σ 2 1σ 2 2C ^第1条,我们不需要[R 完全[R 简单,只需[R 2 × 22ñ×2ñ[R充分σ1个2σ22C[R充分[R简单[R2×2

谢谢。首先,难道(1)中的sigma不应该说从n = 0到N-1吗?(不是从i = 1到n)。
Spacey 2012年

我不确定我是否仍然了解我们为什么这样做。您是在说(1),因为两个向量中的噪声都是彼此完全独立的,所以我们必须使用该方法,从而获得2x2的协方差矩阵,但是在第二种情况(2)中,因为向量中的噪声不是独立的,我们必须将两个向量连接起来,然后计算它们的协方差矩阵?为什么呢 恐怕我仍然不明白这里的动机...
Spacey 2012年

谢谢,我会再读一遍。同样,sigma的下标必须为“ n”,而不是“ i”。
Spacey 2012年

明天我会再记一些问题/评论,但现在,的“正式”名称是什么[R2X2[R充分[R简单

X1个X2
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.