# 検算
import numpy as np
a = np.array([3,4])
b = np.array([1,2])
np.outer(a,b)array([[3, 6],
[4, 8]])Aはのときは零にならないが、で零になる
Source
# 検算
from sympy import symbols, Matrix
a, b, c = symbols('a b c')
A = Matrix([
[0, a, b],
[0, 0, c],
[0, 0, 0],
])# 検算
A @ ALoading...
# 検算
A @ A @ ALoading...
2 つの行列 が可換となるように a の値を求めよ。
Source
# 検算
from sympy import Symbol, Matrix
a = Symbol('a')
A = Matrix([
[1, a],
[0, a**2],
])
B = Matrix([
[a**2, a],
[0, 1],
])
A @ BLoading...
Source
# 検算
import numpy as np
for a in (0, -1, 1):
A = np.array([
[1, a],
[0, a**2]
])
B = np.array([
[a**2, a],
[0, 1],
])
assert (A @ B == B @ A).all()2.5¶
n 次の正方行列 A, B に対して、 とおき、これを A と B の交換子積という。 交換子積に関して、次の 1~3 が成り立つことを示せ。ただし、A, B, C はすべて n 次の正方行列で、O は n 次の零行列である。
(交代性または反対称性)
Source
# 例えば n=2 のとき
from sympy import symbols, Matrix
a11, a12, a21, a22 = symbols("a_11, a_12, a_21, a_22")
b11, b12, b21, b22 = symbols("b_11, b_12, b_21, b_22")
A = Matrix([
[a11, a12],
[a21, a22]
])
B = Matrix([
[b11, b12],
[b21, b22]
])
(A @ B - B @ A) == (- 1 * (B @ A - A @ B))True2.7¶
次の問いに答えよ。
任意の正方行列は対称行列と交代行列の和で一意的に表されることを、次の(ア)、(イ)の手順で示せ。
(ア)正方行列 A が対象行列 X と交代行列 Y の和で、 A = X + Y \tag{2.1} と表されると仮定する。 このとき、式 2.1 の両辺の転置を取ることにより、 A^T = X − Y \tag{2.2} が成り立つ。
式 2.1 と式 2.2 を連立させることにより、X, Y を A および を用いて表せ。
Invalid delimiter type 'ordgroup' at position 46: …= X + Y\\
+\big{̲)̲}̲& A^T = X - Y\\…
\begin{array}{rr}
& A = X + Y\\
+\big{)}& A^T = X - Y\\
\hline
& A + A^T = 2 X\\
& \to X = \frac{1}{2} (A + A^T)
\end{array}Invalid delimiter type 'ordgroup' at position 46: …= X + Y\\
-\big{̲)̲}̲& A^T = X - Y\\…
\begin{array}{rr}
& A = X + Y\\
-\big{)}& A^T = X - Y\\
\hline
& A - A^T = 2 Y\\
& \to Y = \frac{1}{2} (A - A^T)
\end{array}