鞋垫花样图案大全-贺词怎么写

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










论文标题： 椭圆周长和旋转椭球面积近似值的简单算法
作者姓名： 吴欢庆 、吴斯乾（合作）
性别： 男 学校： 北京市通州区第四中学
年级： 高二 指导教师： 曹凤华、李江涛
准考证号： 0311

[内容摘要]：我选这个课 题是因为我父亲在工作中遇到了需要计算椭圆周长而又
不会算的问题，我为了解决我父亲的困难，所以对 这方面做了一些研究。我们用
类比推理的方法，提出了椭圆周长近似的简单算法，并推广到怎样计算旋转 椭球
体的表面积。我们采用简单推导加上用C语言编程计算的方法，让程序运行了约
一个小时， 得到了一组比较精确的数据，一定程度上解决了计算椭圆周长难的问
题。对一些需要计算椭圆周长的工程 应用有一定的帮助。
[关键词]：椭圆周长，椭球表面积，椭圆积分，修正因子








本人郑重声明：所呈交的数 学应用论文是本人在指导教师的指导
下独立进行研究的成果，除文中已经注明引用的内容外，本文不含其
他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要
贡献的个人和集体均已在文 中以明确方式标明。
论文作者签名：吴欢庆 吴斯乾
2013年3月30日

椭圆周长和旋转椭球面积近似值的简单算法

圆是最美的图形，但自然界中的正圆， 少之又少。行星运行的轨道是椭圆的，
地球是一个椭球体，在生物界里很难看到完美的圆。许多的图形都 是不规则的，
椭圆是独特的，与众不同。
众所周知，椭圆积分属于高等微积分的知识，求解椭 圆周长的积分存在着很
大的难度，在实际的生活应用中给人们带来很大的困扰。当你需要求一个椭圆物< br>体的周长时，总是会没有头绪，仿佛不知道它的计算公式（其实本来就没有公式），
没有了这些数 据，做起研究就会不知所措。
我的一个朋友以前做关于种子完整度的研究，用到了圆形度来描述种子的 形
状。他在他的博客中这样写道：
“为了检测表皮破裂种子，我调用IMAQ Vision 里的形状分析函数得到了面
积周长长度宽度等一些形状参数。根据种子图像轮廓，发现用种子的圆形度：
P
2
C
，（P为周长，A为面积）可以比较好地区分完整和破碎种子。后来 我想，
4

A
种子的轮廓更接近椭圆，何不用‘椭圆形度’衡量面积和周长的 关系呢？可是我
记不得椭圆周长的计算方法，百度之，发现我原本就不可能知道。”
椭圆是一 个不怎么完美的图形，椭圆周长的积分不能用像圆的周长和面积那
样简单的公式表达，所以我们只能计算 椭圆周长的近似值。因为没有一个确定的
公式，当我们需要精确计算椭圆周长的时候，麻烦就来了。 < br>我父亲是做机床生意的，有次我和父亲去公司，遇到了一位做管道设备的客
户，他需要把圆形的管 道截成斜面，然后沿着斜面再接上一个椭圆形的管道，这
位客户要求把铝板卷成椭圆的管道，至于用多宽 的铝板，就需要计算椭圆截面的
周长。由于这样一个斜的截面不方便测量，虽然可以做一些试验，来确定 铝板的
宽度，而这样必定会耗费很多时间和精力。如果材料比较贵的话，用试验的方法
就不太可 取了。对于普通的工人来说，没有多少数学知识，又不能设计程序计算，
怎样求椭圆的周长就成了一个问 题。而我，懂得一些微积分，也想帮父亲解决这
个问题，所以我想找出一个计算椭圆周长的简单方法。

§ 椭圆周长的积分
我们在计算曲线图形的面 积和周长通常都要用到微积分，下面我们简单地推
导一下椭圆周长的积分表达式：
我们知道椭 圆的参数方程为
xacos

,ybsin

,
0

2







我们先求四分之一椭圆的周长，取参数

为积 分变量，它的变化区间为

0,



2




相应于

0 ,

上任一小区间

t,tdt

的小弧段的长度
s
近似等于对应的长度

2

(x)
2
( y)
2
，

