矩阵逆码(c, c++, fortran 或 tksolver)

Question

我需要一个代码来使用行运算确定 nxn 矩阵 A 的逆。 我写这段代码的时间很荒谬。一切都有帮助，我是一名研究生，已经很多年没有上过编程课了。谢谢。

Answer 1

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

你没有说你想用什么语言来编写这个应用程序。Java 有 Apache Commons Math；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 代码是免费的，您可以获得大量有关如何最好地进行操作的线索。