我有一个真实的不可对角矩阵,我希望尽可能整齐地分解。我希望将其以Jordan的正常形式表示,但是由于从数字上来说这是有问题的,因此我正在寻找下一个最好的东西。我发现有FORTRAN和MATLAB例程可以对矩阵进行块对角Schur分解。 SLICOT中的FORTRAN实现为MB03RD
,而MATLAB实现为bdschur
(据我所知,它只是MB03RD
的包装)。
我的计算机上没有MATLAB,而生成矩阵的代码在Python中,因此我正在寻找Python中的等效函数。 Python控制系统库indicated的旧文档计划了对bdschur
的仿真,但并未在current docs中的任何地方显示。 Slycot存储库中有一个MB03RD
的FORTRAN file,但是我似乎找不到太多有关Slycot的文档,并且当我导入它时,几乎没有任何函数看起来像是被包装为Python函数。>
我想知道是否有人知道在Python中调用等效例程的方法,或者是否存在在Python中具有实现的其他类似分解。