initialization.m

function path = initialization(transport_time,customer,vehicle,espo,n)
path(espo).infor = []; 
espoo = 1;
while espoo <= espo
    customers = customer;
    vehicles = vehicle;
    
    %车辆的路径初始化

 for i = 1:n 
     select(i).path = 0; 
     select(i).num = i;
     select(i).price = vehicles(4,i); 
 end
    %
 %或者设置成自行输入   
%    for v = 1:n
    % Assign the number for each vehicle
%    select(v).num = v;
    % Ask the user to input the price for each vehicle
%    select(v).price = input(['请输入第',num2str(v),'辆车的价格：']);
%    end
    
    
    cus = false ;
    while cus ==  false
        % 取出需求量最大的需求点的需求量，并赋值为空。
        max_demand = max(customers(2,:)) ;
        [i,j] = find(customers == max_demand);
        customer_i = customers(1,j(1));
        customers(: ,j(1)) = [];
        vel = false;
        %进行成本的计算
        while vel == false
            % 随机选取车辆
            a = floor(rand(1,1)*size(vehicles,2)+1);
            num = select(a).path(1,end);
            %取出车辆到达此地的时间
            t = transport_time(num+1,customer_i + 1);
            % 判断车辆是否满足需求量和运输时间的要求
            if vehicles(2,a) >= max_demand && vehicles(3,a) > t
                % 如果为真，更新路径，更新时间，更新剩余载重量,将 vel = true
                select(a).path(end + 1) = customer_i;
                vehicles(2,a) = vehicles(2,a) - max_demand;
                vehicles(3,a) = vehicles(3,a) - t;
                vel = true;
            end
        end
        if isempty(customers)
            cus = true;
        end
    end
 % 将路径进行提纯，用到的车辆去除
    for pat = 1:size(select,2)
        if select(pat).path(end) == 0
            continue
        else
            for path_num = 1:size(select(pat).path,2)
                path(espoo).infor(1,end+1) = select(pat).num;
                path(espoo).infor(2,end) = select(pat).price;
                path(espoo).infor(3,end) = select(pat).path(1,path_num);
            end
        end
    end
    espoo =espoo + 1;
end
end
% 将计算的路径进行组合