用python写一个最小二乘法的代码

好的，以下是用 Python 实现最小二乘法的代码：

import numpy as np

# 计算最小二乘法的函数
def least_squares(x, y):
    x_mean = np.mean(x)  # 计算 x 的平均值
    y_mean = np.mean(y)  # 计算 y 的平均值
    
    # 计算最小二乘解的系数
    numerator = 0  # 分子
    denominator = 0  # 分母
    
    for i in range(len(x)):
        numerator += (x[i] - x_mean) * (y[i] - y_mean)
        denominator += (x[i] - x_mean) ** 2
    
    beta = numerator / denominator
    alpha = y_mean - beta * x_mean
    
    return alpha, beta

# 测试最小二乘法的函数
if __name__ == "__main__":
    x = np.array([0, 1, 2, 3, 4, 5])  # 自变量 x 的取值
    y = np.array([1, 3, 5, 7, 9, 11])  # 因变量 y 的取值

    alpha, beta = least_squares(x, y)  # 计算最小二乘解的系数
    print("y = " + str(alpha) + " + " + str(beta) + "x")  # 输出最小二乘法的方程式

这个代码的实现是基于 numpy 库完成的，numpy 是 Python 的一个常用科学计算库，提供了大量的数学函数和工具以及支持高维数组和矩阵运算。