基于GA遗传优化的OSPF协议网络最优路由matlab仿真，含仿真操作录像

clc; clear; close all; warning off; addpath(genpath(pwd)); webpoint=30; %节点数目 [D,M]=tsp(webpoint); popsize=20; N=199; %N为迭代次数 pcro=0.85; pmut=0.5; gen=1; minval=inf; pop=initialize(popsize,webpoint); while gen<=N val=adapt(pop,popsize,webpoint,D); fitval=(1000./val).^15; [valmin,minind]=min(val); %保存最优 if valmin<minval minval=valmin; minpath=pop(minind,:); end minf(gen)=minval; %最短路径长度 pop=select(pop,fitval,popsize,webpoint,minpath); %比例选择 pop=crossover(pop,pcro,popsize,webpoint); %OX交叉 pop=mutation(pop,pmut,popsize,webpoint); %变异 drawTSP(M,minpath,minval,gen,0); gen=gen+1; end drawTSP(M,minpath,minval,gen,1); hold on r=0.022; t=0:pi/180:2*pi; x=r*cos(t)+0.03; y=r*sin(t)+0.03; plot(x,y,'b--');hold on; plot([0.03,0.01],[0.03,0.03],'b--');hold on; plot([0.03,0.04],[0.03,0.01],'b--');hold on; plot([0.03,0.02],[0.03,0.01],'b--');hold on; plot([0.03,0.05],[0.03,0.03],'b--');hold on; plot([0.03,0.03],[0.03,0.05],'b--');hold on; plot([0.02,0.04],[0.01,0.01],'b--');hold on; plot([0.04,0.05],[0.01,0.03],'b--');hold on; plot([0.05,0.03],[0.03,0.05],'b--');hold on; plot([0.03,0.01],[0.05,0.03],'b--');hold on; plot([0.01,0.02],[0.03,0.01],'b--');hold on; %================================================ %================================================ hold on t=0:pi/180:2*pi; x=r*cos(t)+0.13; y=r*sin(t)+0.03; plot(x,y,'b--');hold on plot([0.13,0.11],[0.03,0.03],'b--');hold on; plot([0.13,0.14],[0.03,0.01],'b--');hold on; plot([0.13,0.12],[0.03,0.01],'b--');hold on; plot([0.13,0.15],[0.03,0.03],'b--');hold on; plot([0.13,0.13],[0.03,0.05],'b--');hold on; plot([0.12,0.14],[0.01,0.01],'b--');hold on; plot([0.14,0.15],[0.01,0.03],'b--');hold on; plot([0.15,0.13],[0.03,0.05],'b--');hold on; plot([0.13,0.11],[0.05,0.03],'b--');hold on; plot([0.11,0.12],[0.03,0.01],'b--');hold on; %================================================ %================================================ t=0:pi/180:2*pi; x=r*cos(t)+0.23; y=r*sin(t)+0.03; plot(x,y,'b--');hold on plot([0.23,0.21],[0.03,0.03],'b--');hold on; plot([0.23,0.24],[0.03,0.01],'b--');hold on; plot([0.23,0.22],[0.03,0.01],'b--');hold on; plot([0.23,0.25],[0.03,0.03],'b--');hold on; plot([0.23,0.23],[0.03,0.05],'b--');hold on; plot([0.22,0.24],[0.01,0.01],'b--');hold on; plot([0.24,0.25],[0.01,0.03],'b--');hold on; plot([0.25,0.23],[0.03,0.05],'b--');hold on; plot([0.23,0.21],[0.05,0.03],'b--');hold on; plot([0.21,0.22],[0.03,0.01],'b--');hold on; %================================================ %================================================ t=0:pi/180:2*pi; x=r*cos(t)+0.08; y=r*sin(t)+0.08; plot(x,y,'b--');hold on plot([0.08,0.07],[0.08,0.06],'b--');hold on; plot([0.08,0.09],[0.08,0.06],'b--');hold on; plot([0.08,0.06],[0.08,0.08],'b--');hold on; plot([0.08,0.10],[0.08,0.08],'b--');hold on; plot([0.08,0.08],[0.08,0.10],'b--');hold on; plot([0.07,0.09],[0.06,0.06],'b--');hold on; plot([0.09,0.10],[0.06,0.08],'b--');hold on; plot([0.10,0.08],[0.08,0.1],'b--');hold on; plot([0.08,0.06],[0.1,0.08],'b--');hold on; plot([0.06,0.07],[0.08,0.06],'b--');hold on; %================================================ %================================================ t=0:pi/180:2*pi; x=r*cos(t)+0.18; y=r*sin(t)+0.08; plot(x,y,'b--');hold on; plot([0.18,0.17],[0.08,0.06],'b--');hold on; plot([0.18,0.19],[0.08,0.06],'b--');hold on; plot([0.18,0.16],[0.08,0.08],'b--');hold on; plot([0.18,0.20],[0.08,0.08],'b--');hold on; plot([0.18,0.18],[0.08,0.10],'b--');hold on; plot([0.17,0.19],[0.06,0.06],'b--');hold on; plot([0.19,0.20],[0.06,0.08],'b--');hold on; plot([0.20,0.18],[0.08,0.1],'b--');hold on; plot([0.18,0.16],[0.1,0.08],'b--');hold on; plot([0.16,0.17],[0.08,0.06],'b--');hold on; %================================================ %================================================ figure;plot(minf,'-r');xlabel('迭代次数');ylabel('最短路径长度'); disp('最短路径为:');disp(minpath); disp('最短路径长度为:');disp(minval);