Pythonの線形代数でコレスキー分解を返す

コレスキー分解を返すには、numpy.linalg.cholesky（）メソッドを使用します。正方行列aのコレスキー分解L*L.Hを返します。ここで、Lは下三角行列で、.Hは共役転置演算子です。 aはエルミートで正定値でなければなりません。 aがエルミートであるかどうかを確認するためのチェックは実行されません。さらに、aの下三角要素と対角要素のみが使用されます。実際にはLのみが返されます。

次に、パラメーターaは、エルミート（すべての要素が実数の場合は対称）の正定値入力行列です。このメソッドは、aの上三角または下三角のコレスキー因子を返します。 aが行列オブジェクトの場合、行列オブジェクトを返します。

ステップ

まず、必要なライブラリをインポートします-

import numpy as np

numpy.array（）メソッドを使用して2Dnumpy配列を作成する-

arr = np.array([[1,-2j],[2j,5]])

配列を表示する-

print("Our Array...\n",arr)

寸法を確認してください-

print("\nDimensions of our Array...\n",arr.ndim)

データ型を取得-

print("\nDatatype of our Array object...\n",arr.dtype)

形をとる-

print("\nShape of our Array object...\n",arr.shape)

コレスキー分解を返すには、numpy.linalg.cholesky（）メソッドを使用します-

print("\nCholesky decomposition in Linear Algebra...\n",np.linalg.cholesky(arr))

例

import numpy as np

# Creating a 2D numpy array using the numpy.array() method
arr = np.array([[1,-2j],[2j,5]])

# Display the array
print("Our Array...\n",arr)

# Check the Dimensions
print("\nDimensions of our Array...\n",arr.ndim)

# Get the Datatype
print("\nDatatype of our Array object...\n",arr.dtype)

# Get the Shape
print("\nShape of our Array object...\n",arr.shape)

# To return the Cholesky decomposition, use the numpy.linalg.cholesky() method.
print("\nCholesky decomposition in Linear Algebra...\n",np.linalg.cholesky(arr))

出力

Our Array...
[[ 1.+0.j -0.-2.j]
[ 0.+2.j 5.+0.j]]

Dimensions of our Array...
2

Datatype of our Array object...
complex128

Shape of our Array object...
(2, 2)

Cholesky decomposition in Linear Algebra...
[[1.+0.j 0.+0.j]
[0.+2.j 1.+0.j]]