教育教学故事-张信哲的歌曲

Team # 1739
第五届“认证杯”数学中国
数学建模网络挑战赛
承 诺 书

我们仔细阅读了第五届“认证杯”数学中国数学建模网络挑战赛的竞赛规则。
我们完全明白， 在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮
件、网上咨询等）与队外的任何人（包括指导 教师）研究、讨论与赛题有关的问
题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果 引用别人的成果或其他
公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正< br>文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平 性。如有违反
竞赛规则的行为，我们将受到严肃处理。
我们允许数学中国网站()公布论文， 以供网友之间学习交流，
数学中国网站以非商业目的的论文交流不需要提前取得我们的同意。
我们的参赛队号为：1739
参赛队员 ：


队员2：

队员3：

参赛队教练员 ：

参赛队伍组别：


队员1：

Team # 1739
第五届“认证杯”数学中国
数学建模网络挑战赛
编 号 专 用 页



参赛队伍的参赛队号：
# 1739








竞赛统一编号（由竞赛组委会送至评委团前编号）：















竞赛评阅编号（由竞赛评委团评阅前进行编号）：

Team # 1739
2012年第五届“认证杯”数学中国
数学建模网络挑战赛



题 目
蜘蛛网最优模型
关 键 词 蜘蛛网 螺旋线 层次分析法 能量守恒
摘 要：
世界上有很多种类的蜘蛛，而大部分种类的蜘蛛都是通过结网来捕食食物的。
蜘蛛网的结 构主要受外因和内因的影响，本文研究并建立合理的数学模型，确定
出了蜘蛛网的最优结构，得出结果与 实际情况基本相符。
首先，蜘蛛织网是由一根丝完成整个工作，为了构造合理的网状，设蜘蛛只有材料
L
，已知
L
的情况下易知圆的面积最大，即捕食面的面积最大值，得 到最大
圆面积为
nd
1
S
max


d< br>2

rdr
1
，其中
d
1
为捕丝间距。 < br>其次，由以上可以确定
S
面积只与
d
1
有关，根据相关生物学 家的资料查得影响
蜘蛛网捕丝间距的因素，主要有风、湿度、空间、温度、食物状况等。再运用层
次分析的方法计算出各个因素对
d
1
的权重，从而计算出捕丝的间距。然而，以上< br>并没有计算出放射丝的间距。因此本文再利用能量守恒的关系
11
mv
2
E
p
kx
2
22
，即假设食物的最大垂直速度与质量已知， 求解与该蜘蛛丝
所受的最大弹性形变
x
。随后在做合理的假设与计算，就可以得到最 大的放射
丝间距
L
a
。
最后，针对蜘蛛网本身的的螺旋问题，使用 极坐标建立基本模型

d
1

r
1
，
其中

变化则其放射丝的半径变化,由此模拟出捕丝的的螺旋线图形。蜘蛛网的
最外的 边框（框丝、停泊丝）根据各个外点对整个平面图根据空间的具体关系可
以得出，但是要保证整个蜘蛛网 为平面状，那么边框至少为五边形。
参赛队号 1739

参赛密码
所选题目 A题




（由组委会填写）

Team # 1739
The best structure of cobweb
There are many types of spiders in the world, and most species of spiders prey
food by webs. This article study that the structure of spider web is mainly affected by
the external and internal factors，and a rational mathematical model is established. So
we derive the optimal structure of the spider web, which is consistent with the real
situation.
First, the spider web is constructed by a wire. To established a better model for the
spider web, we suppose that under the condition that the length of spider silk is
L
,
we compute the biggest circular area. So we can find out the biggest measure of area
is
nd
1
S
max


d
2

rdr
1
, where
d
1
means the spacing of spider web.
Second, in a word, we know that the
S
just related to
d
1
. According to some
biologists’ reporting, we can see the reasons of spacing of spider web such as,
wind、 humidity、space、temperature、size of food. Then using the method of
Hierarchical analysis, we can find out

