Processing math: 100%

2018年1月28日日曜日

学習環境

線型代数入門(松坂 和夫(著)、岩波書店)の第5章(行列式)、7(積の行列式)、問題1.を取り組んでみる。


  1. (abcd-ba-dc-cda-b-d-cba)(a-b-c-dbad-cc-dabdc-ba)=(a2+b2+c2+d20000a2+b2+c2+d20000a2+b2+c2+d20000a2+b2+c2+d2)
    det(AAT)=(a2+b2+c2+d2)4det(AAT)=detAdetAT=detAdetA=(detA)2

    よって、

    (detA)2=(a2+b2+c2+d2)4detA=(a2+b2+c2+d2)2

    (aの4乘の係数は1。)

コード(Emacs)

Python 3

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

a, b, c, d = symbols('a, b, c, d')
A = Matrix([[a, b, c, d],
            [-b, a, -d, c],
            [-c, d, a, -b],
            [-d, -c, b, a]])

for t in [A, A.T, A * A.T, A.det(), A.det().factor()]:
    pprint(t)
    print()

入出力結果(Terminal, Jupyter(IPython))

$ ./sample1.py
⎡a   b   c   d ⎤
⎢              ⎥
⎢-b  a   -d  c ⎥
⎢              ⎥
⎢-c  d   a   -b⎥
⎢              ⎥
⎣-d  -c  b   a ⎦

⎡a  -b  -c  -d⎤
⎢             ⎥
⎢b  a   d   -c⎥
⎢             ⎥
⎢c  -d  a   b ⎥
⎢             ⎥
⎣d  c   -b  a ⎦

⎡ 2    2    2    2                                                         ⎤
⎢a  + b  + c  + d           0                  0                  0        ⎥
⎢                                                                          ⎥
⎢                    2    2    2    2                                      ⎥
⎢        0          a  + b  + c  + d           0                  0        ⎥
⎢                                                                          ⎥
⎢                                       2    2    2    2                   ⎥
⎢        0                  0          a  + b  + c  + d           0        ⎥
⎢                                                                          ⎥
⎢                                                          2    2    2    2⎥
⎣        0                  0                  0          a  + b  + c  + d ⎦

 4      2  2      2  2      2  2    4      2  2      2  2    4      2  2    4
a  + 2⋅a ⋅b  + 2⋅a ⋅c  + 2⋅a ⋅d  + b  + 2⋅b ⋅c  + 2⋅b ⋅d  + c  + 2⋅c ⋅d  + d 

                   2
⎛ 2    2    2    2⎞ 
⎝a  + b  + c  + d ⎠ 

$

0 コメント:

コメントを投稿