https://www.sciencedirect.com/science/article/pii/S0924013615300595#bib0040，请通过MATLAB求解公式（1）

公式（1）为：


其中，N为总样本数，X为输入变量，Y为输出变量，b为截距，β为斜率，ε为误差。

可以通过MATLAB中的regress函数求解，示例代码如下：

% 生成随机数据
X = randn(100, 1);
Y = 2*X + randn(100, 1);

% 构造回归模型
X = [ones(size(X)), X];
[beta, bint, r, rint, stats] = regress(Y, X);

% 输出结果
disp('回归系数：');
disp(beta');
disp('回归统计信息：');
disp(stats);

其中，X和Y为输入输出变量，regress函数返回的beta为回归系数，stats为回归统计信息。