旋转体表面积和体积的数值计算.pdf

2013 高教社杯全国大学生数学建模竞赛 编 号 专 用 页 赛区评阅编号 （由赛区组委会评阅前进行编号）： 赛区评阅记录 （可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号 （由赛区组委会送交全国前编号）： 全国评阅编号 （由全国组委会评阅前进行编号）： 旋转体表面积和体积的数值计算 摘要 本文针对已给的花瓶侧面投影图，要求在读入图片的基础上，分析其中像素点的特 征，观察其边界曲线的变化轨迹和趋势，利用 matlab 等数学工具，通过积分、微分、 拟合、插值、和误差分析等不同数学模型对得到的数据进行处理，画出该花瓶的三维立 体图形，并计算其表面积和体积。 针对问题一，我们首先通过 matlab 程序实现图片数据的导入，并取出其边界点， 对其边界像素坐标进行分析后，我们发现图片中黑色像素点的每一行的中点完全处于同 一竖直轴线上，即该图形严格对称。于是我们把花瓶的每一个横截面看作圆来处理。分 离出花瓶一侧的曲线，以它的中轴为横轴，以平面上垂直于中轴的方向为纵轴建立坐标 系。采用斜率单变点分析方法，利用斜率变化幅度的二阶差分，找到斜率变化加速度最 大的分界点，将之分成 4 段，然后利用多种曲线拟合模型得到每一段的函数表达式，并 在原图像上取点，使用残差分析的方法，比较各拟合模型的优劣，确定最佳拟合结果。 最后重新确立三维坐标轴，使边界线沿原二维坐标系中图像的对称轴旋转，获得所需旋 转曲面模型的数值表达式，即  2 3 121.21331.2804z 0.0013z  0.00017z (1 z  50)  80.5(51  z  87) 2 2   x  y  2 3 793.0326  20.3933z0.1538z 0.0003959z (88  z 141)   192.3740 0.4033z(142 