数学建模社区-数学中国

标题: 第三届全国研究生数学建模赛优秀论文 Ad Hoc网络问题的数学模型（2006年） [打印本页]

---

作者: xueyuan828    时间: 2009-9-3 19:07
标题: 第三届全国研究生数学建模赛优秀论文 Ad Hoc网络问题的数学模型（2006年）

全国第三届研究生数学建模竞赛

题 目
8 {2 D% x9 ?7 _- L5 O  M
" g/ B; L( y. P" J+ B4 m8 f9 O
Ad Hoc网络问题的数学模型

摘
  D0 h- r8 k" v. w' s要：

本文主要利用计算几何与图论的有关知识，分析和解决了Ad Hoc网络区域覆盖、信道分配、节点划分等问题，并应用仿真手段对各种条件下网络的连通性进行了研究。对前四个问题给出了求解的过程和结果，针对第五问，提出了解决问题的思想。
+ I# Q/ x5 q4 J) D, p8 d8 E问题一考虑到区域的对称性，将圆分别按照正三角形、正方形的排列方式去覆盖给定区域，覆盖方案为：当相交面积不小于圆面积的5％时，应采用正三角形排列，所需圆的个数是45，当相交面积不小于圆面积的18％时，应采用正方形排列，所需圆的个数是61。将信道分配问题转化为图的点着色问题，得出的结果为正三角形排列时分配3个信道，正方形排列时分配2个信道。通过对网络邻接矩阵D的研究，提出了网络连通的一个充要条件： 中的任意元素非零。最后通过仿真研究了网络的抗毁性。
0 N$ |# Y: a1 J  f问题二中首先证明了：在半径和一定时，用大圆比用小圆覆盖的面积大。因此在第一问的结果上对内、外边界上的圆进行局部调整，调整后半径和为4346。
问题三采用了改进的LBG算法对节点聚类，并定义了一个指标用于判断LBG算法执行的效果，选择较优的结果进行局部搜索，得到圆的数量为48，其半径和为4247。在判断连通性时，将每个圆看成一个节点建立邻接矩阵，然后根据问题一指出的充要条件进行判断，由此推出了一个更为直观的网络连通的充要条件。最后还提出了一种基于节点度数的划分方法。
  W# u8 y+ @0 g) I" ]2 c问题四首先通过仿真研究了节点随机运动下网络连通性。另外我们定义了网络的连通强度，深入讨论了节点随机运动条件下网络连通性的变化。
! |/ x  k) r# K* `  t问题五对通信过程中节点能量的变化规律进行了分析，建立了以节点退出网络所用的时间最大为目标的规划模型。
8 z; T# G! a6 u8 g) b: i

