entrama.asv

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

function [sys]=Entrama(serial)

[a b]=size(serial);
h=0;

curBlk=([gcs,'/','Entramador']);
Tipo_Senal=get_param(curBlk,'SSS');

switch (Tipo_Senal)
    case 'STS-1'
        y=b/6480; % 9*720=6480
        trama=zeros(y*9,720);
        s=0;
        q=1;
        for k=1:y
            trama(q,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+1,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+2,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+3,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+4,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+5,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+6,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+7,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            trama(q+8,1:720)=serial(1,((1+s):(720+s)));
            s=s+720;
            q=q+9;
        end
        
    case 'STS-3'
        y=b/19440;  %  9*2160=19440
        trama=zeros(y*9,2160);
        s=0;
        q=1;
        for k=1:y
            trama(q,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+1,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+2,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+3,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+4,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+5,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+6,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+7,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            trama(q+8,1:2160)=serial(1,((1+s):(2160+s)));
            s=s+2160;
            q=q+9;
        end
        
    case 'STS-12'
        y=b/77760;  %  9*8640=77760
        trama=zeros(y*9,8640);
        s=0;
        q=1;
        for k=1:y
            trama(q,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+1,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+2,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+3,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+4,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+5,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+6,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+7,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            trama(q+8,1:8640)=serial(1,((1+s):(8640+s)));
            s=s+8640;
            q=q+9;
        end
        
    case 'STS-48'
        y=b/311040;  %  9*34560=311040
        trama=zeros(y*9,34560);
        s=0;
        q=1;
        for k=1:y
            trama(q,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+1,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+2,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+3,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+4,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+5,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+6,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+7,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            trama(q+8,1:34560)=serial(1,((1+s):(34560+s)));
            s=s+34560;
            q=q+9;
        end
        
    case 'STS-192'
        y=b/1244160;  %  9*138240=1244160
        trama=zeros(y*9,138240);
        s=0;
        q=1;
        for k=1:y
            trama(q,1:138240)=serial(1,((1+s):(1244160+s)));
            s=s+138240;
            trama(q+1,1:138240)=serial(1,((1+s):(1244160+s)));
            s=s+138240;
            trama(q+2,1:138240)=serial(1,((1+s):(1244160+s)));
            s=s+138240;
            trama(q+3,1:1244160)=serial(1,((1+s):(1244160+s)));
            s=s+138240;
            trama(q+4,1:1244160)=serial(1,((1+s):(1244160+s)));
            s=s+1244160;
            trama(q+5,1:1244160)=serial(1,((1+s):(1244160+s)));
            s=s+1244160;
            trama(q+6,1:1244160)=serial(1,((1+s):(1244160+s)));
            s=s+1244160;
            trama(q+7,1:1244160)=serial(1,((1+s):(1244160+s)));
            s=s+1244160;
            trama(q+8,1:1244160)=serial(1,((1+s):(1244160+s)));
            s=s+1244160;
            q=q+9;
        end
        
    case 'STS-768'
        y=b/552960;
        trama=zeros(y*9,4976640);
        s=0;
        q=1;
        for k=1:y
            trama(q,1:4976640)=serial(1,((1+s):(2160+s)));
            s=s+4976640;
            trama(q+1,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+2,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+3,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+4,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+5,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+6,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+7,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            trama(q+8,1:4976640)=serial(1,((1+s):(4976640+s)));
            s=s+4976640;
            q=q+9;
        end
        
end

sys=trama(:,:);

end