initial commit / add all books

4 files changed, 77 insertions, 0 deletions

diff --git a/803/CH6/EX6.1/ex6_1.sce b/803/CH6/EX6.1/ex6_1.sce
new file mode 100644
index 000000000..3e250b662
--- /dev/null
+++ b/803/CH6/EX6.1/ex6_1.sce
@@ -0,0 +1,32 @@
+clc

+f=1.2*(10^6);...............................//Assigning values to parameters

+Yn=56;

+deltalat=37.45;

+x=23.44*sind(0.9856*(Yn-80.7));

+deltah=acosd(tand(x)/tand(deltalat));

+y=acosd((sind(deltalat)*sind(x))+(cosd(deltalat)*cosd(x)*cosd(deltah)));

+disp("degree",x,"Solar declination");

+disp("degree",deltah,"hour angle of the sun");

+disp("degree",y,"Solar zenith angle");

+m=4.8*(10^(-26));................//mean molecular mass of air

+k=1.38*(10^(-23));...............//Boltzmann constant

+g=9.8;...........................//gravitational constant

+e=1.6*(10^(-19));....//electron charge

+me=9.11*(10^(-31));

+epsilon=8.85*(10^(-12));

+hmax=[100 200 300 120 250];

+T=[341 1360 1710 341 1540];

+f=1.2*(10^6);

+Nm=[1.5*10^11 3*10^11 12.5*10^11 0.8*1010 4*1011];

+h=[122 256 335 132 276];

+for i=1:5

+    H(i)=(g * m * [h(i)-hmax(i)])/ (k * T(i));

+    No(i)=(Nm(i)*(secd(x)^0.5));

+    Ne(i)=((No(i))*exp(0.5*(1-(H(i)-secd(x)*exp(-H(i))))));

+    fc(i)=(1/(2*3.14))*(sqrt(((e^2)*(Ne(i)))/(epsilon*me)));

+    n(i)=sqrt(1-((fc(i)^2)/(f^2)));

+    disp(Ne(i),"Electron density of each layer:");

+    disp(fc(i),"Critical frequency of each layer:");

+    disp(n(i),"Refractive index of each layer:");       

+end

+

diff --git a/803/CH6/EX6.2/ex6_2.sce b/803/CH6/EX6.2/ex6_2.sce
new file mode 100644
index 000000000..dc8c19a90
--- /dev/null
+++ b/803/CH6/EX6.2/ex6_2.sce
@@ -0,0 +1,22 @@
+clc

+deltalat=-5;..//geographic latitude

+ag=9;..//apparent elevation angle

+pi=3.14;

+x=6378.4*(10^3);..//equatorial radius of earth

+f=15*(10^6);

+r=[6493 6593 6693 6793];

+fc=[3.04*(10^6) 4.38*(10^6) 5.86*(10^6)];

+for i=1:4

+    j=1:3

+    n(j)=sqrt(1-((fc(j)^2)/f^2));..//refractive index

+    a(i)=acosd((x*cosd(ag))/r(i));..//apparent elevation angles

+    phi(i)=sind((r(i)*cosd(a(i)))/(r(i+1)));

+    theta(i)=(%pi/2)-a(i)-phi(i);

+    R012=sqrt((r(1)^2)+(r(4)^2)-(2*r(1)*r(4)*cosd(sum(theta(i)))));

+    R(i)=[r(i+1)*sind(theta(i))/cosd(a(i))];

+    dela=a(i)-acosd((r(4)/R012)*sind(sum(theta(i))));

+    delR=sum(R(i)/n(j))-R012;

+    disp(n(j),"Refractive index of each layer");

+    disp("degree",dela,"The refraction angle error");

+    disp("km",delR,"The range error");

+end;
\ No newline at end of file
diff --git a/803/CH6/EX6.3/ex6_3.sce b/803/CH6/EX6.3/ex6_3.sce
new file mode 100644
index 000000000..80b2582f6
--- /dev/null
+++ b/803/CH6/EX6.3/ex6_3.sce
@@ -0,0 +1,16 @@
+clc

+f=12*(10^6);..........................//Assigning values to parameters

+c=3*(10^8);

+Vt=85*(10^3);

+alphag=7.6;

+fc=5.6*(10^6);

+y=6377;

+ro=6527*(10^3);

+x=15;................................//orientation of target velocity

+ho=150;

+n=sqrt(1-((fc^2)/f^2));

+b=acosd(y*cosd(alphag)/(n*(ro+ho)));

+delf=((-2*f*Vt*sind(x)*b)/c)

+disp(n,"Refractive index");

+disp(b,"a");

+disp("Hz",delf,"The doppler frequency error is");

diff --git a/803/CH6/EX6.4/ex6_4.sce b/803/CH6/EX6.4/ex6_4.sce
new file mode 100644
index 000000000..5f78705e3
--- /dev/null
+++ b/803/CH6/EX6.4/ex6_4.sce
@@ -0,0 +1,7 @@
+clc

+T=[341 1360 1710 341 1540];

+N=[57.6*10^5 81.5*10^5 166.5*10^5 423 947];

+for i=1:5

+    lambda=69*sqrt(T(i)/N(i));

+    disp("mm",lambda,"The debye length is");

+end
\ No newline at end of file