matlab求非圆齿轮的节曲线,MATLAB样条曲线在非圆齿轮离散节曲线拟合中的应用

第3 期 (总第142 期) 2007 年6 月 机 械 工 程 与 自 动 化 M ECHAN ICAL 　EN GIN EER IN G　&　AU TOMA T ION No3 Jun1文章编号 : 167226413(2007) 030007202 MA TLAB 样条曲线在非圆齿轮离散节曲线拟合中的应用 曲国毅1, 王建华2, 廖　璇1 (1. 北京机械工业学院 机械工程系, 北京　100085; 2. 北京机械工业学院 基础部, 北京　100085) 摘要: 介绍一种只需要知道非圆齿轮齿顶上的一些离散点, 就可以用具有统一表达式的分段三次样条曲线来代替表达式形式各异的非圆齿轮节曲线的方法。以插秧机上的一个非圆齿轮为例, 利用MA TLAB 自带的计算工具箱和绘图功能来说明该方法的实用性。 关键词: 分段三次样条曲线; MA TLAB; 离散点 中图分类号: TP30114　　　文献标识码: A收稿日期 : 20061217; 修回日期: 200701220作者简介 : 曲国毅(19802) , 男, 山东烟台人, 硕士研究生。 0　引言 设计一对非圆齿轮时, 首先要先根据其满足的传动关系和几何关系得到它们的节曲线方程, 再根据选用的加工方法计算出非圆齿轮的加工数据。为了计算这些加工数据, 要根据节曲线方程推导出大量繁琐的公式。同时, 许多非圆齿轮的节曲线是由几段不同曲线组成, 每一段曲线的公式是根据一种特定的非圆齿轮推导出来的, 对于具有不同节曲线的非圆齿轮, 计算公式都要重新推导, 这样使得非圆齿轮加工数据的计算非常复杂并且常常无法通用。 本文采用的是通过非圆齿轮齿顶的测量数据点, 用离散点拟合的方法对非圆齿轮节曲线进行拟合。由于拟合后的方程只有一种形式, 只是系数不同, 这样就会使得非圆齿轮加工数据计算易于通用, 大大简化了非圆齿轮的计算。由于在计算非圆齿轮加工数据过程中要用到节曲线的二阶导数, 同时要保证给出的节曲线上的所有点都在拟合后的曲线上, 所以采用分段三次样条曲线拟合方法, 它可以非常好地满足以上对节曲线拟合方法的要求。 1　样条曲线拟合方法 样条曲线是一条没有误差地经过数据点的光滑曲线, 在数据点的每一个区间可构造一个三次函数, 使得分段曲线与它的一阶导数和二阶导数在区间内连续, 两相邻段在节点处二阶连续。用样条曲线对数据 进行处理的方法在各个领域得到了广泛的应用。 设{ (x k, y k) }N k= 0有N + 1 个点, 其中x 0< x 1< ⋯< xN , 如果存在N 个三次多项式: S (x ) = S k (x ) = sk, 0+ sk, 1 (x - x k) + sk, 2 (x - x k) 2+ sk, 3 (x - x k) 3 。 (1)⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯ 其中, x ∈[x k, x k+ 1 ]　k= 0, 1, ⋯,N - 1, S (x k) = y k　k= 0, 1, ⋯,N , S k (x k+ 1) = S k+ 1 (x k+ 1)　k= 0, 1, ⋯,N - 2, S k′( x k+ 1) = S k+ 1′( x k+ 1)　k= 0, 1, ⋯,N - 2, S k" (x k+ 1) = S k+ 1" (x k+ 1)　k= 0, 1, ⋯,N - 2。则称函数为三次样条函数。根据上面三次样条函数的表达式和性质, 经过整 理可以得到下面的关系式: hk- 1m k- 1+ 2(hk- 1+ hk)m k+ hkm k+ 1= uk　k= 1, 2,