初等模型一.pptVIP

二次插值也称作抛物线插值，几何意义是用抛物线来代替曲线y=f(x), % 按照表1输入原始数据 x=[0:0.2:5, 4.8:-0.2:0]; y=[5 4.71 4.31 3.68 3.05 2.5 2.05 1.69 1.4 1.18 1 0.86 0.74 0.64 0.57 0.5 ... 0.44 0.4 0.36 0.32 0.29 0.26 0.24 0.2 0.15 0 -1.4 -1.96 -2.37 -2.71 ... -3 -3.25 -3.47 -3.67 -3.84 -4 -4.14 -4.27 -4.39 -4.49 -4.58 -4.66 ... -4.74 -4.8 -4.85 -4.9 -4.94 -4.96 -4.98 -4.99 -5]; % 逆时针方向转900，节点（x, y）变为（u, v） v0=x; u0=-y; % 按0.05的间隔在u方向产生插值点 u=-5:0.05:5; % 在v方向计算分段线性插值 v1=interp1(u0,v0,u); % 在v方向计算三次样条插值 v2=spline(u0,v0,u); % 在（x, y）坐标系输出结果 [v1 v2 -u] subplot(1,3,1),plot(x,y),axis([0 5 -5 5]) gtext(原轮廓线,FontSize,12) subplot(1,3,2),plot(v1,-u),axis([0 5 -5 5]) gtext(分段线性插值,FontSize,12) subplot(1,3,3),plot(v2,-u),axis([0 5 -5 5]) gtext(三次样条插值,FontSize,12) 1. 拉格朗日插值:自编程序,如名为 lagr.m 的M文件， 第一行为 function y=lagr(x0,y0,x) 输入:节点x0,y0, 插值点x (均为数组，长度自定义)）； 输出:插值y (与x同长度数组)）。 应用时输入x0,y0,x后,运行 y=lagr1(x0,y0,x) 2. 分段线性插值:已有程序 y=interp1(x0,y0,x) 3. 三次样条插值:已有程序 y=interp1(x0,y0,x,’spline’) 或 y=spline(x0,y0,x) 用MATLAB作插值计算 插值的应用 图1 零件的轮廓线 （x间隔0.2） 表1 x间隔0.2的加工坐标x,y（图1右半部的数据） 例：数控机床加工零件 问题：待加工零件的外形根据工艺要求由一组数据(x,y)给出，用数控机床加工时刀具必须沿这些数据点前进，并且由于刀具只能沿x方向或y方向走非常小的一步，所以需要将已知数据加密，得到加工所要求的步长很小的(x,y)坐标。 ……… 2.6,0.64 2.4,0.74 2.2,0.86 2.0,1.00 1.8,1.18 1.6,1.40 1.4,1.69 1.2,2.05 1.0,2.50 0.8,3.05 0.6,3.68 0.4,4.31 0.2,4.71 0.0,5.00 图1是待加工零件的轮廓线，表1给出了轮廓线上x每间隔0.2的加工坐标(x,y)(顺时针方向为序，由轮廓线的左右对称性，表中只给出右半部数据）。 加工时需要x每改变0.05时的y值，试完成加工所需的加密数据，画出曲线。 建立模型： 若直接利用表1数据加密，为得到经过这些点的曲线，须将零件轮廓线的右半部分分为y≥0和y≤0两部分，分别计算两个（单值）函数在加密点（即差值点）的值，还要设法保证连接点处的光滑性。 图2 逆时针方向转90度的结果 一种简捷的方法：将图1逆时针方向转90度，轮廓线上下对称，只需对上半部计算一个函数在插值点的值。 三、曲 线 拟 合 问 题 的 提 法 已知一组（二维）数据，即平面上 n个点（xi,yi) i=1,…n, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所有数据点最为接近，即曲线拟合得最好。 + + + + + + + + + x y y=f(x) (xi,yi) ?i ?i 为点（xi,yi) 与曲线 y=f(x) 的距离 拟合与插值的关系 函数插值与曲线拟合都是要根据一组数据构造一个函数作为近似，由于近似的要求不同，二者的数学方法上是完全不同的。 问题：给定一批数据点，需确定满足特定要求的曲线或曲面 解决方案： 若不要求曲线（面）通过所有数据点，而是要求它反映对象整体的变化趋势，这就是数据