regenerador.asv

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

function  [output]=Regenerador(u)%snr_in_dB,s)

we=([gcs,'/','Reg']);
snr_in_dB=str2num(get_param(we,'Potencia'));


% En dB
SNR=exp((snr_in_dB)*(log(10)/10));
sigma=(1/sqrt(2*SNR));


[a b]=size(u);
j=1;
r=zeros(1,b);
output=zeros(1,b);
while j<=b
     if u(1,j)==0
           r(1,j) = 0+gngauss(sigma);
        else
           r(1,j) = 1+gngauss(sigma);
        end
    % circuito de decision
        if r(1,j) < 0.5
            output(1,j) = 0;
        else
            output(1,j) = 1;
        end   
    j=j+1;
end
end

%calculo teorico
r=0:0.1:20;
snr=exp(r*log(10)/10)
teorica=Qfunct(sqrt(sr/2));
semilogy(r,teorica);
end