图与网络分析例题讲解.doc

图与网络分析例题讲解例1求救信号的采集问题紧急呼救电话发挥着极其重要的作用，现在的问题是往往在呼救时当事者大多处于紧张或身体状况不佳的状态，难以清晰表达自己所处位置，给救援工作带来极大的困难，对于有线电话来说，定位相对容易，而对于移动设备由于其可移动性，则确定位置相对比较困难。一种可行的办法是依赖通信基站，按照移动设备接收附近几个基站信号强弱进行定位。区域内的某个点接收到各基站的信号强度组成一个向量，该向量唯一标志区域内的一个点。采用这种方法定位就需要采集区域内各点的信号强度，派遣一辆装载信号采集设备和GPS的车辆，从研究所出发，依次到达各主要地点采集信号，最后回到研究所提交数据。考察某大城市的一个特定区域，示意图共5个节点。主要信号采集点在图中已标出（即图中的节点），如何选择一条最短路线，使得信号采集车辆能够顺利地采集信号并返回研究所。图的邻接矩阵为：。解 该问题实际上就是一个TSP（旅行商问题），要求寻找遍历图中所有节点，并返回起点的最短路。TSP属于组合优化的范畴，可以采用组合优化的方法求解TSP。设表示两个城市之间的距离，决策变量是或1（0表示不连接，1表示连接），由组成的邻接矩阵是图的哈密顿圈等价于中每个节点都只有一个入度和一个出度，且去掉任何一个节点将不是圈。此时求解TSP就等价于求解下面0-1规划问题： （1）例2 装备的合理配置问题设有套不同型号的装备要配备给个部队，由于各个部队的基础设施、训练特点等条件的差异，不同的装备在不同的部队所产生的效能是不同的，具体的数据如表1所示。试问如何分配这批装备，保证每个部队都有一套设备，并且使总的效能最大？表1 装备在不同部队效能表 装备部队ABCDEFGHI10.140.170.230.550.470.260.190.170.1220.370.400.490.090.050.530.420.390.1230.590.620.670.220.170.060.030.020.0840.110.120.160.060.030.190.140.120.4650.120.140.190.240.190.460.370.350.1060.100.120.150.060.030.390.330.300.2170.110.140.180.470.390.060.030.020.2580.630.650.730.070.040.220.170.140.0990.290.300.360.050.030.050.020.010.44解 由题意可以知道，这个问题是属于一标准指派问题，即属于组合优化的范畴，在这里我们来建立组合优化模型，并且相应的方法进行求解。将各部队关于各种装备的效能（表1）数据用矩阵表示，即用表示分配装备给部队产生的效能。用表示决策矩阵，为一个0-1矩阵，即表示将装备分配给部队；表示不将装备分配给部队，则此时可以建立如下的优化规划模型： （2）例3 网络的数据传输问题分组交换技术在计算机网络发挥着重要作用，从源节点到目的节点传送文件不再需要固定的一条“虚路径”，而是将文件分割为几个分组，再通过不同的路径传送到目的节点，目的节点在根据分组信息进行重组、还原文件。分组交换技术具有文件传输时不需要始终占用一条线路，不怕单条线路掉线，多路传提高传输速率等优点。现在考虑如图2所示的网络，图中连接两个节点间的数字表示两交换机得可用宽带，此时从节点1到节点9的最大传输宽带是多少？解 将此问题视为一个求网络最大流问题，将分组的传输方式用以下矩阵来刻画：其中表示从节点到节点的实际传输宽带。记容量矩阵为由此可以建立线性模型如下： （3）例4 出租车的最短行驶路线问题某市的出租车公司为了更好地为乘客服务，向乘客承诺：“出租车走最短的行驶路线，方便快捷。”乘客上车后只要告知司机目的地，出租车上电脑就可以计算出到达目的地最短的行驶路线。解 首先将地图视为一个赋权图。function d,DD=dijkstra(D,s)%Dijkstra最短路算法Matlab程序用于求从起始点s到其它各点的最短路%D为赋权邻接矩阵%d为s到其它各点最短路径的长度%DD记载了最短路径生成树m,n=size(D);d=inf.*ones(1,m);d(1,s)=0;dd=zeros(1,m);dd(1,s)=1;y=s;DD=zeros(m,m);DD(y,y)=1;counter=1;while length(find(dd=1)<m     for i=1:m         if dd(i)=0             d(i)=min(d(i),d(y)+D(y,i);         end     end     ddd=inf;     for i=1:m         if dd(i)=0&&d(i)<ddd             ddd=d(i);         end     end     yy=find(d=ddd);     counter=counter+1;     DD(y,yy(1,1)=counter;     DD(yy(1,1),y)=counter;     y=yy(1,1);     dd(1,y)=1;end