数学软件4(实验三理论)程序.pptVIP

实验三 函数式－直接确定型模型 例1：经验函数(机床加工问题） 用程控铣床加工机翼断面的下轮廓线时 每一刀只能沿x方向和y方向走非常小的一步。 表3-1给出了下轮廓线上的部分数据 但工艺要求铣床沿x方向每次只能移动0.1单位. 这时需求出当x坐标每改变0.1单位时的y坐标。 试完成加工所需的数据，画出曲线. 观察插值、拟合的效果 运行观察程序zxy3_1.m 选取一个已知函数作为参考，并将这一函数的图象用虚线显示在图中。 观察程序允许用鼠标选取节点 按鼠标左键选点，按右键选最后一个点 观察不同的选点方式对各种插值和拟合效果的影响 程序注解（inline指令） 定义内联函数：inline指令 g=inline(x^2-x^4); 程序 程序注解（ginput） [x,y,button] = ginput(n) 用鼠标在屏幕选n个点，返回这n个点，存于x,y中。 button 记录了选点时使用的鼠标键方式： 1为左键、2为中间键、 3为右键。 程序注解（插值拟合） xx=linspace(a,b,n); %定义自变量xx ynearest=interp1(x1,y1,xx,nearest); ylinear=interp1(x1,y1,xx,linear); yspline=interp1(x1,y1,xx,spline); [p,c]=polyfit(x1,y1,4); ypolyfit=polyval(p,xx); 程序注解（插值拟合） subplot(2,2,1), h=plot(xx,ynearest,r-);set(h,linewidth,2) subplot(2,2,2), h=plot(xx,ylinear,r-);set(h,linewidth,2); subplot(2,2,3), h=plot(xx,yspline,r-);set(h,linewidth,2) subplot(2,2,4), h=plot(xx,ypolyfit,r-);set(h,linewidth,2) 插值拟合效果观察 沿曲线选取３个节点，保持等间隔。当节点较少时，插值的效果如何？ 加密节点，共８个等距节点，观察插值的效果，如果去掉中间的一个节点，插值效果又会如何？ 有意偏离原来的曲线，如果误差较大，将会怎样呢？ * * 返回 表3-2 某河流水深和流速的观测数据 2.9757 3.0594 3.1266 3.1870 3.2282 3.2516 3.2611 3.2532 3.2299 3.1950 v 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 x 例2：试根据下述测量数据，确定该河流水深和流速之间的 关系，并画出曲线。 返回 一、插值 已知有n +1个节点（xj,yj），j = 0,1,…, n，其中xj互不 相同；节点(xj, yj)可看成由某个函数 y= f（x）产生. f 的解析 表达式可能十分复杂或不存在封闭形式, 也可以是未知的 插值方法： 构造一个相对简单的函数 y=g(x)，使g通过全部节点，即 使g (xj) = yj，j=0,1,…, n ，用g (x)作为函数f ( x )的近似。 插值指令： ★ 一维插值指令interp1 其调用格式 其中两个同维向量 称为节点向量，满足 的函数关系。输出值 是与 对应的函数值。 参数’method’是该插值指令的多种算法，对应四个选项： ‘nearest’ 为近邻插值；‘linear’为线性插值； ‘spline’ 为样条插值； cubic为立方函数插值。 程序： clear,clf, x=[ 0 3 5 7 9 11 12 13 14 15]; y=[0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6]; xi=0:0.1:15; yi=interp1(x,y, xi,spline); plot(x,y,o,xi,yi); legend(插值节点,样条插值) 例1 例3 在区间[0,10]画出 y=sin(x) 的曲线，取插值节点 和节点处的函数值 作分段线性插值， 并画出相应的折线图，将两图绘在一张图上。 程序： clear,clf, x=0:10; y=sin(x); xi=0:0.5:10; yi=interp1(x,y,xi,linear); t=0:0.001:10; z=sin(t); plot(x,y,ro,xi,yi,t,z); legend(插值节点,线性插值,sinx,4) 二、多项式拟合 拟合方法： 对给定的