Cannot create a matrix using pinv.

조회 수: 1 (최근 30일)
은영 차
은영 차 2021년 1월 19일
댓글: Christine Tobler 2021년 1월 20일
I want to find the inverse function of the MxN matrix using the pinv function.
However, a message indicating that the creation of the matrix has failed is displayed, which is as follows."SVD values do not converge"
The matrix is a full rank matrix.
Is this why I can't build a pseudo inverse matrix?Or is there any other reason that I can't create a pseudo inverse matrix?
  댓글 수: 1
Matt J
Matt J 2021년 1월 19일
I suggest that you attach the matrix in a .mat file, so that we can examine it.

댓글을 달려면 로그인하십시오.

답변 (1개)

Christine Tobler
Christine Tobler 2021년 1월 19일
You've hit on a rare matrix that the SVD algorithm isn't able to work with (the SVD being the factorization called withing the PINV function). Please attach the matrix here or send it to Tech Support as a bug report. If you can also tell us your MATLAB release and operating system, that would be very helpful.
  댓글 수: 4
Christine Tobler
Christine Tobler 2021년 1월 20일
Hi,
I just tried to reprocude this and didn't see the error in either 2019a or 2018b. This is likely to be machine specific.
Could you tell me what operating system you are using and also enter the line
>> version -lapack
in MATLAB and send me the printout that MATLAB returns? This should help me to try and reproduce this.
The reason you can get the inverse but not the pseudoinverse is that different decompositions are used to compute these two. The pseudoinverse is the more expensive algorithm. Since the matrix here is very well conditioned, both the inverse and the pseudoinverse will be the same, so you would be save just using INV here.
Christine Tobler
Christine Tobler 2021년 1월 20일
A colleague was able to reproduce this (using Windows on an Intel machine with AVX2 instructions). He could reproduce this in R2019a, but not anymore in R2019b. So this is fixed from R2019b onward.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Linear Algebra에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by