Finding inv(A) for Ax=b system
1 次查看(过去 30 天)
显示 更早的评论
I tried 3 methods to solve the system for inv(A), where A is highly sparse matrix spy(A) is given below:
1) pinv(A)......... Matlab solve it in 190 sec without any warnings.
2) A\b............. Matlab took only 10 sec but gives warning "Matrix is singular to working precision"
3) [L,U]=lu(A); inv(L)*inv(U) .......... Matlab took 50 sec but give the same warning as in 2nd method.
Is there anyway, I can get inv(A) without warnings.
3 个评论
Jan
2017-7-27
编辑:Jan
2017-7-27
@Asif Arshid: What is your question? You observed that the slash operator has a different sensitivity to detect near to singular matrices than lu and inv(L)*inv(U). What is the condition number of the matrix? Do you think that slash is to pessimistic or the lu method too sloppy? Or are you surprised by the speed of the slash operator?
回答(1 个)
Star Strider
2017-7-27
4 个评论
M.Shaarawy
2019-5-20
Is there regularized parameterized trust region sub problem (RPTRS) in MATLAB to solve this kind of problem?
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!