correlation matrix 相關矩陣
一個穩態離散時間隨機過程可以一個時間序列來表示(M×1),而當我們把這個隨機時間序列乘以它的赫密特轉置矩陣(1×M),所得到的期望值(因為時間序列是隨機的)即是這個時間序列的相關矩陣(M×M)。
R = E[ u(n) uH(n) ]
R: correlation matrix
u(n) : stochastic process
uH(n) : Hermitian transpose
E: expectation
Hermitian transpose 赫密特轉置矩陣
將一個矩陣轉置並取共軛複數(complex conjugate)
在Matlab中,可直接以'符號求得Hermitian transpose,或使用ctranspose function
>> a = [1 1+j;2-j -2]
a =
1.0000 + 0.0000i 1.0000 + 1.0000i
2.0000 - 1.0000i -2.0000 + 0.0000i
>> a'
ans =
1.0000 + 0.0000i 2.0000 + 1.0000i
1.0000 - 1.0000i -2.0000 + 0.0000i
Hermitian matrix/self-adjoint matrix 埃爾米特矩陣/厄米特矩陣/自伴隨矩陣
共軛對稱的方塊矩陣,即 A = (AT)*,則A的Hermitian transpose等於A自己,即 AH = A
Example:
>> a = [1 3+j; 3-j -2]
a =
1.0000 + 0.0000i 3.0000 + 1.0000i
3.0000 - 1.0000i -2.0000 + 0.0000i
>> a'
ans =
1.0000 + 0.0000i 3.0000 + 1.0000i
3.0000 - 1.0000i -2.0000 + 0.0000i
>> isequal(a,a')
ans =
logical
1
>>
References:
Haykin, S. S. (2014). Adaptive filter theory. 5th edition, Pearson Education. pp 52-56.
Conjugate transpose (Wikipedia)
Hermitian matrix (Wikipedia) 埃爾米特矩陣(維基百科)
沒有留言:
張貼留言