因为：
xdxasin

d



ydybcos

d


所以：
s
的近似值，即弧长元素为

ds (dx)
2
(dy)
2
a
2
sin
2

(d

)
2
b
2
cos
2

(d

)
2

=
a< br>2
sin
2

b
2
cos
2
< br>d

=
a
2
(1cos
2

) b
2
cos
2

d


=
a
2
a
2
cos
2

b
2
cos
2

d


a
2
b2
cos
2

d

=
a1
2
a

a
2
 b
2
2
其中=
e
（
e
是离心率）
2
a
所以椭圆周长C的积分表达式如下：

C=
4a
2
0
1e
2
cos
2

d



就是上面这个看似简单的积分实际上是很难积出的，如果你觉得你的 微积分
很厉害，不妨也试试。实际上人们把椭圆周长的积分叫做“第二类不完全椭圆积
分”，到 此为止，接下来只能将
1e
2
cos
2

展开成无穷级数 ，再进行积分，过
程稍复杂，这里就不再讨论了。




§ 猜想椭圆周长的公式
我们知道，圆的面积：
S
圆


(rr)
，类似的，椭圆的面积：
S
椭圆


( ab)
；
圆的周长公式是：
C
圆


(rr)
，而椭圆的周长却不等于

乘上半长轴与半短轴
之和。那么椭圆的周长和
(ab)
有什么关系呢？其实很容易想到，当离心率很
小的时候，周长是近似 等于

(ab)
的；离心率越接近1，

(ab)
给出 的结果就
越不精确，我们先猜想
C
椭圆


(ab)，下面来看看这个猜想公式的精度如何。
要想知道精确度，首先要得到精确的周长。我们的方法很 简单，只要借助计
算机将椭圆n等分，把每一个弧长“微元”相加就行了。根据上面推出的积分表
达式，先将四分之一椭圆20万等分。我们采用C语言编写计算不同离心率下椭
圆周长、面积的程序， 程序的原理很简单，其实就是取
d

的值为，这样
200000
10
5

n

22
22
2
4a1ecos 
求
4a

1ecos

d

的近似值 ，就变成了算出

0
2
n1
的值，以下是程序内容：


