2017年11月13日月曜日

学習環境

線型代数入門(松坂 和夫(著)、岩波書店)の第3章(線型写像)、12(連立1次方程式(I))、問題3.を取り組んでみる。


    1. ( - 4 0 - 8 - 4 1 - 5 8 0 16 )
      ( - 4 0 - 8 0 1 3 0 0 0 )
      ( 1 0 2 0 1 3 0 0 0 )

      基本解

      u 1 = - 2 , - 3 , 1

      一般解

      x 1 = - 2 α x 2 = - 3 α x 3 = α

    2. ( 1 2 3 3 0 5 0 - 4 3 )
      ( 3 6 9 3 0 5 0 - 4 3 )
      ( 0 6 4 3 0 5 0 - 4 3 )
      ( 0 12 8 3 0 5 0 - 4 3 )
      ( 0 0 17 3 0 5 0 - 4 3 )
      ( 0 0 17 3 0 0 0 1 0 )

      方程式の解は0のみから成る。いいかえれば、自明な解しかもたない。

      x 1 = 0 x 2 = 0 x 3 = 0

コード(Emacs)

Python 3

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

print('1.')
xs = symbols('x1, x2, x3')

MS = [Matrix([[4, -2, 2],
              [-4, 1, -5],
              [4, 1, 11]]),
      Matrix([1, 2, 3, 3, 0, 5, 0, -4, 3]).reshape(3, 3)]

x = Matrix(xs).reshape(3, 1)
for i, M in enumerate(MS):
    print(f'({chr(ord("a") + i)})')
    for t in [M, solve(M * x, xs)]:
        pprint(t)
        print()
    print()

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

$ ./sample1.py
1.
(a)
⎡4   -2  2 ⎤
⎢          ⎥
⎢-4  1   -5⎥
⎢          ⎥
⎣4   1   11⎦

{x₁: -2⋅x₃, x₂: -3⋅x₃}


(b)
⎡1  2   3⎤
⎢        ⎥
⎢3  0   5⎥
⎢        ⎥
⎣0  -4  3⎦

{x₁: 0, x₂: 0, x₃: 0}


$

0 コメント:

コメントを投稿