Câu hỏi về ma trận hiệp phương sai của 2 tín hiệu không gian


9

Mỗi khi tôi nghĩ rằng tôi đã hiểu được ma trận hiệp phương sai, một người khác sẽ đưa ra một công thức khác nhau.

Tôi hiện đang đọc bài viết này:

J. Benesty, "Thuật toán phân rã eigenvalue thích ứng cho nội địa hóa nguồn âm thụ động" , J. Acoust. Sóc. Là. Tập 107 , Số 1, trang 384-391 (2000)

và tôi đã bắt gặp một công thức mà tôi không hiểu lắm. Ở đây, tác giả đang xây dựng ma trận hiệp phương sai giữa hai tín hiệu x 2 . Hai tín hiệu đó là từ các cảm biến khác nhau.x1x2

Đối với ma trận hiệp phương sai của một tín hiệu, tôi biết rằng chúng ta có thể lấy được nó bằng cách tính ma trận hồi quy, sau đó nhân nó với Hermiti của cùng một ma trận và chia cho , độ dài của vectơ gốc. Kích thước của ma trận hiệp phương sai ở đây có thể tùy ý, với kích thước tối đa là N × N .NN×N

Đối với ma trận hiệp phương sai của hai tín hiệu không gian, nếu chúng ta đặt tín hiệu đầu tiên ở hàng đầu tiên và tín hiệu thứ hai ở hàng thứ hai của ma trận, sau đó nhân với Hermiti của nó và cũng chia cho , thì chúng ta có được 2 × 2 ma trận hiệp phương sai của cả hai tín hiệu không gian.N2×2

Tuy nhiên, trong bài báo này, tác giả tính toán những gì trông giống như bốn ma trận, R 2R11,R12,R21 , và sau đó đặt chúng vào một siêu ma trận và gọi đó là ma trận hiệp phương sai.R22

Tại sao cái này rất? Đây là một hình ảnh của văn bản:

nhập mô tả hình ảnh ở đây

Câu trả lời:


6

x1[n]x2[n]N

  1. n=1Nxi[n]xj[n], i,j{1,2}2×2

    R2×2=[σ12CCσ22].
    x1[n]σ12=R1,1R2,2=σ22x2[n]R1.2=R2,1=Cnn1n+1
  2. x1[n]x1[m]x1[n]x2[m]2N2N×2NRfull=E[XXT]

    X=(x1[1],x1[2],,x1[N],x2[1],x2[2],,x2[N])T=(x1,x2)T
    Rfull=[Rx1,x1Rx1,x2Rx2,x1Rx2,x2]
    Rxi,xj=E[xixjT]Rxi,xj(xi[1],xi[2],,xi[N])(xj[1],xj[2],,xj[N])iji=jn=m
    RfullRsimple=[σ12ICICIσ22I]
    IN×Nσ12,σ22C2N×2NRfull2×2R2×22N×2NRfullσ12,σ22CRfullRsimpleR2×2

Cảm ơn. Đầu tiên, không nên sigma trong (1) nói từ n = 0 đến N-1? (Không phải từ i = 1 đến n).
Spacey

Tôi không chắc là tôi vẫn hiểu những gì / tại sao chúng ta đang làm theo cách này. Bạn có nói rằng với (1), vì các tạp âm trong cả hai vectơ hoàn toàn độc lập với nhau, chúng ta phải sử dụng phương pháp đó và do đó có được ma trận đồng phương 2x2, nhưng trong trường hợp thứ hai (2), vì các tạp âm trong các vectơ không độc lập, chúng ta phải ghép cả hai vectơ và sau đó tính ma trận đồng phương của chúng? Tại sao mặc dù? Tôi sợ tôi vẫn không hiểu động lực ở đây ...
Spacey

Cảm ơn tôi sẽ đọc lại. Ngoài ra, chỉ số cho sigma phải là 'n', không phải 'i'.
Spacey

R2x2,RfullRsimple

x1x2
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.