用python进行qr分解

QR分解是将一个矩阵分解为一个正交矩阵和一个上三角矩阵的乘积，可以用于求解线性方程组、最小二乘问题等。在Python中，可以使用NumPy库中的`qr`函数进行QR分解。下面是一个示例代码：

import numpy as np

# 生成一个随机矩阵
A = np.random.rand(4, 4)

# 对矩阵进行QR分解
Q, R = np.linalg.qr(A)

# 打印分解结果
print("Q = ", Q)
print("R = ", R)

# 验证分解结果
print("Q*R = ", np.dot(Q, R))
print("A = ", A)

运行代码后，会输出分解结果和验证结果。其中，`Q`是一个正交矩阵，`R`是一个上三角矩阵，`Q*R`应该等于原矩阵`A`。