Фон
Ковариационная матрица для вектора случайных величин содержит процедуру вычисления дисперсии любой линейной комбинации этих случайных величин. Правило заключается в том, что для любого вектора коэффициентов , X = ( X 1 , X 2 , … , X n ) ′ λ = ( λ 1 , … , λ n )AX=(X1,X2,…,Xn)′λ=(λ1,…,λn)
Var(λX)=λAλ′.(1)
Другими словами, правила умножения матриц описывают правила дисперсий.
Два свойства являются непосредственными и очевидными:A
Поскольку отклонения являются ожиданиями квадратов, они никогда не могут быть отрицательными. Таким образом, для всех векторов ,Ковариационные матрицы должны быть неотрицательно определенными.0 ≤ Var ( λ X ) = λ A λ ′ .λ
0≤Var(λX)=λAλ′.
Отклонения - это просто числа - или, если вы читаете матричные формулы буквально, они представляют собой матрицы . Таким образом, они не меняются, когда вы их транспонируете. Транспонирование дает Поскольку это верно для всех , должен быть равен его транспонированной : ковариационные матрицы должны быть симметричными.1×1(1)
λAλ′=Var(λX)=Var(λX)′=(λAλ′)′=λA′λ′.
λAA′
Более глубокий результат состоит в том, что любая неотрицательно определенная симметричная матрица является ковариационной матрицей. A Это означает, что на самом деле существует некоторая векторная случайная величина с качестве ее ковариации. Мы можем продемонстрировать это путем явного построения . Один из способов - заметить, что (многомерная) функция плотности со свойством имеет для своей ковариации. (Некоторая деликатность необходима, когда не является обратимым - но это просто техническая деталь.)XAXf(x1,…,xn)
log(f)∝−12(x1,…,xn)A−1(x1,…,xn)′
AA
Решения
Пусть и ковариационные матрицы. Очевидно, они квадратные; и если их сумма имеет какой-то смысл, они должны иметь одинаковые размеры. Нам нужно только проверить два свойства.XY
Сумма.
Я оставляю это как упражнение.
Этот хитрый. Один метод, который я использую, чтобы продумать сложные матричные задачи, состоит в том, чтобы сделать некоторые вычисления с матрицы Есть несколько распространенных, знакомых ковариационных матриц такого размера, таких как с и . Проблема заключается в том, что может не быть определенным: то есть может ли он давать отрицательное значение при вычислении дисперсии? Если это произойдет, то нам лучше иметь некоторые отрицательные коэффициенты в матрице. Это предполагает рассмотрение для . Чтобы получить что-то интересное, мы могли бы изначально стремиться к матрицам2×2
(abba)
a2≥b2a≥0XYX=(a−1−1a)
a≥1Y с различными структурами. На ум приходят диагональные матрицы, такие как с . (Обратите внимание, как мы можем свободно выбирать некоторые коэффициенты, такие как и , потому что мы можем масштабировать все записи в любой ковариационной матрице без изменения ее фундаментальных свойств. Это упрощает поиск интересных примеров.)Y=(b001)
b≥0−11
Я оставляю вам задачу вычислить и проверить, всегда ли она является ковариационной матрицей для любых допустимых значений и .XYab