nrz.txt

对光纤系统中的NRZ时域和频域特性进行了描述

clear all 
T=2;   %时间
t=0:0.01:10;           %  NRZ在时域上的长度，取其中一段进行分析
omiga=2*pi/T;          %  NRZ码频域角频率
tao=1;                 %  初始幅度  
A=1                    %  幅度
f1=0;                  %  第一次时域描述
f2=0;
f3=0;
for k=1:1001          %   循环次数 ，进行1000次抽样
    f1=f1+((A*tao)/T)*((sin(k*omiga*tao/2)/(k*omiga*tao/2))*exp((k*omiga*t+31.4+...
        (tan(k*omiga-omiga))/0.0000005)*i));
end
figure(1)             %   画出NRZ时域形式
plot(t,f1)
xlabel('时间（sec）') %   对横坐标进行标注
ylabel('归一化幅度')  %   对纵坐标标注
title('NRZ码时域特性分析')
fftf1=fft(f1);
f=(0:length(fftf1)-1)*100/length(fftf1);
figure(2)
% plot(f,log10(fftf1.*conj(fftf1)))
% axis([0 500 -2 6])
plot(f,fftf1.*conj(fftf1)) %NRZ码的频域特性描述
%  axis([0 2 -150 350])
xlabel('频率（GHz）')
ylabel('功率（dB）')
title('NRZ码频域特性分析')
 

% 上述代码对光纤通信系统中的常用NRZ码时域和频域特性进行了分析，NRZ适
%  于短距离低速率光纤通信系统中使用，其频谱特性相对于RZ码频谱较窄，
%   但峰值功率较高，对非线性容忍度较小



%==============================================
%  对PMD进行了分析
%================================================
t=0:0.01:13;
 x=sqrt(2*pi)*3.75/4;
y=sqrt(2/pi).*(t.^2./x.^3).*exp((-t.^2)/(2.*x.^2));
figure(1)
plot(t,y.*4)
xlabel('时间')
ylabel('归一化强度')
title('PMD的Maxwell分布')
grid on
pmd=0.05:0.1:2;
L=(3.75./pmd).^2;
L1=(15./pmd).^2;
L2=(60./pmd).^2;
figure(2)
semilogy(pmd,(L),'--',pmd,(L1),'-+',pmd,(L2),'-o')
axis([0.05 2 0 10000000])
 xlabel('PMD')
 ylabel('传输距离')
 title('PMD各种速率传输距离的影响')

 legend('40G bit/s','10G bit/s','2.5G bit/s')