summaryrefslogtreecommitdiff
path: root/3543/CH2/EX2.7/Ex2_7.sce
blob: 9bc4986f12dcccfb665c038b50d83ae2e016bf7e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37

// Example 2.7
// Calculation of attenuation and Rayleigh scattering coefficient for fiber
// Page no 50
clc;
clear;
close;

// Given data
n=1.46;                                      // Refractive index
p=0.286;                                     // Average photoelastic coefficient
B=7.25*10^-11;                               // Isothermal compressibility
k=1.38*10^-23;                               // Boltzmann's constant
T=1350;                                      // Fictive temperature
l1=1*10^-6;                                  // Wavelength 1
l2=1.3*10^-6;                                // Wavelength 2
L=10^3;                                      // Length

// Rayleigh scattering coefficient for length 1
y1=8*(%pi)^3*(n)^8*(p)^2*B*k*T/(3*(l1)^4);   
// Rayleigh scattering coefficient for length 2
y2=8*(%pi)^3*(n)^8*(p)^2*B*k*T/(3*(l2)^4);   
y1=y1;
y2=y2;
//Attenuation 1
T1=exp(-(y1*L));                               
//Attenuation 2
T2=exp(-(y2*L));                               

//Display result on command window
printf("\n First Rayleigh scattering coefficient = %0.6f m^-1 ",y1);
printf("\n Second Rayleigh scattering coefficient = %0.6f m^-1 ",y2);

printf("\n Attenuation (@ Length 1) = %0.2f (dB/km) ",T1);
printf("\n Attenuation (@ Length 2) = %0.2f (dB/km) ",T2);