voltage2current.asv

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

function I=voltage2current(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
[a b]=size(V);
if isempty(find_system(gcs,'Name','Laser'))~=1

    LB=([gcs,'/','Laser']);
    LR=get_param(LB,'sel'); 
    LM=get_param(LB,'frec');
    
    
    switch (LM)
    case '1310'
        switch (LR)
            case 'STS-12/STM-4'
                Ibias=0.012;
                Iop=0.02-Ibias;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibias=0.012;
                Iop=0.03-Ibias;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibias=0.011;
                Iop=0.03-Ibias;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LR)
            case 'STS-12/STM-4'
                Ibias=0.012;
                Iop=0.02-Ibias;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibias=0.012;
                Iop=0.035-Ibias;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibias=0.02;
                Iop=0.065-Ibias;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LB,'Iop',num2str(Iop));
 set_param(LB,'Ibias',num2str(Ibias));
end

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

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

 set_param(LB11,'Iop',num2str(Iop11));
 set_param(LB11,'Ibias',num2str(Ibias11));
end

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

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

 set_param(LB22,'Iop',num2str(Iop22));
 set_param(LB22,'Ibias',num2str(Ibias22));
end

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

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

 set_param(LB33,'Iop',num2str(Iop33));
 set_param(LB33,'Ibias',num2str(Ibias33));
end

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

    LB44=([gcs,'/','Laser4']);
    LR44=get_param(LB44,'sel'); 
    LM=get_param(LB,'frec');
    
    
    switch (LM44)
    case '1310'
        switch (LR44)
            case 'STS-12/STM-4'
                Ibias44=0.012;
                Iop44=0.02-Ibias44;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibias44=0.012;
                Iop44=0.03-Ibias44;  %Las unidades de la resistencia son "Amperes".
            otherwise
                Ibias44=0.011;
                Iop44=0.03-Ibias44;  %Las unidades de la resistencia son "Amperes".
                
        end
        
    case '1550'
        switch (LR44)
            case 'STS-12/STM-4'
                Ibias44=0.012;
                Iop44=0.02-Ibias44;  %Las unidades de la resistencia son "Amperes".
                
            case 'STS-48/STM-16'
                Ibias44=0.012;
                Iop44=0.035-Ibias44;  %Las unidades de la resistencia son "Amperes".
                
            otherwise
                Ibias44=0.02;
                Iop44=0.065-Ibias44;  %Las unidades de la resistencia son "Amperes".
                
        end
    end

 set_param(LB44,'Iop',num2str(Iop44));
 set_param(LB44,'Ibias',num2str(Ibias44));
end

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

    LB55=([gcs,'/','Laser5']);
    LR55=get_param(LB55,'sel'); 
    LM55=get_param(LB55,'frec');
    
    
    switch (LM55)
    case '1310'
        switch (LR55)
            case 'STS-12/STM-4'
                Ibias55=0.012;
                Iop55=0.02-Ibias55;  %Las unidades de la resistencia son "Amperes".
            case 'STS-48/STM-16'
                Ibias55=0.012;