插值和拟合

Posted 2021-03-13 code-fun

一维插值

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)计算函数值

 

经典例题：黄河小浪底调水调沙问题