I have developed Matlab code of an algorithm, but I am unable to run it. Please help me to find error.
    3 次查看(过去 30 天)
  
       显示 更早的评论
    
I have developed Matlab code of  an algorithm, but I am unable to run it. Please help me to ffix the errors. 
Urgent reply will be highly appreciated. Thanks in advnce!
distance_matrix = [0, 29, 20, 21;
                   29, 0, 15, 18;
                   20, 15, 0, 30;
                   21, 18, 30, 0];
tabu_tenure = 5;
max_iterations = 1000;
[best_tour, best_distance] = tabu_search_tsp(distance_matrix, tabu_tenure, max_iterations);
function [best_tour, best_distance] = tabu_search_tsp(distance_matrix, tabu_list_size, max_iterations)
    n = size(distance_matrix, 1);  % Number of cities
    best_tour = randperm(n);       % Initialize the best tour with a random permutation
    best_distance = calculate_tour_distance(best_tour, distance_matrix);
    tabu_list = zeros(tabu_list_size, n); % Initialize the tabu list to store forbidden moves
    tabu_tenure = 5;  % Tabu tenure (number of iterations for which a move is forbidden)
    % Main loop of Tabu Search
    for iter = 1:max_iterations
        initial_solution = nearest_neighbor_heuristic(distance_matrix);
        [current_tour, current_distance] = reversal_heuristic(initial_solution, distance_matrix);
        % Check if the current tour is better than the best tour found so far
        if current_distance < best_distance
            best_tour = current_tour;
            best_distance = current_distance;
        end
        % Add the current tour to the tabu list
        tabu_list = update_tabu_list(tabu_list, current_tour, tabu_tenure);
        % Perform Tabu Search on the current tour
        [current_tour, current_distance] = perform_tabu_search(current_tour, distance_matrix, tabu_list);
        % Check if the current tour is better than the best tour found so far
        if current_distance < best_distance
            best_tour = current_tour;
            best_distance = current_distance;
        end
        % Display the current best distance every 100 iterations (optional)
        if mod(iter, 100) == 0
            fprintf('Iteration %d, Best Distance: %.2f\n', iter, best_distance);
        end
    end
end
function distance = calculate_tour_distance(tour, distance_matrix)
    n = length(tour);
    distance = 0;
    for i = 1:n-1
        distance = distance + distance_matrix(tour(i), tour(i+1));
    end
    distance = distance + distance_matrix(tour(n), tour(1)); % Return to the starting city
end
function tour = nearest_neighbor_heuristic(distance_matrix)
    n = size(distance_matrix, 1);  % Number of cities
    tour = zeros(1, n);
    unvisited = 1:n;
    current_city = randi(n);  % Start from a random city
    for i = 1:n
        tour(i) = current_city;
        unvisited(unvisited == current_city) = [];
        if isempty(unvisited)
            break;
        end
        [~, next_city] = min(distance_matrix(current_city, unvisited));
        current_city = unvisited(next_city);
    end
end
function [tour, distance] = reversal_heuristic(tour, distance_matrix)
    n = length(tour);
    distance = calculate_tour_distance(tour, distance_matrix);
    for i = 1:n-1
        for j = i+1:n
            % Try reversing the segment between city i and city j
            new_tour = apply_reverse(tour, i, j);
            new_distance = calculate_tour_distance(new_tour, distance_matrix);
            % If the new tour is better, update the tour and distance
            if new_distance < distance
                tour = new_tour;
                distance = new_distance;
            end
        end
    end
end
function new_tour = apply_reverse(tour, i, j)
    new_tour = [tour(1:i-1), fliplr(tour(i:j)), tour(j+1:end)];
end
function tabu_list = update_tabu_list(tabu_list, tour, tabu_tenure)
    tabu_list(2:end, :) = tabu_list(1:end-1, :);
    tabu_list(1, :) = tour;
    tabu_list(:, tabu_tenure+1:end) = 0;
end
function [best_tour, best_distance] = perform_tabu_search(initial_tour, distance_matrix, tabu_list)
end
1 个评论
  John D'Errico
      
      
 2023-8-4
				Does your need for an urgent reply make anyone reply more quickly? No. All it does is make somepne like me waste my time, telling you why your question is not in fact urgent. Unless, of course, you are much more important than any of the other people who have asked questions. I am sure they would like a quick response too.
采纳的回答
  Voss
      
      
 2023-8-4
        The function perform_tabu_search doesn't do anything.
function [best_tour, best_distance] = perform_tabu_search(initial_tour, distance_matrix, tabu_list)
end
That's the reason for the error you got. Specifically, perform_tabu_search doesn't assign a value to its output arguments best_tour and best_distance.
更多回答(0 个)
另请参阅
类别
				在 Help Center 和 File Exchange 中查找有关 Get Started with Optimization Toolbox 的更多信息
			
	Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!


