circuito_decision.asv

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

function y = Circuito_Decision(u)

i=0; band=0;
if isempty(find_system(gcs,'Name','Decision'))~=1
    Blks=([gcs,'/','Decision']);
    snr_in_dB=str2num(get_param(Blks,'Potencia'));
    
%SNR=exp((snr_in_dB)*(log(10)/10));
sigma=(1/sqrt(2*SNR));
%==============================================
% Punto de decision
% Id = I(0)+I(1)/2
Id=(max(u(:,:))+min(u(:,:)))/2;
%---------------------------------------------
[a b]=size(u);
j=1;
r=zeros(1,b);
y=zeros(1,b);
%output=zeros(1,b);
while j<=b
        if u(1,j)<Id
           r(1,j) = 0+gngauss(sigma);
        else
           r(1,j) = 1+gngauss(sigma);
        end
        if r(1,j) < 0.5
            r(1,j) =0;
        else
            r(1,j) =1;
        end   
    j=j+1;
end
y=r;
end