最小二乘法进行曲线拟合@TOC
先假设拟合曲线的阶次,再利用最小二乘方法估计待估计参数。
clc;
clear all;
z=[66.6,84.9,88.6,78.0,96.8,105.2,93.2,111.6,88.3,117.0,115.2];
plot(z,‘o–’);
xlabel(‘第k年’);
ylabel(‘钢产量/百万吨’);
一次曲线拟合
clc;
clear all;
z=[66.6,84.9,88.6,78.0,96.8,105.2,93.2,111.6,88.3,117.0,115.2];
plot(z,‘o–’);
xlabel(‘第k年’);
ylabel(‘钢产量/百万吨’);
k=1:11;
Hk=[k;ones(1,11)]’; %计算Hk
estim=inv(Hk’*Hk)*Hk’*z’;%计算最小二乘估计方法式(注z的共轭)
ze=estim(1)k+estim(2);%利用估算的结果重新计算钢产量
%评价使用一次曲线的拟合结果
for i=1:4
ze(11+i)=estim(1)(11+i)+estim(2);%估算未来四年的钢产量
end
hold on;
plot(ze);
参考:[金学波]Kalman滤波器理论与应用——基于MATLAB实现