Matlab solution to Ax=b fails

3 visualizaciones (últimos 30 días)
Daniel Walsh
Daniel Walsh el 15 de Mzo. de 2016
Comentada: Daniel Walsh el 15 de Mzo. de 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.

Respuesta aceptada

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

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Linear Algebra en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by