Loading [MathJax]/jax/output/CommonHTML/jax.js

2019年2月18日月曜日

学習環境

ラング線形代数学(上)(S.ラング (著)、芹沢 正三 (翻訳)、ちくま学芸文庫)の3章(行列)、2(行列の乗法)、練習問題21の解答を求めてみる。



    1. trAB=ni=1(nk=1aikbki)trBA=ni=1(nk=1bikaki)=ni=1(nk=1akibik)=ni=1(nk=1aikbki)

      よって、

      trAB=trBA

    2. B-1=(cij)

      とおく。
      このとき、(i, j)成分 について、

      BB-1=I(nk=1bikckj)=Ii=jnk=1bikckj=nk=1bikcki=1ijnk=1bikckj=0

      また、

      B-1AB=(nk=1cikakj)B=(nl=1(nk=1cikakl)blj)=(nl=1(nk=1bljcikakl))tr(B-1AB)=ni=1(nl=1(nk=1blicikakl))=ni=1(nl=1(nk=1bilclkaki))=ni=1(nl=1bilcliaii)=ni=1(nl=1bilclk)aii=ni=1aii

      よって、行列 B が可逆ならば、

      tr(B-1AB)=trA

      が成り立つ。

      (証明終)

コード

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, randMatrix

print('21.')

for n in range(1, 6):
    A = randMatrix(n)
    B = randMatrix(n)
    try:
        for t in [A, B, (A * B).trace() == (B * A).trace(),
                  (B**-1 * A * B).trace() == A.trace()]:
            pprint(t)
            print()
        print()
    except Exception as err:
        print(type(err), err)

入出力結果(cmd(コマンドプロンプト)、Terminal、Jupyter(IPython))

C:\Users\...> py -3 sample21.py
21.
[66]

[69]

True

True


⎡95  54⎤
⎢      ⎥
⎣35  17⎦

⎡24  42⎤
⎢      ⎥
⎣62  36⎦

True

True


⎡56  80  19⎤
⎢          ⎥
⎢13  29  80⎥
⎢          ⎥
⎣77  18  74⎦

⎡91  13  29⎤
⎢          ⎥
⎢90  66  10⎥
⎢          ⎥
⎣20  44  9 ⎦

True

True


⎡55  60  31  63⎤
⎢              ⎥
⎢31  63  82  3 ⎥
⎢              ⎥
⎢65  87  79  73⎥
⎢              ⎥
⎣54  68  99  56⎦

⎡98  33  73  49⎤
⎢              ⎥
⎢26  96  9   44⎥
⎢              ⎥
⎢81  65  11  86⎥
⎢              ⎥
⎣96  6   66  26⎦

True

True


⎡0   45  99  94  11⎤
⎢                  ⎥
⎢74  2   90  15  1 ⎥
⎢                  ⎥
⎢0   95  44  9   0 ⎥
⎢                  ⎥
⎢52  26  26  85  90⎥
⎢                  ⎥
⎣33  21  59  48  41⎦

⎡67  49  38  18  30⎤
⎢                  ⎥
⎢82  55  91  41  26⎥
⎢                  ⎥
⎢74  15  50  22  98⎥
⎢                  ⎥
⎢64  30  24  4   57⎥
⎢                  ⎥
⎣82  0   68  79  8 ⎦

True

True


C:\Users\...>

0 コメント:

コメントを投稿