2018年8月4日土曜日

学習環境

ラング線形代数学(下)(S.ラング (著)、芹沢 正三 (翻訳)、ちくま学芸文庫)の9章(多項式と行列)、4(特性多項式)、練習問題3-(a).を取り組んでみる。



    1. 特性多項式。

      det t - 1 - i - i t + 2 = t 2 + t - 2 + 1 = t 2 + t - 1

      固有値。

      - 1 ± 1 + 4 2 = - 1 ± 5 2

      国有ベクトル。

      - 3 ± 5 2 - i - i 3 ± 5 2 x y = 0 0 - 3 ± 5 2 x - y i = 0 - x i + 3 ± 5 2 y = 0 x = 2 , y = 3 5 i 2 , 3 5 i

コード(Emacs)

Python 3

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

print('3.')

t, x, y = symbols('t, x, y')

A = Matrix([[1, I],
            [I, -2]])
B = Matrix([[t, 0],
            [0, t]])
C = B - A
D = C.det()
ts = solve(D, t)
X = Matrix([[x],
            [y]])
for s in [A, B, C, D, ts]:
    pprint(s)
    print()

for t0 in ts:
    E = C.subs({t: t0}) * X
    a, b = E
    pprint(solve(a - b))

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

$ ./sample3.py
3.
⎡1  ⅈ ⎤
⎢     ⎥
⎣ⅈ  -2⎦

⎡t  0⎤
⎢    ⎥
⎣0  t⎦

⎡t - 1   -ⅈ  ⎤
⎢            ⎥
⎣ -ⅈ    t + 2⎦

(t - 1)⋅(t + 2) + 1

⎡  1   √5    √5   1⎤
⎢- ─ + ──, - ── - ─⎥
⎣  2   2     2    2⎦

⎡⎧   y⋅(√5 + 3 + 2⋅ⅈ)⎫⎤
⎢⎨x: ────────────────⎬⎥
⎣⎩    -3 + √5 + 2⋅ⅈ  ⎭⎦
⎡⎧   y⋅(-3 + √5 - 2⋅ⅈ)⎫⎤
⎢⎨x: ─────────────────⎬⎥
⎣⎩      √5 + 3 - 2⋅ⅈ  ⎭⎦
$

0 コメント:

コメントを投稿