2015-03-03 75 views
2

在matlab中计算三对角矩阵的逆时有没有办法关闭pivoting?我试图看看我解决三对角线系统时遇到的问题是否来自不旋转,我可以通过解决相同的系统和关闭旋转来在matlab中简单地测试它。任何帮助表示赞赏!matlab反转关闭枢轴

+1

对文档进行粗略浏览并不能揭示任何内容。但是,如果您使用['mldivide'](http://www.mathworks.com/help/matlab/ref/mldivide.html)来解决您的系统问题,MATLAB正在使用[Hessenberg求解器](http:///www.mathworks.com/help/matlab/ref/mldivide.html#bt42oms),无论如何可能无法做到透视。您可以尝试将您的结果与直接[三角网格求解器](http://en.wikipedia.org/wiki/Tridiagonal_matrix_algorithm)进行比较。 – TroyHaskin 2015-03-03 03:40:41

+0

@TroyHaskin:我建议你改变你的评论到一个答案。 – Jonas 2015-03-03 08:34:25

回答

0

mldivide的文档没有列出设置低级选项的任何选项。 我想象那是因为自动旋转不仅仅是期望的,而是现在大多数工具所期望的。

对于一个完整的三对角矩阵,MATLAB将使用它的Hessenberg解算器(我想它类似于this flow),对于稀疏三对角矩阵,将使用tridiagonal solver。在这两种情况下,都可以使用部分旋转来确保系统的精确解决方案。

为了解决MATLAB没有旋转的切换的事实,您可以实现自己的三角对角线求解器(参见上面的链接),而不需要旋转并查看解决方案如何受到影响。