2015-03-03 70 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沒有旋轉的切換的事實,您可以實現自己的三角對角線求解器(參見上面的鏈接),而不需要旋轉並查看解決方案如何受到影響。