第二章 内插法.ppt

第二章 内插法 在航海数值计算中，经常要用到一些专用表册，如：航海表、吃水差表等等，这些表册都是按一定的函数关系编排的，如： y＝f（x） 根据已知的x值，查表可求得y值，但是表内不可能一一列出全部y值，当所求的函数值y正好在两表列数值之间，利用表列数据间的引数求y值的方法称为内插法。 内插法： 利用函数表册，根据任意居间引数查取相应函数的方法。 内插分类： a：按使用目的： 正内插－已知引数求函数； 反内插－已知函数求引数。 b：按引数的个数： 单内插、双内插、三内插 c：按函数的性质： 线性内插、变率内插、高次内插 第一节 比例内插（线性内插） 一．比例单内插（一元函数 y＝f(x) ） 比率内插的几何意义： 用表列引数两点的直线代替曲线进行内插，即以弦代替曲线进行内插。 结论： 1） f（x）为线性函数，求得的y值没有误差 2) f（x）为非线性函数，求得的y值有df误差 ①：只要在误差允许的范围内，均可采用线性内插。 ②：对非线性函数，表间距越小，利用线性内插求得的函数值的误差越小。但是表的篇幅会增大。 例2－1－1：某船静水力性能数据表部分内容如表2-1-2： 2. 比例反内插 内插的逆运算，y＝f（x），已知y求x？ 二. 比例双内插（二元函数） 当函数有两个自变量时，用比例双内插求近似解。 例2－1－1：设物标高h，垂直角α，水平距离D＝h ctgα，利用该式编表2-1-2如下： （1）求α＝4′，h＝13.4m时的D？ （2）求α＝5′，h＝13.4m时的D？ 求h＝13.4m，α＝4?.4时的D？ 比例双内插的简便算法 例2－1－3：求h＝13.4m，α＝4?.4时的D？ 第二节 变率内插 当函数是非线性函数时，如果用比例内插计算将会导致一定的计算误差， 为了尽量减小该误差，则引进了变率内插。 一. 变率单内插（一元函数） 已知比例内插计算公式： （1）?? 用比例内插 y＝5.5 （2）?? 用x＝2变率内插 y＝4＋4（2.3－2）＝5.2 （3）?? 用x＝3变率内插 y＝9＋6（2.3－3）＝4.8 （4）?? 用y＝x2直接计算 y＝5.29 分析：① 比例内插误差大；② x＝2的变率内插较准。 结论： 使用变率内插时，为减小误差，应使用最接近实际引数的表列引数所对应的函数为基准 。 综上所述 查算由非线性函数造的函数表，不论用比例内插还是变率内插都会导致一定的误差。 在精度允许的情况下可采用任意一种方法计算。在航海实际工作中大多数采用比例内插。 * * … … y1 x1 yo x0 函数值 引数 1．比例正内插 已知 x 求 y。 y＝f（x） 比例内插公式： y y y y - - 0 1 0 ＝ x x x x - - 0 1 0 y ＝ y 0 ＋ x x x x - - 0 1 0 （ y 1 － y 0 ） ＝ y 0 ＋ y y x x 1 0 1 0 - - （ x － x 0 ） O x0 y0 x1 y1 x y f(x) x y a c d e b O x0 y0 x1 y1 x y f(x) x y a b c e d f 9600 4035 10040 4475 10500 4935 10960 5395 排水量△（t） 总载重量DW（t） 5.00 5.20 5.40 5.60 型吃水d（m） 已知型吃水5.25，求排水量和总载重量 ＝ △ 0 ＋ △ △ d d 1 0 1 0 - - （ d － d 0 ） ＝10155t △ ＝ △ 0 ＋ △ △ d d 1 0 1 0 - - （ d － d 0 ） ＝4590t DW 6.2 12.3 4.6 9.3 3.7 7.4 10 20 3 4 5 α h 13.4 6.2 D ＝D0＋