www.gusucode.com > 《MATLAB智能算法30个案例》分析源代码 > 《MATLAB智能算法30个案例》分析源代码/code/chapter24/data.m
clc clear h=[2000 800 500 1000 900 700 1100 1000 1300 1200 1500 200 600 600 1200 1500 600 600 1100 1300 1600 1400 700 500 600 500 300 600 600 600 600 1100 1500 300 600 600 600 400 600 800 600 600 600 1400 800 400 600 1700 900 1000 400 1700 1800 900 600 500 700 800 1600 900 1600 2100 1700 900 1200 1100 300 0 900 1500 1000 1400 1900 2100 1200 600 1800 900 1100 1000 1700 1600 1200 1600 900 1600 1100 1300 600 800 700 600 800 1600 900 1100 1000 600 900 900 800 600 900 800 800 800 1100 600 1400 1300 600 900 1100 800 800 600 900 1200 600 1100 600]; for i=1:11 for j=1:11 h1(2*i-1,j)=h(i,j); end end for i=1:10 for j=1:11 h1(2*i,j)=(h1(2*i-1,j)+h1(2*i+1,j))/2; end end for i=1:21 for j=1:11 h2(i,2*j-1)=h1(i,j); end end for i=1:21 for j=1:10 h2(i,2*j)=(h2(i,2*j-1)+h2(i,2*j+1))/2; end end z=h2; % 初始地形 x=1:21; y=1:21; [x1,y1]=meshgrid(x,y); figure(2) mesh(x1,y1,z) axis([1,21,1,21,0,2000]) xlabel('km') ylabel('km') zlabel('m') for i=1:21 information(i,:,:)=ones(21,21); %初始信息素 end HeightData=z; save HeightData HeightData