椭圆周长和旋转椭球面积近似值的简单算法（第十六届北京高中数学知识应用竞赛论文）.doc

椭圆周长和旋转椭球面积近似值的简单算法 ——第十六届北京高中数学知识应用竞赛论文 论文标题： 椭圆周长和旋转椭球面积近似值的简单算法 作者姓名： 吴欢庆 、吴斯乾（合作） 性别： 男 学校： 北京市通州区第四中学 年级： 高二 指导教师： 曹凤华、李江涛 准考证号： 0311 第 第 PAGE \* MERGEFORMAT 5 页 ， 共 17 页 [内容摘要]：我选这个课题是因为我父亲在工作中遇到了需要计算椭圆周长而又不会算的问题，我为了解决我父亲的困难，所以对这方面做了一些研究。我们用类比推理的方法，提出了椭圆周长近似的简单算法，并推广到怎样计算旋转椭球体的表面积。我们采用简单推导加上用C语言编程计算的方法，让程序运行了约一个小时，得到了一组比较精确的数据，一定程度上解决了计算椭圆周长难的问题。对一些需要计算椭圆周长的工程应用有一定的帮助。 [关键词]：椭圆周长，椭球表面积，椭圆积分，修正因子 本人郑重声明：所呈交的数学应用论文是本人在指导教师的指导下独立进行研究的成果，除文中已经注明引用的内容外，本文不含其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。 论文作者签名：吴欢庆 吴斯乾 2013年3月30日 椭圆周长和旋转椭球面积近似值的简单算法 圆是最美的图形，但自然界中的正圆，少之又少。行星运行的轨道是椭圆的，地球是一个椭球体，在生物界里很难看到完美的圆。许多的图形都是不规则的，椭圆是独特的，与众不同。 众所周知，椭圆积分属于高等微积分的知识，求解椭圆周长的积分存在着很大的难度，在实际的生活应用中给人们带来很大的困扰。当你需要求一个椭圆物体的周长时，总是会没有头绪，仿佛不知道它的计算公式（其实本来就没有公式），没有了这些数据，做起研究就会不知所措。 我的一个朋友以前做关于种子完整度的研究，用到了圆形度来描述种子的形状。他在他的博客中这样写道： “为了检测表皮破裂种子，我调用IMAQ Vision里的形状分析函数得到了面积周长长度宽度等一些形状参数。根据种子图像轮廓，发现用种子的圆形度：，（P为周长，A为面积）可以比较好地区分完整和破碎种子。后来我想，种子的轮廓更接近椭圆，何不用‘椭圆形度’衡量面积和周长的关系呢？可是我记不得椭圆周长的计算方法，百度之，发现我原本就不可能知道。” 椭圆是一个不怎么完美的图形，椭圆周长的积分不能用像圆的周长和面积那样简单的公式表达，所以我们只能计算椭圆周长的近似值。因为没有一个确定的公式，当我们需要精确计算椭圆周长的时候，麻烦就来了。 我父亲是做机床生意的，有次我和父亲去公司，遇到了一位做管道设备的客户，他需要把圆形的管道截成斜面，然后沿着斜面再接上一个椭圆形的管道，这位客户要求把铝板卷成椭圆的管道，至于用多宽的铝板，就需要计算椭圆截面的周长。由于这样一个斜的截面不方便测量，虽然可以做一些试验，来确定铝板的宽度，而这样必定会耗费很多时间和精力。如果材料比较贵的话，用试验的方法就不太可取了。对于普通的工人来说，没有多少数学知识，又不能设计程序计算，怎样求椭圆的周长就成了一个问题。而我，懂得一些微积分，也想帮父亲解决这个问题，所以我想找出一个计算椭圆周长的简单方法。 § 椭圆周长的积分 我们在计算曲线图形的面积和周长通常都要用到微积分，下面我们简单地推导一下椭圆周长的积分表达式： 我们知道椭圆的参数方程为 我们先求四分之一椭圆的周长，取参数为积分变量，它的变化区间为 相应于上任一小区间的小弧段的长度近似等于对应的长度， 因为： 所以：的近似值，即弧长元素为 == = = 其中=（是离心率） 所以椭圆周长C的积分表达式如下： C= 就是上面这个看似简单的积分实际上是很难积出的，如果你觉得你的微积分很厉害，不妨也试试。实际上人们把椭圆周长的积分叫做“第二类不完全椭圆积分”，到此为止，接下来只能将展开成无穷级数，再进行积分，过程稍复杂，这里就不再讨论了。 § 猜想椭圆周长的公式 我们知道，圆的面积：，类似的，椭圆的面积：；圆的周长公式是：，而椭圆的周长却不等于乘上半长轴与半短轴之和。那么椭圆的周长和有什么关系呢？其实很容易想到，当离心率很小的时候，周长是近似等于的；离心率越接近1，给出的结果就越不精确，我们先猜想，下面来看看这个猜想公式的精度如何。 要想知道精确度，首先要得到精确