Pythonで疑似逆行列を求める方法を現役エンジニアが解説【初心者向け】
初心者向けにPythonで疑似逆行列を求める方法について現役エンジニアが解説しています。疑似逆行列とは逆行列を持たない行列に対しても疑似的に逆行列を定義したものです。Pythonで線形代数を扱うには、Numpyモジュールのlinalgライブラリを使用します。
テックアカデミーマガジンは受講者数No.1のプログラミングスクール「テックアカデミー」が運営。初心者向けにプロが解説した記事を公開中。現役エンジニアの方はこちらをご覧ください。 ※ アンケートモニター提供元：GMOリサーチ株式会社 調査期間：2021年8月12日～8月16日 調査対象：2020年8月以降にプログラミングスクールを受講した18～80歳の男女1,000名 調査手法：インターネット調査
Pythonで疑似逆行列を求める方法について、TechAcademyのメンター（現役エンジニア）が実際のコードを使用して、初心者向けに解説します。
Pythonについてそもそもよく分からないという方は、Pythonとは何なのか解説した記事を読むとさらに理解が深まります。
なお本記事は、TechAcademyのオンラインブートキャンプ、Python講座の内容をもとに紹介しています。
今回は、Pythonに関する内容だね！
どういう内容でしょうか？
Pythonで疑似逆行列を求める方法について詳しく説明していくね！
お願いします！
目次
疑似逆行列とは
ある行列Bが「AB=BA=E」を満たす場合(Eは単位行列)、Aは正則行列であると言い、BをAの逆行列と呼びます。逆行列は行列の右上に「-1」を付けて表現します。逆行列は連立一次方程式の解を求める際に用いられます。
正則ではない行列は逆行列を持ちません。この場合連立一次方程式の解は不定(一意に定まらない)または不能(解が無い)となってしまいますが、逆行列の概念を一般化することで正則ではない行列に対しても逆行列のような行列を得ることができます。これをムーア・ペンローズの擬似逆行列と呼びます。
疑似逆行列の求め方
Pythonで線形代数を扱うには、Numpyモジュールのlinalgライブラリを使用します。linalgは線形代数(Linear Algebra)の略です。linalgについて詳しくは公式サイトも参考にしてください。
逆行列を求めるにはinv関数を使用します。
A_inv = np.linalg.inv(A)
また擬似逆行列を求めるには pinv 関数を使用します。pinvは擬似逆行列(pseudo-inverse of a matrix)の略です。
A_inv = np.linalg.pinv(A)
実際に書いてみよう
今回のサンプルプログラムでは、正則ではない行列に対して擬似逆行列を求めます。まずはもととなる行列を定義します。
import numpy as np A = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]]) A
実行結果は以下のとおりです。
array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]])
行列式を求めるにはlinalg.det()関数を使用します。
np.linalg.det(A)
実行結果は以下のとおりです。行列式=0のため、正則ではないことが確認できました。そのまま逆行列を求めてみましょう。
A_inv = np.linalg.inv(A) A_inv
実行結果は以下のとおりエラーとなりました。正則ではない(特異行列である)ためにエラーとなっていることが分かります。
LinAlgError: Singular matrix
次に擬似逆行列を求めてみましょう。
A_inv = np.linalg.pinv(A) A_inv
実行結果は以下のとおりです。エラーとならずに擬似逆行列を求めることができました。
array([[-6.38888889e-01, -1.66666667e-01, 3.05555556e-01],
[-5.55555556e-02, 3.36727575e-17, 5.55555556e-02],
[ 5.27777778e-01, 1.66666667e-01, -1.94444444e-01]])
実は私も！？独学で損する人の特徴
「スクールは高いから独学で成功する」という気持ちの方は多いと思います。
もちろんその方が金額は低く抑えられるでしょう。
ただ
独学には向き不向きがあり、実はスクールが向いている人も大勢います。
そんな方のために参考として、
テックアカデミー卒業生がスクールを選んだ理由
をご紹介します。
- ・困って挫折しそうなときに、質問や相談できる相手がいる環境で学んでいきたいなと思った
- ・わかった気になっているだけだったので、自分を追い込む環境に置いた方がいいと感じた
- ・スクールのカリキュラムで市場に求められるスキルを学ぶべきと思った
少しでも当てはまる部分があれば、
スクールが向いているかもしれません。
お試しのつもりで、まずは一度
無料相談
に参加してみませんか？
現役エンジニア・デザイナーに何でも気軽に相談できる30分
を
すべて無料で
できます。
無理な勧誘は一切ない
ので、お気軽にご参加ください。
まとめ
今回の記事ではPythonで疑似逆行列を求める方法を学習しました。
監修してくれたメンター
| 太田和樹（おおたかずき）
ITベンチャー企業のPM兼エンジニア。 普段は主に、Web系アプリケーション開発のプロジェクトマネージャーとプログラミング講師を行っている。守備範囲はフロントエンド、モバイル、サーバサイド、データサイエンティストと幅広い。その幅広い知見を生かして、複数の領域を組み合わせた新しい提案をするのが得意。 開発実績：画像認識技術を活用した駐車場混雑状況把握（実証実験）、音声認識を活用したヘルプデスク支援システム、Pepperを遠隔操作するアプリの開発、大規模基幹系システムの開発・導入マネジメント。 地方在住。仕事のほとんどをリモートオフィスで行う。通勤で消耗する代わりに趣味のDIYや家庭菜園、家族との時間を楽しんでいる。 |
内容分かりやすくて良かったです！
ゆかりちゃんも分からないことがあったら質問してね！
分かりました。ありがとうございます！
TechAcademyでは、初心者でも、Pythonを使った人工知能（AI）や機械学習の基礎を習得できる、オンラインブートキャンプを開催しています。
また、現役エンジニアから学べる無料体験も実施しているので、ぜひ参加してみてください。
プログラミングを独学していて、このように感じた経験はないでしょうか？
- ・調べてもほしい情報が見つからない
- ・独学のスキルが実際の業務で通用するのか不安
- ・目標への学習プランがわからず、迷子になりそう
テックアカデミーでは、このような
学習に不安を抱えている方へ、現役エンジニア講師とマンツーマンで相談できる機会を無料で提供
しています。
30分間、オンラインでどんなことでも質問し放題です。
「受けてよかった」と感じていただけるよう
厳しい試験を通過した講師
があなたの相談に真摯に向き合います。
「ただ気になることを相談したい」
「漠然としているがプロの話を聞いてみたい」
こんな気持ちでも大丈夫です。
無理な勧誘は一切ありません
ので、まずはお気軽にご参加ください。
※体験用のカリキュラムも無料で配布いたします。（1週間限定）
