2018年11月16日金曜日

学習環境

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



    1. x , y , z = 2 , 1 , 1 + t 1 , - 2 , 0 + s 2 , 0 , - 2 x = 2 + t + 2 s y = 1 - 2 t z = 1 - 2 s t = 1 - y 2 s = 1 - z 2 x = 2 + 1 - y 2 + 1 - z 2 x = 4 + 1 - y + 2 - 2 z 2 x + y + 2 z = 7

    2. x , y , z = - 2 , 3 , - 1 + s 4 , - 1 , 4 + t - 2 , - 4 , 2 x = - 2 + 4 s - 2 t y = 3 - s - 4 t z = - 1 + 4 s + 2 t x + z = - 3 + 8 s 2 x - y = - 7 + 9 s 9 x + 9 z = - 27 + 725 16 x - 8 y = - 56 + 72 s 7 x - 8 y - 9 z = - 29

    3. x , y , z = - 5 , - 1 , 2 + s 6 , 3 , - 3 + t 8 , 0 , 0 x = - 5 + 6 s + 8 t y = - 1 + 3 s z = 2 - 3 s y + z = 1

コード(Emacs)

Python 3

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

ps = [((2, 1, 1), (3, -1, 1), (4, 1, -1)),
      ((-2, 3, -1), (2, 2, 3), (-4, -1, 1)),
      ((-5, -1, 2), (1, 2, -1), (3, -1, 2))]

x = Matrix(symbols('x, y, z'))
s, t = symbols('s, t')
for i, (u, v, w) in enumerate(ps):
    print(f'({chr(ord("a") + i)})')
    a = Matrix(v) - Matrix(u)
    b = Matrix(w) - Matrix(u)
    eq = Matrix(u) + s * a + t * b - x
    st = solve(eq, dict=True)[0]
    s0, t0 = st[s], st[t]
    pprint(eq.subs({s: s0, t: t0}))

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

$ ./sample8.py
(a)
⎡     y       7⎤
⎢-x - ─ - z + ─⎥
⎢     2       2⎥
⎢              ⎥
⎢      0       ⎥
⎢              ⎥
⎣      0       ⎦
(b)
⎡     8⋅y   9⋅z   29⎤
⎢-x + ─── + ─── - ──⎥
⎢      7     7    7 ⎥
⎢                   ⎥
⎢         0         ⎥
⎢                   ⎥
⎣         0         ⎦
(c)
⎡    0     ⎤
⎢          ⎥
⎢-y - z + 1⎥
⎢          ⎥
⎣    0     ⎦
$

0 コメント:

コメントを投稿

関連コンテンツ