Test of whether matrix is Symmetric Positive Definite is giving wrong result when matrix is not symmetric

3 次查看（过去 30 天）
Bill Tubbs
Bill Tubbs 2022-11-8
编辑： Bruno Luong 2022-11-9
I'm computing a multivariate normal probability density with an estimated covariance matrix as follows:
% Update the conditional likelihood given the data
p_yk_g_seq_Ykm1(j) = mvnpdf(yk, ykp1_est, Sk);
and getting the following error
Error using mvnpdf (line 127)
SIGMA must be a square, symmetric, positive definite matrix.
Here are some checks I did in the debugger when this error occurred:
K>> Sk
Sk =
0.0540 -0.0001
-0.0001 0.0540
K>> issymmetric(Sk)
ans =
logical
0
Clearly it is not symmetric.
But when I tried to check if Sk is symmetic positive definite using the method described here in the documentaion:
K>> try chol(Sk)
disp('Matrix is symmetric positive definite.')
catch ME
disp('Matrix is not symmetric positive definite')
end
ans =
0.2323 -0.0002
0 0.2323
Matrix is symmetric positive definite.
This confused me. What is going on here? Is it non-symmetric or not PD or both?
Looking in the code for mvnpdf it actually uses this check:
[R,err] = cholcov(Sigma,0);
So maybe the method from the documentation page linked above is wrong or out of date?
Or maybe it tests for positive semi-definite but not symmetric positive semi-definite.

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

采纳的回答

Steven Lord
Steven Lord 2022-11-8
From the chol documentation page: "If A is nonsymmetric , then chol treats the matrix as symmetric and uses only the diagonal and upper triangle of A." Is the symmetric matrix generated using the diagonal and upper triangle of your matrix SPD?
I'll ask the documentation staff to take a look at the page to which you linked.
  4 个评论
显示 2更早的评论
Matt J
Matt J 2022-11-9
The definition of PD and PSD can be extended to non-symmetric matrices, but all of the interesting theorems connecting eigenvalues and determinants to positive-definiteness hold only for symmetric\Hermitian matrices.

请先登录，再进行评论。

更多回答（1 个）

Matt J
Matt J 2022-11-8
编辑：Matt J 2022-11-9
But when I tried to check if Sk is symmetic positive definite using the method described here in the documentaion:
All of the tests at this link are predicated on the assumption that Sk is already symmetric. In other words, they are really intended as tests of postive or non-negative definiteness, not of symmetry.
Therefore, you should just symmetrize the matrix and be done with it,
Sk=(Sk+Sk.')/2;
This does not change the status of Sk as a PSD or PD matrix.
  3 个评论
显示 1更早的评论
Bill Tubbs
Bill Tubbs 2022-11-9
Alternatively, the page could be amended to provide the following test which seems to check both symmetric and positive definite (I assume so since this is what raised the error in my case):
[~,err] = cholcov(Pk,0)
assert(err == 0)

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by