%TSP問題的統(tǒng)計與畫圖
iterfile = fopen( 'c:\saitersfile.txt', 'rt' );
cityfile = fopen( 'c:\sacitysfile.txt', 'rt' );

if( iterfile == -1 || cityfile == -1 )
    display( 'File not exist!!' );
    exit;
end

iters = fscanf( iterfile, '%d %f',[ 2,inf] );
citys = fscanf( cityfile, '%f %f',[ 2,inf] );

fclose( iterfile );
fclose( cityfile );

subplot(1,2,1);
plot( iters(1,:), iters(2,:) ), xlabel('迭代次數(shù)'), ylabel('最短路徑路程'), title('路徑優(yōu)化過程');
subplot(1,2,2);
plot( citys(1,:), citys(2,:),'-b' ), xlabel('Ｘ坐標'), ylabel('Ｙ坐標'), title('求解路徑');
hold on;

plot( citys(1,1), citys(2,1), 'r*', 'MarkerSize', 12 );
for i = 2:length(citys)
    if( mod( i, 2 ) == 1 )
        plot( citys(1,i), citys(2,i), 'r*', 'MarkerSize', 4 );
    else
        plot( citys(1,i), citys(2,i), 'bs', 'MarkerSize', 4 );
    end
end