function [so,sm,sv] = heuristic_Population()
%% 功能：启发式初始化
global AGV_num ps job_number machine_number operations_number operations_vector operations_start_vector feasible_machines_table time_table trans_time_table operations_number_vector
so1=zeros(ps,operations_number);
sm1=zeros(ps,operations_number); % 存机器
sv1=zeros(ps,operations_number); % 存机器
for i=1:ps
%% 工序：初始化
indices=randperm(operations_number); % 将300以内的数随机排列
so1(i,:)=operations_vector(indices); % 得到工件号排序的工序
%% 机器：选择完工时间最早的机器
ML=zeros(1,machine_number); % 记录全局机器时间
TT=zeros(1,machine_number); % 用来临时比较
% 生成机器 
Number=operations_number_vector;
WPNumberTemp=Number;
for j=1:operations_number
    val=so1(i,j);
    WPNumberTemp(val)=WPNumberTemp(val)-1;
    temp=[val,operations_number_vector(val)-WPNumberTemp(val)]; % 第一个数表示工件号，第二个数表示工件内对应的工序号
    op_position=sum(operations_number_vector(1:val-1))+temp(2); % 总工序表中第几个 如O13表示3
    machines_list = 1:machine_number;
    Temp1 = machines_list(logical(feasible_machines_table(op_position,:))); % 工序Oij可选机器
    time1 = time_table(op_position,:); % 例[3 7]
    time1(find(time1==0))=[];
    Temp2=time1;
    Temp=[Temp1;Temp2];
    SizeTemp=size(Temp,2);%查看该工序有几个可用机器
    if SizeTemp == 1
        sm1(i,op_position)=Temp(1,1); % 如果只有一个可选机器，就直接存入
        ML(Temp(1))=Temp(2,1); % 更新ML表
    elseif SizeTemp>=2
        % TT=ML+time 选全局最小机器
        for ii=1:SizeTemp
            job=Temp(1,ii); % 取可选机器集内第ii个机器
            time=Temp(2,ii); % 取对应时间
            TT(job)=ML(job)+time; % 加入TT时间表用来比较选机器
        end
        [~,indice]=min(TT(Temp(1,:)));
        TT=zeros(1,machine_number); % 清空TT表，以备下一个工序选机器
        ML(Temp(1,indice))=Temp(2,indice); % 更新ML
        sm1(i,op_position)=Temp(1,indice);
    end
end
    %% AGV初始化
%     AGV_current_machine=zeros(AGV_num,operations_number*2); % 存AGV当前在哪个机器
    machineoRagV=zeros(AGV_num,machine_number+1); % 用0-1存当前机器是否有AGV
    machineoRagV(:,1)=1; % 首先对装卸站位置赋值（所有的AGV初始时刻均在装卸站）
    job_current_op=ones(1,job_number);
    for j=1:operations_number
%         j
        job=so1(i,j); operation=job_current_op(1,job);
        index=operations_start_vector(1,job)+operation; % 该工序索引
        job_current_op(1,job)=job_current_op(1,job)+1; % 更新当前加工工序
        if operation==1  % 判断是否是首道工序
%             sv1(1,index)=randi([1,AGV_num],1,1); % 首道工序随机选择AGV
%         else
            if ismember(1,machineoRagV(:,1)) % 判断起点是否有AGV
               tempm=sm1(i,index); % 当前工序对应机器编号
               tempm=tempm+1;
               startindex=find(machineoRagV(:,1)==1); % 装卸站可用的AGV编号
               AGV_code=startindex(randi([1 length(startindex)])); % 选择在装卸站的AGV
               machineoRagV(AGV_code,1)=0; machineoRagV(AGV_code,tempm)=1; % 更新AGV当前所在的位置
               sv1(i,index)=AGV_code; % 将选择的AGV的编号放入初始编码中
            else % 装卸站没有AGV
                tempm=sm1(i,index); % 当前工序对应机器编号
                tempm=tempm+1;
                distance=zeros(1,AGV_num); % 初始距离矩阵
                agvmac=[];
                for k=1:AGV_num % 循环所有的AGV
                    % 计算距离
                    nowAgvIndex=find(machineoRagV(k,:)==1); % 当前AGV所在的位置
                    agvmac=[agvmac;nowAgvIndex]; % 存一下agv所在机器
                    distance(1,k)=trans_time_table(nowAgvIndex,1)+trans_time_table(1,tempm);
                end
                % 选择距离最小的AGV
                 [~,sv1(i,index)]=min(distance);
                  % 更新一下当前AGV的位置
               machineoRagV(sv1(i,index),agvmac(sv1(i,index),1))=0; machineoRagV(sv1(i,index),tempm)=1; % 更新AGV当前所在的位置
            agvmac=[];
            end
        else % 非首道工序，只需要计算与上一个工序所在机器
            tempm=sm1(i,index); % 当前工序对应机器编号
            premac=sm1(i,index-1); % 上一个工序所在机器编号
            premac1=premac+1;
            tempm=tempm+1;
            distance=zeros(1,AGV_num); % 初始距离矩阵
            agvmac=[];
            for k=1:AGV_num % 循环所有的AGV
                % 计算距离
                nowAgvIndex=find(machineoRagV(k,:)==1); % 当前AGV所在的位置
%                 nowAgvIndex
                agvmac=[agvmac;nowAgvIndex]; % 存一下agv所在机器
                distance(1,k)=trans_time_table(premac1,tempm)+trans_time_table(nowAgvIndex,premac1);
            end
            % 选择距离最小的AGV
            [~,sv1(i,index)]=min(distance);
            % 更新一下当前AGV的位置
            machineoRagV(sv1(i,index),agvmac(sv1(i,index),1))=0; machineoRagV(sv1(i,index),tempm)=1; % 更新AGV当前所在的位置
            agvmac=[];
        end
    end

 so=so1;sm=sm1;sv=sv1;
end