参赛队号
* d: D3 `, A8 }' H
! R; r; g: Z: F% N9 N
目
录一、问题重述... 2

; U* I  n& c, N1 F6 ]! _" W一、问题重述... 2

2 U; Y4 G; m* H( g# Z% Q9 w二、问题分析... 2

三、基本假设... 3

四、符号说明... 4

五、问题求解... 4

9 ?6 r7 E9 ?( h5．1 问题一的求解... 4

* f2 J6 L! [" C2 Q/ Q5．1．1 区域的覆盖问题... 4
6 i  f$ D+ n0 e* m4 L4 |) O
% z- m# |' |+ n! K- b3 \4 L7 X5．1．2 着色问题... 7
4 H% X4 q) w) H3 R- i
6 Z/ m5 }# M( t; d/ B5．1．3 抗毁性讨论... 8

5．2、问题二的求解... 9

5．3、问题三的求解... 11
- X& E7 N+ C. y& l2 M
5．3．1 基于胞腔划分的分簇方式... 11
! \: a2 N0 Q' q- O
5．3．2 抗毁性讨论... 14
) n- n, p& B$ m  N* ]# f- W" A+ M  w" L
  b" s0 Q8 \4 t5．3．3 分簇方法二：基于权值的划分方法... 15
$ ~; r4 A% V9 s2 i1 Y
5．4、问题四的求解... 16
1 C$ w  h5 j) K0 L% R/ G) v: Q
( m+ M: k/ \4 l& v, Q' S5．4．1 节点移动过程的仿真... 16

8 k! Q! _' y3 q- h; u9 W5．4．2 连通强度... 16

5．5、问题五的讨论... 18

六、模型评价... 21

3 j& x& r; u& }9 S$ N3 K# l七、参考文献... 21
' y! p$ @. W# N. C1 }4 j


# Q  f& D: m" |  f一、问题重述
7 U8 @% X- @" D3 {! _9 F
, q3 h( i3 b/ t* ]; t2 t) l
问题一：用若干个半径为100的圆完全覆盖一个边长1000的正方形区域，在相邻两个圆的公共面积分别不小于一个圆面积的5%和18%情况下，求圆的最小个数和信道分配方案；并讨论网络的抗毁性。
问题二：设正方形区域中有一椭圆形湖泊，节点仅能设置在地面上，研究使全部圆半径之和为最小的区域分划和信道分配方案。
问题三：将节点分簇，以完全覆盖某一簇内所有节点、且半径不大于100的圆作为一个一跳覆盖区。在满足相邻一跳覆盖区的公共面积不小于大圆面积的5%和网络连通等条件下，研究使全部一跳覆盖区半径之和为最小的一跳覆盖区划分和信道分配方案。找出区域连通的充分、必要条件。并讨论网络的抗毁性。
问题四：假设前10个节点作随机折线运动，其他节点不移动。节点到达正方形区域边界后只可能向区域内运动。讨论一定时间后Ad Hoc网络的连通性。
问题五：考虑Ad Hoc网络的节能的要求。请按照问题三中给出的办法（无湖的情况），找到比较节能的区域分划方式，使出现第一个退出网络的节点的时间尽量长。通过对该网络的运行状况进行分析，提出对组网方式的改进意见。
7 v5 i4 b7 u- ]问题六： 假设Ad Hoc网络中通信实行先到先服务，对Ad Hoc网络的通信质量进行定量评价。
  f+ }: k7 x$ J1 ~; R; y$ r二、问题分析
9 ^+ a1 n% j. G9 |0 q) d
问题一：
由于正方形区域具有对称性，若用一系列的圆覆盖该区域，那么这些圆的排列应该是规则的，圆心的排列自然也按照某种规则的方式排列，如正三角形、正方形、正六边形等。因此我们只需要根据问题要求对圆心的排列方式进行研究即可。
给每个圆分配一个信道，使相邻圆拥有不同的信道，本质上是一个着色问题，利用图论的有关知识可以解决。
对网络抗毁性的讨论实际上是对图的连通性的检验，可以把网络转化为图结构，根据图的邻接矩阵的性质来判断图的连通性。
+ ?1 f- s/ ?) _, Z7 m# N5 G问题二：
正方形区域中有一椭圆形湖泊，相当于区域中有一个内边界。考虑到问题一中的排列方式是一种较优的结果，可以在区域内沿用问题一的排列方式，只在边界（包括内、外边界）附近进行局部调整。
问题三：
与第一问相比，该问题有以下几点不同：
( O5 M. V+ c) c. O' T; U) B2 i/ N+ E1．要求用一些圆覆盖所有的点而不是整个区域；
& F1 i) w6 T- k% e+ |2．圆的半径可变；
3．所求目标是圆的半径之和最小而不是圆的个数最少。
这时圆的排列不再具有规则性，各圆的半径也不相同，可以先按照一定的原则将所有的节点分成一些小区域（簇），然后对各区域用尽量小的圆覆盖。
对网络抗毁性的讨论可以使用借鉴问题一判断连通性的的方法。
% \3 g5 D1 a  D: U0 [问题四：
. c' U. _* T4 J1 q) H9 w节点在做随机游动以后将处于一个新的位置，新位置可以用仿真的方法得到。同样用上面的方法来判断网络的连通性。为了进一步分析此网络的连通性能，可以另外建立一个指标——连通度，再用仿真的方法研究节点的随机移动对网络的连通度的影响。
# H. n* X* o; R2 ~9 Z" O问题五：
+ j$ R6 }8 l+ {考虑电池能耗的情形，需要我们对区域进行较优的划分，使得第一个因电池能量耗尽而退出网络的节点的时间尽量晚。各个节点的地位应该都是均等的，可以设想最佳的划分区域一定是半径相等或近似相等的圆，并且每个节点的最大通讯半径也应该相同，边界的区域和节点的半径可以不同。将电池能量转化为节点待机状态下的工作时间来描述，可以考虑每隔一段时间后电池的剩余能量，以电池能量消耗至0所用时间最大为目标，建立规划模型求出合适的圆半径，再在第三问的方法的基础上，提出合适的指标作为权值，结合算出的半径进行聚类分区。
) r8 Z# c  I+ p* ^* l三、基本假设

$ t) m5 l% I! v9 d7 ]1．假设相切的两个圆不是相邻的圆，即它们可以分配相同的信道；
2．三、四、五问中，两个圆的相交区域没有节点时，其相交面积大小没有要求；
3、假设网络具有好的分布式算法（信道接入、路由等）来协调节点间的通信。
四、符号说明

* @5 r/ m& }, j1 t4 |问题一：d--圆心之间的距离
6 b8 |# z) s  Y' @  ^8 bR--圆的半径
S--两个圆的公共面积   
1 N" c% j! Q$ c: o% k5 o--邻接矩阵
问题三： --节点集合
( E: O# M8 p3 g3 `/ y) _5 O--圆心集合
9 _, |* ^0 j, `0 @5 n! v' F--第i个胞腔
4 i6 E# o5 f! P5 {问题四： --连通强度
& W  v6 F- j) Z3 t  _8 k0 O问题五：W(ai) --节点ai的权值

E(ai) --节点ai在一段时间通讯后电池的剩余的能量
        r --节点在每段时间内平均转发的次数

En--节点在第n段时间后电池剩余的能量
' G6 W8 A3 Q- |  E; s! j+ F! q8 O; ^五、问题求解
* o5 K* N' p' b; d( v
5．1 问题一的求解5．1．1 区域的覆盖问题根据问题的分析，应该考虑规则的排列方式。首先考虑圆心按等边三角形排列，以三个圆为例，当相邻圆的圆心距离不同时，有如下三种情况（图1）：

a 3 }3 M2 h& y# L: i" K

b

c # C2 e' B  h6 f- P+ B. ?- X9 u 6 P4 N8 T8 y. l. Q/ ~1 Q: c

6 d6 ?& ?# Y0 W% m

图1 三角形排列的圆心的三种情况

如果不考虑相邻圆的公共面积的要求，显然要使得覆盖方形区域用的圆数最少，最优的选择是c（三个圆相交于同一点）。因为 a
0 j5 i7 @5 q, V6 f' P0 H5 G$ G结构中间有空隙，要覆盖方形区域，还需要在中间补一个圆，导致圆的个数增加；b
3 Z5 s8 Q' G; m  H$ ?2 [中三个圆有相交区域，显然比 c
# s& Y; W) x& K2 P: [: d+ O/ p结构的覆盖面积小。
但此题中两个相邻圆的公共面积是有一定要求的，不能简单的选用最后一种方式。因此首先计算情况 c
$ z7 H0 B! B) c( b3 A+ e* p的两个圆心距离以及公共面积。
2 s0 p. N; C$ w  h( k2 S设圆心距离为 d
4 a9 W3 \4 t/ {$ X0 t, k，圆的半径为 R
，夹角为a，公共部分的面积为：S    = p      -       =1811.722, 占圆面积的5.77%。这个值的意义是，若两个相邻圆的公共面积与圆面积之比小于5.77%时，符合情况a，但这种情况并不让人满意。因此我们将这种情况下的圆心都按照 c 来排列；当公共面积大于5.77%时，符合情况b。
题中要求公共面积不小于一个圆面积的5%，所以我们用c结构去覆盖方形。由于正三角形的边长确定，显然只需确定一个圆心就能确定整个正三角形网，因此问题转化为确定一个圆心的位置，使得方形区域内圆的个数最少。也可以认为是需要确定正方形的位置，使它覆盖到的圆的个数最少。需要注意的是，有些圆心并不在正方形内，但圆的一部分落在正方形内，该圆也需要计入总数，这类圆的圆心与正方形边界的距离满足一定的关系，如图2所示，实线表示正方形边界，当圆心落在实线与点划线之间时，该圆需要计入。
8 r8 D& Y) T$ ]* W7 J2 K! M

……

R

d/2 * [) x6 C4 V3 w 7 k. C4 n) M. v3 A! G9 B* [- h& _

图2 边界圆确定的示意图

综上，算法的步骤如下：
1 [% B( G2 r. S9 H1）以一个圆的圆心为原点，计算周围足够多个圆的圆心位置；
2）设正方形中心坐标为 ，其变化范围是：
3 U5 n. R, D4 r0 p0 Y

，

且 （中心坐标在一个圆内进行搜索）；
) p) `1 |  m( G/ Y' n8 ~
, m4 \4 _6 v# X5 I% U( `
3）统计 ， 范围内圆心的个数，即覆盖正方形所需的圆的个数。
- [3 h" d5 ?8 U5 T4 U/ z$ h  L. L4）对正方形中心坐标进行遍历搜索，并记录其所覆盖的圆的个数，其中最小值对应的 就是正方形的中心。
5）坐标平移：将原点移到正方形左下的顶点，计算每个圆心的坐标。
" x8 q5 u' y6 A* E# T4 h经计算，覆盖正方形区域最少需要45个圆，他们的位置关系如图3。
: {; @$ D$ }$ T3 j' j2 s

图3 公共面积比例不低于5％时的覆盖方法和信道划分

当公共面积比率不低于18%时，显然是情况 b。先计算圆心之间的距离 ，然后，通过搜索得到最少圆的个数是64，位置关系如图4。
3 W5 u, M* y2 u' q' O

  

图4 公共面积比率不低于18％时的覆盖方法和信道划分

同样，也可以考虑圆心按四边形结构排列，根据圆心之间距离的不同，也有三种情况（图5），显然， f
是一种边界状态。它的公共面积为：S    = 2(p      - )，占圆面积的18.17％。
# n5 r) \* o0 d% s

f

d   U% s0 L6 H7 }- M7 R+ H 6 x- A0 q8 \, F) i

