牛顿插值法

题目：牛顿插值法在凸轮修正设计中

的应用

算法：Newton插值法

组号：6

组员：赵冬冬闫鹏田二方李婵娟张帅军郑亚军刘洋郭洋波

牛顿插值法在凸轮修正设计中的应用

赵冬冬，闫鹏，田二方，李婵娟，郭洋波，张帅军，郑亚军，刘洋（河南理工大学机械与动力工程学院,河南焦作 454000）

摘要：本文利用牛顿插值法，提出了一种简单实用的凸轮工作轮廓线的修正方法。首先对要进行修正的的曲线附近的一些离散点的数据进行分析处理，确定插值多项式的阶次以满足高精度和低运算量的要求。然后利用Matlab编程计算出插值点的值，并进行误差分析，实现对凸轮的局部工作廓线进行修正。

关键词：凸轮轮廓线；牛顿插值；修正

Interpolation method Newton inthe design of CAM fixed application ZHAO Dongdong,YAN Peng,TIAN Erfang,LI Chanjuan,,GUO Yangbo,ZHANG

Shuaijun,ZHENG Yajun,LIU Yang

(School of Machinery and power engineering Henan polytechnic

uiversity ,Jiaozuo 454000)

Abstract: Based on the Newton interpolation method, we put forward a simple but practical solution to the work of the cam contour correction. Firstly,we rehandle the discrete data nearby the premodifying curve and get the order of the polynomial to meet the demand of high precision and low computation.Then The Newton interpolation and error analysis are realized by matlab programming. SO far ,we’ve resolved the problem of the cam contour correction .

Key words: Newton interpolation; cam contour;correction

0．问题背景

在自动包装机或包装线中，为保证各个机械间歇运动的快捷与准确，常常采用凸轮机构来实现。包装材料、产品和包装地间歇输送、翻转或转移、工作转台的间歇转位，工作机构带停留段的往复运动，有特定位移、速度或加速度要求的动作等，均属于简谐运动范围，正确设计或选用简谐运动机构，对包装机的运行性能具有关键性的作用。凸轮机构在高速包装机械设备中应用更广泛，是一种不可缺少和替代的重要机构。

1．问题分析及模型

高速包装机械中凸轮工作廓线的设计多采用解析法，这样既保证了凸轮的运动特性，又便于对凸轮机构进行运动学和动力学分析，因此这就使得在不同工况下，凸轮设计的解析方程式往往是不相同的。这样虽然能保证凸轮的精度，但同时也对凸轮在实际使用中的修正提高了难度，因为只有建立新的解析方程式才能对凸轮进行修正，尤其是只需对凸轮局部曲线进行修正时，也要建立相应的解析方程，这样就使曲线修正的工作量大增，工作效率降低[1]。

基于减小修正量和提高工作效率的考虑，所以考虑插值法中的牛顿插值，提出了一种简单、实用的凸轮工作廓线的修正设计方法，这种方法不必对原有的解析方程进行修改计算，只需通过对要进行修正的曲线附近的一些离散点的数据进行处理，就能对现有凸轮工作廓线进行修正，特别适合凸轮曲线在实际使用中的局部修正设计[2]。

已知两个变量X,Y 之间的一个离散的函数关系式Yi=F(Xi) i=0，l ，2，…，n ，即给出一 个数据表例如： _________________

作者简介：赵冬冬（1988—），男，河南济源人，硕士研究生，主要从事先进制造方向研究。

邮 箱：451896301@

表1

找到一个与原函数近似的简单的多项式函数，使其在一部分离散点 i x ，0x ，

1x ，……n x 的值分别与()i x f ，()0x f ，()1x f ……()n x f 相等。根据以上要求，

可以使用插值多项式来逼近原解析函数。 2. 算法的数学原理

Lagrange 插值公式结构紧凑，便于理论分析。利用插值基函数也容易得到插值多项式的，Lagrange 插值公式的缺点是，当插值节点增加，或其位置变化时，全部插值基函数均要随之变化，从而整个插值公式的结构也发生变化，这在实际计算中是非常不利的。下面引入的Newton 插值公式也可以克服这个缺点。

Newton 插值多项式可以灵活的增加插值节点进行递推计算。该公式形式对称，结构紧凑，因而容易编写计算程序。 Newton 插值多项式为[3]：

()()[]()x x x x f x f x N k N

k k n ω∑=+=1

,,1,00...

其中：

[][][]0

1101010,...,,,...,,,...,,x x x x x f x x x f x x x f k k k k --=

-

i x 0x 1x ... n x ()i x f ()0x f ()1x f ... ()n x f

称之为()x f 在k x x x ...,,10上得k 阶均差。实际计算中，可以由插商表计算。如图所示：

表2

k x ()k x f 一阶均差 二阶均差 三阶均差 四阶均差

0x ()0x f

1x ()1x f []10,x x f

2x ()2x f []21,x x f []210,,x x x f

3x ()3x f []32,x x f []321,,x x x f []3210,,,x x x x f 4x

()4x f

[]43,x x f []432,,x x x f []4321,,,x x x x f

[]43210,,,,x x x x x f

3. 算法的MATLAB 实现 3.1实验数据

高速包装机上有一凸轮，其工作廓线共有分A 、B 、C 三段，在实际的使用中发现A 段和c 段的行程符合设计要求，而B 段的行程须进行修正设计。 已知凸轮A 段曲线数据，如表3所示。

表3

已知凸轮C 段曲线数据，如表4所示。

表4

一般来说，使用的已测数据点越多，可获得的信息量越大，对未知点的估计越精确。在牛顿插值算法中，使用已测数据点越多，会导致插值多项式的次数过高，出现龙格现象，曲线也会有较多折点。而凸轮的工作多为翻转，需要其工作廓线有很好的圆滑性。

由于距离插值点越近的数据对插值点数值的精度影响越大，本文使用

i x 239 240 241 242 243 ()i x f 14.227 14.03 13.854 13.681 13.526 i x 249 250 251 252 253 ()i x f 13.098 13.095 13.085 13.067 13.039