この記事では、Python言語とNumPyを用いて疑似逆行列を求める方法をソースコード付きで解説します。
疑似逆行列
正則でない行列の逆行列は計算できませんが、擬似逆行列を求めることは出来ます。
Pythonでは、NumPyのnumpy.linalg.pinvを利用することで、疑似逆行列を求めることが出来ます。
ソースコード
サンプルプログラムのソースコードです。
# -*- coding: utf-8 -*- import numpy as np def main(): A = np.array([[1.,1.] # 行列Aの生成 ,[1.,1.]]) pinvA = np.linalg.pinv(A) # Aの擬似逆行列 print( "invA=\n" + str(pinvA) ) # 計算結果の表示 if __name__ == '__main__': main()
実行結果
サンプルプログラムの実行結果です。
invA= [[ 0.25. 0,25 ] [ 0.25 0.25]]
コメント