一般用曲线拟合33数据的插值计算一函数插值.PPT

* 在生产和科学实验中，自变量x与因变量y之间的函数关系式有时不能直接写出表达式，而只能得到函数在若干个点的函数值或导数值. 当要求知道观测点之外的函数值时，需要估计函数在该点的数值. 这就要根据观测点的值，构造一个比较简单的函数y=φ(x)，使函数在观测点的值等于已知的数值或导数值，寻找这样的函数φ(x)，办法是很多的. 根据测量数据的类型有如下两种处理观测数据的方法： ① 测量值是准确的，没有误差，一般用插值. ② 测量值与真实值有误差，一般用曲线拟合. 3.3 数据的插值计算 一.函数插值 1.一维插值 格式一： YI = interp1(X,Y,XI,method) X,Y为原始数据，XI, YI为插值出的数值， method是插值所用的方法. . Available methods are: nearest - nearest neighbor interpolation linear - linear interpolation spline - piecewise cubic spline interpolation (SPLINE) pchip - piecewise cubic Hermite interpolation (PCHIP) cubic - same as pchip 例1.对 y=cosx的数据进行插值，比较各种插值方法 x=[-2*pi:0.5*pi:2*pi];y=cos(x);xi=-2*pi:0.3*pi:2*pi; y_nearest=interp1(x,y, xi,nearest); y_linear= interp1(x,y,xi); y_spline= interp1(x,y,xi, spline ); y_cubic= interp1(x,y,xi, cubic ); plot(x,y,o,xi,y_nearest,-,xi,y_linear, r* , xi,y_spline,k:,xi,y_cubic,k -); legend (original data,nearest,linear,spline,cubic) (图1) 例2. 在某处测得海洋不同深处水温 如下： 2.13 2.54 3.40 4.28 7.04 水温(℃) 1634 1422 950 714 446 深度(m) 试求在深度500米、1000米、1500米处的水温 解法一：x=[-2*pi:0.5*pi:2*pi];y=cos(x);xi=-2*pi:0.3*pi:2*pi; y_nearest=interp1(x,y, xi,nearest); y_linear= interp1(x,y,xi); y_spline= interp1(x,y,xi, spline ); y_cubic= interp1(x,y,xi, cubic ); 例2. 在某处测得海洋不同深处水温 如下： 2.13 2.54 3.40 4.28 7.04 水温(℃) 1634 1422 950 714 446 深度(m) 试求在深度500米、1000米、1500米处的水温 解法一：x=[446,714,950,1422,1634]; y=[7.04,4.28,3.40,2.54,2.13]; xi=[500,1000,1500]; yi1= interp1(x,y,xi);%线性插值算法 yi2= interp1(x,y,xi, spline );%三次样条插值. yi3= interp1(x,y,xi, cubic );%分段三次Hermite多项式插值 Y=[yi1;yi2;yi3] 例2. 在某处测得海洋不同深处水温 如下： 2.13 2.54 3.40 4.28 7.04 水温(℃) 1634 1422 950 714 446 深度(m) 试求在深度500米、1000米、1500米处的水温 结果为： Y =6.4839 3.3089 2.3892 6.2569 3.2853 2.4025 6.3212 3.2803 2.3914 用不同的方法，插值的结果是不一样的 解法二：h=[446,714,950,1422,1634]; w=[7.04,4.28,3.40,2.54,2.13]; hi=[500,1000,1500]; p=spline(h,w); ci=ppval(p,hi) plot(h,w, or) ,hold on,plot(hi,ci,*) legend (original data,spline) 我们将原始数据与插值结果在同一个坐标系