用matlab写一段小二乘法代码

以下是一个简单的 Matlab 程序，实现了最小二乘法的计算：

function [b] = least_squares(X,y)
% 最小二乘法求解函数
% X: 自变量矩阵
% y: 因变量向量

% 求解系数矩阵b
b = (X' * X) \ (X' * y);

end

使用方法：将自变量矩阵 X 和因变量向量 y 作为参数输入该函数中，该函数会返回一个系数矩阵 b，即最小二乘法的解。

例如，我们可以使用以下代码来测试该函数：

% 生成测试数据
X = [1 2; 2 3; 3 4; 4 5];
y = [2; 4; 6; 8];

% 计算最小二乘法解
b = least_squares(X, y);

% 输出结果
disp(b);

该程序的输出结果应该为：

1.0000
1.0000

这意味着最小二乘法的解为 y = 1x1 + 1x2，即通过线性回归得到的直线方程为 y = x + 1。