交通巡警调度

que1

file='cmm.xls';
AbiaoX=xlsread(file,1,'B2:B93');%A区路口节点标号横坐标
AbiaoY=xlsread(file,1,'C2:C93');%A区路口节点标号纵坐标
Anfalv=xlsread(file,1,'E2:E93');%案发率表格
% plot(AbiaoX,AbioaY,'g*');
S=xlsread(file,2,'A2:A144');%路线起点标号
T=xlsread(file,2,'B2:B144');%路线终点标号

%去除A区以外的点
ex=find(T>92);
S(ex)=[];T(ex)=[];

%画相邻点的连线
figure(1);
x1=[];y1=[];x2=[];y2=[];%起终点坐标
bS=[];bT=[];bW=[];%构造稀疏矩阵需要的起终点标号及权值
for bh=1:140;
%     x1=[x1 AbiaoX(S(bh))];y1=[y1 AbiaoY(S(bh))];
%     x2=[x2 AbiaoX(T(bh))];y2=[y2 AbiaoY(T(bh))];
     x1=AbiaoX(S(bh)); y1=AbiaoY(S(bh));
     x2=AbiaoX(T(bh)); y2=AbiaoY(T(bh));
     bS=[bS S(bh)];bT=[bT T(bh)];
     dis=sqrt((x1-x2)^2+(y1-y2)^2);bW=[bW dis];
     h=plot([x1,x2],[y1,y2],'b-d','LineWidth',2,'MarkerSize',3);
