www.gusucode.com > 《MATLAB智能算法30个案例》分析源代码 > 《MATLAB智能算法30个案例》分析源代码/code/chapter24/data1.m

    clc
clear

h=[1800 1600 1900 1700 2300 2100 2500 2400 2700 2600 2900
   1600 1700 2000 2500 1900 2900 2000 2500 2700 3000 2800
   2100 1900 2500 1900 1700 2400 2700 3000 2300 2500 2900
   1700 2000 2400 2000 1800 2300 2500 2000 2500 2000 2800
   2200 1800 3500 3100 2300 2400 1800 3100 3200 2300 2000
   1900 2100 2600 2600 2300 3000 3500 3100 2300 2600 2500
   1700 1400 2300 2900 2400 2800 1800 3500 2600 2000 3200
   2300 2500 2400 3100 3000 2600 3000 2300 3000 2500 2700
   2000 1800 2600 2000 2200 3000 2300 2500 2400 2000 2300
   2300 1500 2000 2700 2800 2700 3000 2500 2000 2800 2700
   2000 2300 2500 1500 2500 2000 2300 2600 2000 2500 2000];

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);
mesh(x1,y1,z)

for i=1:21
    information(i,:,:)=ones(21,21);     %初始信息素
end

save data z information