covariance QR 矩阵

covariance matrix是对称矩阵，用QR分解，逆序相乘可以很快求得近似解
对于很大的矩阵没有必要通过解线性方程得到精确解，可以用收敛比较快的方法得到相
对精确的近似解，比如用matlab的QR函数

a=[.....];

for i=1:10
[q,r]=qr(a);
a=r*q;
end

target = qr(a)


太大的矩阵不知道有效果没。。。