习题2.7（2）

import numpy as np  # 定义系数矩阵A和常数项向量b  
A = np.array([[2, 3, 1],  [1, -2, 4],  [3, 8, -2],  [4, -1, 9]])  
b = np.array([4, -5, 13, -6])  # 使用numpy的lstsq函数求解最小二乘解  
# 对于这个特定的问题，由于方程数和未知数数量相同，且没有矛盾，lstsq将给出唯一解  
x, residuals, rank, s = np.linalg.lstsq(A, b, rcond=None)  # 输出解  
print("解 x:", x)  # 验证解是否正确（可选）  
print("验证结果:", np.dot(A, x))  # 检查是否精确等于b（对于精确解，这应该非常接近）  
print("与b的误差:", np.linalg.norm(np.dot(A, x) - b))  # 计算系数矩阵的秩（可选，以确认方程组是否有唯一解）  
print("系数矩阵的秩:", np.linalg.matrix_rank(A))  # 由于秩等于未知数数量，且没有矛盾，我们可以确信有一个唯一解print("学号：3008")

结果如下