第五节几何设计.ppt

* * 第五节 几何设计 曲线按几何性质可分为两类: 一类为规则曲线,如圆、抛物线、渐开线、圆锥曲线等,这 类曲线都可用方程式表示 一类是不规则曲线,像三次参数样条曲线、贝塞尔曲线等 就属于这类曲线。 本节主要介绍最基本的方法, 拉格朗日插值法和最小二乘法. 在包装设计过程中,除经常遇到一些规则曲线,还要用到一些 不规则曲线。不规则的曲线的绘制,通常采用拟合方法来实 现的. 不规则曲线的拟合方法大致分成两类: 一类是插值法,要求拟合曲线通过全部数据点 一类是逼近法,由于数据点本身存在一定的误差, 拟合的曲线不要求通过全部数据点 一、多项式插值 一个函数,无法用确切的函数式表达时,可利用由实验或测量等手段得到的一组离散的数据点yi = f ( x i ) i= 0,1,2,…….n 构造一个通过这组数据的函数 P(x) , 近似地代替未知函数 f ( x ),从而能求出任一点 x*∈[ x0 , xn ]处函数f( x )的近似值 P (x* ) , 这就是函数的插值, P(X) 叫插值函数 函数插值有多种方法,其中多项式插值就是一种形式简单, 计算方便,出现最早的插值方法。 什么叫插值函数： 1、线性插值 最简单的插值是两点插值, 即用一个一次多项式 y = p 1( x ) 作为插值多项式, 使两个插值点满足此式. 从几何图形上看,此一次多项式y = p 1( x ) 表示通过两点 (x1 , y1 ), (x2 , y2 ) 的直线,因此两点插值也称为线性插值. 由解析几何可写出p 1( x )的两种表达式: 线性插值 点斜式 对称式 若记： 则式可以表为 p 1( x ) = A1 (x) y1 + A2 (x) y2 A0(x)和A1（x）称为以X1 , X2 为结点的基本插值多项式. 由上式可见,一次插值多项式p 1( x )可由两个基本插值多项式 A1 (x)和A2 (x) 通过线性组合的方法构成。 线性插值只利用两个结点上的信息,精度很低。为改善精度,常用多于两点的插值。 2、抛物线插值 用通过三点（x0 、y0）、 （x1 、y1）、（x2 、y2）所作的 抛物线来近似表示函数曲线y=f(x) 若给定含有三个结点的函数表 X X0 X1 X2 Y Y0 Y1 Y2 作二次多项式y = p 2( x ) 使它满足下面条件 p 2( x0 ) = y0 , p 2( x 1) = y1 , p 2( x 2) = y2 作二次多项式y = p 2( x ) 使它满足下面条件p 2( x0 ) = y0 , p 2( x 1) = y1 , p 2( x 2) = y2 p 2( x ) = A0 (x) y0 + A1 (x) y1 + A2(x) y2 三个基本插值多项式可用以下方法求得： 为满足结点函数值要求，须满足以下条件： 当x = x0时 y=y0, 即要求 A0(x0)=1 A1(x0)=0 A2(x0)=0 当x = x1时 y=y1, 即要求 A0(x1)=0 A1(x1)=1 A2(x1)=0 当x = x2时 y=y2, 即要求 A0(x2)=0 A1(x2)=0 A2(x2)=1 可以概括地写成 ： 由以上条件就可以求得基本插值多项式。 以A0(x)为例： （1）由条件A0（x1）= A0(x2)=0 ,可知以A0(x)中含有（x-x1）和 （x-x2）两个因子，令A0(x)＝λ（x-x1）（x-x2） （2）用条件A0(x0)=1 确定系数λ 得到 ： 同理得到 ： 于是插值多项式为： 3、拉格朗日插值公式 给定n+1个结点的函数表 x x0 x1 x2 ………… xn y Y0 y1 y2 ……… yn 作n次多项式 pn(x) 使满足： pn(xi) = yi (I = 0,1,2,3…n) 同样仿照上述方法, pn(x)为n+1个基本插值多项式线性组合 要求： 由此就可以求得基本插值多项式，以Ai(x) 为例说明其求法。 由此就可以求得基本插值多项式，以Ai(x) 为例说明其求法。 （1）除Ai(xi)以外其他所有结点Ai(x)都为零，因此Ai（x）多项 式中含有下述因子： （2）由Ai(xi)=1 确定λ 最后得n+1个结点的n阶插值多项式为：