clc
clear
close all
L1=50;
L2=100;
L3=200;
c=3e8;
lambda=[1:0.01:2];
N=length(lambda);
del_lambda=1;
n=1.45;
n=sqrt(1+((0.961663*(lambda).^2)/(lambda.^2-(0.684043)^2))+...
((0.4076426.*lambda.^2)/(lambda.^2-(0.1162414)^2))+...
((0.8974794.*lambda.^2)/(lambda.^2-(9.896161)^2)));
for k=1:N
Dm(k)=(-lambda(k)/c)*((6.51593e-9)-(1.94069e-6.*lambda(k).^2)-(0.00164847.*lambda(k).^4)+...
(1.04025.*lambda(k).^6)-(249.523.*lambda(k).^8)+32646.1.*lambda(k).^10-2.51246e6.*lambda(k).^12+...
1.07969e8.*lambda(k).^14-1.99892e9.*lambda(k).^16+5.6188e7.*lambda(k).^18+...
6.91612e8.*lambda(k).^20+7.16339e6.*lambda(k).^22-455368.*lambda(k).^24+...
3166.17.*lambda(k).^26-2.66595e-15.*lambda(k).^28)/((4*(-0.00619187+...
1.7816.*lambda(k).^2-97.9522.*lambda(k).^4+lambda(k).^6).^5)*(n.^3));
BR1(k)=1/(4*abs(Dm(k))*L1*del_lambda);
BR2(k)=1/(4*abs(Dm(k))*L2*del_lambda);
BR3(k)=1/(4*abs(Dm(k))*L3*del_lambda);
end
figure(1)
plot(Dm,BR1,'r',Dm,BR2,'b',Dm,BR3,'g')
legend('光纤长度：50km','光纤长度：100km','光纤长度：200km')
xlabel('色散（ps/nm.km）')
ylabel('比特率（Bps）')
title('色散和比特率之间的关系')
grid on
figure(2)
plot(lambda,Dm)
xlabel('波长（um）')
ylabel('色散（ps/km nm）')
title('光纤中的色散')
grid on
figure(3)
plot3(Dm,[1:101]-[1:101]+50,BR1,'r--',Dm,[100:200]-[100:200]+100,BR2,'b-',Dm,[300:400]-[300:400]+200,BR3,'r-.')
grid on
% axis([-3e-11 5e-11 0 2.5e10])
legend('传输50km','传输100km','传输200km')
- 1
- 2
前往页