分配管辖范围程序附件.doc

%matlab第一问1小问-分配管辖范围clcclearload shuju11zb=shuju11;jds=size(zb,1);a=zeros(jds,jds);a1=zeros(jds,jds);for i=1:jds for j=1:jds a(i,j)=sqrt(zb(i,2)-zb(j,2)2+(zb(i,3)-zb(j,3)2); endendload shuju12m=size(shuju12,1);lianj=shuju12;index=1:m;j=0;for i=1:m if lianj(i,2)92 j=j+1; index1(j)=i; endendlianj(index1,:)=;jiaot=lianj;lianj(:,2:-1:1);1:92;1:92; n=size(jiaot,1);for i=1:n a1(jiaot(i,1),jiaot(i,2)=1;end for i=1:jds for j=1:jds if a1(i,j)=0 a1(i,j)=inf; end endend a0=a1.*a;aa=a0/10;d,r=floyd(aa);%r?sjd=d(1:20,21:end); m n=size(sjd);for i=1:m% k=0; for j=1:n if sjd(i,j)=3 sjd(i,j)=inf;% k=k+1;% gld(i,k)=j+20; end endend for i=1:n % k=0; for j=1:m if sjd(j,i)min(sjd(:,i) sjd(j,i)=inf; % k=k+1; % gld(i,k)=j+20; end endendgld=zeros(m,n);for i=1:m k=0; for j=1:n if sjd(i,j)=inf k=k+1; gld(i,k)=j+20; end endendgld1=1:20,gld;%gld0=zeros(m,n);%for i=1:m k=0; for j=1:n if sjd(i,j)=inf k=k+1; gld0(i,k)=sjd(i,j); end endendgld2=1:20.*0,gld0;%matlabload shuju11x=shuju11(:,2);y=shuju11(:,3);plot(x(1:20),y(1:20),ro);hold onplot(x(7,8,9,30,31,33,34,35,36,45,46,47,48),y(7,8,9,30,31,33,34,35,36,45,46,47,48),r*);for k=1:92 text(x(k)+1,y(k)+1,sprintf(%d,k);endD=a1;for i=1:92 for j=1:92 if D(i,j)=1 line(x(i),x(j),y(i),y(j) % end endendlegend(,3,Location,best);xlabel();ylabel();title(A);! 第一问2小问lingo程序model: sets:si/1.13/;sj/1.20/;sij(si,sj):c,x;endsetsdata:c=22.23615 20.46392 18.35227 21.99738 17.62819 17.65878 14.91494 14.09251 13.01071 7.586585 3.791353 0 5.977002 11.95028 17.02961 14.54326 21.89211 24.24718 22.54653 26.9458 16.02847 14.12972 12.76723 15.00851 12.96963 13.00021 10.90122 9.433943 8.274202 12.77566 8.337298 11.95028 5.97328 0 13.29808 6.741662 14.90324 18.51448 16.96148 21.21311 9.286812 7.388063 6.025566 8.266853 6.227967 6.258552 4.159559 2.692282 1.53254 6.95667 11.39503 14.54326 12.71494 6.741662 6.556421 0 8.16158 11.77282 10.21982 14.47145 19.29344 17.39469 16.03219 18.27348 16.23459 16.26518 14.16619 12.69891 11.53917 9.510693 5.072332 8.685316 2.708314 3.264966 16.56305 10.00663 18.16821 21.77945 20.22645 24.47808 21.09621 19.19747 17.83497 20.07626 17.74952 17.78011 15.03627 14.21384 13.13205 7.707918 3.269557 6.882541 0.905539 5.067741 17.15094 11.8094 19.97098 23.58223 22.02923 26.28085 22.50175 20.603 19.24051 21.48179 19.15506 19.18565 16.44181 15.61938 14.53759 9.113456 4.675095 6.477002 0.5 6.47328 18.55648 13.21494 21.37652 24.98776 23.43476 27.68639 22.8932 21.12097 19.00932 22.65443 18.28524 18.31583 15.57199 14.74956 13.66776 8.243635 3.805274 3.59163 2.385372 8.358652 17.68666 15.10031 22.54916 24.90423 23.20358 27.60286 19.00116 17.22893 15.11728 16.22691 11.30687 11.33745 8.570218 10.22803 9.775722 14.19486 18.63323 21.78145 22.80832 18.04992 4.751842 11.30826 18.65712 21.01219 19.31153 23.01082 19.51581 17.74358 15.63192 15.53534 10.61529 10.64588 8.015457 10.49318 10.72441 15.14355 19.58191 22.73013 23.757 18.91668 5.700525 12.17502 19.52387 21.52683 19.82618 22.31925 12.08344 10.31121 8.19956 8.102976 3.182933 3.213518 0.583095 3.06082 3.492304 7.911446 12.34981 15.49803 16.5249 11.48432 4.401472 4.742655 12.09151 14.09447 12.39382 14.88689 5.880935 3.982186 6.09384 4.860976 9.421119 9.451704 7.352711 5.885434 4.725692 10.14982 14.58818 17.73641 16.12082 10.14754 9.749573 3.405877 4.755703 8.366946 7.639281 11.06557 11.85011 10.30954 8.197884 7.395869 2.475826 2.506411 1.290202 3.099467 4.19941 8.618553 13.05691 16.20514 17.23201 12.19142 5.108579 5.449761 12.79862 13.69926 11.99861 14.17978 4.885217 6.035068 4.393385 0.35 5.255075 5.337332 7.991722 8.677283 9.336668 14.7608 19.19916 22.34738 21.33179 15.35851 11.8101 8.616853 7.820525 6.734362 5.033709 6.44888 ;!ole(D:date1.xls,date1);enddatamin=max(si(I):(sum(sj(j):c(i,j)*x(i,j);for(sij:bin(x);for(sj(j):sum(si(i):x(i,j)=1);for(si(i):sum(sj(j):x(i,j)=1);End！第一问3小问lingo程序model: sets:si/1.92/:c;sj/1.92/:y;!交巡平台;sij(si,sj):d,x;endsetsdata:d=ole(d:date1.xls,date1);c=ole(d:date3.xls,date3);enddata!目标函数;min=max(sj(j):(sum(si(i):c(i)*x(i,j);!for(si(I):(sum(sj(j):d(i,j)*x(i,j)=3);for(sij:d*x=2;sum(sj(j)|(j#ge#21)#and#(j#le#92):y(j)=5;for(si(i):for(sj(j):x(i,j)=y(j);End%第二问1小问-评价合理性.matlab程序clcclearload shuju2zb=shuju2;sizezb=size(zb,1);jds=size(zb,1);a=zeros(jds,jds);for i=1:jds for j=1:jds a(i,j)=sqrt(zb(i,2)-zb(j,2)2+(zb(i,3)-zb(j,3)2); endendload shuju3lianj=shuju3;jiaot=lianj;lianj(:,2:-1:1);1:sizezb;1:sizezb;a1=zeros(jds,jds);n=size(jiaot,1);for i=1:n a1(jiaot(i,1),jiaot(i,2)=1;endfor i=1:jds for j=1:jds if a1(i,j)=0 a1(i,j)=inf; end endend a0=a1.*a;aa=a0/10;d,r=floyd2(aa);%r?lj(1:57)=r(1,582,1:57);ds1=d(1:20,1:92);ds2=d(93:100,93:165);ds3=d(166:182,166:319);ds4=d(320:328,320:371);ds5=d(372:386,372:474);ds6=d(475:485,475:582);tmin1=min(ds1);tmin2=min(ds2);tmin3=min(ds3);tmin4=min(ds4);tmin5=min(ds5);tmin6=min(ds6);dbl0= dbl(tmin1) dbl(tmin2) dbl(tmin3) dbl(tmin4) dbl(tmin5) dbl(tmin6);%pjcs0= pjcs(tmin1) pjcs(tmin2) pjcs(tmin3) pjcs(tmin4) pjcs(tmin5) pjcs(tmin6);%fal0=sum(fal(1:92)/size(fal(1:92),1) sum(fal(93:165)/size(fal(93:165),1) sum(fal(166:319)/size(fal(166:319),1) sum(fal(320:371)/size(fal(320:371),1) sum(fal(372:474)/size(fal(372:474),1) sum(fal(475:582)/size(fal(475:582),1);%zhb=dbl0;pjcs0;fal0;g=shang(zhb);pdx1=1 1/5 1/4;5 1 3;4 1/3 1;w1 cr1=cenquan(pdx1);w01=zqz(g,w1);zhb1=gyh(zhb);pj=zhb1*w01;%第二问2小问求解即仿真.matlab程序clcclearload shuju2zb=shuju2;sizezb=size(zb,1);jds=size(zb,1);a=zeros(jds,jds);for i=1:jds for j=1:jds a(i,j)=sqrt(zb(i,2)-zb(j,2)2+(zb(i,3)-zb(j,3)2); endendload shuju3lianj=shuju3;jiaot=lianj;lianj(:,2:-1:1);1:sizezb;1:sizezb;a1=zeros(jds,jds);n=size(jiaot,1);for i=1:n a1(jiaot(i,1),jiaot(i,2)=1;endfor i=1:jds for j=1:jds if a1(i,j)=0 a1(i,j)=inf; end endenda0=a1.*a;aa=a0/10;%32qwd=1:31 33:582;i=0;for j=qwd; tpt,tpp=graphshortestpath(a,32,j,Directed,0);%3217 if tpt3 i=i+1; tplx(i)=tpp; tpsj(i)=tpt; zd3(i)=j; endendzd3;tplx1,1:13;% a=sparse(aa);i=0;csk=151,153,177,202,203,264,317,325,328,332,362,387,418,483,541,572,578;path1=cell(1,size(csk,2);d1=zeros(1,size(csk,2);for j=csk i=i+1;d,path=graphshortestpath(a,32,j,Directed,0);%3217path1(i)=path;d1(i)=d;end%g