2017年10月8日日曜日

学習環境

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


  1. 2平面の法線ベクトルはそれぞれ、(2, -1, 1)、(3, 1, 1)。この2つのベクトルと垂直なベクトルを求める。

    ( 2,1,1 )·( x,y,z )=0 ( 3,1,1 )·( x,y,z )=0 2xy+z=0 3x+y+z=0 x+2y=0 x=2y 4yy+z=0 z=5y ( 2y,y,5y ) =y( 2,1,5 )

    よって、求める問題の2平面の交線に平行なベクトルはt(-2, 1, 5)。


  2. 2平面の法線ベクトルはそれぞれ、(2, 1, 5)、(3, -2, 1)。この2つのベクトルと垂直なベクトルを求める。

    ( 2,1,5 )·( x,y,z )=0 ( 3,2,1 )·( x,y,z )=0 2x+y+5z=0 3x2y+z=0 7x+11z=0 z= 7 11 x y=2x5z =2x+ 35 11 x = 13 11 x ( x, 13 11 x, 7 11 x ) = x 11 ( 11,13,7 )

    よって求める問題の2平面の交線に平行なベクトルは。t(11, 13, -7)。

コード(Emacs)

Python 3

#!/usr/bin/env python3

from sympy import pprint, symbols, Matrix, solve

a = [((2, -1, 1), (3, 1, 1)),
     ((2, 1, 5), (3, -2, 1))]

x1, x2, x3 = symbols('x1 x2 x3')
X = Matrix([x1, x2, x3])

for i, (A, B) in enumerate(a, 12):
    print(f'{i}.')
    A = Matrix(A)
    B = Matrix(B)
    pprint(solve((A.dot(X), B.dot(X)), dict=True))
    print()

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

$ ./sample12.py
12.
⎡⎧    -2⋅x₃       x₃⎫⎤
⎢⎨x₁: ──────, x₂: ──⎬⎥
⎣⎩      5         5 ⎭⎦

13.
⎡⎧    -11⋅x₃       -13⋅x₃ ⎫⎤
⎢⎨x₁: ───────, x₂: ───────⎬⎥
⎣⎩       7            7   ⎭⎦

$

0 コメント:

コメントを投稿