invert nearly singular matrix

4 次查看(过去 30 天)
Peter
Peter 2012-5-22
I am trying to solve the poisson equation in cylindrical coordinates using a finite difference which essentially narrows down to Ax = B where A is a finite difference matrix, x and B are column vectors. I am solving for x and unfortunately A is nearly singular. So when I do x = A\B i get inaccurate answers. Is there a way to accurately solve this issue?

回答(1 个)

Jan
Jan 2012-5-22
Some instable systems caused by nearly singular matrices can be solved sufficiently using quadrupel precision. See INTLAB of Rump, http://www.ti3.tu-harburg.de/rump/Research/topics.php.

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by