matlab spapi,第20章MATLAB多项式回归与数据插值重点.ppt

* 华文行楷，字号88，颜色红，位于每一章的第一节的最前面。　不使用动画 * 【例20.3-1续】针对例20.3-1中的数据，调用interp1函数作一维插值。 >> ylin = interp1(x0,y0,x); >> yspl = interp1(x0,y0,x,'spline'); >> plot(x0,y0,'ko'); >> hold on; >> plot(x,f,'k', 'linewidth',2); >> plot(x,ylin,':', 'linewidth',2); >> plot(x,yspl,'r-.', 'linewidth',2); >> xlabel('X'); >> ylabel('$$f(x)=\frac{1}{1+25x^2}$$','Interpreter','latex'); >> legend('插值节点','原函数图像','分段线性插值','三次样条插值') * 【例20.3-2】在加工机翼的过程中，已有机翼断面轮廓线上的20组坐标点数据，如表20.3-2所列，其中(x, y1) 和(x, y2) 分别对应轮廓线的上下线。假设需要得到 x坐标每改变0.1时的y 坐标，试通过插值方法计算加工所需的全部数据，并绘制机翼断面轮廓线，求加工断面的面积。 0 3 5 7 9 11 12 13 14 15 0 1.8 2.2 2.7 3.0 3.1 2.9 2.5 2.0 1.6 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 * x0 = [0,3,5,7,9,11,12,13,14,15]; y01 = [0,1.8,2.2,2.7,3.0,3.1,2.9,2.5,2.0,1.6]; y02 = [0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6]; x = 0:0.1:15; ysp1 = interp1(x0,y01,x,'spline'); pp = interp1(x0,y02,'spline','pp'); ysp2 = ppval(pp,x); xx = [x,fliplr(x)]; ysp = [ysp1,fliplr(ysp2)]; plot([x0,x0],[y01,y02],'o'); hold on; plot(xx,ysp,'r','linewidth',2); xlabel('X'); ylabel('Y'); legend('插值节点','三次样条插值','location','northwest'); % 截面面积： S1 = trapz(x,ysp1)-trapz(x,ysp2) S2 = trapz(xx,ysp) * 三、spline函数 yy = spline(x, Y, xx) 自变量观测值向量 因变量观测值向量 用户指定的插值点横坐标 1. spline函数的用法 计算得到的近似函数值 * 四、csape和csapi函数 pp = csape( x, y, conds) 自变量观测值向量 因变量观测值向量 边界条件参数 1. csape函数的用法 分段多项式形式的插值结果 conds参数取值 说 明 'complete'或 'clamped' 给定端点处的一阶导数值，默认为拉格朗日边界条件 'not-a-knot' 非纽结边界条件，csapi函数使用的就是这种边界条件 'periodic' 周期边界条件 'second' 给定端点处的二阶导数值。默认为[0, 0]，同'variational'情形 'variational' 设定端点处的二阶导数值为0 * values = csapi(x, y, xx) 2. csapi函数的用法 自变量观测值向量 因变量观测值向量 用户指定的插值点横坐标 计算得到的近似函数值 * * fun = @(x)sin(pi*x/2).*(x>=-1&x<1) + x.*exp(1-x.^2).*(x>=1 | x