Loading [MathJax]/jax/output/CommonHTML/jax.js

2019年5月14日火曜日

学習環境

ラング線形代数学(上)(S.ラング (著)、芹沢 正三 (翻訳)、ちくま学芸文庫)の4章(線形写像と行列)、3(線形写像の合成)、練習問題3の解答を求めてみる。



    1. id(1,1,0)=a11(2,1,1)+a21(0,0,1)+a31(-1,1,1)id(-1,1,1)=a12(2,1,1)+a22(0,0,1)+a32(-1,1,1)id(0,1,2)=a13(2,1,1)+a23(0,0,1)+a33(-1,1,1)2a11-a31=1a11+a31=1a11+a21+a31=0a11=23a31=1-23=13a21=-23-13=-12a12-a32=-1a12+a32=1a12+a22+a32=1a12=0a32=1a22=02a13-a33=0a13+a33=1a13+a23+a33=2a13=13a33=23a23=1

      よって、 求める問題の2つの基底と恒等を後に対応する行列は、

      Mββ,(id)=[23013-10113123]

    2. a11-a21+2a31=3a11+2a21-a31=24a21+a31=1a11+6a21=3a11=3-6a21a31=3-6a21+2a21-2=1-4a213-6a21-a21+2-8a21=3a21=215a11=3-6·215=3315a31=1-815=715a12-a22+2a32=0a12+2a22-a32=-24a22+a32=5a12+6a22=3a12=3-6a223-6a22-a22+2a32=0a32=7a22-324a22+7a22-32=5152a22=132a22=1315a12=3-6·1315=-3315a32=7·1315-32=2315a13-a23+2a33=1a13+2a23-a33=14a23+a33=2-3a23+3a33=0a23=a33a33=25a13-25+45=1a13=35

      よって求める行列は、

      115[33-33921367236]

コード

Python 3

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


print('3.')

t = [(Matrix([[1, -1, 0],
              [1, 1, 1],
              [0, 1, 2]]),
      Matrix([[2, 0, -1],
              [1, 0, 1],
              [1, 1, 1]])),
     (Matrix([[3, 0, 1],
              [2, -2, 1],
              [1, 5, 2]]),
      Matrix([[1, -1, 2],
              [1, 2, -1],
              [0, 4, 1]]))]

X = Matrix([[symbols(f'a{i + 1}{j + 1}') for j in range(3)]
            for i in range(3)])

for i, (A, B) in enumerate(t):
    print(f'({chr(ord("a") + i)})')
    pprint(solve(B * X - A))
    print()

入出力結果(cmd(コマンドプロンプト)、Terminal、Jupyter(IPython))

C:\Users\...>py sample3.py
3.
(a)
{a₁₁: 2/3, a₁₂: 0, a₁₃: 1/3, a₂₁: -1, a₂₂: 0, a₂₃: 1, a₃₁: 1/3, a₃₂: 1, a₃₃: 2
/3}

(b)
⎧                                                 13                          
⎨a₁₁: 11/5, a₁₂: -11/5, a₁₃: 3/5, a₂₁: 2/15, a₂₂: ──, a₂₃: 2/5, a₃₁: 7/15, a₃₂
⎩                                                 15                          

  23          ⎫
: ──, a₃₃: 2/5⎬
  15          ⎭


C:\Users\...>

0 コメント:

コメントを投稿