第一讲 插值法.pptVIP

根据已知的x值，查表可求得y值，但是表内不可能一一列出全部y值，当所求的函数值y正好在两表列数值之间，利用表列数据间的引数求y值的方法称为内插法。 内插法： 利用函数表册，根据任意居间引数查取相应函数的方法。 内插分类： a：按使用目的： b：按引数的个数： 单内插、双内插、三内插 c：按函数的性质： 线性内插、变率内插、高次内插 一．比例单内插 (proportional single interpolation) 比例内插的几何意义： 用表列引数两点的直线代替曲线进行内插，即以弦代替曲线进行内插。 结论： 1）f(x)为线性函数，求得的y值没有误差 2) f(x)为非线性函数，求得的y值有df误差 ②：对非线性函数，表间距越小，利用线性内插求得的函数值的误差越小。但是表的篇幅会增大。 （1）求α＝4′，h＝13.4m时的D*？ 内插的逆运算，，已知求？ ? 当函数有两个自变量时，用比例双内插求近似解。 比例双内插是比例单内插的自然推广。 比例单内插是比例双内插的特殊情形。 比例双内插可按照比例单内插的计算方式进行计算。 例2－1－2：由例2－1－1的计算结果，求h＝13.4m，α＝4?.4时的D？ 例2－1－2：由例2－1－1的计算结果，求h＝13.4m，α＝4?.4时的D？ 例2－1－2：由例2－1－1的计算结果，求h＝13.4m，α＝4?.4时的D？ 例2－1－3：由例2－1－1的计算结果，求h＝13.5m，α＝3?.5时的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 结论： 使用变率内插时，为减小误差，应使用最接近实际引数的表列引数所对应的函数为基准 。 例2－2－2：求h＝13.4m，α＝4?.4时的D？ 查算由非线性函数造的函数表，不论用比例内插还是变率内插都会导致一定的误差。 在精度允许的情况下可采用任意一种方法计算。在航海实际工作中大多数采用比例内插。 航海数学 Nautical Mathematics September 7, 2015 武汉理工大学 航运学院 刘文 第一章 航海数值内插法 插值法是一种古老的数学方法，早在一千多年前的隋唐时期定制历法时就广泛应用了二次插值。刘焯将等距节点的二次插值应用于天文计算。 插值理论却是在17世纪微积分产生后才逐步发展起来的，Newton插值公式理论是当时的重要成果。 由于计算机的使用以及航空、造船、精密仪器的加工，插值法在理论和实践上都得到进一步发展，获得了广泛的应用。 在航海数值计算中，经常要用到一些专用表册，如：航海表、吃水差表等等，这些表册都是按一定的函数关系编排的，如： ? 数值内插在航海工程中的意义 正内插－已知引数求函数； 反内插－已知函数求引数。 第一节 比例内插 (Proportional Interpolation) 引数 函数值 x0 yo x1 y1 … … 1．比例正内插 已知 求 。 ? 比例内插公式： O x0 y0 x1 y1 x y f(x) x y a b c e d ? ? ? O x0 y0 x1 y1 x y f(x) x y a b c e d f O x0 y0 x1 y1 x y f(x) x y a b c e d f ①：只要在误差允许的范围内，均可采用线性内插。 例2－1－1：设物标高h，垂直角α，水平距离D＝h ctgα，利用该式编表如下： α h 10 20 3 4 5 6.2 12.3 4.6 9.3 3.7 7.4 13.4 6.2 ? （2）求α＝5′，h＝13.4m时的D**？ α h 10 20 3 4 5 6.2 12.3 4.6 9.3 3.7 7.4 13.4 5.0 ? 2. 比例反内插 (inverse propo