遗传算法染色体实数编码

最近在做毕业论文，论文是关于车辆路径的，但是卡在染色体编码这个问题上很久了。

我的染色体编码的思路是：

假设在某车辆调度问题中存在N个客户点K辆车进行作业。为避免将不同路线组合在一起时无法区分的问题，在构造染色体时添加虚拟节点0。在编码时首先生成N个客户点的一个全排列，然后将（K+1）个0随机插入此全排列中，并保证首尾为0，可编成长度为（N+K+1）的染色体。举个例子，如果某一个解的编码是：0258607430190，那么它表示有3条路线：2—5—8—6；7—4—3；1—9。

但问题就是，生成N个客户点的一个全排列后，怎么编码将（K+1）个0随机插入此全排列中。下面是我的部分代码：

%% 初始化种群(染色体长度为N+K+1，)
%输入：
% NIND：种群大小
% N：   客户点的数量
% K：   牵引车的数目

%输出：
%初始种群
function Chrom=InitPop(NIND,N,K)
Chrom=zeros(NIND,N,K);%用于存储种群
for i=1:NIND
    Chrom(i,:)=randperm(N);%随机生成染色体
  ???（接下来该怎么写呢？） 
    
end
end

有大佬能为我提供点什么思路吗？