数学 - Python - 解析学 - 行列式 - 行列式の他の性質(転置行列、積)

解析入門〈3〉

学習環境

• Surface 3 (4G LTE)、Surface 3 タイプ カバー、Surface ペン(端末)
• Windows 10 Pro (OS)
• Nebo(Windows アプリ)
• iPad Pro + Apple Pencil
• MyScript Nebo(iPad アプリ)
• 参考書籍
  • Pythonからはじめる数学入門
  • JavaScriptリファレンス 第6版
  • インタラクティブ・データビジュアライゼーション

解析入門〈3〉(松坂 和夫(著)、岩波書店)の第16章(行列式)、16.2(行列式の他の性質)、問題2.を取り組んでみる。

2. 
   
   $A^T=\left(\begin{array}{cccc}a& -b& -c& -d\\ b& a& d& -c\\ c& -d& a& b\\ d& c& -b& a\end{array}\right)$
   $A{A}^T=\left(a^2+b^2+c^2+d^2\right)\left(\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right)$
   $\begin{array}{}{\left(\det A\right)}^2={\left(a^2+b^2+c^2+d^2\right)}^4\\ \det A={\left(-1\right)}^2{\left(a^2+b^2+c^2+d^2\right)}^2\\ ={\left(a^2+b^2+c^2+d^2\right)}^2\end{array}$

コード(Emacs)

Python 3

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

init_printing(order='rev-lex')

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().factor(), A.T.det().factor()]:
    pprint(t)
    print()

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

$ ./sample3.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                                                         ⎤
⎢d  + c  + b  + a           0                  0                  0        ⎥
⎢                                                                          ⎥
⎢                    2    2    2    2                                      ⎥
⎢        0          d  + c  + b  + a           0                  0        ⎥
⎢                                                                          ⎥
⎢                                       2    2    2    2                   ⎥
⎢        0                  0          d  + c  + b  + a           0        ⎥
⎢                                                                          ⎥
⎢                                                          2    2    2    2⎥
⎣        0                  0                  0          d  + c  + b  + a ⎦

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

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

$