d
1
0.07 86x
1
0.2075x
2
0.2312x
3
0.4 782x
4
0.8175x
5
.
Hence, the answer for our mathematical model is gotten.
Third, above all, we haven’t found out the spacing of Radial wire, so this article
11
mv
2
E
p
kx
2
2
uses The Energy Conservation which is
2
, and we suppose
the fastest speed and quality are known, the purpose of us is to find out the biggest
deformation
x
.Moreover when we do some reasonable hypotheses and calculations,
we can get the biggest spacing of Radial wire
L
a
.

Fourth, in the light of spider web’s spiral question, with some known results, we
study the basic model is

d
1

r
1
. So we can use Matlab to get the figure of
Spiral curve. However, in order to make the plan of web be supported, we need to
know that what a space state is. When we think of the gravity, we can draw a
conclusion that the polygon at least is a pentagon.

Team # 1739
1、 问题重述
蜘蛛目共有105个，可约37000种，世界各地均有分布。蜘蛛主要有结 网、
游猎和穴居3种生活类型。蛛网的进化经历了拌丝、偏网和圆网阶段，并在圆网
的基础上， 继续进化形成了其他类型的网。其中圆网在蛛网进化的地位较为特殊，
且其结构较其它种类比较简单、规 则，因此，本文在圆网的基础上分析、改进、
求解。现在建立合理的数学模型，对蜘蛛网的最大面积、捕 丝和放射丝间距合理
性以及整个蜘蛛网的稳定性的角度分析，并在有外来事物的冲击下，怎样的蜘蛛网才是最合适的。

图1
2、 符号说明

S
：捕丝的面积。
n
：捕丝的圈数。
d
：捕丝间距（网眼）。
‘
L
：捕丝最外圈的周长。
L
a
：放射丝的最大间距。
L
a
'
：捕丝间每一个小面积的切向长度
。
a
：放射丝的条数（半径的条数）

：蜘蛛从起始点螺旋运动到极心绕的
n
圈的 总共角度求和。

：相邻放射线的夹角。
k
：捕丝的弹性系数。
x
：捕丝的伸长量。

：螺旋线半径。
- 1 -

Team # 1739
L
：材料总长度。
m
：外来食物的质量。
v
：外来食物垂直冲到蜘蛛网上。

3、 问题分析

图2蜘蛛网的结构图

1、当一只蜘蛛的材料一定时，假设该材料总长为
L
，我们需要建立一个捕丝
结构有利于节约材料，且达到最大捕食效果，那么我们易知长度一定时 ，其围成
一个圆的面积都比其他图形面积大，由此：
由简单到难，先假设其捕丝间的距离为< br>d
1
，当捕丝圈数为
n
圈时，那么其最
大的半径可以确定为< br>nd
1
。那么此种情况选的
S
max
可以求出，为一个与d< br>1
有关的
式子。
2、根据相关生物学家的调查可知，
其一：蜘蛛丝 具有多级结构的蛋白质纤维，外观呈金黄色，透明，横截面为
圆形，具有皮芯层结构，芯层内含有数十根 纳米级的微纤维。蜘蛛丝具有高强
度、高韧性、高弹性和良好的耐热性能，被称为“生物钢”。而且蜘蛛 为变温
动物，温度过低时蜘蛛新陈代谢较慢，只有温度达到一定了之后，蜘蛛新城代
谢快且蜘蛛 丝的黏性最好。
其二：当捕丝的间距越小时，结构更牢固，捕获的食物更大。但是不同的织
网 的蜘蛛种类不一样，因此其间距合适才能不浪费材料的同时，捕获自身需要
的足够食物。
其三：捕丝的搭建直接与当时的空间有关系。只有空间允许条件下才能搭出
合理的捕丝结构。
其四：空气里有一定湿度，特别是凌晨时湿度最大，会影响到捕丝的黏性情
况。
其五：在没有外来食物的情况，蜘蛛网还要受到空气的流动的制约。
因此我们可以得出捕丝的 间距
d
1
主要由食物的状况、温度、空间、湿度、空气
流动，那么就需要找出 这五项因素与
d
1
的关系。本文利用层次分析的方法求解出
- 2 -

