插值法
1.介绍
介绍比较常见的插值法优缺点,但我们在建模中主要用到两类——埃尔米和样条,其他不用可略。
1.1多项式插值
n阶多项式有n+1个未知参数的值,代入n+1个互异点,得到n+1个方程解出未知数,使得:
因为n+1个方程确定n+1个未知数,解存在;同时每个方程不同,所以解唯一。
注:但是一般选取二到三阶,阶数过高会出现过拟合
1.2拉格朗日插值
一种不用计算参数的多项式函数,但是他的公式已经注定他是n阶函数
记公式即可(推导略),因为分别带入x0,x1等点可得对应的y0、y1,所以这样设计。
1.3分段插值
为了防止拉格朗日n阶函数,阶数过高,我们只选取这个点附近的三个点来
1.4牛顿插值
与拉格朗日插值法有些相似,但是具有继承性,新添加的点只需要在原来的函数上添加新项。
1.5三次埃尔米插值(常用)
不仅确定了点,还确定这个点的导数,更加精确
1.6三次样条插值(常用)
不仅确定了点,还确定这个点的二阶可微,更加精确
2.代码
matlab中自带函数
p1 = pchip(x,y,new_x); %分段三次埃尔米特插值
p2 = spline(x,y,new_x); %三次样条插值
拟合法
前言:
高中学过的一元线性回归
与插值法的区别:不用经过所有的点
1.步骤:
首先确定拟合曲线,然后求参
利用最小二乘法求解参数
评估拟合的好坏,可以直接用最小二乘法,算出误差SSE,或者用R方
2.操作:
利用matlab工具箱