icon
icon

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で疑似逆行列を求める方法について詳しく説明していくね!

大石ゆかり

お願いします!

 

目次

  1. 疑似逆行列とは
  2. 疑似逆行列の求め方
  3. 実際に書いてみよう
  4. まとめ

 

疑似逆行列とは

ある行列Bが「AB=BA=E」を満たす場合(Eは単位行列)、Aは正則行列であると言い、BをAの逆行列と呼びます。逆行列は行列の右上に「-1」を付けて表現します。逆行列は連立一次方程式の解を求める際に用いられます。

正則ではない行列は逆行列を持ちません。この場合連立一次方程式の解は不定(一意に定まらない)または不能(解が無い)となってしまいますが、逆行列の概念を一般化することで正則ではない行列に対しても逆行列のような行列を得ることができます。これをムーア・ペンローズの擬似逆行列と呼びます。

[PR] Pythonで挫折しない学習方法を動画で公開中

疑似逆行列の求め方

Pythonで線形代数を扱うには、Numpyモジュールのlinalgライブラリを使用します。linalgは線形代数(Linear Algebra)の略です。linalgについて詳しくは公式サイトも参考にしてください。

Linear algebra (numpy.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分か60分お好きな時間が選べて、かつ3回まで すべて無料で ご利用できます。
無理な勧誘は一切ない ので、お気軽にご参加ください。

今なら相談した方限定の割引・参加特典付き! 無料相談はこちら

まとめ

今回の記事ではPythonで疑似逆行列を求める方法を学習しました。

監修してくれたメンター

太田和樹(おおたかずき)

ITベンチャー企業のPM兼エンジニア。

普段は主に、Web系アプリケーション開発のプロジェクトマネージャーとプログラミング講師を行っている。守備範囲はフロントエンド、モバイル、サーバサイド、データサイエンティストと幅広い。その幅広い知見を生かして、複数の領域を組み合わせた新しい提案をするのが得意。

開発実績:画像認識技術を活用した駐車場混雑状況把握(実証実験)、音声認識を活用したヘルプデスク支援システム、Pepperを遠隔操作するアプリの開発、大規模基幹系システムの開発・導入マネジメント。

地方在住。仕事のほとんどをリモートオフィスで行う。通勤で消耗する代わりに趣味のDIYや家庭菜園、家族との時間を楽しんでいる。

 

大石ゆかり

内容分かりやすくて良かったです!

田島悠介

ゆかりちゃんも分からないことがあったら質問してね!

大石ゆかり

分かりました。ありがとうございます!

 

TechAcademyでは、初心者でも、Pythonを使った人工知能(AI)や機械学習の基礎を習得できる、オンラインブートキャンプを開催しています。

また、現役エンジニアから学べる無料体験も実施しているので、ぜひ参加してみてください。

初心者・未経験でもできる。まずはテックアカデミーに相談しよう

プログラミングを独学で学習していて、このように感じた経験はないでしょうか?

  • ・調べてもほしい情報が見つからない
  • ・独学のスキルが実際の業務で通用するのか不安
  • ・目標への学習プランがわからず、迷子になりそう

テックアカデミーでは、このような 学習に不安を抱えている方へ、マンツーマンで相談できる機会を無料で提供 しています。
30分間、オンラインでどんなことでも質問し放題です。

「受けてよかった」と感じていただけるよう カウンセラーやエンジニア・デザイナー があなたの相談に真摯に向き合います。

「自分に合っているか診断してほしい」
「漠然としているが話を聞いてみたい」

こんなささいな悩みでも大丈夫です。

無理な勧誘は一切ありません ので、まずはお気軽にご参加ください。
※体験用のカリキュラムも無料で配布いたします。(1週間限定)

今なら参加者限定の割引特典付き! 無料相談を予約する