diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3547/CH10/EX10.6/EX10_6.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3547/CH10/EX10.6/EX10_6.sce')
| -rw-r--r-- | 3547/CH10/EX10.6/EX10_6.sce | 41 |
1 files changed, 41 insertions, 0 deletions
diff --git a/3547/CH10/EX10.6/EX10_6.sce b/3547/CH10/EX10.6/EX10_6.sce new file mode 100644 index 000000000..1d641faf7 --- /dev/null +++ b/3547/CH10/EX10.6/EX10_6.sce @@ -0,0 +1,41 @@ +// Example 10.6
+// to find the nonlinear phase shift at the center of the pulse. Compare the exact results with those obtained using first and second-order perturbation theory
+// Page no 469
+
+clc;
+clear;
+close;
+
+//Given data
+P=6*10^(-3); // The peak power of rectangular pulse
+L=40*10^3; // Fiber of length
+Floss=0.2; // The fiber loss (dB/Km)
+gamm=1.1*10^(-3);
+
+alpha=Floss/4.343; // Attenuation coefficient
+Zeff=(1-exp(-alpha*10^(-3)*L))/alpha*10^3;
+
+// The nonlinear phase shift at the center of the pulse
+phi=gamm*P*Zeff; // Nonlinear phase shift
+
+//Displaying results in the command window
+printf("\n The nonlinear phase shift at the center of the pulse = %0.4f rad ",phi);
+
+
+// Results using first order
+B01=sqrt(1+gamm^2*P^2*(Zeff)^2); // Amplitude shift
+thet1=atan(gamm*P*Zeff); // Non-linear phase shift
+
+//Displaying results in the command window
+printf("\n\n Amplitude shift using first order = %0.3f ",B01);
+printf("\n Non-linear shift using first order = %0.5f rad",thet1);
+
+// Results using second order
+x=1-((gamm)^2/2*P^2*Zeff^2);
+y=gamm*P*Zeff;
+thet2=atan(y/x); // Nonlinear phase shift
+B02=x/cos(thet2); // Amplitude shift
+
+//Displaying results in the command window
+printf("\n\n Amplitude shift using second order = %0.5f ",B02); // Answer is varying due to round-off error
+printf("\n Non-linear shift using second order = %0.5f rad",thet2); // Answer is varying due to round-off error
|