matrix - 行列逆コード (c、c++、fortran または tksolver)

Question

行操作を使用して nxn 行列 A の逆数を決定するコードが必要です。私はこのコードを書くのにばかげた時間を過ごしています。私は大学院生で、何年もプログラミングの授業を受けていません。ありがとう。

Answer 1

あなたは本当に逆を望んでいないかもしれません。連立方程式を解こうとしている場合は、LU 分解を使用したほうがよいでしょう。

このアプリをどの言語で書きたいかは言いません。Java には Apache Commons Math があります。Python には NumPy があります。FORTRAN には LinPack があります。言語を選択してライブラリを使用します。自分で書かないでください。

Answer 2

何年もプログラミングクラスを受講していない場合は、行列演算にScilabまたは MATLAB をお勧めします。

1 回限りの計算の場合は、WolframAlpha を見てください。

ライブラリを使用することが C++ のオプションである場合は、Armadilloを参照してください。

Answer 3

  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 4

Numerical Recipies Onlineを見てください。Fortran 77、Ansi C、および Fortran 90 のコードは無料で、最適な方法についての手がかりがたくさん得られます。