Team # 1739
各项因素的权重，从而可以求解出
d
1
。
3、根据科学生物现象， 知道蜘蛛是用一根丝即不会断的情况下完成整个织网
过程。因为要考虑其捕丝面积大，蜘蛛网为一根丝完 成，可以得出其构建的是一
个螺旋线的形状，且螺旋线的间距为
d
1
。再者根 据蜘蛛丝的本身的高韧性可知，
捕丝的伸长率为原来的4到5倍，因此可以利用能量守恒定律
1 1
mv
2
E
p
kx
2
求解出放射丝的间距< br>L
a
和放射丝的条数
a
。此处已完成整
22
个网的基 本框架。
4、根据空间的关系以及为了能把蜘蛛网能稳定的支撑起来，确定最后的框丝
和停泊线满足：
瑞伊尔说：“当昆虫碰撞入网，蛛网必须承受住碰撞力，而对称网的优势可能
在于它可以使这种力均匀地 分布在全网以减少某一处的受力，这样可以尽可能地
避免网被撕破。”所以无论是停泊丝还是边框丝都为 对称形态。蜘蛛网停泊丝与
边框丝的形态主要受三个方面的影响，一是为了保证蜘蛛网在受昆虫撞击力或 风
力的影响时结构变化较小，二是为了把蜘蛛网撑平，有效利用捕食面积，三是受
周围空间形状 与大小的影响。
具体情况，本文考虑蜘蛛网为一个平面状，其目的则是支撑它。那么本身其受
到重力的作用，因此通过受力分析，知道其边框至少为四的边。但是四边形不具
有很好的稳定性。因此我 们可以得出此种蜘蛛网的的边框为五边。

4、 模型假设
假设1：蜘蛛正常健康，吐的丝一定。
假设2：此期间蜘蛛处于一个相对安全状态，无天敌情况。
假设3：天气状况良好。
假设4：本文只考虑在空间织平面状的一类蜘蛛，不考虑再空穴里会织网的蜘蛛。
假设5：捕丝间距一定。

5、 模型建立

5.1 最大圆环捕丝模型

一张蜘蛛网的有效期为一天，每年因蜘蛛的生物钟的原因，蜘蛛会在 每天一定
的时间吐出新的蜘蛛网。因此我们可合理假设已知一只蜘蛛在一天的蜘蛛丝的总
共有L 。那么为了获得最大面积，我们可以建立圆面积最大模型。
设捕丝间距为
d
1
，总共有
n
圈的。
那么总材料长
L
L满足的方程式为：
L2


d
1
nd
1

n
2
，
- 3 -

Team # 1739
则最大面积：
nd
1
S
max


d
2

rdr

1
此模型可以求解出最大的捕环形捕丝面积与捕丝的间距的关系。

5．2 层次分析得出捕丝间距
d
1

通过问题的分析我们知道影响捕丝间距< br>d
1
主要有食物状况、温度、空间、湿
度、风。建立层次分析模型如下：

影响捕丝间距
d
1
的因素
A







风

x
2





B
1

湿度

x
2





B
2

空间

x
3





B
3

温度

x
4





B
4

食物

x
5


B
5


确定准则
B
i

i1,2,3,4,5

对目标
A
的权重向量：
两两比较
B
i

i1,2,3,4,5

对目标矩阵
A
的重要程度，构造准则
B< br>i
对目标
A
的成对比较判断矩阵：


1151315

2325

51

A

332113

1

5523
< br>7632


17

10.20.3330.20.143


16

510.6670.40.167
13



