线性最小二乘拟合
1、polyfit
例:a= polyfit(x0,y0,m)
其中,输入参数x0,y0为要拟合的数据,m为拟合多项式的次数(一般不超过3次),输出参数a为拟合多项式的次数(从0次开始)
%例:
clc;clear
x=[19,25,31,38,44];
y=[19.0,32.3,49.0,73.3,97.8];
ab=polyfit(x,y,2)
x0=19:0.1:44;
y0=ab(3)+ab(2)*x0+ab(1)*x0.^2;
plot(x,y,'o',x0,y0,'r')
2、polyval
例:y=polyval(a,x)
其中,输入参数a为拟合多项式的次数,输入参数x为所要求y值的x值,输出参数y为所要求的值。
例:某乡镇企业1990-1996年的生产利润如下
年份 1990 1991 1992 1993 1994 1995 1996
利润(万元) 70 122 144 152 174 196 202
试预测1997和1998年的利润。
%代码如下
x0=[1990 1991 1992 1993 1994 1995 1996