codeooc4.m.txt

this is the code for orthogonal code divison multiple access

%/////////////////////////////////////////
%  COPYRIGHT 
% Dr C. Litchfield 
% The University of Kent at Canterbury
% [e-mail]
% optical orthogonal code (OOC)
%////////////////////////////////////////
function [code]=codeooc4()

%clear all

F=100;
w=3;
%code=zeros(N,F);

%--------------------------------------------------------------------------
%-------------------------------------------------------------------------

code{1}=[1,zeros(1,33),1, zeros(1,25) 1,zeros(1,36)];
code{2}=[zeros(1,1),1 zeros(1,22),1, zeros(1,44),1 ,zeros(1,27)];
code{3}=[zeros(1,3),1 zeros(1,38),1, zeros(1,6),1 ,zeros(1,47)];
code{4}=[zeros(1,2),1 zeros(1,4),1, zeros(1,77),1 ,zeros(1,11)];
code{5}=[zeros(1,5),1 zeros(1,9),1, zeros(1,58),1 ,zeros(1,22)];
code{6}=[zeros(1,11),1 zeros(1,19),1, zeros(1,20),1 ,zeros(1,44)];
code{7}=[zeros(1,8),1 zeros(1,14),1, zeros(1,39),1 ,zeros(1,33)];
code{8}=[zeros(1,17),1 zeros(1,12),1, zeros(1,16),1 ,zeros(1,49)];
code{9}=[zeros(1,35),1 zeros(1,25),1, zeros(1,33),1 ,zeros(1,1)];
code{10}=[zeros(1,26),1 zeros(1,44),1, zeros(1,22),1 ,zeros(1,2)];
code{11}=[zeros(1,46),1 zeros(1,6),1, zeros(1,38),1 ,zeros(1,4)];
code{12}=[zeros(1,10),1 zeros(1,77),1, zeros(1,4),1 ,zeros(1,3)];
code{13}=[zeros(1,21),1 zeros(1,58),1, zeros(1,9),1 ,zeros(1,6)];
code{14}=[zeros(1,43),1 zeros(1,20),1, zeros(1,19),1 ,zeros(1,12)];
code{15}=[zeros(1,32),1 zeros(1,39),1, zeros(1,14),1 ,zeros(1,9)];
%code{16}=[zeros(1,48),1 zeros(1,20),1, zeros(1,18),1 ,zeros(1,18)];

L = 15;  % L = 16
for i=1:L
    code1(i,:)=code{i};
end
code=code1;

%--------------------------------------------------------------------------
%--------------------------------------------------------------------------
l1=length(code);
for t=1:L
    tr=randperm(l1);
    t0=tr(2);
    index1=[t0 0 0 ];
    index0=find(code(t,:));
    for i=1:w-1
       retard(i)=index0(i+1)-index0(i);
    end
    retard(w)=length(code)+index0(1)-index0(3);
  
for r=1:w
    if r<w
    retard1(r)=retard(r)+retard(r+1);
    else 
     retard1(r)=retard(r)+ retard(1) ;
    end
end

for k=2:w
        index1(k)=mod((index1(k-1)+retard1(k-1)),length(code)); 
        if index1(k)==0
            index1(k)=1;
        else
        end         
end
code3=zeros(1,length(code));
code3(index1)=1;
code(t,:)=code3;
end