#include
#include
#define pi 3.97932384626
int main()
{
printf(不同离心率下椭圆的周长和面积--（a=1）n离心率 计算周
长 猜想周长 计算面积 理论面积n
double e,x1,x2,y1,y2;
for(e=0.05;e<1.00;e+=0.05)
{







for(x1=0,y1=0;x1{

}
for(x2=0,y2=0;x2<1;x2+=1.0100000)


{
y2+=4*sqrt((1-e*e)*(1-x2*x2))*1100000;
}
y1+=4*sqrt(1-e*e*cos(x1)*cos(x1))*pi200000;

printf(
sqrt(1-e*e)),y2,pi*sqrt(1-e*e));
}
}

（注：此程序在C-Free 5.0中测试通过，不同的编译器可能需要稍做修改。）

为了粗略地估计程序给出的椭圆周 长的精度，计算周长的同时，程序还将椭
圆分割成10万个矩形求面积。因为

S
椭圆


(ab)
是一个精 确公式，对比计算机得出的面积和公式给出的面积，
可以估计出计算机计算的误差。
下面是运行结果：
--不同离心率下椭圆的周长和面积--（a=1）
离心率 计算周长 猜想周长 计算面积 理论面积
0.05 6.279319266 6.279255859 3.137683143 3.137663205
0.10 6.267510442 6.267437876 3.125865086 3.125845223
0.15 6.247754348 6.247641317 3.106068401 3.106048664
0.20 6.219936382 6.219712246 3.078139152 3.078119592
0.25 6.183890859 6.183426661 3.041853336 3.041834007
0.30 6.139395245 6.138481041 2.996907431 2.996888387
0.35 6.086161920 6.084478629 2.942904675 2.942885975
0.40 6.023826659 6.020909881 2.879335524 2.879317228
0.45 5.951932559 5.947124601 2.805549775 2.805531948
0.50 5.869907460 5.862291700 2.720716335 2.720699046
0.55 5.777031606 5.765339796 2.623763815 2.623747143
0.60 5.672390126 5.654866776 2.513290093 2.513274123
0.65 5.554800713 5.528996378 2.387418895 2.387403724
0.70 5.422698394 5.385138562 2.243560165 2.243545909
0.75 5.273940628 5.219560874 2.077981425 2.077968221
0.80 5.105450038 5.026548246 1.884967570 1.884955592
0.85 4.912481219 4.796529275 1.654947138 1.654936622
0.90 4.686833321 4.510981143 1.369397192 1.369388490
0.95 4.410927819 4.122554645 0.980968225 0.980961992

对比计算面积和理论面积，可以知道计算周长的误差大约 在0.00002以内，
直到离心率为0.5时，

(ab)
给出的值的误 差还控制在0.01以内。
根据以上结果，我们提出椭圆周长近似值的公式：
C
椭圆
p

(ab)
，其中
p
是一个修正因子，为了精确地得 到修正因子和离心率之间的关系，我们将椭圆2
亿等分，程序运行结果如下：
离心率 计算周长 猜想周长 计算面积 理论面积 修正因子
0.05 6.279256474 6.279255859 3.137663220 3.137663205 1.000000098
0.10 6.267447768 6.267437876 3.125845237 3.125845223 1.000001578
0.15 6.247691871 6.247641317 3.106048678 3.106048664 1.000008092
0.20 6.219874184 6.219712246 3.078119607 3.078119592 1.000026036
0.25 6.183829023 6.183426661 3.041834021 3.041834007 1.000065071
0.30 6.139333858 6.138481041 2.996888401 2.996888387 1.000138930
0.35 6.086101074 6.084478629 2.942885989 2.942885975 1.000266653
0.40 6.023766447 6.020909881 2.879317241 2.879317228 1.000474441

0.45 5.951873084 5.947124601 2.805531961 2.805531948 1.000798450
0.50 5.869848833 5.862291700 2.720699059 2.720699046 1.001289109
0.55 5.776973948 5.765339796 2.623747155 2.623747143 1.002017947
0.60 5.672333572 5.654866776 2.513274135 2.513274123 1.003088808
0.65 5.554745415 5.528996378 2.387403735 2.387403724 1.004657091
0.70 5.422644533 5.385138562 2.243545919 2.243545909 1.006964718
0.75 5.273888422 5.219560874 2.077968230 2.077968221 1.010408452
0.80 5.105399760 5.026548246 1.884955601 1.884955592 1.015687010
0.85 4.912433239 4.796529275 1.654936629 1.654936622 1.024164132
0.90 4.686788194 4.510981143 1.369388496 1.369388490 1.038973129
0.95 4.410886572 4.122554645 0.980961996 0.980961992 1.069940110
（由于运算数据较长，这里就只给出了上面这么多数据，更详细的结果在文章最后的附
录A中）
根据上面这些比较精确的数据，可以看出当离心率很小的时候

(ab)
给 出
的结果是很精确的，所以引入修正因子是个不错的选择。
修正因子与离心率的关系图像如下：

根据计算机给出的数据和图像可以得出修正因子和椭圆离心率的函数近似
关系：


1,e(0,0.66)


1.0 1,e0.66,0.80



p

1.02,e 0.80,0.86


0.31e0.76,e0.86,0.90



0.60e0.50,e0.90,1






§ 猜想公式和其它近似公式的对比
查阅资料，我们 找到了一个常见的椭圆周长的近似公式：
C2

b4(ab)
，
我们取了一些不同的离心率值，计算后和我们的猜想公式给出的数值对比，猜想
公式中的修正因子的值 是根据上面给的近似函数关系算出的，对比结果如下：（其
中
a1
）

离心率
e

C2

b4(ab)

Cp

(ab)

两亿等分周长的计算结果

0.05 6.280329540 6.279255859 6.279256474
0.10 6.271740697 6.267437876 6.267447768
0.15 6.257353341 6.247641317 6.247691871
0.20 6.237055596 6.219712246 6.219874184
0.25 6.210684668 6.183426661 6.183829023
0.30 6.178019969 6.138481041 6.139333858
0.35 6.138773151 6.084478629 6.086101074
0.40 6.092573899 6.020909881 6.023766447
0.45 6.038949676 5.947124601 5.869848833
0.50 5.977296478 5.862291700 5.862291700
0.55 5.906835668 5.765339796 5.776973948
0.60 5.826548246 5.654866776 5.672333572
0.65 5.735070617 5.528996378 5.554745415
0.70 5.630520446 5.438989948 5.422644533
0.75 5.510185130 5.271756483 5.273888422
0.80 5.369911184 5.127079211 5.105399760
0.85 5.202742493 4.892459861 4.796529275
0.90 4.995217402 4.691420389 4.510981143
0.95 4.712924384 4.411133470 4.122554645

可以看出我们的猜想公式引入修正因子后比近似公式更接近标准值，而且
一直到
e
为 0.4时，
Cp

(ab)
给出的值误差还在0.005之内。这里的修 正因

子取的只是近似值，如果要更精确的结果，可以查看文章最后的

< br>“附录A”，我们给出了0.01到0.99的共99个离心率下修正因子的较精确的值，
并且精 确到了小数点后9位，足以满足对精度要求较高的计算。
至于前面提到的怎么计算管道截面的 周长，只需用卡尺量出椭圆的长半轴、
短半轴，在附录A中找到对应的修正因子值，就可以很精确地算出 来了。
§ 推广：计算旋转椭球体的表面积
生活中需要计算椭球体的例子并不多，但我们还 是推广一下。举个例子，比
如我们要给一个椭球状的建筑物刷漆，要用多少油漆，就需要计算表面积。而 椭
球表面积的积分也是不完全椭圆积分，也没有具体的公式。
我们一般通过将椭圆绕坐标轴旋 转得到旋转椭球体，再沿一个轴拉伸得到一
般的椭球。因为我们能力有限，就只讨论到旋转椭球的表面积 。
以焦点在x轴的椭圆为例，可以绕焦点所在轴旋转，也可以绕y轴旋转，这
样就得到了两种 椭球，我们称前者为长椭球，后者为扁椭球。以绕x轴旋转得到
长椭圆为例，我们简单推导一下它的面积 积分。

前面已经推导出了椭圆的弧微分元素
dsa1e
2
c os
2

，将每一个弧元素
ds
绕x轴旋转，得到面积元素
dA
：
dA2

yds2

absin

1e
2
cos
2

d


< br>积分可以得到半个椭球面积的积分表达式：
2


2
absi n

1e
2
cos
2

d


0

所以整个椭球的面积
A4


2
a bsin

1e
2
cos
2

d

0

这个积分比椭圆周长的积分更难积出。而且对于怎样计算椭球的表面积，我
们查阅资料，并没有找到近似公式。

和前面类比推出椭圆 周长的近似公式一样，对比球的表面积是
4

r
2
，有了
上面求椭圆周长的经验，当离心率很小时，我们推测旋转椭球的面积近似等于
4

(< br>ab
2
ab
2
)
，而对于较大的离心率，它的面积
Aq4

()
，其中
q
也是一
2
2
个修正因子。
为了得到较精确的面积值，我们还是采用C语言编写程序计算，其实只要对
前面 的程序稍做修改就行了，程序内容如下：

#include
#include
#define pi 3.97932384626
int main()
{
printf(不同离心率下由椭圆绕x轴旋转所得椭球的表面积--（a=1）n
离心率 计算面积 猜想面积 修正因子n
double e,x1,y1;
for(e=0.01;e<1.00;e+=0.01)
{
for(x1=0,y1=0;x1 {

y1+=4*pi*sqrt(1-e*e)*sin(x1)*sqrt(1-e*e*cos(x 1)*cos(x1))*pi20000
0000;
}

pri ntf(
e*e))*(1+sqrt(1-e*e)),y1(pi*(1+sqrt(1-e*e) )*(1+sqrt(1-e*e))));
}
}

运行的结果详见“附录B”
我们得到了修正因子
q
和用来旋转的椭圆的离心率之间的图像关系：

§ 总结
生活中、生产中涉及椭圆的计算也许不多，但 是要是遇到了，那就是一个坎
儿。椭圆积分是积分学中重要的一个分支，它最初出现在有关椭圆弧长的问 题中，
关于椭圆弧长的积分，最早的研究者正是大名鼎鼎的数学巨匠——欧拉。而本文
的内容主 要是简单推导和巨量的计算，其中这些巨量计算和数据分析都是由计算
机完成的，计算机的发展使得我们 ，21世纪的中学生能够得到几百年前的大数
学家算不出来的数据。
数学推动了计算机的发展 ，计算机也必将推动数学的发展，事实即是如此，
本文就是一个很好的例子。今天的计算机有着强大的数 据处理能力，并且体积更
小，更便携，这让今天的应用数学更离不开计算机。
也许这篇文章理 论性的东西很多，我想说的是，理论和应用是紧密相连的，
应用需要理论指导，而理论必然要应用到实践 中去，它们都是科学发展中必不可
少的一部分。

附录A
--不同离心率下椭圆的周长和面积（a=1）--
离心率 计算周长 猜想周长 计算面积 理论面积
0.01 6.283028225 6.283028224 3.141435585 3.141435570 1.000000000
0.02 6.282556941 6.282556926 3.140964287 3.140964272 1.000000002
0.03 6.281771352 6.281771272 3.140178633 3.140178619 1.000000013
0.04 6.280671279 6.280671027 3.139078388 3.139078373 1.000000040
0.05 6.279256474 6.279255859 3.137663220 3.137663205 1.000000098
0.06 6.277526618 6.277525342 3.135932703 3.135932688 1.000000203
0.07 6.275481319 6.275478953 3.133886314 3.133886300 1.000000377
0.08 6.273120115 6.273116074 3.131523435 3.131523420 1.000000644
0.09 6.270442468 6.270435987 3.128843348 3.128843333 1.000001033
0.10 6.267447768 6.267437876 3.125845237 3.125845223 1.000001578
0.11 6.264135331 6.264120826 3.122528187 3.122528172 1.000002316
0.12 6.260504398 6.260483818 3.118891179 3.118891165 1.000003287
0.13 6.256554132 6.256525733 3.114933094 3.114933079 1.000004539
0.14 6.252283620 6.252245343 3.110652704 3.110652690 1.000006122
0.15 6.247691871 6.247641317 3.106048678 3.106048664 1.000008092
0.16 6.242777813 6.242712214 3.101119575 3.101119560 1.000010508
0.17 6.237540293 6.237456480 3.095863841 3.095863827 1.000013437
0.18 6.231978076 6.231872450 3.090279811 3.090279796 1.000016949
0.19 6.226089842 6.225958339 3.084365700 3.084365686 1.000021122
0.20 6.219874184 6.219712246 3.078119607 3.078119592 1.000026036
0.21 6.213329608 6.213132145 3.071539506 3.071539491 1.000031782
0.22 6.206454530 6.206215885 3.064623245 3.064623231 1.000038453
0.23 6.199247270 6.198961184 3.057368544 3.057368530 1.000046151
0.24 6.191706058 6.191365627 3.049772987 3.049772973 1.000054985
0.25 6.183829023 6.183426661 3.041834021 3.041834007 1.000065071
0.26 6.175614195 6.175141589 3.033548950 3.033548935 1.000076534

修正因子

0.27 6.167059500 6.166507568 3.024914929 3.024914915 1.000089505
0.28 6.158162758 6.157521601 3.015928962 3.015928947 1.000104126
0.29 6.148921678 6.148180532 3.006587893 3.006587878 1.000120547
0.30 6.139333858 6.138481041 2.996888401 2.996888387 1.000138930
0.31 6.129396777 6.128419635 2.986826996 2.986826982 1.000159444
0.32 6.119107791 6.117992645 2.976400006 2.976399992 1.000182273
0.33 6.108464134 6.107196215 2.965603575 2.965603561 1.000207611
0.34 6.097462906 6.096026293 2.954433654 2.954433640 1.000235664
0.35 6.086101074 6.084478629 2.942885989 2.942885975 1.000266653
0.36 6.074375461 6.072548755 2.930956115 2.930956102 1.000300814
0.37 6.062282747 6.060231987 2.918639347 2.918639333 1.000338396
0.38 6.049819457 6.047523404 2.905930764 2.905930750 1.000379668
0.39 6.036981957 6.034417842 2.892825202 2.892825189 1.000424915
0.40 6.023766447 6.020909881 2.879317241 2.879317228 1.000474441
0.41 6.010168952 6.006993830 2.865401190 2.865401176 1.000528571
0.42 5.996185316 5.992663713 2.851071073 2.851071059 1.000587652
0.43 5.981811192 5.977913254 2.836320613 2.836320600 1.000652057
0.44 5.967042034 5.962735860 2.821143219 2.821143206 1.000722181
0.45 5.951873084 5.947124601 2.805531961 2.805531948 1.000798450
0.46 5.936299365 5.931072195 2.789479554 2.789479541 1.000881320
0.47 5.920315666 5.914570978 2.772978338 2.772978325 1.000971277
0.48 5.903916532 5.897612891 2.756020250 2.756020237 1.001068846
0.49 5.887096250 5.880189445 2.738596804 2.738596791 1.001174589
0.50 5.869848833 5.862291700 2.720699059 2.720699046 1.001289109
0.51 5.852168005 5.843910232 2.702317591 2.702317579 1.001413056
0.52 5.834047184 5.825035101 2.683442460 2.683442448 1.001547129
0.53 5.815479463 5.805655814 2.664063173 2.664063161 1.001692083
0.54 5.796457591 5.785761286 2.644168645 2.644168633 1.001848729
0.55 5.776973948 5.765339796 2.623747155 2.623747143 1.002017947
0.56 5.757020524 5.744378942 2.602786300 2.602786288 1.002200687
0.57 5.736588892 5.722865582 2.581272941 2.581272928 1.002397979
0.58 5.715670176 5.700785784 2.559193143 2.559193131 1.002610937
0.59 5.694255024 5.678124756 2.536532115 2.536532103 1.002840774
0.60 5.672333572 5.654866776 2.513274135 2.513274123 1.003088808
0.61 5.649895401 5.630995113 2.489402471 2.489402459 1.003356474
0.62 5.626929503 5.606491935 2.464899293 2.464899282 1.003645340
0.63 5.603424226 5.581338217 2.439745575 2.439745563 1.003957117
0.64 5.579367228 5.555513621 2.413920979 2.413920968 1.004293682
0.65 5.554745415 5.528996378 2.387403735 2.387403724 1.004657091
0.66 5.529544880 5.501763140 2.360170497 2.360170486 1.005049607
0.67 5.503750828 5.473788826 2.332196183 2.332196172 1.005473723
0.68 5.477347495 5.445046436 2.303453793 2.303453782 1.005932192
0.69 5.450318058 5.415506844 2.273914201 2.273914190 1.006428062
0.70 5.422644533 5.385138562 2.243545919 2.243545909 1.006964718

0.71 5.394307654 5.353907464 2.212314821 2.212314811 1.007545926
0.72 5.365286741 5.321776471 2.180183828 2.180183818 1.008175892
0.73 5.335559548 5.288705183 2.147112540 2.147112530 1.008859326
0.74 5.305102086 5.254649451 2.113056808 2.113056798 1.009601522
0.75 5.273888422 5.219560874 2.077968230 2.077968221 1.010408452
0.76 5.241890439 5.183386203 2.041793559 2.041793550 1.011286876
0.77 5.209077564 5.146066637 2.004473993 2.004473984 1.012244483
0.78 5.175416443 5.107536977 1.965944333 1.965944324 1.013290059
0.79 5.140870556 5.067724605 1.926131960 1.926131951 1.014433687
0.80 5.105399760 5.026548246 1.884955601 1.884955592 1.015687010
0.81 5.068959742 4.983916445 1.842323800 1.842323792 1.017063548
0.82 5.031501348 4.939725686 1.798133041 1.798133032 1.018579101
0.83 4.992969770 4.893858024 1.752265379 1.752265371 1.020252272
0.84 4.953303535 4.846178098 1.704585452 1.704585444 1.022105138
0.85 4.912433239 4.796529275 1.654936629 1.654936622 1.024164132
0.86 4.870279947 4.744728638 1.603135992 1.603135985 1.026461220
0.87 4.826753116 4.690560319 1.548967673 1.548967666 1.029035507
0.88 4.781747892 4.633766488 1.492173842 1.492173835 1.031935447
0.89 4.735141481 4.574034888 1.432442241 1.432442234 1.035221986
0.90 4.686788194 4.510981143 1.369388496 1.369388490 1.038973129
0.91 4.636512481 4.444122882 1.302530235 1.302530229 1.043290792
0.92 4.584098836 4.372840491 1.231247843 1.231247837 1.048311469
0.93 4.529276567 4.296314938 1.154722290 1.154722285 1.054223592
0.94 4.471695717 4.213423775 1.071831126 1.071831121 1.061297405
0.95 4.410886572 4.122554645 0.980961996 0.980961992 1.069940110
0.96 4.346185831 4.021238597 0.879645947 0.879645943 1.080807748
0.97 4.276585853 3.905329271 0.763736621 0.763736617 1.095064092
0.98 4.200369021 3.766761698 0.625169047 0.625169045 1.115114084
0.99 4.113903209 3.584768835 0.443176183 0.443176181 1.147606275



附录B
--不同离心率下由椭圆绕x轴旋转所得椭球的表面积（a=1）--
离心率 计算面积 猜想面积 修正因子
0.01 12.565532750 12.565742288 0.999983325
0.02 12.563019350 12.563857215 0.999933312
0.03 12.558830015 12.560715111 0.999849921
0.04 12.552964241 12.556315506 0.999733101
0.05 12.545421322 12.550657736 0.999582778
0.06 12.536200348 12.543740950 0.999398855
0.07 12.525300202 12.535564103 0.999181217

0.08 12.512719562 12.526125955 0.998929725
0.09 12.498456895 12.515425074 0.998644219
0.10 12.482510455 12.503459826 0.998324514
0.11 12.464878283 12.490228381 0.997970406
0.12 12.445558199 12.475728703 0.997581664
0.13 12.424547799 12.459958549 0.997158036
0.14 12.401844456 12.442915470 0.996699245
0.15 12.377445307 12.424596800 0.996204988
0.16 12.351347255 12.404999656 0.995674937
0.17 12.323546958 12.384120933 0.995108738
0.18 12.294040826 12.361957298 0.994506010
0.19 12.262825013 12.338505183 0.993866342
0.20 12.229895410 12.313760786 0.993189296
0.21 12.195247638 12.287720054 0.992474404
0.22 12.158877036 12.260378685 0.991721165
0.23 12.120778658 12.231732116 0.990929048
0.24 12.080947256 12.201775517 0.990097485
0.25 12.039377275 12.170503780 0.989225877
0.26 11.996062837 12.137911515 0.988313584
0.27 11.950997733 12.103993032 0.987359932
0.28 11.904175407 12.068742338 0.986364202
0.29 11.855588942 12.032153122 0.985325637
0.30 11.805231046 11.994218743 0.984243434
0.31 11.753094034 11.954932216 0.983116744
0.32 11.699169810 11.914286203 0.981944668
0.33 11.643449851 11.872272989 0.980726257
0.34 11.585925184 11.828884476 0.979460507
0.35 11.526586364 11.784112157 0.978146356
0.36 11.465423452 11.737947103 0.976782682
0.37 11.402425989 11.690379940 0.975368298
0.38 11.337582965 11.641400829 0.973901950
0.39 11.270882797 11.590999442 0.972382309
0.40 11.202313290 11.539164938 0.970807970
0.41 11.131861606 11.485885935 0.969177447
0.42 11.059514225 11.431150482 0.967489164
0.43 10.985256907 11.374946026 0.965741453
0.44 10.909074645 11.317259382 0.963932546
0.45 10.830951621 11.258076690 0.962060565
0.46 10.750871153 11.197383384 0.960123520
0.47 10.668815637 11.135164139 0.958119297
0.48 10.584766490 11.071402834 0.956045647
0.49 10.498704080 11.006082493 0.953900181
0.50 10.410607654 10.939185236 0.951680352
0.51 10.320455261 10.870692215 0.949383448

0.52 10.228223662 10.800583549 0.947006578
0.53 10.133888237 10.728838252 0.944546651
0.54 10.037422877 10.655434155 0.942000366
0.55 9.938799871 10.580347815 0.939364192
0.56 9.837989779 10.503554427 0.936634341
0.57 9.734961286 10.425027711 0.933806754
0.58 9.629681053 10.344739800 0.930877068
0.59 9.522113537 10.262661110 0.927840590
0.60 9.412220801 10.178760198 0.924692263
0.61 9.299962300 10.093003599 0.921426631
0.62 9.185294634 10.005355655 0.918037794
0.63 9.068171282 9.915778310 0.914519365
0.64 8.948542298 9.824230892 0.910864412
0.65 8.826353962 9.730669859 0.907065401
0.66 8.701548397 9.635048520 0.903114123
0.67 8.574063123 9.537316709 0.899001615
0.68 8.443830553 9.437420428 0.894718066
0.69 8.310777416 9.335301426 0.890252713
0.70 8.174824098 9.230896724 0.885593712
0.71 8.035883868 9.124138072 0.880727999
0.72 7.893862000 9.014951311 0.875641113
0.73 7.748654732 8.903255642 0.870316999
0.74 7.600148059 8.788962766 0.864737770
0.75 7.448216306 8.671975881 0.858883420
0.76 7.292720436 8.552188490 0.852731490
0.77 7.133506039 8.429482990 0.846256650
0.78 6.970400908 8.303728984 0.839430203
0.79 6.803212112 8.174781235 0.832219471
0.80 6.631722419 8.042477193 0.824587035
0.81 6.455685901 7.906633951 0.816489791
0.82 6.274822452 7.767044471 0.807877755
0.83 6.088810918 7.623472869 0.798692541
0.84 5.897280331 7.475648419 0.788865393
0.85 5.699798624 7.323257858 0.778314616
0.86 5.495857837 7.165935350 0.766942146
0.87 5.284854428 7.003249159 0.754628931
0.88 5.066062526 6.834683626 0.741228534
0.89 4.838596803 6.659614235 0.726558121
0.90 4.601359578 6.477272238 0.710385392
0.91 4.352963142 6.286692888 0.692409065
0.92 4.091611510 6.086636959 0.672228611
0.93 3.814912339 5.875466390 0.649295237
0.94 3.519561058 5.650936281 0.622827950
0.95 3.200772184 5.409821921 0.591659436

0.96 2.851156450 5.147185404 0.553925345
0.97 2.458196630 4.854734013 0.506350425
0.98 1.997324521 4.516337812 0.442244270
0.99 1.404656777 4.090462710 0.343398015









参考文献：
1.同济大学数学系编，《高等数学 第六版》，高等教育出版社
2.刘里鹏 著，《从割圆术走向无穷小——揭秘微积分》，湖南科学技术出版社
3.[美]Jesse Liberty,Siddhartha Rao,Bradley Jones著，袁国忠 陈秋萍译，
《21天学通C++（第6版）》，人民邮电出版社



合作经历：
这个课题我们用了两个星期左右的时期来完成。我们一个是理科生（吴欢 庆）
一个是文科生（吴斯乾），文中的一些积分式子是吴欢庆完成的；大量的数据是
吴斯乾处理 的。我们两人都没有多少编程的基础，关于程序我们还请教了安徽工
业大学的一位电器学老师。翻阅大量 的资料，我们俩人合作完成了这篇论文。虽
然这篇文章写得不是很有说服力，一些观点都是个人的，而文 中的猜想也只是来
自对实际情况的分析，几乎没有什么理论依据，但是经过我们两周不断地计算，
得出来的结果还是比较接近实际的。在这里还要感谢支持我们的数学老师，曹凤
华老师在内容以及这些 方法的实际应用上给了我们很大的启发；在论文格式上，
李江涛老师辅导我们完成写作，并给我们提了一 些建议。他们耐心地帮我们检查

了这么长的一篇论文，在此，我们特别致谢！

（另附.cpp文件和编译的计算程序已上传到“其他附件”中。）




（注：本资料素材和资料部分来自网络，仅供参考。请预览后才下载，期待你
的好评与关注！）