function [PPMTHseq,THseq] = ...
  cp0201_2PPM_TH(seq,fc,Tc,Ts,dPPM,THcode)
dt = 1 ./ fc;                   
framesamples = floor(Ts./dt);  
chipsamples = floor (Tc./dt);  
PPMsamples = floor (dPPM./dt); 
THp = length(THcode);          
totlength = framesamples*length(seq);
PPMTHseq=zeros(1,totlength);
THseq=zeros(1,totlength);
for k = 1 : length(seq)
    index = 1 + (k-1)*framesamples;
    kTH = THcode(1+mod(k-1,THp));
    index = index + kTH*chipsamples;
    THseq(index) = 1;
    index = index + PPMsamples*seq(k);
    PPMTHseq(index) = 1;
        
end