k1_k2.asv

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

function [h]=K1_K2(ent)

[aa bb]=size(ent);
%%%%%%%%%%%%%%
u=zeros(aa,720);
punt=zeros(1,10);

%%%%%%%%%%%%%%
u(:,:)=ent(:,1:720);
punt(:,:)=ent(1,721:730);
%%%%%%%%%%%%%%
~
for x=1:aa/9
    m2=(x-1)*9;
    
    %Estos vectores pertenecen a los denominados punteros H1 y H2 son
    %complemento de los mismos.
                  %New Rate Flag (NDF). Este valor es en "Operacion Normal".
                   u((4+m2),1:4)= [0 1 1 0];
    
                  %Spare Field
                   u((4+m2),5:6)= [0 0];
                   
                   %Coloca el valor del retraso de la trama introducido por el usuario.
                   u((4+m2),7:16)= punt(1,1:10);
    
    %K1
    u((5+m2),9:16) = [0 0 0 0 1 0 0 1];
    
    %K2
    u((5+m2),17:24) = [0 0 0 1 1 1 0 0];
    
    %D4 - D12
    u((6+m2),1:24)= [0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1];
    u((7+m2),1:24)= [0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1];
    u((8+m2),1:24)= [0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1];
     
     %S1
     u((9+m2),1:8)=[0 0 0 0 0 0 0 0];
     
     %M0
     u((9+m2),9:16)=[0 0 0 0 1 1 1 1];
     
     %E2
     u((9+m2),17:24)=[0 0 0 0 0 0 0 0];
 end
     
     %Salida
     h=u;