31.510.3330.333


12

52.5310.5


1

321


76

A
利用MATLAB软件求得 判断矩阵
A
的最大特征值

max
和相应的特征向量
WA
，求得
d
1
与
x
i
（
i1,2, 3,4,5
）的关系式子。具体结果见后文求解。

5.3 蜘蛛网的螺旋模型
- 4 -

Team # 1739

首先，蜘蛛由一根蜘蛛丝织成一个网状，其先固定好4个点，再利用空间条件
和 稳定条件，构建
a
条放射丝。最后由最外圈做等间距的向极心的螺旋运动。由
此可建立 模型为：

d
1

r
1
，其中
r< br>1
为捕丝面与中枢区的距离。
上式可知：
d
1
以求出，但是

和
r
1
未知。
1
其次，第一步可以知道n
与
d
1
的的关系：
n
4L
-1

d
1
2
。由此可以知：
（1）从
r
n
 nd
1
的最外圈上任意一点出发，向极心螺旋运动。
（2）螺旋线的角度的知识可以 知道

的范围为：

n-2

2




n-1

2

。
再者，一般 的湿性捕丝具有极高的弹性，根据生物学家的相关证明知道捕丝的
最大伸长率为原长的
4~5< br>倍，即：
4
x
5
，
L
a
此处x
为伸长量，
L
a
为放射丝的间距。
根据能量守恒，设外来 食物重
m
以
v
速度垂直冲击带蜘蛛网上，极限假设食
物的全部动能全 部转换为捕丝的势能得：
11
mv
2
E
p
kx
2

22
可以求解出La的大小范围。
当最外圈捕丝的周长
L
与
L
a
做商时可以求解出放射丝（半径）的条数：
aL
'
L
a
。易求出放射线平分圆的角度

2

。
a< br>一旦求出
L
a
之后就可求出
r
1
的大小。主要求解思 想是根据中枢区与捕食面距
离的大小主要材料节省的方面分析。
前文已求出放射线之间的夹角

与捕丝间距
d
1
，如果靠近中心点的位置处捕
丝间距仍为
d
1
，那么靠近中心位置处捕丝排列得将非常密，会比较浪费材料，所
以本文 适当减少中间位置所用的材料，即产生了中枢区与捕食面之间无捕丝的现
象，但是中枢区与捕食面距离满 足一定比值关系：
捕丝间距
d
1
，即网眼最短边为
d
1< br>，靠近中心点区域放射线切向间距较小、
- 5 -

Team # 1739
变化也较小，所以我们令放射线切向间距为
d
1
，则径向长为
r1

d
1
，也就是蜘蛛

在距中心
r
1
处停止织网。这样中枢区与捕食面之间放射线的切向距离近似于捕食
面中的网眼，中枢区与捕 食面之间的距离近似于捕食面中的
L
a
'
，符合之前捕丝
间距
d
1
的影响因素最后，由以上可以得出

d
1


形的
S
max
以及环形最外圈的周长
L
’
进行 修正：

n-1

2

S
’
max


n-2

2

1

2

L
''




n
n2< br>
2

d
1
。利用此结果再对前面环

1< br>
d
1

r
1

2
d


2

d
1
r
1

2
d
1
2
d



6、 模型求解域误差分析

6.1 层次分析的求解结果和层次单排序及一致性检验
确定准则
Bi

i1,2,3,4,5

对目标
A
的权重向量：
两两比较
B
i

i1,2,3,4,5
对目标
A
的重要程度，构造准则
B
i
对目标
A
的对
比较判断矩阵：

115131517

10.20.33 30.20.143


51232516510.6670.40.167< br>
10.3330.333


A

33211313



31.5

552 311252.5310.5


7

6321
< br>321


76

A
利用MATLAB软 件求得判断矩阵
A
的最大特征值

max
和相应的特征向量
W
A
为：
A
5.2171,W
A


