demapeocept4c4.m

Matlab实现SONET/SDH光通信系统simulink仿真程序源代码

function y =DemapeoCEPT4C4(u)

% Demapeo de los bytes Z
[a b]=size(u);
Nciclo=a/9;%(N de tramas)
p=17408*Nciclo;
q=zeros(1,p);
s=0;
f=1;%(Filas)
contador=0;
for j=1:a
    q(1,1+s:104+s)=u(f,1:104);
    q(1,105+s:200+s)=u(f,113:208);
    q(1,201+s:296+s)=u(f,217:312);
    q(1,297+s:392+s)=u(f,321:416);
    q(1,393+s:488+s)=u(f,425:520);
    q(1,489+s:584+s)=u(f,529:624);
    q(1,585+s:680+s)=u(f,633:728);
    q(1,681+s:776+s)=u(f,737:832);
    q(1,777+s:872+s)=u(f,841:936);
    q(1,873+s:968+s)=u(f,945:1040);
    q(1,969+s:1064+s)=u(f,1049:1144);
    q(1,1065+s:1160+s)=u(f,1153:1248);
    q(1,1161+s:1256+s)=u(f,1257:1352);
    q(1,1257+s:1352+s)=u(f,1361:1456);
    q(1,1353+s:1448+s)=u(f,1465:1560);
    q(1,1449+s:1544+s)=u(f,1569:1664);
    q(1,1545+s:1640+s)=u(f,1673:1768);
    q(1,1641+s:1736+s)=u(f,1777:1872);
    q(1,1737+s:1838+s)=u(f,1881:1982);
    q(1,1839+s:1934+s)=u(f,1985:2080);
   
    contador = contador+1;
   if contador==9
       q(1,1935+s)=u(f-8,1983);% Primera fila de la trama en el bit 'S' se coloca un bit de dato 
       q(1,1936+s)=u(f-7,1983);% Segunda fila de la trama en el bit 'S' se coloca un bit de dato
       s=s+1936;
       contador=0;
   else
       s=s+1934;
   end    
   f=f+1;
end
y=q;
end