Matlab協方差矩陣的計算原理

2022-07-13 18:30:17 字數 1028 閱讀 8717

a =-

112-

2314

03fori=1

:size(a,2) 

forj=1

:size(a,2) 

c(i,j)

=sum((a(:,i)

-mean(a(:,i))).

*(a(:,j)

-mean(a(:,j))))

/(size(a,1)

-1);end 

endc =

10.3333   -4.1667    3.0000

-4.1667    2.3333   -1.5000

3.0000   -1.5000    1.0000 

c為求得的協方差矩陣,在matlab以矩陣a的每一列為變數,對應的每一行為樣本。這樣在矩陣a中就有3個列變數分別為a(:,1), a(:,2), a(:,3)。

在協方差矩陣c中,每乙個元素c(i,j)為對第i列與第j列的協方差,例如c(1,2) = -4.1667為第一列與第二列的協方差。

拿c(1,2)的求解過程來說

=sum((a(:,1)

-mean(a(:,1))).

*(a(:,2)

-mean(a(:,2))))

/(size(a,1)

-1);-

mean(a(:,1)),第一列的元素減去該列的均值得到

-1.3333

-2.3333

3.6667

-mean(a(:,2)),第二列的元素減去該列的均值得到

-0.3333

1.6667

-1.3333

3, 再將第一步與第二部的結果相乘

-1.3333        -0.3333           0.4444

-2.3333  .*     1.6667  =     -3.8889

3.6667         -1.3333          -4.8889

4, 再將結果求和/size(a,1)-1 得 -4.1667,該值即為c(1,2)的值。

協方差矩陣計算

cov x,y exy ex ey 協方差的定義,ex為隨機變數x的數學期望,同理,exy是xy的數學期望 舉例 xi 1.1 1.9 3 yi 5.0 10.4 14.6 e x 1.1 1.9 3 3 2 e y 5.0 10.4 14.6 3 10 e xy 1.1 5.0 1.9 10.4 ...

詳解協方差與協方差矩陣計算

協方差的定義 對於一般的分布,直接代入e x 之類的就可以計算出來了,但真給你乙個具體數值的分布,要計算協方差矩陣,根據這個公式來計算,還真不容易反應過來。網上值得參考的資料也不多,這裡用乙個例子說明協方差矩陣是怎麼計算出來的吧。記住,x y 是乙個列向量,它表示了每種情況下每個樣本可能出現的數。比...

協方差 協方差矩陣

期望 離散型隨機變數的一切可能的取值xi與對應的概率pi xi 之積的和稱為該離散型隨機變數的數學期望 設級數絕對收斂 記為 e x 隨機變數最基本的數學特徵之一。它反映隨機變數平均取值的大小。又稱期望或均值。求法 設離散型隨機變數x的取值為 方差 方差是各個資料與平均數之差的平方的平均數。在概率論...