voltage2current2.m

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

function I=voltage2current2(V)
%==============================================
% voltage2current= Esta funcion permite transformar el valor del voltage de la entrada
%                           a corriente.
%==============================================
%
%	AUTORES : Z. AGUI%O & D. STECKLER
%                       Universidad de Carabobo
%                       Valencia edo. Carabobo, VENEZUELA
%
%	FECHA     : Diciembre 2003
%	VERSION : 1.0
%
%==============================================
%	Realizado bajo MATLAB 6.5 Release 13/ Diciemdre 2003 (DZ)
%==============================================
% NOTA : SOLO VALIDO PARA 10 BLOQUES POR SIMULACION

if isempty(find_system(gcs,'Name','Regenerator'))~=1

    LBa=([gcs,'/','Regenerator']);
    LRa=get_param(LBa,'sel'); 
    LMa=get_param(LBa,'frec');
    [a b]=size(V);
    
    switch (LMa)
    case '1310'
        switch (LRa)
            case 'STS-12/STM-4'
                Ibiasa=0.012;
                Iopa=0.02-Ibiasa;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibiasa=0.012;
                Iopa=0.03-Ibiasa;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibiasa=0.011;
                Iopa=0.03-Ibiasa;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LRa)
            case 'STS-12/STM-4'
                Ibiasa=0.012;
                Iopa=0.02-Ibiasa;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibiasa=0.012;
                Iopa=0.035-Ibiasa;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibiasa=0.02;
                Iopa=0.065-Ibiasa;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LBa,'Iopa',num2str(Iopa));
 set_param(LBa,'Ibiasa',num2str(Ibiasa));
end

if isempty(find_system(gcs,'Name','Regenerator1'))~=1

    LB11=([gcs,'/','Regenerator1']);
    LR11=get_param(LB11,'sel'); 
    LM11=get_param(LB11,'frec');
    [a b]=size(V);
    
    switch (LM11)
    case '1310'
        switch (LR11)
            case 'STS-12/STM-4'
                Ibiasa11=0.012;
                Iopa11=0.02-Ibiasa11;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibiasa11=0.012;
                Iopa11=0.03-Ibiasa11;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibiasa11=0.011;
                Iopa11=0.03-Ibiasa11;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LR11)
            case 'STS-12/STM-4'
                Ibiasa11=0.012;
                Iopa11=0.02-Ibiasa11;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibiasa11=0.012;
                Iopa11=0.035-Ibiasa11;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibiasa11=0.02;
                Iopa11=0.065-Ibiasa11;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LB11,'Iopa',num2str(Iopa11));
 set_param(LB11,'Ibiasa',num2str(Ibiasa11));
end

if isempty(find_system(gcs,'Name','Regenerator2'))~=1

    LB22=([gcs,'/','Regenerator2']);
    LR22=get_param(LB22,'sel'); 
    LM22=get_param(LB22,'frec');
    [a b]=size(V);
    
    switch (LM22)
    case '1310'
        switch (LR22)
            case 'STS-12/STM-4'
                Ibiasa22=0.012;
                Iopa22=0.02-Ibiasa22;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibiasa22=0.012;
                Iopa22=0.03-Ibiasa22;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibiasa22=0.011;
                Iopa22=0.03-Ibiasa22;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LR22)
            case 'STS-12/STM-4'
                Ibiasa22=0.012;
                Iopa22=0.02-Ibiasa22;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibiasa22=0.012;
                Iopa22=0.035-Ibiasa22;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibiasa22=0.02;
                Iopa22=0.065-Ibiasa22;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LB22,'Iopa',num2str(Iopa22));
 set_param(LB22,'Ibiasa',num2str(Ibiasa22));
end

if isempty(find_system(gcs,'Name','Regenerator3'))~=1

    LB33=([gcs,'/','Regenerator3']);
    LR33=get_param(LB33,'sel'); 
    LM33=get_param(LB33,'frec');
    [a b]=size(V);
    
    switch (LM33)
    case '1310'
        switch (LR33)
            case 'STS-12/STM-4'
                Ibiasa33=0.012;
                Iopa33=0.02-Ibiasa33;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibiasa33=0.012;
                Iopa33=0.03-Ibiasa33;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibiasa33=0.011;
                Iopa33=0.03-Ibiasa33;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LR33)
            case 'STS-12/STM-4'
                Ibiasa33=0.012;
                Iopa33=0.02-Ibiasa33;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibiasa33=0.012;
                Iopa33=0.035-Ibiasa33;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibiasa33=0.02;
                Iopa33=0.065-Ibiasa33;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LB33,'Iopa',num2str(Iopa33));
 set_param(LB33,'Ibiasa',num2str(Ibiasa33));
end

if isempty(find_system(gcs,'Name','Regenerator4'))~=1

    LB44=([gcs,'/','Regenerator4']);
    LR44=get_param(LB44,'sel'); 
    LM44=get_param(LB44,'frec');
    [a b]=size(V);
    
    switch (LM44)
    case '1310'
        switch (LR44)
            case 'STS-12/STM-4'
                Ibiasa44=0.012;
                Iopa44=0.02-Ibiasa44;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibiasa44=0.012;
                Iopa44=0.03-Ibiasa44;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibiasa44=0.011;
                Iopa44=0.03-Ibiasa44;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LR44)
            case 'STS-12/STM-4'
                Ibiasa44=0.012;
                Iopa44=0.02-Ibiasa44;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibiasa44=0.012;
                Iopa44=0.035-Ibiasa44;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibiasa44=0.02;
                Iopa44=0.065-Ibiasa44;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LB44,'Iopa',num2str(Iopa44));
 set_param(LB44,'Ibiasa',num2str(Ibiasa44));
end

if isempty(find_system(gcs,'Name','Regenerator5'))~=1

    LB55=([gcs,'/','Regenerator5']);
    LR55=get_param(LB55,'sel'); 
    LM55=get_param(LB55,'frec');
    [a b]=size(V);
    
    switch (LM55)
    case '1310'
        switch (LR55)
            case 'STS-12/STM-4'
                Ibiasa55=0.012;
                Iopa55=0.02-Ibiasa55;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibiasa55=0.012;