Why doesnt chol work for hermitian positive semi-definite matrix?

Question

Nikhil Challa 2022-9-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1811895-why-doesnt-chol-work-for-hermitian-positive-semi-definite-matrix

回答： Pravarthana P 2022-9-28

Hello,

I am trying to perform factorization of form A = R'*R using matlab function 'chol' but A is not positive definite, rather its hermitian positive semi-definite. As per theory from https://en.wikipedia.org/wiki/Hermitian_matrix , I should be able to perform factorization but matlab expects A to be positive definite matrix. Is there an alternative function in matlab?

Thanks,

3 个评论
显示 1更早的评论隐藏 1更早的评论

Nikhil Challa 2022-9-25

@Geoff Hayes , the https://en.wikipedia.org/wiki/Cholesky_decomposition link also talks about decomposition for hermitian positive semi-definite matrix except the solution may not be unique. Matlab has lot of functions that provide solution even for non-unique cases, so I would expect matlab to do the same for 'chol' function also.

Let me know if my understanding is incorrect.

Torsten 2022-9-25

https://de.mathworks.com/matlabcentral/answers/477584-problem-with-cholesky-s-decomposition-of-a-positive-semi-definite-matrix

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Pravarthana P 2022-9-28

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1811895-why-doesnt-chol-work-for-hermitian-positive-semi-definite-matrix#answer_1062295

在 MATLAB Online 中打开

Hi Nikil Challa,

I understand that you are trying to perform factorization using “chol” function and are facing an issue with positive semi-definite matrix.

The “chol” function expects a symmetric positive definite matrix as input as mentioned in the documentation link.

The following can likely be a workaround to factorize positive semi-definite matrix:

Compute the LDL decomposition: A lower-triangular matrix L, a block-diagonal matrix D (1-by-1 and 2-by-2 blocks), and a permutation matrix P, such that A is P*L*D*L'*P' :

[L, D, P] = ldl(A);

2.Compute the eigenvalue decomposition, set its negative eigenvalues to zero, and use QR to get two triangular factors for this modified matrix:

[U, d] = eig(A, 'vector'); % A == U*diag(d)*U'
d(d < 0) = 0;  % Set negative eigenvalues of A to zero
[~, R] = qr(diag(sqrt(d))*U');  
% Q*R == sqrt(D)*U', so A == U*diag(d)*U'== R'*Q'*Q*R == R'*R
% In this case, check that d(d<0) are all nearly zero.