0.0786,0.2075,0.2321,0.4782,0.8175



max
T
一致性检验：
A

max
 n
5.217250.2172
CI=0.0543
n1514

RI1.12

CR
0.0543
0.04850.1
1.12
- 6 -

Team # 1739
故通过检验。故准则
B
i

i1,2,3,4,5
对目标
A
的权重向量为：

W
A

0.0786,0.2075,0.2321,0.4782,0.8175



d
1
0.0786x
1
0.2075x
2
0.2321x
3
0.4782x
4
0.8175 x
5


6.2 螺旋线图形求解结果和图形比较
（1）捕丝圈数
T
1
n
4L
-1

d
1
2< br>，
其中
d
1
0.0786x
1
0.2075x
2
0.2321x
3
0.4782x
4
0.8175 x
5

（2）由前文模型建立中所列的公式：
x

4 5

L
a

1

1
2
mvk x
2

2


2

L
’
2

nd
1

d
1
*L
’


r
1

2

*L
a

可以求解出：
放射线的间距
L
a
：
mv
2
5k
放射线数
a
为：
L
a

mv
2
4k

8
nd
1
k
mv
中枢区与捕丝面积
r
1
： 2
a
10

nd
1
k
mv
2
r
1

d
1
，

mv
2
mv
2
2


其中可计算出

即
a
5nd
1
k4nd
1
k
因此可以求解出

d
1

r
1
的具体图形。
- 7 -

Team # 1739
根据以上结果，带入实际捕丝间距
d
1
0.5cm
时，
r
1
20
，运用matlab求解出
圆面积和螺旋面积



当
d
1
0.5cm
，
r
1
50










- 8 -

Team # 1739



当
d
1
0.8cm
，
r
1
20



当
d
10.8cm
，
r
1
50



6.2 圆最大面积与螺旋最大面积的数据为：
根据上一步骤，可见当
d
1
0.5cmr
1
20
时，与实际效果最想符合，因此带入
此数据 ，通过网络积分计数器求解出：

n

2
3
4 nd
1
S
max


d
2

rdr

1

n-1

2

S
’
max



n-2

2

0. 376991
1.00531
1.88496
- 9 -
1

d
1

r
1

2
d

2
1.76819
2.16997
2.55776

Team # 1739
5
6
7
8
9
10
3.01593
4.39823
6.03186
7.81681
10.0831
12.4407
2.95284
3.34733
3.74211
4.13689
4.83168
4.97646
其中
n
为螺旋圈数，
S
max
为最大圆环面积，
S
’
max为螺旋线的最大面积。
由此数据表可以看出，以圆为基础求得的捕丝面积变化远远大于螺旋线的面积
变化。分析原因：
如果蜘蛛网为同心圆形，会发生半径突变，但是因为蜘蛛是一根丝织完整
个网，不可能会是同心 圆的织法。如果为同心圆形，那么在确定半径
nd
的
情况下求得的圆面积会变化更大。 但是如果是蜘蛛是以螺旋线的方式运动
到中心，那么不会存在突变的问题，而且每增加一圈的螺距等间距 ，其螺
旋线的半径变化受螺旋线的的总角度来决定的。
综上所述，我们可以得出螺旋线更加能符合现实的蜘蛛网的捕丝的形状和结
构。
7、 模型评价与改进
7.1 模型评价
本文从节约材料的基础上，建立最大的捕丝面积的圆模 型。在此基础上，针
对蜘蛛本身的织网特性，建立切合实际的螺旋线的模型。
求解捕丝间距< br>d
1
时，根据生物学家的研究得出捕丝的间距主要由外界因素食
物状况、温度等 的影响。因此本文通过层次分析，找出各个因素影响捕丝间距的
权重。此步骤，起到了关键性的作用，求 解出了蜘蛛网网眼的大小，满足此种织
网结构的为了获得更多更大的捕食量的目的。
以上，只 是得出了在圆基础上改进的螺旋线的蜘蛛网的结构，并不能得出复
杂的蜘蛛网的结构，例如六边形的组合 的蜘蛛网、蜘蛛网边框的结构等。这也正
是模型改进所需要作的工作。
7.2 模型改进
改进（1）
针对蜘蛛网自身和蜘蛛的重力，利用受力平衡的方式，求解蜘蛛网的边框数。
假设此编织蜘蛛网的空间为墙角，此时蜘蛛先至出一个Y型的主轴。再在此基础
上加放射线。
对蜘蛛网受力分析：

