demapeocept3c3.asv

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

function y = DemapeoCEPT3C3(u)

[a b]=size(u);
Nciclo=b/9;
Nfila=Nciclo*9;
Ntotalbits=(Nciclo)*4296;
q=zeros(1,Ntotalbits);
for j=1:Nciclo
    for i=1:3 % (T1 T2 T3)
        % fila 1
        q(1,1+H:24+H)=   u(fila+1,9:32); H=H+24;
        q(1,1+H:24+H)=   u(fila+1,41:64); H=H+24;
        q(1,1+H:24+H)= u(fila+1,73:96); H=H+24;
        q(1,1+H:24+H)=u(fila+1,105:128); H=H+24;
        q(1,1+H:24+H)=u(fila+1,145:168); H=H+24;
        q(1,1+H:24+H)=u(fila+1,177:200); H=H+24;
        q(1,1+H:24+H)=u(fila+1,209:232); H=H+24;
        q(1,1+H:24+H)=u(fila+1,241:264); H=H+24;
        q(1,1+H:24+H)=u(fila+1,273:296); H=H+24;
        q(1,1+H:24+H)=u(fila+1,313:336); H=H+24;
        q(1,1+H:24+H)=u(fila+1,345:368); H=H+24;
        q(1,1+H:24+H)=u(fila+1,377:400); H=H+24;
        q(1,1+H:24+H)=u(fila+1,409:432); H=H+24;
        q(1,1+H:24+H)=u(fila+1,441:464); H=H+24;
        q(1,1+H:24+H)=u(fila+1,481:504); H=H+24;
        q(1,1+H:24+H)=u(fila+1,513:536); H=H+24;
        q(1,1+H:24+H)=u(fila+1,545:568); H=H+24;
        q(1,1+H:24+H)=u(fila+1,577:600); H=H+24;
        q(1,1+H:24+H)=u(fila+1,609:632); H=H+24;
        q(1,1+H:24+H)=u(fila+1,649:672); H=H+24;
        % fila 2
        q(1,1+H:24+H)=  u(fila+2,9:32); H=H+24;
        q(1,1+H:24+H)= u(fila+2,41:64); H=H+24;
        q(1,1+H:24+H)= u(fila+2,73:96); H=H+24;
        q(1,1+H:24+H)=u(fila+2,105:128); H=H+24;
        q(1,1+H:24+H)=u(fila+2,145:168); H=H+24;
        q(1,1+H:24+H)=u(fila+2,177:200); H=H+24;
        q(1,1+H:24+H)=u(fila+2,209:232); H=H+24;
        q(1,1+H:24+H)=u(fila+2,241:264); H=H+24;
        q(1,1+H:24+H)=u(fila+2,273:296); H=H+24;
        q(1,1+H:24+H)=u(fila+2,313:336); H=H+24;
        q(1,1+H:24+H)=u(fila+2,345:368); H=H+24;
        q(1,1+H:24+H)=u(fila+2,377:400); H=H+24;
        q(1,1+H:24+H)=u(fila+2,409:432); H=H+24;
        q(1,1+H:24+H)=u(fila+2,441:464); H=H+24;
        q(1,1+H:24+H)=u(fila+2,481:504); H=H+24;
        q(1,1+H:24+H)=u(fila+2,513:536); H=H+24;
        q(1,1+H:24+H)=u(fila+2,545:568); H=H+24;
        q(1,1+H:24+H)=u(fila+2,577:600); H=H+24;
        q(1,1+H:24+H)=u(fila+2,609:632); H=H+24;
        q(1,1+H:24+H)=u(fila+2,649:672); H=H+24;
        % fila 3
        q(1,1+H:24+H)=  u(fila+3,9:32); H=H+24;
        q(1,1+H:24+H)= u(fila+3,41:64); H=H+24;
        q(1,1+H:24+H)= u(fila+3,73:96); H=H+24;
        q(1,1+H:24+H)(fila+3,105:128)=u; H=H+24;
        q(fila+3,145:168)=u(1,1+H:24+H); H=H+24;
        q(fila+3,177:200)=u(1,1+H:24+H); H=H+24;
        q(fila+3,209:232)=u(1,1+H:24+H); H=H+24;
        q(fila+3,241:264)=u(1,1+H:24+H); H=H+24;
        q(fila+3,273:296)=u(1,1+H:24+H); H=H+24;
        q(fila+3,313:336)=u(1,1+H:24+H); H=H+24;
        q(fila+3,345:368)=u(1,1+H:24+H); H=H+24;
        q(fila+3,377:400)=u(1,1+H:24+H); H=H+24;
        q(fila+3,409:432)=u(1,1+H:24+H); H=H+24;
        q(fila+3,441:464)=u(1,1+H:24+H); H=H+24;
        q(fila+3,481:504)=u(1,1+H:24+H); H=H+24;
        q(fila+3,513:536)=u(1,1+H:24+H); H=H+24;
        q(fila+3,545:568)=u(1,1+H:24+H); H=H+24;
        q(fila+3,577:600)=u(1,1+H:24+H); H=H+24;
        q(fila+3,609:632)=u(1,1+H:24+H); H=H+24;
        q(fila+3,657:672)=u(1,1+H:16+H); H=H+16;
        fila = fila+3;
    end
end