import numpy as np

# 定义矩阵 A 和向量 b
A = np.array([[1, 1], [1, 2], [1, 3]])  # 设计矩阵（3x2）
b = np.array([1, 2, 2])                # 观测值（3x1）

# 求解最小二乘问题
x, residuals, rank, singular_values = np.linalg.lstsq(A, b, rcond=None)

print("最优解 x:", x)
print("残差平方和:", residuals)