interbyte64stm1.m

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

              %**********************************************
              %* Multiplexa byte a byte 64 STM-1 en un STM-64 * 
              %**********************************************

%   1 -----)  \              
%   2 -----)   \
%   3 -----)    \
%   4 -----)     \
%   5 -----)      \
%   6 -----)       \
%  .. -----)        \----\ columnas...12 11 10 9 8 7 6 5 4 3 2 1                
%  .. -----)        /----/ 
%  .. -----)       /
%  .. -----)      /
%  .. -----)     /
%  .. -----)    /
%  .. -----)   /
%  .. -----)  /
%  63 -----) /
%  64 -----)/

function [e]=interbyte64STM1(u)

 [r,t]=size(u);
 e=zeros(r,t);
 s=1; i=1; s1=1; s2=1; s3=1; s4=1; s5=1; s6=1; s7=1; s8=1; s9=1; s10=1; s11=1; S12=1; S13=1; S14=1; S15=1; s16=1; s17=1; s18=1; s19=1;
 s20=1; s21=1; s22=1; s23=1; s24=1; s25=1; s26=1; S27=1; S28=1; S29=1; S30=1; s31=1; s32=1; s33=1; s34=1; s35=1; s36=1; s37=1; s38=1;
 s39=1; s40=1; s41=1; S42=1; S43=1; S44=1; S45=1; s46=1; s47=1; s48=1; s49=1; s50=1; s51=1; s52=1; s53=1; s54=1; s55=1; s56=1; S57=1;
 S58=1; S59=1; S60=1; s61=1; s62=1; s63=1; s64=1;
 while i<=(t)
   Q=1;W=1;H=1;E=1;O=1;D=1;G=1;F=1;L=1;Z=1;C=1;V=1;QQ=1;QW=1;QH=1;QE=1; 
   Q1=1;W1=1;H1=1;E1=1;O1=1;D1=1;G1=1;F1=1;L1=1;Z1=1;C1=1;V1=1;QQ1=1;QW1=1;QH1=1;QE1=1; 
   Q2=1;W2=1;H2=1;E2=1;O2=1;D2=1;G2=1;F2=1;L2=1;Z2=1;C2=1;V2=1;QQ2=1;QW2=1;QH2=1;QE2=1;
   Q3=1;W3=1;H3=1;E3=1;O3=1;D3=1;G3=1;F3=1;L3=1;Z3=1;C3=1;V3=1;QQ3=1;QW3=1;QH3=1;QE3=1;
   while H<=8
   e(:,i)=u(:,s);       % columna 1
   i=i+1;
   s=s+1;
   H=H+1;
   end
   while Q<=8
   e(:,i)=u(:,s1+2160);  % columna 2
   i=i+1;
   s1=s1+1;
   Q=Q+1;
   end
   while W<=8
   e(:,i)=u(:,s2+4320);  % columna 3
   i=i+1;
   s2=s2+1;
   W=W+1;
   end
   while E<=8
   e(:,i)=u(:,s3+6480);  % columna 4
   i=i+1;
   s3=s3+1;
   E=E+1; 
   end
   while O<=8
   e(:,i)=u(:,s4+8640);   % columna 5
   i=i+1;
   s4=s4+1;
   O=O+1;
   end
   while D<=8
   e(:,i)=u(:,s5+10800);  % columna 6
   i=i+1;
   s5=s5+1;
   D=D+1;
   end
   while F<=8
   e(:,i)=u(:,s6+12960);  % columna 7
   i=i+1;
   s6=s6+1;
   F=F+1;
   end
   while G<=8
   e(:,i)=u(:,s7+15120);  % columna 8
   i=i+1;
   s7=s7+1;
   G=G+1;
   end
   while L<=8
   e(:,i)=u(:,s8+17280);  % columna 9
   i=i+1;
   s8=s8+1;
   L=L+1;
   end
   while Z<=8
   e(:,i)=u(:,s9+19440);  % columna 10
   i=i+1;
   s9=s9+1;
   Z=Z+1;
   end
   while C<=8  
   e(:,i)=u(:,s10+21600);  % columna 11
   i=i+1;
   s10=s10+1;
   C=C+1;
   end
   while V<=8
   e(:,i)=u(:,s11+23760);  % columna 12
   i=i+1;
   s11=s11+1;
   V=V+1;
   end
   while QQ<=8
   e(:,i)=u(:,s12+25920);  % columna 13
   i=i+1;
   s12=s12+1;
   Q1=Q1+1;
   end
    while QW<=8
   e(:,i)=u(:,s13+28080);  % columna 14
   i=i+1;
   s13=s13+1;
   Q2=Q2+1;
   end
    while QH<=8
   e(:,i)=u(:,s14+30240);  % columna 15
   i=i+1;
   s14=s14+1;
   Q3=Q3+1;
   end
    while QE<=8
   e(:,i)=u(:,s15+32400);  % columna 16
   i=i+1;
   s15=s15+1;
   Q4=Q4+1;
   end
   
   
   while H1<=8
   e(:,i)=u(:,s16+34560);       % columna 1
   i=i+1;
   s16=s16+1;
   H1=H1+1;
   end
   while Q1<=8
   e(:,i)=u(:,s17+36720);  % columna 2
   i=i+1;
   s17=s17+1;
   Q1=Q1+1;
   end
   while W1<=8
   e(:,i)=u(:,s18+38880);  % columna 3
   i=i+1;
   s18=s18+1;
   W1=W1+1;
   end
   while E1<=8
   e(:,i)=u(:,s19+41040);  % columna 4
   i=i+1;
   s19=s19+1;
   E1=E1+1; 
   end
   while O1<=8
   e(:,i)=u(:,s20+43200);   % columna 5
   i=i+1;
   s20=s20+1;
   O1=O1+1;
   end
   while D1<=8
   e(:,i)=u(:,s21+45360);  % columna 6
   i=i+1;
   s21=s21+1;
   D1=D1+1;
   end
   while F1<=8
   e(:,i)=u(:,s22+47520);  % columna 7
   i=i+1;
   s22=s22+1;
   F1=F1+1;
   end
   while G1<=8
   e(:,i)=u(:,s23+49680);  % columna 8
   i=i+1;
   s23=s23+1;
   G1=G1+1;
   end
   while L1<=8
   e(:,i)=u(:,s24+51840);  % columna 9
   i=i+1;
   s24=s24+1;
   L1=L1+1;
   end
   while Z1<=8
   e(:,i)=u(:,s25+54000);  % columna 10
   i=i+1;
   s25=s25+1;
   Z1=Z1+1;
   end
   while C1<=8  
   e(:,i)=u(:,s26+56160);  % columna 11
   i=i+1;
   s26=s26+1;
   C1=C1+1;
   end
   while V1<=8
   e(:,i)=u(:,s27+58320);  % columna 12
   i=i+1;
   s27=s27+1;
   V1=V1+1;
   end
   while QQ1<=8
   e(:,i)=u(:,s28+60480);  % columna 13
   i=i+1;
   s28=s28+1;
   QQ1=QQ1+1;
   end
    while QW1<=8
   e(:,i)=u(:,s29+62640);  % columna 14
   i=i+1;
   s29=s29+1;
   QW1=QW1+1;
   end
    while QH1<=8
   e(:,i)=u(:,s30+64800);  % columna 15
   i=i+1;
   s30=s30+1;
   QH1=QH1+1;
   end
    while QE1<=8
   e(:,i)=u(:,s31+66960);  % columna 16
   i=i+1;
   s31=s31+1;
   QE1=QE1+1;
   end
   
   while H2<=8
   e(:,i)=u(:,s32+69120);       % columna 1
   i=i+1;
   s32=s32+1;
   H2=H2+1;
   end
   while Q2<=8
   e(:,i)=u(:,s33+71280);  % columna 2
   i=i+1;
   s33=s33+1;
   Q2=Q2+1;
   end
   while W2<=8
   e(:,i)=u(:,s34+73440);  % columna 3
   i=i+1;
   s34=s34+1;
   W2=W2+1;
   end
   while E2<=8
   e(:,i)=u(:,s35+75600);  % columna 4
   i=i+1;
   s35=s35+1;
   E2=E2+1; 
   end
   while O2<=8
   e(:,i)=u(:,s36+77760);   % columna 5
   i=i+1;
   s36=s36+1;
   O2=O2+1;
   end
   while D2<=8
   e(:,i)=u(:,s37+79920);  % columna 6
   i=i+1;
   s37=s37+1;
   D2=D2+1;
   end
   while F2<=8
   e(:,i)=u(:,s38+82080);  % columna 7
   i=i+1;
   s38=s38+1;
   F2=F2+1;
   end
   while G2<=8
   e(:,i)=u(:,s39+84240);  % columna 8
   i=i+1;
   s39=s39+1;
   G2=G2+1;
   end
   while L2<=8
   e(:,i)=u(:,s40+86400);  % columna 9
   i=i+1;
   s40=s40+1;
   L2=L2+1;
   end
   while Z2<=8
   e(:,i)=u(:,s41+88560);  % columna 10
   i=i+1;
   s41=s41+1;
   Z2=Z2+1;
   end
   while C2<=8  
   e(:,i)=u(:,s42+90720);  % columna 11
   i=i+1;
   s42=s42+1;
   C2=C2+1;
   end
   while V2<=8
   e(:,i)=u(:,s43+92880);  % columna 12
   i=i+1;
   s43=s43+1;
   V2=V2+1;
   end
   while QQ2<=8
   e(:,i)=u(:,s44+95040);  % columna 13
   i=i+1;
   s44=s44+1;
   QQ2=QQ2+1;
   end
    while QW2<=8
   e(:,i)=u(:,s45+97200);  % columna 14
   i=i+1;
   s45=s45+1;
   QW2=QW2+1;
   end
    while QH2<=8
   e(:,i)=u(:,s46+99360);  % columna 15
   i=i+1;
   s46=s46+1;
   QH2=QH2+1;
   end
    while QE2<=8
   e(:,i)=u(:,s47+101520);  % columna 16
   i=i+1;
   s47=s47+1;
   QE2=QE2+1;
   end
   
   while H3<=8
   e(:,i)=u(:,s48+103680);       % columna 1
   i=i+1;
   s48=s48+1;
   H3=H3+1;
   end
   while Q3<=8
   e(:,i)=u(:,s49+105840);  % columna 2
   i=i+1;
   s49=s49+1;
   Q3=Q3+1;
   end
   while W3<=8
   e(:,i)=u(:,s50+108000);  % columna 3
   i=i+1;
   s50=s50+1;
   W3=W3+1;
   end
   while E3<=8
   e(:,i)=u(:,s51+110160);  % columna 4
   i=i+1;
   s51=s51+1;
   E3=E3+1; 
   end
   while O3<=8
   e(:,i)=u(:,s52+112320);   % columna 5
   i=i+1;
   s52=s52+1;
   O3=O3+1;
   end
   while D3<=8
   e(:,i)=u(:,s53+114480);  % columna 6
   i=i+1;
   s53=s53+1;
   D3=D3+1;
   end
   while F3<=8
   e(:,i)=u(:,s54+116640);  % columna 7
   i=i+1;
   s54=s54+1;
   F3=F3+1;
   end
   while G3<=8
   e(:,i)=u(:,s55+118800);  % columna 8
   i=i+1;
   s55=s55+1;
   G3=G3+1;
   end
   while L3<=8
   e(:,i)=u(:,s56+120960);  % columna 9
   i=i+1;
   s56=s56+1;
   L3=L3+1;
   end
   while Z3<=8
   e(:,i)=u(:,s57+123120);  % columna 10
   i=i+1;
   s57=s57+1;
   Z3=Z3+1;
   end
   while C3<=8  
   e(:,i)=u(:,s58+125280);  % columna 11
   i=i+1;
   s58=s58+1;
   C3=C3+1;
   end
   while V3<=8
   e(:,i)=u(:,s59+127440);  % columna 12
   i=i+1;
   s59=s59+1;
   V3=V3+1;
   end
   while QQ3<=8
   e(:,i)=u(:,s60+129600);  % columna 13
   i=i+1;
   s60=s60+1;
   QQ3=QQ3+1;
   end
    while QW3<=8
   e(:,i)=u(:,s61+131760);  % columna 14
   i=i+1;
   s61=s61+1;
   QW3=QW3+1;
   end
    while QH3<=8
   e(:,i)=u(:,s62+133920);  % columna 15
   i=i+1;
   s62=s62+1;
   QH3=QH3+1;
   end
    while QE3<=8
   e(:,i)=u(:,s63+136080);  % columna 16
   i=i+1;
   s63=s63+1;
   QE3=QE3+1;
   end
end
end