matrix - 矩阵逆码（c、c++、fortran 或 tksolver）

Question

我需要一个代码来确定使用行操作的 nxn 矩阵 A 的逆矩阵。我写这段代码的时间很可笑。有什么帮助，我是一名研究生，多年来没有上过编程课。谢谢。

Answer 1

您可能并不真正想要相反的情况。如果您尝试求解方程组，最好使用 LU 分解。

你没有说你想用什么语言编写这个应用程序。Java 有 Apache Commons Math；Python 有 NumPy；FORTRAN 有 LinPack。选择一种语言并使用一个库；不要自己写。

Answer 2

  SUBROUTINE MATINV(A,N)
  DIMENSION A(N,N)
  DO 1 I=1,N
     Z=A(I,I)
     A(I,I)=1.0
     DO 2 J=1,N
2       A(I,J)=A(1,J)/Z
     DO 1 K=1,N
        IF (K-I) 3,1,3
3          Z=A(K,I)
        A(K,I)=0.0
        DO 4 J=1,N
4          A(K,J)=A(K,J)-Z*A(I,J)
1 CONTINUE
  RETURN
  END

Answer 3

如果您多年没有上过编程课，我建议您使用Scilab或 MATLAB 进行矩阵运算。

如果是一次性计算，请查看 WolframAlpha。

如果使用库是 C++ 的一个选项，请查看Armadillo。

Answer 4

在线查看数字食谱。Fortran 77、Ansi C 和 Fortran 90 代码是免费的，您可以获得关于如何最好地进行的大量线索。