丁伯伦 张力娟 周金倩 朱贝贝
(扬州工业职业技术学院,江苏 扬州 225127)
以2020年全国大学生数学建模竞赛D题第一问为例,根据某一工件1在水平状态下的测量数据(数字信号),来标注出被测工件表面轮廓线的各项参数值。

图1 某工件1在水平状态下测量的轮廓线
假设被测工件的轮廓线是由直线和圆弧构成的光滑曲线,首先根据工件1的水平测量数据(数字信号),运用Matlab软件可得工件x轴与z轴数字信号的平面图,如图2所示:

图2 工件1在水平状态下的数字信号坐标图
可设定z轴初始点为-1.770,那么近似水平状态轮廓线的表达式为z0=-1.770,再以第一段倾斜线为例,运用Excel进行数据拟合,可得轮廓线D1的表达式为z1=-2.7468x+134.95。采用类似的方法可得每段倾斜线的表达式,如表1所示:

表1 每段倾斜线的表达式
结合“人字形”两边轮廓线的表达式,求出相应交点坐标,即可得到“人字形”轮廓线的长度。如表1,已知“人字形”两边的轮廓线D7、D8的表达式,可求出轮廓线D7、D8的交点坐标为(76.810,-0.782),运用Excel进行数据拟合,可求得轮廓线D7、D8的长度,D7=5.235,D8=5.147,结合之前水平线z0=-1.770,得到人字线性的高度为0.998。

图3 工件1在水平状态下的数字信号坐标图的局部图
运用Excel软件拟合出工件1数字信号的局部图,利用水平线段的表达公式:z0=-1.770与轮廓线D1的表达公式z1=-2.7468x+134.95求出交点A1的横坐标为49.772,同理利用水平线段的表达公式z0=-1.770与轮廓线D2的表达公式z2=2.779x-148.05,求出A2的横坐标为52.676,再用A2的横坐标减去A1的横坐标,即可求出槽口宽度x1=2.904;利用水平线段的表达公式z0=-1.770与轮廓线D3的表达式z3=-3.7092x+212.1,求出交点A3的横坐标为57.653,再利用A3的横坐标减去A2的横坐标,即可求出水平线段长度为x2=4.977。……