e ( u0 I- F! m3 e

图5 四边形排列的圆心的三种情况

当公共面积仅要求大于圆的面积的5％时，不会采用四边形排列，因为这个比例与18.17％相差较远。当要求公共部分面积大于18％时，与18.17％比较接近，因此用 f
排列方式来覆盖圆，具体的实现步骤与三角形排列相同，经计算，至少也需要64个圆。
! F' t: _" y/ u' \5 p  s6 v然而我们发现，此时边界圆处于正方形外面的部分较大，并没有充分利用，将f 旋转45度，得到一种新的排列方式（图6），经计算，这种方式只需要61个圆（图7）。

图6 旋转后的排列方式
2 K2 K' \& L1 I" c! E
图 7公共面积比例为18％时四边形排列圆的覆盖方法
5．1．2 着色问题构造图G(V,E)，V是覆盖方形区域的圆的圆心构成的顶点集，eij Î E当且仅当圆i和j相交。给每个圆分配一个信道，使得有公共部分的圆拥有不同的信道，实际上是对图G(V,E)的顶点进行染色使相邻的两点颜色不同。根据图论知识已知以下定理：
定理1. 一个简单图是1可着色的当且仅当它是空图。
# h# S- I" r# ^0 E# [5 |1 m, Z定理2. 一个简单图是2可着色的当且仅当它是偶图。
当圆心按照三角形排列时，图G(V,E)是正三角形网的一部分，顶点的最大次数为6，它既不是空图也不是偶图，而六个相邻的正三角形构成的图是3色图，G(V,E)正是这样的图的扩展，所以G(V,E)是3可着色的，即只需分配3个不同信道即可。着色方案标注在图3和图4中。
5 o6 }; v, d" I& z& Q当圆心按照四边形排列时，图G(V,E)是正四边形网的一部分时，G(V,E)是一个偶图，所以它是2可着色的，从而只需分配2个信道。着色方案见图7。
! B$ ?! G$ Z7 Z& v" ?; ?5．1．3 抗毁性讨论此问的关键问题是判断图的连通性，即判断是否任意两点之间都有通路。设方形区域内有N个节点，以这N个节点为顶点构造图G(V,E)，eij Î E当且仅当节点i和j在同一圆内。设其邻接矩阵为 ，其中

