遥测遥感网部分解答.docx

1、遥测遥感网部分解答遥测遥感网（部分解答）大气污染所引起的地球气候异常，导致大面积严重森林大火的频频发生，给人民的生命财产造成巨大损失。因此，不少国家政府都在研究有效的森林防火措施。在容易出现高森林火险的重点地区放置高科技的监视装置，建立遥测遥感网，使人们能准确而及时地掌握险情的发展情况，为有效地防止火灾发生或在酿成严重灾害之前将其扑灭创造条件。科技的迅速发展使人们可以制造不太昂贵且具有收发报通讯功能的监视装置。放置在同一监视区域内的这种监视装置（以下简称为装置）构成一个Ad Hoc无线网络，即通常所说的遥测遥感网。如果监视区域的每一点都处于放置在该区域内某个装置的监视范围内，则称这些装置能覆盖

2、该监视区域。研究能确保有效（即按一定概率）覆盖且数量最少的装置系统的随机放置问题显然具有重要意义。第一个问题涉及能覆盖给定监视区域的装置数目及分配问题，具体如下：A1 设监视区域为边长b=100(长度单位)的正方形，每个装置的监视半径均为r=10(长度单位)。请参考蜂窝网格的特性讨论覆盖该区域所需装置的最少数量。根据文献（凡志刚，郭文生，桑楠，一种基于蜂窝网格的传感器节点部署算法，传感器与微系统，2008年，(27)4：15-17.）结论，传感器节点覆盖圆盘与其邻居节点覆盖圆盘的交点构成正六边形蜂窝时，有效覆盖面积最大。建立直角坐标系，设监视区域为边长b=100(长度单位)的正方形，每个装置的

3、监视半径（即正六边形的边长）均为r=10(长度单位) ，如上图，设表示第m层沿x轴正方向的第n个正六边形。当m为奇数时，的中心坐标为，n取不小于（这里）的最小正整数（这里取n=6）。当m为偶数时，的中心坐标为，如果（这里） ，则（这里取n=7），否则，。经分析，第m层所能覆盖的最远距离为，由确定层数m取不小于（这里）的最小正整数（这里取m=7）。再结合上图可知，覆盖该区域所需装置的最少数量为45（=4层偶数层6+3层奇数层7）。利用matlab可以画出图像，代码如下：%a1.close all;b=100;%正方形边长r=10;%圆的半径Ar=linspace(0,pi*2,200); % 圆

