project2.m

fiber material dispertion with different dopping

x=0;
c=3*10^8;
syms f
a1=0.6961603+0.110700*x;
a2=0.4079426+0.310321588*x;
a3=0.8974794-0.04331091*x;
b1=(0.0684043+0.00568306*x);
b2=(0.1162414+0.03772465*x);
b3=(9.896161+1.94577*x);
q=(a1.*(f.^2))./((f.^2)-b1.^2);
w=(a2.*(f.^2))./((f.^2)-b2.^2);
e=(a3.*(f.^2))./((f.^2)-b3.^2);
n=sqrt(1+q+w+e);
D=-diff(n,2)*10^12;
D1=(D.*f.*f)/(c.*f);
ezplot(D1,[1.1 1.7]);
set(findobj('Type','line'),'Color','r')
hold on
x=0.05;
c=3*10^8;
syms f
a1=0.6961603+0.110700*x;
a2=0.4079426+0.310321588*x;
a3=0.8974794-0.04331091*x;
b1=(0.0684043+0.00568306*x);
b2=(0.1162414+0.03772465*x);
b3=(9.896161+1.94577*x);
q=(a1.*(f.^2))./((f.^2)-b1.^2);
w=(a2.*(f.^2))./((f.^2)-b2.^2);
e=(a3.*(f.^2))./((f.^2)-b3.^2);
n=sqrt(1+q+w+e);
D=-diff(n,2)*10^12;
D1=(D.*f.*f)/(c.*f);
ezplot(D1,[1.1 1.7]);
set(findobj('Type','line','Color','b'),'Color','g')
hold on
x=0.1;
c=3*10^8;
syms f
a1=0.6961603+0.110700*x;
a2=0.4079426+0.310321588*x;
a3=0.8974794-0.04331091*x;
b1=(0.0684043+0.00568306*x);
b2=(0.1162414+0.03772465*x);
b3=(9.896161+1.94577*x);
q=(a1.*(f.^2))./((f.^2)-b1.^2);
w=(a2.*(f.^2))./((f.^2)-b2.^2);
e=(a3.*(f.^2))./((f.^2)-b3.^2);
n=sqrt(1+q+w+e);
D=-diff(n,2)*10^12;
D1=(D.*f.*f)/(c.*f);
ezplot(D1,[1.1 1.7]);
set(findobj('Type','line','Color','b'),'Color','c')
hold on
x=0.15;
c=3*10^8;
syms f
a1=0.6961603+0.110700*x;
a2=0.4079426+0.310321588*x;
a3=0.8974794-0.04331091*x;
b1=(0.0684043+0.00568306*x);
b2=(0.1162414+0.03772465*x);
b3=(9.896161+1.94577*x);
q=(a1.*(f.^2))./((f.^2)-b1.^2);
w=(a2.*(f.^2))./((f.^2)-b2.^2);
e=(a3.*(f.^2))./((f.^2)-b3.^2);
n=sqrt(1+q+w+e);
D=-diff(n,2)*10^12;
D1=(D.*f.*f)/(c.*f);
ezplot(D1,[1.1 1.7]);
set(findobj('Type','line','Color','b'),'Color','m')
hold on 
x=0.2;
c=3*10^8;
syms f
a1=0.6961603+0.110700*x;
a2=0.4079426+0.310321588*x;
a3=0.8974794-0.04331091*x;
b1=(0.0684043+0.00568306*x);
b2=(0.1162414+0.03772465*x);
b3=(9.896161+1.94577*x);
q=(a1.*(f.^2))./((f.^2)-b1.^2);
w=(a2.*(f.^2))./((f.^2)-b2.^2);
e=(a3.*(f.^2))./((f.^2)-b3.^2);
n=sqrt(1+q+w+e);
D=-diff(n,2)*10^12;
D1=(D.*f.*f)/(c.*f);
ezplot(D1,[1.1 1.7]);
xlabel('wavelength \lambda \mum');
ylabel('Material Dispersion,ps/nm.km');
title('');
grid;
hold off;
legend('x=0','x=0.05','x=0.1','x=0.15','x=0.2');
axis([1.1 1.7 -35 35]);