论文主要程序附录问题的主要程序1、问题一的主要程序1.1问题一、问题三模型程序的main.m源文件%%若要解决问题三时给pingtai加一些相应的点如293988...clear;%%建立所有可行的路线sj2;loaddata;[m,n]=size(link);link1=sort(link,2);k=1;fori=1:mif(link1(i,2)<=92)luxian(k,:)=link1(i,:);k=k+1;endendsaveluxianluxian;%%%解决问题三时给pingtai加一些相应的点如293988.。。pingtai=[1:20];num=size(pingtai,2);savepingtaipingtai;%%初始运算loadfananlv2;loadju;%%找出路得属性并作图jieg=luu(ju,luxian,pingtai);loadjuli;loadzuobiao;juli1=juli(1:92,1:92);zuotu(juli1,x,y,jieg,pingtai);biaodian2%%增添中点，给属性，作图，给出新路的犯案率jieg1=zhongdian(jieg);loadjuli2;loadzuobiao2;loadju2;jieg2=luu(ju2,jieg1(:,[12]),pingtai);fananlv=fanan1(jieg2);figure(2)zuotu(juli2,x,y,jieg2,pingtai);biaodian;%%0-1线性规划得到分配方案jieg3=juliw(jieg2,ju2,pingtai);%%给一些分配了路线的付初值fanan=zeros(num,1);[mm,nn]=size(jieg3);fori=1:mmfanan(jieg3(i,1),1)=[fanan(jieg3(i,1),1)+fananlv(jieg3(i,4),4)];endsavefananfanan%%退火算法第二次规划[jieg4fananSum]=tuih(jieg2,pingtai,1.5);%%做出最终的结果图像jieguofigure(3)zuitu(jieg4,jieg3,pingtai)1.2模拟退火算法在第一问中的求解二次规划问题function[jiegqfananSum]=tuih(jieg2,pingtai,fazhi)Sum=inf;loadfanan;loadfananlv1;loadju2;%loadaaaa[mm,nn]=size(jieg2);mm1=size(pingtai,2);k=1;fori=1:mmifjieg2(i,3)==2luxian(k,:)=[jieg2(i,[12])iju2(jieg2(i,1),jieg2(i,2))];%luxian(k,4)=exp((luxian(k,4)-3000)/500);fananlv2(k,:)=fananlv1(i,4);k=k+1;endendsavefananlv2fananlv2mm2=k-1;fori=1:mm1forj=1:mm2ss(i,j)=min(ju2(i,luxian(j,1)),ju2(i,luxian(j,2)))+ju2(luxian(j,1),luxian(j,2));endend%t=zeros(140,20);forj=1:mm2k(j)=1;fori=1:mm1ifmax(ju2(pingtai(i),luxian(j,1)),ju2(pingtai(i),luxian(j,2)))<3000t(j,k(j))=i;k(j)=k(j)+1;endendendpp1=mm2;%%设置初始状态bclearxuan0;whileSum==infxuanz=zeros(1,mm1*mm2);fori=1:mm2aa=randperm(k(i)-1);a(i)=t(i,aa(1));xuanz((a(i)-1)*mm2+i)=1;endifmanyi(xuanz)<=fazhi%%以满意函数为目标多次得到的经验值temp=mubiao(xuanz,ss);iftemp
