Hello world

这是我在博客园的第一篇日志，我喜欢Matlab，这里就分享一个用最小二乘法拟合数据的Matlab脚本吧。

clc; clear;
t = 1:10;
R = [1 2 3 5 6 4 7 8 9 10];
figure;
plot(t,R,'*');
hold on;
N = numel(t);

%最小二乘法拟合线性模型参数
den = Nsum(t.^2) - sum(t)^2;
a = (sum(R)sum(t.^2) - sum(R.t)sum(t))/den;
b = (Nsum(R.t) - sum(R)*sum(t))/den;

%绘制拟合直线
R1 = a + b*t;
plot(t, R1);

%计算t2 = 70对应的R2值
t2 = 70;
R2 = a + b*t2;

分享完毕，谢谢关注。