10 - -

Team # 1739
mg

受力之后变成了四边形，但是实际情况中，四边形并不能稳定，因此可以得出
其边框在空间的影响下，一般只是应满足五边形。达到稳定的效果。
改进（2）
蜘 蛛网的捕丝结构的网眼不只是有四边形，还有六边形和三角形的捕丝面。但是
没有五边形的捕丝网眼。因 为只有矩形、三角形、六边形才能填满捕丝面积减小
空隙的存在。如下图所示：

三角形

六边形

五边形

11 - -

Team # 1739
间隙


因此还应该考虑捕丝的网眼为三角形和六边形的捕食状况。此种分 析，还应
结合蜘蛛本身需要什么样的食物来确定。

参考文献
[1] 司守奎 孙玺菁，《数学建模算法与应用》国防工业出版社，2011.
[2] 桌春晖 蒋平等，蛛 网结构性能及适应性，四川大学生命科学学院生物资源
生态环境教育部重点实验室，1000-7083 (2006) 04-0898-05.
[3] 潘自娟，异力学性能的结构机理及其模化，州大学纺织工程博士论文
[4] http:
[5]宋大祥，蜘蛛的进化，中国科学院动物研究所，1997.























12 - -

Team # 1739
附录：
1、
>> A=[1 15 13 15 17;5 1 23 25 16;3 32 1 13 13;5 52 3 1 12;7 6 3 2 1]
A =
1.0000 0.2000 0.3333 0.2000 0.1429
5.0000 1.0000 0.6667 0.4000 0.1667
3.0000 1.5000 1.0000 0.3333 0.3333
5.0000 2.5000 3.0000 1.0000 0.5000
7.0000 6.0000 3.0000 2.0000 1.0000
>> [V,D]=eig(A)
V =
0.0786
0.0029i
0.2075
0.0602i
0.2321
0.2483i
0.4782
0.3802i
0.8175
0.7475
D =
5.2172
0
0
0
0
0
0
0
0
-0.0170 - 0.4763i
>> h= 5.2172
h =
5.2172
>> ci=(h-5)4
ci =
0.0543
>> cr=ci1.12
cr =
0.0485




-0.0736 - 0.0893i -0.0736 + 0.0893i 0.0049 - 0.0029i 0.0049 +
-0.1965 + 0.2964i -0.1965 - 0.2964i 0.0696 + 0.0602i 0.0696 -
0.1413 + 0.0164i 0.1413 - 0.0164i -0.0929 - 0.2483i -0.0929 +
0.1519 - 0.1835i 0.1519 + 0.1835i -0.4668 + 0.3802i -0.4668 -
0.8850 0.8850 0.7475
0 0 0
-0.0916 + 0.9440i 0 0
0 -0.0916 - 0.9440i 0
0 0 -0.0170 + 0.4763i
0 0 0
-

13 -

Team # 1739

2、t=0:0.1:30*pi;
d=0.5;r=20;
r=d*t+r;
polar(t,r)

t=0:0.1:30*pi;
d=0.5;r=50;
r=d*t+r;
polar(t,r)


t=0:0.1:30*pi;
d=0.8;r=20;
r=d*t+r;
polar(t,r)


t=0:0.1:30*pi;
d=0.8;r=50;
r=d*t+r;
polar(t,r)

-

14 -