2017年10月10日火曜日

学習環境

ラング線形代数学(上)(S.ラング (著)、芹沢 正三 (翻訳)、ちくま学芸文庫)の1章(R^n におけるベクトル)、5(直線と平面)、練習問題15.を取り組んでみる。


    1. ( 1,1,1 )·( 1,1,1 )= 1+1+1 1+1+1 cosθ cosθ= 111 3 = 1 3

    2. ( 2,3,1 )·( 1,1,1 )= 4+9+1 1+1+1 cosθ cosθ= 231 14 3 = 2 42

    3. ( 1,2,1 )·( 1,3,1 )= 1+4+1 1+9+1 cosθ cosθ= 1+61 6 · 11 = 4 66

    4. ( 2,1,1 )·( 1,1,1 )= 4+1+1 1+1+1 cosθ cosθ= 21+1 6 3 = 2 3 2 = 2 3

コード(Emacs)

Python 3

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

print('15.')
x = symbols('x')
t = [((1, 1, 1), (1, -1, -1)),
     ((2, 3, -1), (1, -1, 1)),
     ((1, 2, -1), (-1, 3, 1)),
     ((2, 1, 1), (-1, -1, 1))]

for i, (a, b) in enumerate(t):
    print(f'({chr(ord("a") + i)}).')
    A = Matrix(a)
    B = Matrix(b)
    eq = A.dot(B)  - A.norm() * B.norm() * x
    for t in [A.T, B.T, eq, solve(eq)]:
        pprint(t)
        print()
    print()

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

$ ./sample15.py
15.
(a).
[1  1  1]

[1  -1  -1]

-3⋅x - 1

[-1/3]


(b).
[2  3  -1]

[1  -1  1]

-√42⋅x - 2

⎡-√42 ⎤
⎢─────⎥
⎣  21 ⎦


(c).
[1  2  -1]

[-1  3  1]

-√66⋅x + 4

⎡2⋅√66⎤
⎢─────⎥
⎣  33 ⎦


(d).
[2  1  1]

[-1  -1  1]

-3⋅√2⋅x - 2

⎡-√2 ⎤
⎢────⎥
⎣ 3  ⎦


$

0 コメント:

コメントを投稿