untitled.m

求任意一点到任意点的距离的最小值！ 对于求解一些最短距离的问题可以提供一些帮助！

for i=1:66
    for j=1:66
        if i==j
            a0(i,j)=0;
        else
            a0(i,j)=inf;
            p0(i,j)=i;
        end
    end
    r(i)=0;
end
a0(1,2)=2*4;a0(2,3)=2*4;a0(3,4)=3*4;a0(4,5)=0*4;a0(5,6)=2*4;a0(6,7)=3*4;a0(7,8)=4*4;a0(8,9)=5*4;a0(9,10)=2*4;a0(10,11)=2*4;
a0(12,13)=2*4;a0(13,14)=2*4;a0(14,15)=3*4;a0(15,16)=0*4;a0(16,17)=1*4;a0(17,18)=2*4;a0(18,19)=3*4;a0(19,20)=3*4;a0(20,21)=2*4;a0(21,22)=2*4;
a0(23,24)=1*4;a0(24,25)=2*4;a0(25,26)=2*4;a0(26,27)=inf*4;a0(27,28)=1*4;a0(28,29)=1*4;a0(29,30)=3*4;a0(30,31)=3*4;a0(31,32)=5*4;a0(32,33)=4*4;
a0(34,35)=1*4;a0(35,36)=2*4;a0(36,37)=2*4;a0(37,38)=2*4;a0(38,39)=2*4;a0(39,40)=1*4;a0(40,41)=4*4;a0(41,42)=2*4;a0(42,43)=5*4;a0(43,44)=1*4;
a0(45,46)=0*4;a0(46,47)=2*4;a0(47,48)=2*4;a0(48,49)=2*4;a0(49,50)=0*4;a0(50,51)=1*4;a0(51,52)=2*4;a0(52,53)=2*4;a0(53,54)=1*4;a0(54,55)=1*4;
a0(56,57)=0*4;a0(57,58)=3*4;a0(58,59)=3*4;a0(59,60)=4*4;a0(60,61)=0*4;a0(61,62)=5*4;a0(62,63)=2*4;a0(63,64)=3*4;a0(64,65)=1*4;a0(65,66)=3*4;

a0(1,12)=2*3;a0(12,23)=4*3;a0(23,34)=1*3;a0(34,45)=3*3;a0(45,56)=0*3;a0(2,13)=2*3;a0(13,24)=3*3;a0(24,35)=1*3;a0(35,46)=2*3;a0(46,57)=0*3;
a0(3,14)=2*3;a0(14,25)=3*3;a0(25,36)=2*3;a0(36,47)=2*3;a0(47,58)=3*3;a0(4,15)=0*3;a0(15,26)=inf*3;a0(26,37)=2*3;a0(37,48)=2*3;a0(48,59)=3*3;
a0(5,16)=0*3;a0(16,27)=1*3;a0(27,38)=2*3;a0(38,49)=2*3;a0(49,60)=0*3;a0(6,17)=2*3;a0(17,28)=1*3;a0(28,39)=1*3;a0(39,50)=1*3;a0(50,61)=0*3;
a0(7,18)=3*3;a0(18,29)=2*3;a0(29,40)=1*3;a0(40,51)=1*3;a0(51,62)=2*3;a0(8,19)=4*3;a0(19,30)=3*3;a0(30,41)=5*3;a0(41,52)=2*3;a0(52,63)=2*3;
a0(9,20)=2*3;a0(20,31)=3*3;a0(31,42)=inf*3;a0(42,53)=2*3;a0(53,64)=1*3;a0(10,21)=2*3;a0(21,32)=4*3;a0(32,43)=5*3;a0(43,54)=1*3;a0(54,65)=1*3;
a0(11,22)=2*3;a0(22,33)=3*3;a0(33,44)=4*3;a0(44,55)=1*3;a0(55,66)=3*3;

