『Python数値計算ノート』ではアフィリエイトプログラムを利用して商品を紹介しています。

ユニタリ行列

≪【前の記事】エルミート行列

ユニタリ行列

実空間 Rn直交行列 Q は、複素空間 Cn において ユニタリ行列 U として拡張定義されます。ユニタリ行列は
 (1)UU=UU=I
を満たす正方行列です。U が実行列であれば U=UT なので、(1) は直交行列 Q の定義を包含します。たとえば、前回記事で扱ったエルミート行列
 (2)H=[21i1+i3]
を対角化する行列
 (3)[63336(1+i)63(1+i)3]
はノルム 1 の列ベクトル同士が直交するので ユニタリ行列 です。
 
列ベクトル同士が直交する任意の行列は、正規化して各々のベクトルの ノルム1 に揃えることで、ユニタリ行列にすることができます。たとえば、
 (4)A=[211+i1+i]
は列ベクトル同士が直交する行列です。これを正規化してユニタリ行列にしてみましょう。

# python_unitary

# In[1]

import numpy as np
from scipy import linalg
np.set_printoptions(precision=3)

# 共役転置関数
def hermitian(arr):
    return np.conjugate(arr.T)

# 正規化関数
def normalize(arr, axis=0):
    norm = linalg.norm(arr, axis=axis)
    return arr / norm

A = np.array([[-2, 1],
              [1+1j, 1+1j]])

# Aを正規化してユニタリ行列にする
U = normalize(A)

# ユニタリ行列が満たすべき性質を確認
X1 = hermitian(U) @ U
X2 = U @ hermitian(U)

print("U:\n{}\n".format(U))
print("(U^†)U:\n{}\n".format(X1))
print("U(U^†):\n{}".format(X2))

'''
U:
[[-0.816+0.j     0.577+0.j   ]
 [ 0.408+0.408j  0.577+0.577j]]

(U^†)U:
[[1.+0.j 0.+0.j]
 [0.+0.j 1.+0.j]]

U(U^†):
[[ 1.000e+00+0.000e+00j -2.946e-18+2.946e-18j]
 [-2.946e-18-2.946e-18j  1.000e+00+0.000e+00j]]
'''

ユニタリ行列の逆行列

UU=UU=I を満たすので、ユニタリ行列 U の随伴行列 UU の逆行列です。
 (5)U=U1
(4) を正規化して作ったユニタリ行列の随伴行列と逆行列を計算して、両者が一致することを確認してみましょう。

# In[2]

U_H = hermitian(U)
U_inv = linalg.inv(U)

print("Uの随伴行列:\n{}\n".format(U_H))
print("Uの逆行列:\n{}".format(U_inv))

'''
Uの随伴行列:
[[-0.816-0.j     0.408-0.408j]
 [ 0.577-0.j     0.577-0.577j]]

Uの逆行列:
[[-0.816+0.j     0.408-0.408j]
 [ 0.577+0.j     0.577-0.577j]]
'''

応用がみえる線形代数 (Iwanami Mathematics)

新品価格
¥2,640から
(2024/3/30 22:35時点)

ユニタリ変換

任意のベクトル x にユニタリ行列 U をかけたとき、Ux のノルムは
 (6)(Ux,Ux)=(Ux)Ux=xUUx=xx=(x,x)
となって、ユニタリ行列による線形変換 (ユニタリ変換) がベクトルのノルムを保持する等長変換であることがわかります (実数空間の直交行列 Q と同じです)。すなわち、
 (7)Ux∥=∥x
が成り立ちます。一例として、複素ベクトル
 (8)z=[3+5i2i]
について、自身のノルムと Uz のノルムが一致することを確認してみます。

# In[3]

# 複素ベクトルを定義
z = np.array([3 + 5j, 2-1j])

# zのノルムを計算
z_norm = linalg.norm(z)

# zにユニタリ行列Uを作用させる
Uz = U @ z

# Uzのノルムを計算
Uz_norm = linalg.norm(Uz)

print(z_norm)
print(Uz_norm)

# 6.244997998398398
# 6.244997998398399

ユニタリ行列の固有値

ユニタリ行列 U の固有値が λ、すなわち Ux=λx であるとき、λx∥=x より、λ∥=1 となります。

# In[4]

# ユニタリ行列Uの固有値
eigvals = linalg.eig(U)[0]

# 固有値の絶対値
norm = np.abs(eigvals)

print("Uの固有値:\n{}\n".format(eigvals))
print("固有値の絶対値:\n{}".format(norm))

# Uの固有値:
# [-0.997-0.075j  0.758+0.652j]

# 固有値の絶対値:
# [1. 1.]

≫【次の記事】ケーリー・ハミルトンの定理

コメント

  1. HNaito より:

    最初は1/4くらいで挫折しましたが、「高校数学でわかる線形代数」(竹内淳著 ブルーバックス) を読んでから再挑戦し、ようやく最後までたどり着きました。ストラング:線形代数イントロダクションの中古本も見つけたので注文しました。Ax=b, 0, λx の奥深い世界をのぞけたらいいと思います。これからはPythonを復習しながら「機械学習」のほうで勉強させていただきますので、よろしくお願いいたします。

    • あとりえこばと より:

      最後までお読みいただき、ありがとうございます。散見される誤植のせいで、読みにくいところもあったかと思いますが、その都度ご指摘していただいて本当に感謝しています。

      『ストラング:線形代数イントロダクション』は名著で、私も記事を書くために参照しています。ただ、扱う内容がとても幅広く、演習問題の数も膨大なので、読み切るのはなかなか大変です(私も全部は目を通していません)。焦らずゆっくりと読み進めると良いと思います。このサイトで扱っていない内容としては、実行列の複素数固有値、正定値行列、特異値分解 (SVD) などがあります。いずれもとても興味深い内容なので、ぜひ挑戦してください。

      「機械学習」の勉強もぜひ楽しんでください。なるべくライブラリに頼らず、数学をたくさん用いて機械学習を根本から理解できるような構成になっています。これまで学んだ線形代数の知識も生きてくると思います。よろしくお願いします。

      最近は他の仕事などで忙しく、このサイトの記事をなかなか更新できていませんが、久しぶりに線形代数の新しい記事の執筆にとりかかりました。これまで取りこぼしていた内容をいくつか付け加えたり、さらに進んだ内容も扱ってみたいと思いますので、時々お立ち寄りください。次回は「ケーリー・ハミルトンの定理」を予定しています。