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的取值為 方差 方差是各個資料與平均數之差的平方的平均數。在概率論...