2011年数学建模B题交巡警服务平台的设置与调度代码

ShapeX=[]ShapeY=[]N=length(ShapeX);fori=1:Nforj=1:NDistance(i,j)=sqrt((ShapeX(i)-ShapeX(j))^2+(ShapeY(i)-ShapeY(j))^2);endendDistanceA=zeros(N);Max_Value=zeros(N);fork=1:N[max_line,column]=max(Distance(k,:));A(k,column)=max_line;endMax_Value(k,column)=max(max(A))[I,J]=find(Max_Value)point_start=[ShapeX(I)ShapeY(I)]point_end=[ShapeX(J)ShapeY(J)]fori=1:140fork=1:20;text(x,y,'str(k)');enddata1=[413359403343383.5351381377.5339376335383317362334.5353.5333342282325247301219316225270280292290335337328415335432371418374444394251277234271225265212290227300256301250.5306243328246337314367315351326355327350328342.5336339336334331335371330371333388.5330.5411327.5419344411343394346342342342348325372315374342372345382348.5380.5351377348369370363371353354374363382.5357387351382369388335395381381391375392366395361398362401359405360410355408350415351418347422354418.5356405.5364.5405368409370417364420370424372438368438.5373434376438385440392447392448381444.5383441385440.5381.5445380444360];data2=[1751782443453654394635495506597327478984793510341122112612251421157153116141638174017421781188118831979208621222213231324132425251126272610271228292815293030730483132313432333334338349354536353637361636393773839384139404024117419242434324372443454646846554748476475486149504953505151525159525653525354545554635535657575857605745859606261606246285636464656476656666676676674467686869687569706971691702704371727174727373747318741748075767677777877197879798080188182828829083848485852086878688878887928889889189208984899090919192];x=[413403383.5381339335317334.5333282247219225280290337415432418444251234225212227256250.5243246314315326327328336336331371371388.5411419411394342342325315342345348.5351348370371354363357351369335381391392395398401405410408415418422418.5405.5405409417420424438438.5434438440447448444.5441440.5445444];y=[359343351377.5376383362353.5342325301316270292335328335371374394277271265290300301306328337367351355350342.5339334335330333330.5327.5344343346342348372374372382380.5377369363353374382.5387382388395381375366361362359360355350351347354356364.5368370364370372368373376385392392381383385381.5380360];c=[413403383.5381339335317334.5333282247219225280290337415432418444];d=[359343351377.5376383362353.5342325301316270292335328335371374394];x1=data1(data2(i,1),1);y1=data1(data2(i,1),2);x2=data1(data2(i,2),1);y2=data1(data2(i,2),2);a=[x1,x2];b=[y1,y2];plot(x,y,'*',c,d,'r*',a,b);holdonendShapeX=[]ShapeY=[]N=length(ShapeX);fori=1:Nforj=1:NDistance(i,j)=sqrt((ShapeX(i)-ShapeX(j))^2+(ShapeY(i)-ShapeY(j))^2);endendDistanceA=zeros(N);Max_Value=zeros(N);fork=1:N[max_line,column]=max(Distance(k,:));A(k,column)=max_line;endMax_Value(k,column)=max(max(A))[I,J]=find(Max_Value)point_start=[ShapeX(I)ShapeY(I)]point_end=[ShapeX(J)ShapeY(J)]Matlab源代码为functionFloyd(w,router_direction,MAX)%x为此图的距离矩阵%router_direction为路由类型：0为前向路由；非0为回溯路由%MAX是数据输入时的∞的实际值len=length(w);flag=zeros(1,len);%根据路由类型初始化路由表R=zeros(len,len);fori=1:lenifrouter_direction==0%前向路由R(:,i)=ones(len,1)*i;else%回溯路由R(i,:)=ones(len,1)*i;endR(i,i)=0;enddisp('');disp('w(0)');dispit(w,0);disp('R(0)');dispit(R,1);%处理端点有权的问题fori=1:lentmp=w(i,i)/2;iftmp~=0w(i,:)=w(i,:)+tmp;w(:,i)=w(:,i)+tmp;flag(i)=1;w(i,i)=0;endend%Floyd算法具体实现过程fori=1:lenforj=1:lenifj==i||w(j,i)==MAXcontinue;endfork=1:lenifk==i||w(j,i)==MAXcontinue;endifw(j,i)+w(i,k)w(j,k)%Floyd算法核心代码w(j,k)=w(j,i)+w(i,k);ifrouter_direction==0%前向路由R(j,k)=R(j,i);else%回溯路由R(j,k)=R(i,k);endendendend%显示每次的计算结果disp(['w(',num2str(i),')'])dispit(w,0);disp(['R(',num2str(i),')'])dispit(R,1);end%中心和中点的确定[Center,index]=min(max(w'));disp(['中心是V',num2str(index)]);[Middle,index]=min(sum(w'));disp(['中点是V',num2str(index)]);endfunctiondispit(x,flag)%x：需要显示的矩阵%flag：为0时表示显示w矩阵，非0时表示显示R矩阵len=length(x);s=[];forj=1:lenifflag==0s=[ssprintf('%5.2f\t',x(j,:))];elses=[ssprintf('%d\t',x(j,:))];ends=[ssprintf('\n')];enddisp(s);disp('---------------------------------------------------');end%选择后按Ctrl+t取消注释号%%%示例：%a=[%0,100,100,1.2,9.2,100,0.5;%100,0,100,5,100,3.1,2;%100,100,0,100,100,4,1.5;%1.2,5,100,0,6.7,100,100;%9.2,100,100,6.7,0,15.6,100;%100,3.1,4,100,15.6,0,100;%0.5,2,1.5,100,100,100,0%];%%b=[%0,9.2,1.1,3.5,100,100;%1.3,0,4.7,100,7.2,100;%2.5,100,0,100,1.8,100;%100,100,5.3,0,2.4,7.5;%100,6.4,2.2,8.9,0,5.1;%7.7,100,2.7,100,2.1,0%];%%Floyd(a,1,100)%Floyd(b,1,100)