matlab插值曲线拟合课件.ppt

matlab插值曲线拟合;插值就是已知一组离散的数据点集，在集合内部某两个点之间预测函数值的方法。;Matlab常用数据插值函数及功能;四种插值方法比较;函数;函数;注意： （1）只对已知数据点集内部的点进行的插值运算称为内插，可比较准确的估测插值点上的函数值。 （2）当插值点落在已知数据集的外部时的插值称为外插，要估计外插函数值很难。 MATLAB对已知数据集外部点上函数值的预测都返回NaN，但可通过为interp1函数添加‘extrap’参数指明也用于外插。 MATLAB的外插结果偏差较大。;例1 对 在[-1, 1]上, 用n=20的等距分点进行线性插值, 绘制 f(x)及插值函数的图形.;解 在命令窗口输入:;例2. 在普通V带设计中，带轮的包角α与包角系数ka之间的关系如表所示。求α=133.5°时的包角系数ka。;a1=[90,100,110,120,125,130,135,140,145,150,155,160,165,170,175,180]; a2=[0.69,0.74,0.78,0.82,0.84,0.86,0.88,0.89,0.91,0.92,0.93,0.95,0.96,0.98,0.99,1]; ka=interp1(a1,a2,133.5) ka=0.8740;例3. 已知实验数据如表。; x=[0 0.25 0.50 0.75 1.00]; y=[0.9162 0.8109 0.6931 0.5596 0.4055]; x0=0.6; y01=interp1(x,y,x0); y02=interp1(x,y,x0,nearest); y03=interp1(x,y,x0,pchip); y04=interp1(x,y,x0,spline); y01,y02,y03,y04;例 4 对 在[-5, 5]上, 用n=11个等距分点作分段线性插值和三次样条插值, 用m=21个插值点作图,比较结果.; 0 1.0000 1.0000 1.0000 0.5000 0.8000 0.7500 0.8205 1.0000 0.5000 0.5000 0.5000 1.5000 0.3077 0.3500 0.2973 2.0000 0.2000 0.2000 0.2000 2.5000 0.1379 0.1500 0.1401 3.0000 0.1000 0.1000 0.1000 3.5000 0.0755 0.0794 0.0745 4.0000 0.0588 0.0588 0.0588 4.5000 0.0471 0.0486 0.0484 5.0000 0.0385 0.0385 0.0385;解 在命令窗口输入:;例 6 在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼上缘轮廓线的部分数据如下:;例 6 在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼???缘轮廓线的部分数据如下:;例7 天文学家在1914年8月份的7次观测中, 测得地球与金星之间距离(单位: m), 并取其常用对数值与日期的一组历史数据如下所示, 试推断何时金星与地球的距离(单位: m)的对数值为 9.9352.;解 由于对数值 9.9352 位于 24 和 26 两天所对应的对数值之间, 所以对上述数据用三次样条插值加细为步长为1的数据:; x=linspace(0,2*pi,11); y=sin(x).*exp(-x/5); xi=linspace(0,2*pi,21); yi=interpft(y,21); plot(x,y,o,xi,yi); legend(Original,Curve by interpft);Lagrange插值;function yi=lagrange(x,y,xi) yi=zeros(size(xi)); np=length(y); for i=1:np z=ones(size(xi)); for j=1:np if i~=j z=z.*(xi-x(j))/(x(i)-x(j)); end end yi=yi