插值和拟合
一维插值
y=interp1(x0,y0,x,'method')【x0单调的】
method:默认为线性插值
'nearest'最近项插值
‘linear'线性插值
’spline'立方样条插值
'cubic'立方插值
三次样条插值(光滑,它的曲率也是连续的)
pp=csape(x0,y0,conds)conds是边界条件默认为拉格朗日边界条件
y=fnval(pp,x)获得函数值
二维插值(如为了画出精确的等高线图,就要插入更多的点)
1)插值节点为网格节点
y=interp2(x0,y0,z0,x,y,'method') pp=csape({x0,y0},z,conds)
2)插值接点为散乱节点
z=griddata(x,y,z,x1,y1)
拟合
a=polyfit(x0,y0,m) m为最高次数
y=polyval(a,x)计算函数值
经典例题:黄河小浪底调水调沙问题