sustraetrama12a.m

来自「Matlab实现SONET/SDH光通信系统simulink仿真程序源代码」· M 代码 · 共 135 行

function y = sustraetrama12a(u)

if isempty(find_system(gcs,'Name','ADD_DROP_(STS12)'))~=1
   Blk2=([gcs,'/','ADD_DROP_(STS12)']);
   r=get_param(Blk2,'s');
            
[a b]=size(u);
tr=zeros(a,7568);
switch (r)
    case '12'    
      tr=u(:,1:7568);
    case '11'
      tr(:,1:6880)=u(:,1:6880);
      tr(:,6881:7568)=u(:,7569:8256);
    case '10'    
      tr(:,1:6192)=u(:,1:6192);
      tr(:,6193:7568)=u(:,6881:8256);
    case '9'    
      tr(:,1:5504)=u(:,1:5504);
      tr(:,5505:7568)=u(:,6193:8256);
    case '8'   
      tr(:,1:4816)=u(:,1:4816);
      tr(:,4817:7568)=u(:,5505:8256);
    case '7'   
      tr(:,1:4128)=u(:,1:4128);
      tr(:,4129:7568)=u(:,4817:8256);
    case '6'   
      tr(:,1:3440)=u(:,1:3440);
      tr(:,3441:7568)=u(:,4129:8256);
    case '5'    
      tr(:,1:2752)=u(:,1:2752);
      tr(:,2753:7568)=u(:,3441:8256);
    case '4'    
      tr(:,1:2064)=u(:,1:2064);
      tr(:,2065:7568)=u(:,2753:8256);
    case '3'    
      tr(:,1:1376)=u(:,1:1376);
      tr(:,1377:7568)=u(:,2065:8256);
    case '2'
      tr(:,1:688)=u(:,1:688);
      tr(:,687:7568)=u(:,1377:8256);
    case '1'    
      tr(:,1:7568)=u(:,687:8256);
end
y=tr;
end
if isempty(find_system(gcs,'Name','ADD_DROP_(STS12)1'))~=1
   Blk2111=([gcs,'/','ADD_DROP_(STS12)1']);
   r111=get_param(Blk2111,'s');
            
[a b]=size(u);
tr111=zeros(a,7568);
switch (r111)
    case '12'    
      tr111=u(:,1:7568);
    case '11'
      tr111(:,1:6880)=u(:,1:6880);
      tr111(:,6881:7568)=u(:,7569:8256);
    case '10'    
      tr111(:,1:6192)=u(:,1:6192);
      tr111(:,6193:7568)=u(:,6881:8256);
    case '9'    
      tr111(:,1:5504)=u(:,1:5504);
      tr111(:,5505:7568)=u(:,6193:8256);
    case '8'   
      tr111(:,1:4816)=u(:,1:4816);
      tr111(:,4817:7568)=u(:,5505:8256);
    case '7'   
      tr111(:,1:4128)=u(:,1:4128);
      tr111(:,4129:7568)=u(:,4817:8256);
    case '6'   
      tr111(:,1:3440)=u(:,1:3440);
      tr111(:,3441:7568)=u(:,4129:8256);
    case '5'    
      tr111(:,1:2752)=u(:,1:2752);
      tr111(:,2753:7568)=u(:,3441:8256);
    case '4'    
      tr111(:,1:2064)=u(:,1:2064);
      tr111(:,2065:7568)=u(:,2753:8256);
    case '3'    
      tr111(:,1:1376)=u(:,1:1376);
      tr111(:,1377:7568)=u(:,2065:8256);
    case '2'
      tr111(:,1:688)=u(:,1:688);
      tr111(:,687:7568)=u(:,1377:8256);
    case '1'    
      tr111(:,1:7568)=u(:,687:8256);
end
y=tr111;
end
if isempty(find_system(gcs,'Name','ADD_DROP_(STS12)2'))~=1
   Blk2222=([gcs,'/','ADD_DROP_(STS12)2']);
   r222=get_param(Blk2222,'s');
            
[a b]=size(u);
tr222=zeros(a,7568);
switch (r222)
    case '12'    
      tr222=u(:,1:7568);
    case '11'
      tr222(:,1:6880)=u(:,1:6880);
      tr222(:,6881:7568)=u(:,7569:8256);
    case '10'    
      tr222(:,1:6192)=u(:,1:6192);
      tr222(:,6193:7568)=u(:,6881:8256);
    case '9'    
      tr222(:,1:5504)=u(:,1:5504);
      tr222(:,5505:7568)=u(:,6193:8256);
    case '8'   
      tr222(:,1:4816)=u(:,1:4816);
      tr222(:,4817:7568)=u(:,5505:8256);
    case '7'   
      tr222(:,1:4128)=u(:,1:4128);
      tr222(:,4129:7568)=u(:,4817:8256);
    case '6'   
      tr222(:,1:3440)=u(:,1:3440);
      tr222(:,3441:7568)=u(:,4129:8256);
    case '5'    
      tr222(:,1:2752)=u(:,1:2752);
      tr222(:,2753:7568)=u(:,3441:8256);
    case '4'    
      tr222(:,1:2064)=u(:,1:2064);
      tr222(:,2065:7568)=u(:,2753:8256);
    case '3'    
      tr222(:,1:1376)=u(:,1:1376);
      tr222(:,1377:7568)=u(:,2065:8256);
    case '2'
      tr222(:,1:688)=u(:,1:688);
      tr222(:,687:7568)=u(:,1377:8256);
    case '1'    
      tr222(:,1:7568)=u(:,687:8256);
end
y=tr222;
end
end