4、周角度cengshu=ceil(2*b/(3*r)+1/3);%层数zx=;zy=;for m=1:cengshu if mod(m,2)=1%表明为奇数层 for n=1:ceil(b/(sqrt(3)*r) zx=zx,sqrt(3)*r/2+sqrt(3)*r*(n-1); zy=zy,r/2+3*r*(m-1)/2; end else %表明为偶数层 if bceil(b/(sqrt(3)*r)*sqrt(3)*r-sqrt(3)*r/2 for n=1:ceil(b/(sqrt(3)*r)+1 zx=zx,0+sqrt(3)*r*(n-1); zy=zy,2*r+3*r*(m/2-1

5、); end else for n=1:ceil(b/(sqrt(3)*r) zx=zx,0+sqrt(3)*r*(n-1); zy=zy,2*r+3*r*(m/2-1); end end endendaxis(0,b,0,b);axis square;% plot(zx,zy,*)hold onfor j=1:length(zx) x=zx(j);y=zy(j);% fill(x+r*cos(Ar),y+r*sin(Ar),g);%填充圆 plot(x+r*cos(Ar),y+r*sin(Ar),-,LineWidth,2); % 画圆 plot(x,y,o,MarkerSize,10,Ma

6、rkerFaceColor,k); % 画圆心endhold offA2 在设计遥测遥感网时，首先需要知道对给定监视区域在一定的覆盖保证下应放置装置的最佳（越少越佳）数量，并且常假设装置在监视区域内是均匀地随机放置的。请在上述假设下建立数学模型，利用随机模拟实验回答：对于A1中给定的监视区域及监视半径，至少需要随机放置多少个装置，才能使得成功覆盖整个区域的概率在95%以上？并给出一个均匀随机放置装置的分布图。采用随机模拟，对装置个数依次取311396个，可以看到成功覆盖整个区域的概率都小于95%，最大83%，因此，取间距更大(每隔10)的个数400590进行试验，得结果表如下，表明大约至少需要

7、540个装置，才能使得成功覆盖整个区域的概率在95%以上。编号12345678910个数311312313314315316317318319320成功率41394037564841415544编号11121314151617181920个数321322323324325326327328329330成功率50535149465349465345编号21222324252627282930个数331332333334335336337338339340成功率57465147545752545357编号31323334353637383940个数3413423433443453463473483

8、49350成功率60585852534654556362编号41424344454647484950个数351352353354355356357358359360成功率61625964646153665961编号51525354555657585960个数361362363364365366367368369370成功率65745662656060686871编号61626364656667686970个数371372373374375376377378379380成功率72716669646468806666编号71727374757677787980个数38138238338438538

9、6387388389390成功率75726476727683687669编号81828384858687888990个数391392393394395396400410420430成功率72717683807975829383编号919293949596979899100个数440450460470480490500510520530成功率85868385889488919293编号101102103104105106107108109110个数540550560570580590成功率959496959795试验代码如下：clc;clear;Maxs=100;%模拟次数L=100; % 正方

10、形区域边长R=10; % 圆半径m,n=meshgrid(1:L);N0=400:10:600;% N=346; % 统计圆的数目for k=1:length(N0) N=N0(k); ss=0; for i=1:Maxs M=zeros(L); % 覆盖状态 scale=0; % 覆盖面积比例 x=L*rand(1,N); % 随机位置坐标 y=L*rand(1,N); % 随机位置坐标 for j=1:N D=sqrt(m-x(j).2+(n-y(j).2); % 计算坐标点到圆心的距离 m0,n0=find(D=1; % 检测是否满足覆盖比例 ss=ss+1; % 结束循环 end en

11、d fenbu(k)=ss;%存储满足条件的装置个数endresult=N0(1:length(fenbu);fenbu;A3 对一般矩形以及多边形的监视区域进一步探讨以上问题。由于监视旱情的遥测遥感网地处边远地区，它的每个（除极少数例外）装置都只能以电池为能源，电池用尽装置即报废。因此，如何节省电池能耗是设计此类网络运行方案的头等大事。常用的一个很有效的节能措施是：让大多数装置“休眠”只保留尽可能少的装置“值班”。对同时选出的这些值班装置的全体，必须要求它们整体具有与遥测遥感网的每个装置都能联系的功能，从而保证当任何休眠装置定时“苏醒”后若发现“险情”，都能及时向值班者之一传递险情信息。遥测

12、遥感网的若干装置组成的子集S称为一个支配集，如果该遥测遥感网中不属于S的任一装置必位于S中某个装置的通讯范围之内（即二者可互相交换信息）。不言而喻，上述同时“值班”的装置的集合必须要求是该遥测遥感网络的一个支配集。从实际应用的角度来说，这种支配集的优劣以其包含的装置个数来衡量（越少越优）。此外，如果把考虑的遥测遥感网视为一个无向图（每个装置是它的顶点，二顶点相邻接当且仅当二点间的距离小于公共的通讯半径R）。支配集按图论意义是连通者更为可取，因为通过仅在支配集内部传递信息的手段可以让它的每个装置共享任一装置所得到的信息，这样的支配集自然称为连通支配集。第二个问题涉及求元素尽可能少的支配集和连通支

13、配集的问题，具体如下：B1 设监视区域为边长b=100（长度单位）的正方形，每个装置的通讯半径均为R=10（长度单位）。已知在该监视区域内放置了120个装置，它们位置的横、纵坐标依次是：x=57,95,34,31,52,30,15,75,75,65,55,41,36,72,16,85,86,75,32,5,16,25,72,68,61,37,48,81,23,35,6,85,64,22,69,80,76,88,25,62,70,45,35,75,35,56,27,92,25,44,5,17,90,25,58,95,87,68,30,9,32,47,50,56,56,47,80,10,12,63

14、,39,81,43,17,80,45,92,78,89,51,40,65,76,30,26,28,25,29,40,4,74,41,39,95,72,79,78,10,8,15,45,70,90,84,20,40,55,5,73,22,17,50,55,87,72,55,7,85,35,10.y=58,74,12,68,67,4,75,52,30,28,63,61,20,24,10,49,90,90,20,92,35,66,4,33,35,78,46,31,90,66,33,9,37,13,43,83,13,94,95,45,70,42,9,41,91,30,92,90,58,52,80,33

15、,5,74,47,2,72,88,28,9,95,71,43,43,25,25,64,96,33,70,9,89,14,25,55,61,40,22,45,51,90,49,7,98,34,99,8,63,83,11,44,25,21,51,76,8,44,80,89,95,90,82,78,78,70,71,70,95,18,28,80,10,20,22,98,79,2,20,50,68.请建立数学模型找出一个较好的支配集；画出该120个装置的分配图，并在此图上标出所找到的支配集。解答：1、画图2、建模设表示第个装置属于支配集，否则，。建立如下优化模型。其中表示第个装置的坐标。目标就是使得支

16、配集的元素个数最少，约束条件是：第个点要么是支配点，此时，要么不是支配点，即，则必需与支配点（）相连。3、贪婪算法（参见文献：廖飞雄，马良，禁忌遗传算法求解最小支配集，计算机工程与应用，2007，（24）43，81-84.）如果把考虑的遥测遥感网视为一个无向图（每个装置是它的顶点，二顶点相邻接当且仅当二点间的距离小于公共的通讯半径R），将图的顶点记为，其对应的度数（即邻接顶点数）记为。贪婪算法可描述为：Step1 将V 中的顶点按顶点度数从大到小逆序排列成点集V，并全部顶点设置为未标号；Step2 取V中第一个顶点，若该点已经标号，并在V删除该点, 转至步骤3；否则，将该点标号为1，并将与之相

17、关联且未标号的顶点标号为0, 在V删除该点；Step3 若V为空, 转至步骤4；否则，转至步骤2；Step4 取标号为1 的顶点作为支配点，把这些点组成的点集为最小支配集。逆序启发式算法求出的结果是图的一个极小支配集。利用matlab编程求得的极小支配集共有36个装置，他们依次是第3，78， 10，35，47，64，99，4，18，50，52， 93，107，7，17，32，81，8，15，16，23，57， 59，1，26， 41， 42，49， 90，113， 31， 56， 67， 77， 94，112个装置，相应的坐标为X=34.0, 78.0, 65.0, 69.0, 27.0, 5

18、6.0, 8.0, 31.0, 75.0, 44.0, 17.0, 39.0, 55.0, 15.0, 86.0, 85.0, 40.0, 75.0, 16.0, 85.0, 72.0, 87.0, 30.0, 57.0, 37.0, 70.0, 45.0, 25.0, 4.0, 55.0, 6.0, 95.0, 80.0, 92.0, 95.0, 50.0；Y=12.0, 22.0, 28.0, 43.0, 92.0, 43.0, 89.0, 68.0, 90.0, 52.0, 33.0, 21.0, 70.0, 75.0, 90.0, 9.0, 90.0, 52.0, 10.0, 49.0

19、, 4.0, 72.0, 28.0, 58.0, 78.0, 70.0, 42.0, 58.0, 11.0, 20.0, 33.0, 2.0, 64.0, 40.0, 51.0, 10.0；他们的覆盖图形如下 他的代码如下：%b1x=57,95,34,31,52,30,15,75,75,65,55,41,36,72,16,85,86,75,32,5,16,25,72,68,61,37,48,81,23,35,6,85,64,22,69,80,76,88,25,62,70,45,35,75,35,56,27,92,25,44,5,17,90,25,58,95,87,68,30,9,32,47,5

20、0,56,56,47,80,10,12,63,39,81,43,17,80,45,92,78,89,51,40,65,76,30,26,28,25,29,40,4,74,41,39,95,72,79,78,10,8,15,45,70,90,84,20,40,55,5,73,22,17,50,55,87,72,55,7,85,35,10;y=58,74,12,68,67,4,75,52,30,28,63,61,20,24,10,49,90,90,20,92,35,66,4,33,35,78,46,31,90,66,33,9,37,13,43,83,13,94,95,45,70,42,9,41,9

21、1,30,92,90,58,52,80,33,5,74,47,2,72,88,28,9,95,71,43,43,25,25,64,96,33,70,9,89,14,25,55,61,40,22,45,51,90,49,7,98,34,99,8,63,83,11,44,25,21,51,76,8,44,80,89,95,90,82,78,78,70,71,70,95,18,28,80,10,20,22,98,79,2,20,50,68;% plot(x,y,*)%55555贪婪算法r=10;%通讯半径b=100;%正方形边长for i=1:length(x) A(i,:)=(x(i)-x).2+

22、(y(i)-y).2r2;%计算邻接矩阵 A(i,i)=0; dushu(i)=sum(A(i,:);%计算第i个节点的度数endyuanxu=1:length(x);%原始顺序,iv=sort(dushu,descend);%按度数从大到小排序s=;%支配集while isempty(iv) s=s,iv(1);%支配集 liandian=yuanxu(A(iv(1),:)=1);%与支配点相连的点 iv(1)=;%删除已经是支配集的点 for i=1:length(liandian) iv(iv=liandian(i)=;%删除与支配点相连的点 endendclose all;Ar=lin

23、space(0,pi*2,200); % 圆周角度axis(0,b,0,b);axis square;hold onfor j=1:length(s) zx=x(s(j);zy=y(s(j);%支配点的坐标% fill(x+r*cos(Ar),y+r*sin(Ar),g);%填充圆 plot(zx+r*cos(Ar),zy+r*sin(Ar),-,LineWidth,2); % 画圆 plot(zx,zy,o,MarkerSize,10,MarkerFaceColor,k); % 画圆心endplot(x,y,*)%画所有装置点hold off4、遗传算法（参见文献：廖飞雄，马良，禁忌遗传算法

24、求解最小支配集，计算机工程与应用，2007，（24）43，81-84.）一个支配集解对应为，S是一个由0,1组成的序列，图的最小支配集问题转化为一个无约束规划模型在式中n为顶点规模，C(S)表示S序列对应的非支配点，且不与支配点相关联的顶点集，|C(S)|表示C(S)中顶点的个数。采用遗传算法，模拟T=5000代，群体规模N=80，交叉概率取pm=0.05，变异概率pc=0.8，得到的最小支配集数目为：34他们依次是第 10 11 12 13 15 18 22 29 31 35 37 43 47 48 53 57 63 70 75 78 79 90 98 99 101 102 103 106

25、107 110 112 113 119 120个装置。他们的坐标是：x = 65 55 41 36 16 75 25 23 6 69 76 35 27 92 90 87 50 63 80 78 89 4 10 8 45 70 90 40 55 22 50 55 35 10y= 28 63 61 20 10 90 66 90 33 43 13 9 92 90 5 72 43 70 55 22 45 11 80 89 90 82 78 71 70 28 10 20 50 68图形如下Matlab代码为：clc;clear all;% 初始化参数%A表示是否相连的0-1矩阵L=120;%=lengt

26、h(x);%装置个数yuanxu=1:L;%原始顺序x=57,95,34,31,52,30,15,75,75,65,55,41,36,72,16,85,86,75,32,5,16,25,72,68,61,37,48,81,23,35,6,85,64,22,69,80,76,88,25,62,70,45,35,75,35,56,27,92,25,44,5,17,90,25,58,95,87,68,30,9,32,47,50,56,56,47,80,10,12,63,39,81,43,17,80,45,92,78,89,51,40,65,76,30,26,28,25,29,40,4,74,41,3

27、9,95,72,79,78,10,8,15,45,70,90,84,20,40,55,5,73,22,17,50,55,87,72,55,7,85,35,10;y=58,74,12,68,67,4,75,52,30,28,63,61,20,24,10,49,90,90,20,92,35,66,4,33,35,78,46,31,90,66,33,9,37,13,43,83,13,94,95,45,70,42,9,41,91,30,92,90,58,52,80,33,5,74,47,2,72,88,28,9,95,71,43,43,25,25,64,96,33,70,9,89,14,25,55,6

28、1,40,22,45,51,90,49,7,98,34,99,8,63,83,11,44,25,21,51,76,8,44,80,89,95,90,82,78,78,70,71,70,95,18,28,80,10,20,22,98,79,2,20,50,68;r=10;%通讯半径b=100;%正方形边长for i=1:length(x) A(i,:)=(x(i)-x).2+(y(i)-y).2r2;%计算邻接矩阵 A(i,i)=0;end%5T=5000;%仿真代数N=80;% 群体规模pm=0.05;pc=0.8;%交叉变异概率bval=round(rand(N,L);%初始种群bestv=-inf;%最优适应度初值%迭代开始for ii=1:T %计算适应度 for i=1:N S=bval(i,:);%取出一个0-1序列 CS=sum(all(A(yuanxu(1-S=