%      if bh<21%设置交巡警点标号
%          set(h,'Marker','o','MarkerSize',8,'MarkerEdgeColor','r');
%      end
%      if ismember(bh,churuA)
%          set(h,'Marker','*',MarkerSize',6,'MarkerEdgeColor','r');
%      end
     hold on;
end

%标号交巡警位置
c='A';tc=[];
for i=1:20
    tc=num2str(i);
    tc=[c tc];
    h=plot([AbiaoX(i)],[AbiaoY(i)]);
    set(h,'Marker','o','MarkerSize',8,'MarkerEdgeColor','r');
    text(AbiaoX(i),AbiaoY(i),tc,'Color','m','horizontalalignment','left');
%     alpha=linspace(0,2*pi);
%     rx=30*cos(alpha)+AbiaoX(i);
%     ry=30*sin(alpha)+AbiaoY(i);
%     plot(rx,ry,'g.','LineWidth',0.2);
end

for i=21:92
        ttc=num2str(i);
        text(AbiaoX(i),AbiaoY(i),ttc);
end

churuA=[12 14 16 21 22 23 24 28 29 30 38 48 62];%标号出入A区的路口
for i=12:62
    if ismember(i,churuA)
        ttc=num2str(i);
        h=plot([AbiaoX(i)],[AbiaoY(i)]);
        set(h,'Marker','*','MarkerSize',7,'MarkerEdgeColor','r');
    end
end

%求任意A区两点间最短路
DG=sparse(bS,bT,bW,92,92);
UG=tril(DG+DG');
[dist]=graphallshortestpaths(UG,'Directed',false);
save afile.txt -ascii dist;
fid=fopen('afile.txt','wt'); %写入的文件
a=dist;
[m,n]=size(a);
 for i=1:1:m
    for j=1:1:n
       if j==n
         fprintf(fid,'%g
',a(i,j));
       else
        fprintf(fid,'%g	',a(i,j));
       end
    end
 end
fclose(fid);


pdist=zeros(20,13);
for i=1:20
    for j=1:13
        pdist(i,j)=dist(i,churuA(j));
    end
end
save bfile.txt -ascii pdist;
fid=fopen('bfile.txt','wt'); %写入的文件
a=pdist;
[m,n]=size(a);
 for i=1:1:m
    for j=1:1:n
       if j==n
         fprintf(fid,'%g
',a(i,j));
       else
        fprintf(fid,'%g	',a(i,j));
       end
    end
 end
fclose(fid);

cango=[];
for i=21:92
    for j=1:20
        if dist(i,j)<30
            cango=[cango i];
            break;
        end
    end
end
ungo=setdiff(21:92,cango);%规定时间不能到达的点


d32=[];
for i=1:92
   if dist(32,i)<=30
       d32=[d32 i];
   end
end


g1=[];g2=[];g3=[];g4=[];g5=[];g6=[];g7=[];g8=[];g9=[];g10=[];
g11=[];g12=[];g13=[];g14=[];g15=[];g16=[];g17=[];g18=[];g19=[];g20=[];

for j=21:90
        if dist(1,j)<30 g1=[g1 j]; end;
end;
for j=21:90
        if dist(2,j)<30 g2=[g2 j]; end;
end;
for j=21:90
        if dist(3,j)<30 g3=[g3 j]; end;
end;
for j=21:90
        if dist(4,j)<30 g4=[g4 j]; end;
end;
for j=21:90
        if dist(5,j)<30 g5=[g5 j]; end;
end;
for j=21:90
        if dist(6,j)<30 g6=[g6 j]; end;
end;
for j=21:90
        if dist(7,j)<30 g7=[g7 j]; end;
end
for j=21:90
        if dist(8,j)<30 g8=[g8 j]; end;
end
for j=21:90
        if dist(9,j)<30 g9=[g9 j]; end;
end
for j=21:90
        if dist(10,j)<30 g10=[g10 j]; end;
end
for j=21:90
        if dist(11,j)<30 g11=[g11 j]; end;
end
for j=21:90
        if dist(12,j)<30 g12=[g12 j]; end;
end
for j=21:90
        if dist(13,j)<30 g13=[g13 j]; end;
end
for j=21:90
        if dist(14,j)<30 g14=[g14 j]; end;
end
for j=21:90
        if dist(15,j)<30 g15=[g15 j]; end;
end
for j=21:90
        if dist(16,j)<30 g16=[g16 j]; end;
end
for j=21:90
        if dist(17,j)<30 g17=[g17 j]; end;
end
for j=21:90
        if dist(18,j)<30 g18=[g18 j]; end;
end
for j=21:90
        if dist(19,j)<30 g19=[g19 j]; end;
end
for j=21:90
        if dist(20,j)<30 g20=[g20 j]; end;
end
            
% l1=Anfalv(1);l2=Anfalv(2);l3=Anfalv(3);l4=Anfalv(4);l5=Anfalv(5);
% l6=Anfalv(6);l7=Anfalv(7);l8=Anfalv(8);l9=Anfalv(9);l10=Anfalv(10);
% l11=Anfalv(11);l12=Anfalv(12);l13=Anfalv(13);l14=Anfalv(14);
% l15=Anfalv(15);l16=Anfalv(16);l17=Anfalv(17);l18=Anfalv(18);
% l19=Anfalv(19);l20=Anfalv(20);
%  
% for i=1:size(g1,2)
%     l1=l1+Anfalv(g1(i));
% end
% for i=1:size(g2,2)
%     l2=l2+Anfalv(g2(i));
% end
% for i=1:size(g3,2)
%     l3=l3+Anfalv(g3(i));
% end
% for i=1:size(g4,2)
%     l4=l4+Anfalv(g4(i));
% end
% for i=1:size(g5,2)
%     l5=l5+Anfalv(g5(i));
% end
% for i=1:size(g6,2)
%     l6=l6+Anfalv(g6(i));
% end
% for i=1:size(g7,2)
%     l7=l7+Anfalv(g7(i));
% end
% for i=1:size(g8,2)
%     l8=l8+Anfalv(g8(i));
% end
% for i=1:size(g9,2)
%     l9=l9+Anfalv(g9(i));
% end
% for i=1:size(g10,2)
%     l10=l10+Anfalv(g10(i));
% end
% for i=1:size(g11,2)
%     l11=l11+Anfalv(g11(i));
% end
% for i=1:size(g12,2)
%     l12=l12+Anfalv(g12(i));
% end
% for i=1:size(g13,2)
%     l13=l13+Anfalv(g13(i));
% end
% for i=1:size(g14,2)
%     l14=l14+Anfalv(g14(i));
% end
% for i=1:size(g15,2)
%     l15=l15+Anfalv(g15(i));
% end
% for i=1:size(g16,2)
%     l16=l16+Anfalv(g16(i));
% end
% for i=1:size(g17,2)
%     l17=l17+Anfalv(g17(i));
% end
% for i=1:size(g18,2)
%     l18=l18+Anfalv(g18(i));
% end
% for i=1:size(g19,2)
%     l19=l19+Anfalv(g19(i));
% end
% for i=1:size(g20,2)
%     l20=l20+Anfalv(g20(i));
% end

%  [l1;l2;l3;l4;l5;l6;l7;l8;l9;l10;l11;l12;l13;l14;l15;l16;l17;l18;l19;l20]
   
Lv=[6.20 7.7  7.6   6.6  7.3 6.3 5.1 7 7.5 2.8 3.4 4 7.1 3.9 6.4 5.2 7.5 8.6 6.1 7.1];

　　graph2_diaodu

file='cmm.xls';
AbiaoX=xlsread(file,1,'B2:B93');%A区路口节点标号横坐标
AbiaoY=xlsread(file,1,'C2:C93');%A区路口节点标号纵坐标
% plot(AbiaoX,AbioaY,'g*');
S=xlsread(file,2,'A2:A144');%路线起点标号
T=xlsread(file,2,'B2:B144');%路线终点标号

%去除A区以外的点
ex=find(T>92);
S(ex)=[];T(ex)=[];

%画相邻点的连线
x1=[];y1=[];x2=[];y2=[];%起终点坐标
bS=[];bT=[];bW=[];%构造稀疏矩阵需要的起终点标号及权值
for bh=1:140;
%     x1=[x1 AbiaoX(S(bh))];y1=[y1 AbiaoY(S(bh))];
%     x2=[x2 AbiaoX(T(bh))];y2=[y2 AbiaoY(T(bh))];
     x1=AbiaoX(S(bh)); y1=AbiaoY(S(bh));
     x2=AbiaoX(T(bh)); y2=AbiaoY(T(bh));
     bS=[bS S(bh)];bT=[bT T(bh)];
     dis=sqrt((x1-x2)^2+(y1-y2)^2);bW=[bW dis];
     h=plot([x1,x2],[y1,y2],'b-d','LineWidth',2,'MarkerSize',3);
%      if bh<21%设置交巡警点标号
%          set(h,'Marker','o','MarkerSize',8,'MarkerEdgeColor','r');
%      end
%      if ismember(bh,churuA)
%          set(h,'Marker','*',MarkerSize',6,'MarkerEdgeColor','r');
%      end
     hold on;
end

%标号交巡警位置
c='A';tc=[];
for i=1:20
    tc=num2str(i);
    tc=[c tc];
    h=plot([AbiaoX(i)],[AbiaoY(i)]);
    set(h,'Marker','o','MarkerSize',4,'MarkerEdgeColor','r');
    text(AbiaoX(i),AbiaoY(i),tc,'Color','m','horizontalalignment','left');
%     alpha=linspace(0,2*pi);
%     rx=30*cos(alpha)+AbiaoX(i);
%     ry=30*sin(alpha)+AbiaoY(i);
%     plot(rx,ry,'g.','LineWidth',0.2);
end

churuA=[12 14 16 21 22 23 24 28 29 30 28 48 62];%标号出入A区的路口

for i=21:92
        ttc=num2str(i);
        if ismember(i,churuA)
            text(AbiaoX(i),AbiaoY(i),ttc,'Color',[1 0 0]);
        else
             text(AbiaoX(i),AbiaoY(i),ttc);
        end
end


for i=12:62
    if ismember(i,churuA)
        ttc=num2str(i);
        h=plot([AbiaoX(i)],[AbiaoY(i)]);
        set(h,'Marker','*','MarkerSize',10,'MarkerEdgeColor',[1 0 0]);
    end
end

%求任意A区两点间最短路
DG=sparse(bS,bT,bW,92,92);
UG=tril(DG+DG');
[dist]=graphallshortestpaths(UG,'Directed',false);
save afile.txt -ascii dist;
fid=fopen('afile.txt','wt'); %写入的文件
a=dist;
[m,n]=size(a);
 for i=1:1:m
    for j=1:1:n
       if j==n
         fprintf(fid,'%g
',a(i,j));
       else
        fprintf(fid,'%g	',a(i,j));
       end
    end
 end
fclose(fid);

cango=[];
for i=21:92
    for j=1:20
        if dist(i,j)<30
            cango=[cango i];
            break;
        end
    end
end
ungo=setdiff(21:92,cango);%规定时间不能到达的点



g1=[];g2=[];g3=[];g4=[];g5=[];g6=[];g7=[];g8=[];g9=[];g10=[];
g11=[];g12=[];g13=[];g14=[];g15=[];g16=[];g17=[];g18=[];g19=[];g20=[];

for j=21:90
        if dist(1,j)<30 g1=[g1 j]; end;
end;
for j=21:90
        if dist(2,j)<30 g2=[g2 j]; end;
end;
for j=21:90
        if dist(3,j)<30 g3=[g3 j]; end;
end;
for j=21:90
        if dist(4,j)<30 g4=[g4 j]; end;
end;
for j=21:90
        if dist(5,j)<30 g5=[g5 j]; end;
end;
for j=21:90
        if dist(6,j)<30 g6=[g6 j]; end;
end;
for j=21:90
        if dist(7,j)<30 g7=[g7 j]; end;
end
for j=21:90
        if dist(8,j)<30 g8=[g8 j]; end;
end
for j=21:90
        if dist(9,j)<30 g9=[g9 j]; end;
end
for j=21:90
        if dist(10,j)<30 g10=[g10 j]; end;
end
for j=21:90
        if dist(11,j)<30 g11=[g11 j]; end;
end
for j=21:90
        if dist(12,j)<30 g12=[g12 j]; end;
end
for j=21:90
        if dist(13,j)<30 g13=[g13 j]; end;
end
for j=21:90
        if dist(14,j)<30 g14=[g14 j]; end;
end
for j=21:90
        if dist(15,j)<30 g15=[g15 j]; end;
end
for j=21:90
        if dist(16,j)<30 g16=[g16 j]; end;
end
for j=21:90
        if dist(17,j)<30 g17=[g17 j]; end;
end
for j=21:90
        if dist(18,j)<30 g18=[g18 j]; end;
end
for j=21:90
        if dist(19,j)<30 g19=[g19 j]; end;
end
for j=21:90
        if dist(20,j)<30 g20=[g20 j]; end;
end
            

　graph3

clc;
file='cmm.xls';
allX=xlsread(file,1,'B2:B583');%路口节点标号横坐标
allY=xlsread(file,1,'C2:C583');%路口节点标号纵坐标
allS=xlsread(file,2,'A2:B929');
allT=xlsread(file,2,'B2:B929');
Jiaoxun=xlsread(file,3,'B2:B81');
% allJiao=xlsread(file,3,'A2:A81');

figure(2);
long=size(allS,1);
for i=1:long
    px1=allX(allS(i));px2=allX(allT(i));
    py1=allY(allS(i));py2=allY(allT(i));
    h=plot([px1,px2],[py1,py2],'b-d','LineWidth',2);
    hold on;
end

long=size(allX,1);
for i=1:long
    px=allX(i);py=allY(i);
    c=num2str(i);
    if ismember(i,Jiaoxun)
        plot(px,py,'m.','MarkerSize',19);
    end
    text(px,py,c,'Color','r');
end


x1=[];y1=[];x2=[];y2=[];%起终点坐标
allbS=[];allbT=[];allbW=[];%构造稀疏矩阵需要的起终点标号及权值
for bh=1:928;
     x1=allX(allS(bh)); y1=allY(allS(bh));
     x2=allX(allT(bh)); y2=allY(allT(bh));
     allbS=[allbS allS(bh)];allbT=[allbT allT(bh)];
     alldis=sqrt((x1-x2)^2+(y1-y2)^2);allbW=[allbW alldis];
end
%求任意A区两点间最短路
DG=sparse(allbS,allbT,allbW,582,582);
UG=tril(DG+DG');
[distt]=graphallshortestpaths(UG,'Directed',false);
dis=distt;

　　LINGO程序

model:
!求解围堵 13 条要道的方案 ;
sets:
AA/1..20/;
cross/1..13/;
links(AA,cross): dis, x;
Endsets
!数据的定义部分 ;
data: 
dis = @FILE(F:file.txt);
enddata
!目标函数 ;
min=@max(links(i,j):x(i,j)*dis(i,j));
!需求约束 ;
@for(cross(j):@sum(AA(i): x(i,j))=1);
@for(AA(i):@sum(cross(j): x(i,j))<=1);
!整数约束 ;
@for(links(i,j):@bin(x(i,j)));