atenuacion_fibra.m

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

function Atotal=atenuacion_fibra(u)
%==============================================
% atenuacion_fibra= Esta funcion permite calcular el valor de la Atenuacion Total de una
%                             Seccion Elemental de Fibra Optica:
%                                                            
%           Tx ____x______________x____ Rx
%         S                               R         
%              |--Seccion Elemental--|
%
%La Atenuacion sigue la siguiente ecuacion:
%
% Atotal=(1)*10^[-(alfa_n*long_total+alfa_s*Ntotal_emp+alfa_c*Ntotal_conec)/20]
%
%==============================================
%
%	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)
%==============================================
%
[a b]=size(u);
curBlk=([gcs]);
distancia=get_param(curBlk,'DistMaxRep');
long_onda=get_param(curBlk,'lambda');
Ntotal_conec=str2num(get_param(curBlk,'Nconec'));
long_total=str2num(get_param(curBlk,'LongTramo'));
%
%==============================================
%
%Selecciona el valor de la distancia maxima entre repetidores y el
%coeficiente de atenuacion
%
switch (distancia)
    case 'Short_Haul'
        switch (long_onda)
            case '1310'
                distancia_max=15;
                alfa_n=0.35;
            case '1550'
                distancia_max=15;
                alfa_n=0.20;
        end
    case 'Long_Haul'
        switch (long_onda)
            case '1310'
                distancia_max=40;
                alfa_n=0.35;
            case '1550'
                distancia_max=80;
                alfa_n=0.20;
        end
end
%
%Calculo del N