a0(2,1)=2*4;a0(3,2)=2*4;a0(4,3)=3*4;a0(5,4)=0*4;a0(6,5)=2*4;a0(7,6)=3*4;a0(8,7)=4*4;a0(9,8)=5*4;a0(10,9)=2*4;a0(11,10)=2*4;
a0(13,12)=2*4;a0(14,13)=2*4;a0(15,14)=3*4;a0(16,15)=0*4;a0(17,16)=1*4;a0(18,17)=2*4;a0(19,18)=3*4;a0(20,19)=3*4;a0(21,20)=2*4;a0(22,21)=2*4;
a0(24,23)=1*4;a0(25,24)=2*4;a0(26,25)=2*4;a0(27,26)=inf*4;a0(28,27)=1*4;a0(29,28)=1*4;a0(30,29)=3*4;a0(31,30)=3*4;a0(32,31)=5*4;a0(33,32)=4*4;
a0(35,34)=1*4;a0(36,35)=2*4;a0(37,36)=2*4;a0(38,37)=2*4;a0(39,38)=2*4;a0(40,39)=1*4;a0(41,40)=4*4;a0(42,41)=2*4;a0(43,42)=5*4;a0(44,43)=1*4;
a0(46,45)=0*4;a0(47,46)=2*4;a0(48,47)=2*4;a0(49,48)=2*4;a0(50,49)=0*4;a0(51,50)=1*4;a0(52,51)=2*4;a0(53,52)=2*4;a0(54,53)=1*4;a0(55,54)=1*4;
a0(57,56)=0*4;a0(58,57)=3*4;a0(59,58)=3*4;a0(60,59)=4*4;a0(61,60)=0*4;a0(62,61)=5*4;a0(63,62)=2*4;a0(64,63)=3*4;a0(65,64)=1*4;a0(66,65)=3*4;

a0(12,1)=2*3;a0(23,12)=4*3;a0(34,23)=1*3;a0(45,34)=3*3;a0(65,45)=0*3;a0(13,2)=2*3;a0(24,13)=3*3;a0(35,24)=1*3;a0(46,35)=2*3;a0(57,46)=0*3;
a0(14,3)=2*3;a0(25,14)=3*3;a0(36,25)=2*3;a0(47,36)=2*3;a0(58,47)=3*3;a0(15,4)=0*3;a0(26,15)=inf*3;a0(37,26)=2*3;a0(48,37)=2*3;a0(59,48)=3*3;
a0(16,5)=0*3;a0(27,16)=1*3;a0(38,27)=2*3;a0(49,38)=2*3;a0(60,49)=0*3;a0(17,6)=2*3;a0(28,17)=1*3;a0(39,28)=1*3;a0(50,39)=1*3;a0(61,50)=0*3;
a0(18,7)=3*3;a0(29,18)=2*3;a0(40,29)=1*3;a0(51,40)=1*3;a0(62,51)=2*3;a0(19,8)=4*3;a0(30,19)=3*3;a0(41,30)=5*3;a0(52,41)=2*3;a0(63,52)=2*3;
a0(20,9)=2*3;a0(31,20)=3*3;a0(42,31)=inf*3;a0(53,42)=2*3;a0(64,53)=1*3;a0(21,10)=2*3;a0(32,21)=4*3;a0(43,32)=5*3;a0(54,43)=1*3;a0(65,54)=1*3;
a0(22,11)=2*3;a0(33,22)=3*3;a0(44,33)=4*3;a0(55,44)=1*3;a0(66,55)=3*3;


k=0;
while k<66
    for i=1:66
        for j=1:66
            if a0(i,j)>a0(i,k+1)+a0(k+1,j)
                a1(i,j)=a0(i,k+1)+a0(k+1,j);
                p1(i,j)=p0(k+1,j);
            else
                a1(i,j)=a0(i,j);
                p1(i,j)=p0(i,j);
            end
        end
    end
    a0=a1; p0=p1;k=k+1;
end
for j=1:66
    s=0;
    for i=1:66
        s=s+a0(i,j);
    end
    r(j)=s;
end
r;
mina0=min(r)
for i=1:66
    if r(i)==mina0
        i
    end
end
p=0;
for i=1:66
    for j=i+1:66
        p=p+1;
        Q(p)=r(i)+r(j);
    end
end
minQ=min(Q)
p=0;
for i=1:2145
    if Q(i)==minQ
        p=p+1;
        W(p)=i;
    end
end
for s=1:p
    u=0;
    for l=1:66
        for m=l+1:66
            u=u+1;
            if u==W(s)
                L=l
                M=m
            end
        end
    end
end