Cannot create a matrix using pinv.

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?

[Matt J]

I suggest that you attach the matrix in a .mat file, so that we can examine it.

Answer 1

[Christine Tobler]

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.

[은영 차]

Hello.

This is matrix .mat file direct rink -> H_sum.mat.

If it doesn't work, you can use "https://drive.google.com/file/d/19ImNfXJoTIPxQPE_PYG18Fsxyv4Q7lrC/view?usp=sharing"

I can get inverse matrix of that matrix, I don't know why I can't get the pseudo-inverse one.

Thank you

[은영 차]

It occurs '2019a' also '2018b'

[Christine Tobler]

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]

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.