插值和拟合

一维插值

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

 

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

posted @ 2020-01-26 20:51  小帆敲代码  阅读(279)  评论(0编辑  收藏  举报