对下列数据进行回归：
clear
clc
%读人口数据
Y = [33815 33981 34004 34165 34212 34327 34344 34458 34498 34476 34483 34488 34513 34497 34511 34520 34507 34509 34521 34513 34515 34517 34519 34519 34521 34521 34523 34525 34525 34527]
%读时间变量
T = [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30]
%线性处理
for t = 1:30,
x(t) = exp(-t);
y(t) = 1/Y(t);
end
%计算回归系数B
c = zeros(30,1)+1;%产生一个30X1的矩阵全部为0的矩阵再加1
X = [c,x'];
B = inv(X'*X)*X'*y'%inv是求逆矩阵,最小二乘解的矩阵公式对于问题XB = Y
for i = 1:30
z(i) = B(1,1)+B(2,1)*x(i);%回归拟合值
s(i) = y(i)-sum(y)/30;%离差
w(i) = z(i)-y(i);%误差
end
S = s*s';%离差平方和
Q = w*w';%误差平方和
U = S-Q;%回归平方和
F = 28*U/Q
for j = 1:30,
Y(j) = 1/(B(1,1)+B(2,1)*exp(- j));
end
plot(T,Y)
注：我们所需要输出B和F值以及图像，所以后面不加分号。
结果：