用lingo求解数学规划模型实例答案.ppt

附图1：A区的交通网络与平台设置的示意图 说明：（1）图中实线表示市区道路 （2）实圆点“·”表示交叉路口的节点， 没有实圆点的交叉线为道路立体相交； （3）星号“*”表示出入城区的路口节点； （4）圆圈“○”表示现有交巡警服务平台的设置点； （5）圆圈加星号“○* ”表示在出入城区的路口处设置了交巡警服务平台； 为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。 设立目标：各交警平台任务量尽量均衡 各平台到所管辖的路口的用时总量最小 约束：遇突发事件时，尽量能在3分钟内有警车到达事发地。 每个路口只被一个平台管辖。 决策变量： fij ：第i个平台管辖第j个路口 已知量：第j个路口的案发率 各路口的坐标 各条道路的起点与终点标号 由已知可求出各条道路的长度，从而得到各路口的邻接矩阵， 首先提取附件2中：全市交通路口节点数据、全市交通路口路线的数据， 将这两组数据保存在文件‘jiaojing.mat’文件中， 该文件包含两个矩阵，矩阵A是交通路口节点数据（582×5），矩阵B是交通路口线路（928×2）。 建立A区道路的赋权邻接矩阵D1 load jiaojing %jiaojing中矩阵A为路口节点矩阵，为1列编号，2列横坐标，3列纵坐标，4案发率，B为路口线路，1列路口起点标号，2列终点标号 xi=B(:,1);yi=B(:,2);m=length(xi); P=sparse(xi,yi,ones(1,m)); P0=P+P; %全市道路邻接矩阵，有道路连接为1，其余为0 PA=P0(1:92,1:92); %A区道路邻接矩阵,有道路连接为1,其余为0 PA(find(PA~=1))=inf; for i=1:92; PA(i,i)=0; end xx=A(:,2);yy=A(:,3); for i=1:582 %d为全市各路口直线距离矩阵 for j=1:582 d(i,j)=sqrt((xx(i)-xx(j))^2+(yy(i)-yy(j))^2); end end d1=d(1:92,1:92); D1=PA.*d1; %D1为A区道路赋权邻接矩阵 利用floyd法求出A区的最短路矩阵D D=D1; n=length(D);path=zeros(n); for k=1:n for i=1:n for j=1:n if D(i,j)D(i,k)+D(k,j) D(i,j)=D(i,k)+D(k,j); path(i,j)=k; end end end end D, 按尽量在3分钟之内有交警到达事发路口的要求，确定各交巡警平台可以管辖的路口： T=zeros(92,20); c=D(:,1:20)/10; %c为A区各路口到各平台的警车最小用时矩阵 for i=1:92 for j=1:20 if c(i,j)=3i20; T(i,j)=1; end end end T=sparse(T); %T为在3分钟内有警车可达的A区各路口矩阵，T(i,j)=1表示第j平台警车在3分钟内可达第i路口。 TT={}; %TT的i列为第为第i个平台3分钟内可达的路口 for i=1:20 ff=T(:,i); td=find(ff); TT(i)={td}; end kk=1:92;U=[]; for i=1:20 %该循环将第i个平台3分钟内可达的路口显示出来 AA=TT{i}’; [i,AA] kk=setdiff(kk,AA); end kk=setdiff(kk,1:20) %kk存3分钟内没有警车可达的路口 平台编号 可在3分钟内到达路口 1 42 43 44 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 2 40