[求助]急求一道矩阵特征值分解问题
对于特征值分解的做法，Matlab里面=eig(A), A=V*D*V' ，但是如果M=P^T*D*P, P是M的特征矩阵，D是由其对应特征向量组成的对角矩阵。
怎么求得M=P^T*D*P？
= LDL(A) returns unit lower triangular matrix L, block diagonal D,
and permutation matrix P so that P'*A*P = L*D*L'. This is equivalent
to = LDL(A,'matrix').
试试这个函数。 不太对啊，实对称的矩阵对角化，V^-1==V',成立啊。分解后的是一个正交阵。
A=P'*D*P，则A'=(P'*D*P)'=P'*D*P=A，那么A一定是对称矩阵，在实数范围内，可以这样分解的只能实对称矩阵，而实对称矩阵的P是一个正交矩阵：P'=P^-1。所以原来的函数一样适用啊。
页:
[1]