求矩阵的平方根

注意
sqrtm( ) 和 sqrt( )是两个不同的函数
sqrtm( ) 对矩阵整体开平方
sqrt( ) 对矩阵中每个位置上的值开平方

import numpy as np
from scipy.linalg import sqrtm

import numpy as np
from scipy.linalg import sqrtm
a = np.mat(np.array([[1,4],[9,16]]))
'''
#矩阵a
 (1,	 4)
 (9,	 16)
'''

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8

对矩阵a整体开平方

# 对矩阵a整体开平方
b = sqrtm(a)   
c = b.dot(b)
'''
#矩阵b的值
(0.4+0.9j,	 0.9+-0.2j)
(1.9+-0.5j,  3.8+0.1j)
#矩阵c的值
 (1,	 4)
 (9,	 16)
'''

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11

对矩阵a中每个位置上的值开平方

# 对矩阵a中每个位置上的值开平方 
d = np.sqrt(a) 
e = d.dot(d)
'''
#矩阵d的值
1	2
3	4
#矩阵e的值
7	10
15	22
'''

转载自白十月