Loading [MathJax]/jax/output/HTML-CSS/jax.js

2017年9月10日日曜日

学習環境

線型代数入門(松坂 和夫(著)、岩波書店)の第2章(ベクトル空間)、8(基底と次元(II))、問7.を取り組んでみる。


  1. a=(a1,a2,a3)b=(b1,b2,b3)c=(c1,c2,c3)0a+1b+1c=e11a+0b+1c=e21a+1b+0c=e3(b1+c1,b2+c2,b3+c3)=(1,0,0)(a1+c1,a2+c2,a3+c3)=(0,1,0)(a1+b1,a2+b2,a3+b3)=(0,0,1)b1+c1=1b2+c2=0b3+c3=0a1+c1=0a2+c2=1a3+c3=0a1+b1=0a2+b2=0a3+b3=1c1=1b1c2=b2c3=b3c1=a1c2=1a2c3=a3b1=a1b2=a2b3=1a3c1=1+a1c2=a2c3=1+a3a1=1+a1a1=121a2=a2a2=12a3=1+a3a3=12b1=12b2=12b3=112=12c1=12c2=12c3=1+12=12a=(12,12,12)b=(12,12,12)c=(12,12,12)

コード(Emacs)

Python 3

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

from sympy import pprint, symbols, solve, Matrix

print('5.')
a = Matrix(symbols('a1 a2 a3'))
b = Matrix(symbols('b1 b2 b3'))
c = Matrix(symbols('c1 c2 c3'))
e1 = Matrix([1, 0, 0])
e2 = Matrix([0, 1, 0])
e3 = Matrix([0, 0, 1])

s = solve([b + c - e1, a + c - e2, a + b - e3], dict=True)
for s0 in s:
    pprint(s0)

入出力結果(Terminal, IPython)

$ ./sample7.py
5.
{a₁: -1/2, a₂: 1/2, a₃: 1/2, b₁: 1/2, b₂: -1/2, b₃: 1/2, c₁: 1/2, c₂: 1/2, c₃:
 -1/2}
$

0 コメント:

コメントを投稿