Matlab solution to Ax=b fails

2 views (last 30 days)
Daniel Walsh
Daniel Walsh on 15 Mar 2016
Commented: Daniel Walsh on 15 Mar 2016
I have encountered the following problem: I need to solve a large linear system Ax=b, with complex, square A. I have imported A and b into Matlab from a file, and they are in Matlab's sparse matrix format. I then write x = A/b; and Matlab produces the solution for x. Now I check that the answer x satisfies the original equation: norm(A*x-transpose(b)); The transpose() has been applied to b because as far as I can tell Matlab returns x opposite to the mathematical convention (i.e. a row vector rather than a column vector).
When I run this check, the residual (norm of the difference) is on the order of 10^3, indicating that the solution is completely incorrect. Why is Matlab giving the wrong answer?
Any help would be appreciated, thank you.

Accepted Answer

Steven Lord
Steven Lord on 15 Mar 2016
A/b ( MRDIVIDE ) and A\b ( MLDIVIDE ) are two different operations.
  1 Comment
Daniel Walsh
Daniel Walsh on 15 Mar 2016
Yes, this was exactly the problem. Thank you for pointing that out.

Sign in to comment.

More Answers (0)

Categories

Find more on Linear Algebra in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!