2017年9月15日金曜日

学習環境

解析入門〈2〉(松坂 和夫(著)、岩波書店)の第10章(n次元空間)、10.2(ベクトル空間)、問題1.を取り組んでみる。



    • c 1 a+ c 2 b+ c 3 c=0 c 1 ( 1,1,0 )+ c 2 ( 1,0,1 )+ c 3 ( 0,2,1 )=( 0,0,0 ) c 1 + c 2 =0 c 1 +2 c 3 =0 c 2 + c 3 =0 c 3 = c 2 c 2 = c 1 c 3 =( c 1 )= c 1 c 1 +2 c 1 =0 c 1 =0 c 2 =0 c 3 =0

      よって、a、b、cは1次独立なのでR^3の基底をなす。

    • c 1 a+ c 2 b+ c 3 c=( 2,4,1 ) c 1 ( 1,1,0 )+ c 2 ( 1,0,1 )+ c 3 ( 0,2,1 )=( 2,4,1 ) c 1 + c 2 =2 c 1 +2 c 3 =4 c 2 + c 3 =1 c 1 =2 c 2 2+ c 2 +2 c 3 =4 c 2 +2 c 3 =6 c 3 =7 c 2 =8 c 1 =10

      よって{a、b、c}を基底とするベクトルuの座標は(10, -8, 7)。

コード(Emacs)

Python 3

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

from sympy import pprint, symbols, Matrix, solve

print('1.')
a = Matrix([1, -1, 0])
b = Matrix([1, 0, 1])
c = Matrix([0, 2, 1])
z = Matrix([0, 0, 0])
u = Matrix([2, 4, -1])
x1, x2, x3 = symbols('x1 x2 x3')
eq = x1 * a + x2 * b + x3 * c
pprint(eq.T)

for v in [z, u]:
    pprint(v.T)
    pprint(solve(eq - v))

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

$ ./sample1.py
1.
[x₁ + x₂  -x₁ + 2⋅x₃  x₂ + x₃]
[0  0  0]
{x₁: 0, x₂: 0, x₃: 0}
[2  4  -1]
{x₁: 10, x₂: -8, x₃: 7}
$

0 コメント:

コメントを投稿