使用sympy实现奇异值分解（SVD）

import sympy
A=sympy.Matrix([[1/sympy.sqrt(3),1/sympy.sqrt(3)],[1/sympy.sqrt(3),0]])
U,D,V=A.singular_value_decomposition()
U.simplify()
D.simplify()
V.simplify()
print(U)
print(D)
print(V)
B=U@D@V.H
B.simplify()
print(B)
print(A==B)
print(A.equals(B))