Cholesky矩陣實現相關變數和不相關變數的轉換
阿新 • • 發佈:2019-02-02
給定一個協方差矩陣
其中
幾何上講,Cholesky矩陣將無相關的變數轉換為由
在高斯過程的例子中,我們通過下面的形式形成高斯分佈資料點:
n1 = 80;
n2 = 40;
S1 = eye(2);
S2 = [1 0.95 ;0.95 1];% 兩個協方差矩陣
m1 = [0.75;0];
m2 = [-0.75;0]; %兩個均值
x1 = bsxfun (@plus,chol(S1)'*gpml_randn(0.2,2,n1),m1);
x2 = bsxfun(@plus,chol(S2)'*gpml_randn(0.3,2,n2),m2);
x = [x1 x2]';
y = [-ones(1,n1) ones(1,n2)]';
plot(x1(1,:),x1(2,:),'b+');
hold on;
plot(x2(1,:),x2(2,:),'r+');
% check
cov1 = cov(x1');
cov2 = cov(x2');