令 ，表示矩阵 的 次方幂， ，我们可以得到图G(V,E)的连通性的充要条件：
定理3. 图 G(V,E)不连通当且仅当矩阵 中存在元素 ， .
" W+ Q. n( E' n证明：我们先来证明这样一个结论： 中 ，当且仅当i经过m跳不能够到达j  .（ ， ）
时：
9 ]: M. n. ?' B0 O, t  R  g& g

，

由 的定义可知， 当且仅当 i
经过一跳不可以到达 j
；
时：
8 J& q( p5 N, B+ }

，

表示i经过两跳可以到达j
（中间经过1个节点k），如果 ，则对任意的k
， 都为0，即i
& f' U5 K" ]' X8 E8 T2 i经过两跳不能到达 j
6 b+ o6 a0 {: t# E3 P8 U  H；
% j5 B2 D8 S, C$ Z* q! `假设 时成立：
( D: Z' Z  N6 A- r) I当且仅当i经过 n
6 [/ M3 c, Q, V# i; b3 A跳不能到达j
；
时：

，

则 当且仅当对所有的k=1, 2 , , N有 。若 ,由归纳假设知i经过 n
跳不可以到达k点，所以i不可能经过 n
跳到达k再一跳到j. 若 ，则 ，即i能经过 n
% A8 K$ o8 ~) N; t) ^跳到达k，但k不能一跳到j. 所以i经过 n +1跳不可能到达j.
) @& a- E) }( ^因为 ， 即 ，当且仅当 ，其意义是i
8 s3 [/ {4 [# o! L经过1跳不能到达 j ，i
6 ^+ V" L$ \+ f经过2跳不能到达 j ，……，i
经过 跳不能到达 j，而图G(V,E)中有路的任意两点的路长不超过 ，所以i
( B. v! Z) h2 v2 G! m/ Z和 j 之间不存在通路，即图G(V,E)不是连通的，定理得证。
- I" w/ m% f6 j5 P  i/ L+ ^根据这一原则，随机去掉一定比例的节点后，考察网络的连通性。公共面积比例f分别为5%和18%时，求得在2%、5%、10%、15%时网络不连通的概率如图8所示。
8 d2 n5 N% m- l

图8 公共面积比例分别为5%和18%时网络的不连通

概率与节点毁坏率的关系图

5．2、问题二的求解首先给出一个定理：
定理4：在不考虑边界的前提下，如果存在两组大小不同的按照情况 c
8 a8 o8 S# G+ L; @' A排列的圆，当两组圆的半径和相等时，大圆的覆盖面积大于小圆的覆盖面积。
证明：如图9所示的排列方式，一个圆覆盖的有效面积等于它的内接正六边形的面积。
6 T5 X2 Y% ~: m

图9 圆等效面积示意图

设两组圆的半径分别为 、 ，则其内接正六边形的边长分别为 、 ，因此每个六边形的面积分别为：

，

设每组圆的个数分别为 、 ，每组圆覆盖的有效面积分别为 、 ，则

， ，

当 时，

，

因此当 时， ，即大圆的覆盖面积大于小圆的覆盖面积。
9 |2 j+ Z% D. Y1 ^1 G' r: a$ m这说明，为了使半径和最小，在区域内部，应该用尽量大的圆去覆盖。
5 k1 v! b6 W4 d. ^. \0 P因此，当正方形区域内有一个湖时，只需要对边界（包括湖的边界）的圆进行调整，而内部的排列方式不必改变。调整的基本原则是保证公共面积的比例要求下，调整圆心和半径，使圆的半径尽量小。
先考虑内边界，即湖的边界。首先对圆心的位置在可行解的范围内进行微调，使得椭圆边界经过尽可能多的圆的交点，这样做的目的是尽量减少需要调整的圆的个数。调整后，由椭圆与圆的位置关系可以看出，受到影响较大的有三个圆，其中一个圆的绝大部分都在椭圆内，如果去掉该圆，会有一个较小的区域不能被覆盖，如果要覆盖，需要在该处加入一个半径为75的圆，但该小区域的面积很小，在工程应用中，我们认为这样的小区域不被覆盖是允许的，因此可以不再加入圆。对另外两个圆，通过在一定的范围内，对圆心和半径进行搜索，检查是否可以减小半径，结果显示，这两个圆的半径并不能减小。
* V8 t! n9 @) k* i再考虑外边界，即正方形边界。通过观察圆的排列可以发现，左右边界上都有可能进行调整的圆，对这些边界圆的圆心、半径进行搜索发现：左边有三个圆的半径可以减小，最小值是82。
) j' j3 b( b- P( C- L4 s5 r9 M图10显示了边界调整后的结果，全部圆的半径之和为4346。

图10 有湖时的区域分划调整前后的对比

$ m% p0 `5 N& u+ A2 L
, b4 r2 H% j7 L! M+ @3 B# j% e调整后的结果对信道分配并没有影响，该图仍是一个3可着色的，即需分配3个不同信道，不同的是不必给被去掉的圆分配信道，其他圆的信道与问题一中的信道相同。
5．3、问题三的求解5．3．1 基于胞腔划分的分簇方式矢量量化中胞腔的相关理论：
设样本矢量集合 ， 是样本个数， 是一个矢量， 是再生矢量（码本）集合， 是码本的个数， 是与 同维的矢量。
# _% g( I8 }2 h1 l! q3 N对样本矢量集合进行分割：
- F- ?: n1 D) J3 P， ，
$ Z& I, i  b1 c& o就称为胞腔， 表示将 划分到第i
, u* P7 L( A, s' R0 |个胞腔所引起的失真，可以用均方误差等来度量。因此矢量量化就是找到一组码本和胞腔的划分，使平均失真最小。
针对本问题，将节点看成样本矢量，圆心的集合作为码本集合，每一个胞腔中的节点就组成一簇，并用距离来衡量失真，则这样的划分可以使得平均距离最小，相应的圆的半径也较小。
在具体实现时，使用现有的LBG算法，它的算法流程如图11。
! q* v+ X1 x( {7 Y6 L+ `$ Q. c

开始

; T9 \/ m/ }6 ]1 m5 j7 M

划分胞腔，计算平均失真D0 9 d+ Y  w8 p% L4 q8 P ) O8 o* \3 o4 P- Y' S8 k3 ?, f

4 x8 i/ ~' E0 Q$ ^+ h8 r

重新划分胞腔，计算平均失真Dm 2 _/ V; F* D/ W: _3 V4 P / {9 n0 b3 `$ ^

0 k; v3 l3 X( [' P8 c3 ? 0 d0 I/ ~' K" U+ ~. I, Q3 ]  s

结束

k=k+1

图11 LBG算法实现胞腔划分的流程图

在计算时，具体问题的处理如下：
2 J. F  R; H9 L1．码本个数的选择：通过观察已给的数据点看出他们是比较均匀的分布在整个正方形内，因此覆盖这些点与覆盖这个区域所用的圆的个数应该相差不多，因此可以预先 设定码本个数的范围是40到50之间；
  L, u, c$ \  i5 x$ \$ a! _2．初始码本的生成：通常有随机产生法和分裂法，随机产生法是在正方形区域内随机产生 m 个点作为初始码本；分裂法是先计算整个区域的概率中心 ，将 乘以一个常数值 ，将 和 作为码本，重新划分胞腔，计算最佳码本，完成一次分裂，下一次分裂之前要判断每个胞腔中节点和码本距离的最大值，如果大于100，继续分裂，否则不分裂，依次类推，直到完成分裂。后者不需要事先设定码本的个数，但计算复杂，所以这里采用随机产生法；
3．划分胞腔的依据：

  V9 f# @) R2 G5 s; i. C  v

4．平均失真的计算：
) _; Z" ?4 E) q1 W+ x

5．第 n
次码本的产生：第 n
次循环的第 j
个码本是由第 次循环得到的第j
3 _  ^# c* D$ n9 K个胞腔内的所有点的概率中心求得；
3 ]& o7 y& |. g6．结束条件的判断：当前后两次计算的平均失真的相对变化 小于门限 时，结束运算。
经过以上的计算得到了m
: d+ I) f( e; o" P个码本，即圆心的位置，以该胞腔内的点与码本的最大距离为半径画圆，如果最大距离大于90，取半径为90，取90的目的是为后面的调整留有余地，这样就得到了一种覆盖方式。但由于限制了半径的最大值，可能导致部分节点不能被覆盖，另外，这种方式仅仅考虑了半径的因素，对区域是否连通，相邻圆的公共面积是否达到5％并没有加以限制，因此结果可能并不是可行解。为了综合考虑这几个条件的限制因素，定义一个评价指标，用来衡量结果的优劣。该指标要综合反映点的覆盖情况、公共面积的比例要求这三个因素，令
5 N* Q* G( @- a+ m

其中 表示没有被覆盖的点的个数， 表示公共面积不满足5％的要求的个数（题中“满足有转发任务的相邻一跳覆盖区的公共面积不小于较大一跳覆盖区面积的5%”，我们认为，如果两个相邻一跳覆盖区的公共部分没有交点，对该公共面积也不作要求）。
在不同的码本个数下，反复执行LBG算法n次，并分别用指标h进行评价。在所有结果中，若存在 的方案，选择其中半径和最小的一个的作为最终的结果；如果不存在 的方案，则选择指标 最小的方案，然后对不符合要求的圆的圆心进行局部范围内的遍历搜索，并辅助扩大圆半径的方式得到一个可行解。
# Q9 D# o# I3 u. ]$ u8 I' a) b* R' }5 z区域划分好之后，构造 ，其中 是将每个圆看成一个点所形成的顶点集，eij Î E当且仅当圆i和j相交，参考文献[1]中图的顶点染色算法容易求得 的染色方案，即相应的信道分配方案。
+ }8 ], [/ b9 I" J: k8 K; G图12为最终确定的区域划分的结果以及信道分配方案，此时半径和为4247，每个圆的圆心及半径见附录。
, J- s: y- Q/ K/ U- W

      

图12 基于节点的区域划分结果及信道分配方案

5．3．2 抗毁性讨论在问题一中，我们建立了所有节点之间的连接矩阵，通过矩阵运算判断网络的连通性，但在此问中，共有926个节点，如果也采用类似的方法，计算量是很大的，因此必须对节点进行预处理。
2 e/ q+ f/ g) J/ \建立一个新图 ，其中 是将每个圆看成一个点所形成的顶点集，eij Î E当且仅当圆i和j相交且公共部分有节点，显然，图 的连通性与网络的连通性是等价的。这是因为如果 中两点 、 之间有通路，则意味着圆i和j中所有点之间都有通路，因此当 中任意两点都有通路（连通）时，任意两圆中的所有节点都有通路，即网络是连通的。
设邻接矩阵 ，其中
  [; O5 J* }1 j8 i. \

不仅能反映两个圆是否有公共节点，也反映了公共节点的个数。同5.1.3中的讨论，当且仅当矩阵 中， ， ，都有 时，网络是连通的。
1 B7 e) P$ r1 x; k为研究网络的抗毁性，随机去掉一定比例的节点，重新构造矩阵D
" l% g5 d! v1 A$ n6 O7 U，检验网络的连通性，得到网络性能如图13。

图13 在2%、5%、10%、15%毁坏概率下网络的不连通的概率

5．3．3 分簇方法二：基于权值的划分方法我们对节点用圆进行分区覆盖，使得圆半径和最小，这样的圆的圆心节点周围应该较密，所以我们提出一种能反映这种性质的指标作为节点优先作为圆心的权值，我们以距某节点ui 的距离不大于某一定值（如100）的节点数目作为该节点的权值，称为它的度d(ui)。区域划分的步骤如下：
Step1: 输入节点集 及每个节点的坐标；
Step2: 计算每个节点的度d(ui)，并排序 ， 是某个节点 ，记录与 距离小于100（可以调节）的节点集 ， .
1 O" u* g6 P$ d* o9 Q1 nStep3: 比较 与 ， ， 找到最小的正整数l，使得 ，把 作为初步聚类，把 作为圆心 .
Step4: 设与圆心 对应的半径为 ，以半径和最小为目标，利用类似方法一中的算法搜索合适的半径。
! _$ W- T; F/ o* M2 R# l利用此种方法聚类出来的圆心数固定，只须搜索半径即可。按方法一中的方法容易分析网络的连通性和抗毁性。
4 R' ]; J9 i' `7 v& [9 [0 a/ c5．4、问题四的求解5．4．1 节点移动过程的仿真为了研究节点移动后网络的连通性，首先要研究节点的运动，根据题中要求，前10个用户只作折线运动，每30个单位时间可能改变一次运动的方向和速度，运动的方向角、速度是分别服从在[0，2p] 、[0，2]上均匀分布的随机变量，其他节点不移动。节点到达正方形区域边界后只可能向区域内运动。
据此对节点的运动进行仿真，随机选择运动方向和速度，图14是10个点的运动轨迹，图15是其中一个节点运动轨迹的放大图。
0 m" Q  `* J* G: N  R. {( r5 w/ i5．4．2 连通强度前10个节点可以移动，相当于改变了10个节点的位置，在第三问中已经解决了网络连通性的判断问题，因此只需要重新构造连接矩阵，就能判断节点移动后网络是否连通。

图14 10个节点的运动轨迹

图15
其中一个节点的运动轨迹的放大
由于总节点的个数是926，10个节点仅占总结点数的1.08％，如果是去掉这10个节点，网络连通的概率与上一问中随机去掉1.08％的节点后网络连通的概率应该相差不多，但这里并不是去除点，而是将这些点重新放置，因此连通的概率要大于上一问中随机去掉1.08％的节点后网络连通的概率，可见，这时是网络的性能是比较好的，大多情况都能保持连通。为了进一步分析连通的性能，在网络连通的前提下，给出了连通强度的概念。
文献[3]中，连通强度的定义是：
! I) {* ~! w/ Y9 \: t( ?/ U$ ~

其中，E 是网络中存在的链路数， 网络全连通时的链路数，即任意两个节点之间都有一条边。
  ^  m& e, U6 m" l) _9 f& x在本题的网络模型中，我们注意到，处于圆内部和相邻圆公共部分的节点，对网络连通性的贡献是不同，处于公共部分的节点越多，网络的连通性越好，因此定义：
" C% d3 d% f: C5 b! ]0 t

其中， 是节点总数， 表示处于公共部分的节点总数，显然，连通强度越大，网络的抗毁性越好。与文献中的定义方式相比，它的优点是：突出了公共部分的点对连通性的重要性；统计点的个数比边的个数算法更简单、更易实现。
根据连接矩阵的定义， 的值就是处于圆i
、j的公共区域的节点的个数，那么矩阵中所有元素值之和的一半就是处于公共区域的节点数。
* N9 G5 X# \" h6 d

图16 200次仿真时网络的连通强度图

对10个节点的运动多次仿真，在网络连通时，计算连通强度，得到200组仿真结果，将连通强度绘于图16（节点未移动时，连通强度为0.1031）。
从图中可以看出，当节点随机移动时，网络的连通强度在0.1031附近上下浮动，最大可达0.1059，最小为0.102。可见，节点的移动虽未引起网络是否连通的改变，但也影响了网络的抗毁性，例如，如果原节点 处于两圆的公共部分，当它移动到某个圆的内部时，虽然网络仍然是连通的，但其抗毁性是有所降低的，由此也可以看出，节点原来所处的位置对连通强度的变化是有影响的，如果10个节点大部分位于公共区域，那么移动后网络的抗毁性的总趋势是下降的。
5．5、问题五的讨论 本问要求在第三问的基础上，不考虑节点的移动，确定给定节点的较节能的区域划分，即考虑电池在发射、接收和备用等状态下的能量损耗，进行合适的区域划分使得第一个因电池能量耗尽退出网络的节点的时间尽量晚。从附件1给出的数据看这些节点分布较均匀，而目标是使第一个退出的节点时间晚，所以各个节点的地位都是均等的，可以设想最佳的划分区域一定是半径相同的圆，并且每个节点电池的最大通讯半径也应该相同，边界的区域和节点可以不同。根据本题要求，按照第三问的方法对这些节点进行分簇必须提出一个合理的划分指标，既要使节点的寿命更长又不能划分的区域太多，我们提出以节点度数和电池的能量剩余的加权和作为节点的权值。定义节点ai的权值为：

; A* C  U' L$ b7 J% D+ Q7 TW(ai) = w1 d(ai) + w2 E(ai).
(1)

其中w1
' V3 ]6 g" U6 _( f5 x+ w2 =1，d(ai)表示与ai的距离不超过某一定值（如100）的节点个数，反映的是点ai周围节点的稠密程度；E(ai) 表示与ai在一段时间通讯后电池的剩余的能量，它能够反映该节点在网络中的寿命长短。w1和w2分别是节点度数和电池的能量剩余的权重，可根据实际需要设定。本问以考虑节能为主，所以w2要比w1大些。
设电池的最大通讯半径为R , 划分区域的圆的半径为r , 考虑到同一区域的两点可以直接通讯不需中转，所以R ³ 2r ,由于发射功率近似地与最大传输距离的三次方成正比，电池在覆盖半径为100发送状态下的工作总时间是400个时间单位，根据电池总能量相等容易计算电池在覆盖半径为R时发送状态下的工作总时间，并转化为待机备用状态下的工作总时间

                         .
(2)

我们将电池的能量都转化为待机工作时间来描述，最佳的分区一定是使得在一段通讯时间内每个节点进行的发射、接收和转发的任务量分别几乎相同。因为每个节点在1200个时间单位内平均产生25次呼出，通讯是随机的，所以同样应该有25次接受，我们以每48个时间单位作为一个研究的时间段，平均每一段每个节点都有一次呼出、一次接受和一些转发任务，用En表示某一节点在第n段时间后电池能量的剩余，
) V$ a& Q7 W; o( }+ f. o

En=En-1 - (4´11 + 4´10 + 40 + 44r ),
/ {7 r3 u/ s9 Q! b: |1£ n £ 25.
(3)

其中 r 表示节点在每段时间内平均转发的次数。再利用(2)式可得到，

  X1 w" b4 D6 L, ^/ MEn= - (124+ 44r ) n.

3 `: u; J4 h3 l+ ~8 _% V(4)

令En=0，则
  w# z8 y7 q& v0 w

3 ?  X4 m7 O$ o4 [5 an=                               （5）

要节点退出的时间晚即要使得En=0时的n值尽量大。
7 r) j; E  W# c% n3 Z' x

S * v& ]$ F( T. c6 N" Q) [( N# \  P9 J

A1 ! E  ]; P2 W5 _' D: F" t

A2 9 m$ {- |" c) c' I- b$ {, z9 l 8 J. _! q5 N# g' ~4 g6 ^" ]! _+ w

B1 0 t+ Q" l: K/ _   C7 J4 ^7 Q* J% L; E

B2Ï

区域I 0 `; E' ~' T! ~' f# [ 2 T5 p% B# w$ h0 g4 d* O' y

区域II

C - h& T3 _3 O; `9 M+ Y+ P

图17 两个相邻区域的通信示意图

5 s6 r. W' l7 |9 n下面我们来讨论 r 的计算，令
, B- k- [* @4 ?8 c

) a7 ?4 I# n; ?9 u/ `d=
2 u1 d- c) ^% v1 Rmax{ d(ai) , ai 跑遍所有节点}
+ ]7 H( \5 V$ {, e
(6)

用d近似表示区域内节点的个数，如图所示区域I和区域II的半径为r ，均包含d个节点，设相交区的节点数为h，将每个圆相交区外的节点粗略地分成两类A={节点到S的距离接近 }，B={节点到S的距离接近 }，其中l是圆心到S的距离， 他们的个数分别约为 和 ；记 表示在一段时间内 中节点到 中节点通讯的次数，类似地定义 .则

= , =

= , =

则S的转发次数：

r = + + +

记角∠SA2C=a ,则

l = , .

于是可以建立规划模型：
7 m! V  b  h4 S4 H3 Z7 u

我们可以限定R为2r，求解此规划，设使得目标最小的r和a 值为 ， 。取(6)式中d= d(ai)，类似计算ri
= + + + ，进而计算
4 B6 e8 w5 Y  C9 V9 i# L: U% [E(ai)=E0 - (124+ 44ri ),
W(ai) = w1 d(ai) + w2 E(ai)

最后以权值W(ai)为指标，按照第三问中提出的方法对节点进行分簇，确定圆心。区域圆的半径为 ，圆心与覆盖区两端连线成的角为 ，覆盖区面积占原面积的百分比为：

。

网络的改进建议
; O% [% ?0 l5 |% g9 G$ J该网络中个节点的地位近乎相等，呼叫点到目的点的通讯存在多条路径，当网络规模扩大时路由维护的能量开销指数增加，且消耗有限的带宽，我们可以在区域圆相交的区域内选取一节点作为簇头，使得簇头与同簇内节点通讯用一个频率，簇头之间的通讯用另一频率，且对簇头可进一步分簇，确定更高级的簇头，构成分级结构，低级网络的频率可低些，通讯范围在区域内，高级网络频率可高些，通讯范围大些，这样将节点划分级别，不再一律平等，这样整体的能耗将会减少。
( k6 Q/ Q" b# V# p六、模型评价
% ?5 P& \. n) i" c7 V
' C2 R. t0 }2 v! p本文分析解决了Ad Hoc网络的几个有关问题，主要的优点有：
, s; T! |, T# {) e* [" d1．从简单的对称排列入手，较好的解决了区域的覆盖问题。
2．讨论网络的抗毁性时，建立连接矩阵，并给出了网络连通的充要条件，通过矩阵运算判断连通性，避免了求最小生成树等复杂的算法。
3．基于节点的个数定义了连通强度，更进一步的分析了网络的连通性能。
存在的问题是：
1．对问题二、三的求解时，给出了一种我们认为较好的方法，但对该方法的分析、以及解的可行性没有进行深入的研究；
2．由于时间原因，第五问仅给出了思路，并未实现，第六问未考虑。

---

作者: xueyuan828    时间: 2009-9-3 19:11
图标有点问题

---

作者: lingongchun    时间: 2009-9-3 19:12
2009/9/3 签到贴！大家加油哦~~ 吼吼~~~

---

作者: 我想发飙    时间: 2009-9-3 20:44
不愧是研究生的题目啊

---

作者: Pepsi09    时间: 2009-11-25 15:14
LZ很好心，帮我们省钱，但是有些符号看不见啊！

---

作者: 为你奋斗    时间: 2009-11-30 17:31
lz不知道发压缩版本吗？？？！！！！！！！！！！！！！

---

作者: csucsm    时间: 2009-12-1 13:39
不愧是研究生的题目啊
5 X! S5 N. Y9 r1 w- q+ s六楼，看看加分不

---

作者: wanghuizan2004    时间: 2009-12-2 22:27
呵呵 我当时也做这道题 获得了全国二等奖

---

作者: fanarsenal    时间: 2010-9-2 22:29
这个还是压缩一下大家下载吧

---

作者: hbdkfk2    时间: 2012-8-30 13:58
谢谢分享！！！！！！！！！！！！！！

---

作者: tbfuyun    时间: 2013-5-19 11:17

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5